Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")

gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPndData")

gSystem	Load ("libPassive")

geoFace	setMediaFile ("/home/dkhaneft/Programs/PANDAroot/trunk/geometry/media_pnd.geo")

geoFace	readMedia ()

geoFace	print ()

gGeoMan	SetTopVolume (top)

SubunitVol9	AddNode (BoxVol5, b, trrotBox)

	if (b==1)

CrystalVol	SetLineColor (4)

	for (Int_t k=0;k< kNumOfCrystals;k++)

trrotQuar	SetName (name)

trrotQuar	RegisterYourself ()

trrotQuar	ReflectX (1)

top	AddNode (QuarterVol, 1, trrotQuar)

	for (Int_t q=1;q<=3;q++)

gGeoMan	CloseGeometry ()

top	Write ()

fi	Close ()

gGeoManager	Export (outfile)

top	Draw ()

gGeoManager	SetVisLevel (10)

Variables
const Double_t	kHeight = 2.44

const Double_t	kLength = 10.0

const Double_t	kSizeOfCrys = 2.44

const Double_t	kAlveoleThickness = 0.02

const Double_t	kGluedAlveoleThickness = 0.04

const Double_t	kCrystalCrystal = 0.06

const Double_t	kAlveoleAlveole = 0.08

const Double_t	kCrystalAlveole = 0.024

const Int_t	kNumOfQuarters = 4

const Int_t	kNumOfSubunit = 4

const Int_t	kNumOfBoxes4InSubunit = 4

const Int_t	kNumOfBoxes4 = 12

const Int_t	kNumOfBoxes2 = 10

const Int_t	kNumOfBoxes1 = 3

const Int_t	kNumOfCrystalsInBox4 = 4

const Int_t	kNumOfCrystalsInBox2 = 2

Double_t	sizeOfBox1 = kWidth + 2kAlveoleThickness + 2kCrystalAlveole

Double_t	sizeOfBox2 = 2kWidth + 2kAlveoleThickness + 2*kCrystalAlveole + kCrystalCrystal

Double_t	sizeOfBox4 = sizeOfBox2

Double_t	sizeOfSubunit = 2*sizeOfBox4

const Int_t	kNumOfBoxes = 4

const Int_t	kNumOfCrystals = 4

TString	outfile = "./emc_module4_StraightGeo24.4.root"

TFile *	fi = new TFile(outfile,"RECREATE")

FairGeoLoader *	geoLoad = new FairGeoLoader("TGeo","FairGeoLoader")

FairGeoInterface *	geoFace = geoLoad->getGeoInterface()

FairGeoMedia *	Media = geoFace->getMedia()

FairGeoBuilder *	geobuild =geoLoad->getGeoBuilder()

FairGeoMedium *	CbmMediumAir = Media->getMedium("air")

FairGeoMedium *	CbmMediumPWO = Media->getMedium("PWO")

FairGeoMedium *	CbmMediumCarbon = Media->getMedium("carbon")

FairGeoMedium *	CbmMediumAluminium = Media->getMedium("aluminium")

Int_t	nmed =geobuild->createMedium(CbmMediumAir)

TGeoManager *	gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom")

Double_t	sizeOfQuar = 37.644

TGeoVolume *	top = new TGeoVolumeAssembly("Emc4")

TGeoTranslation	trQuar = new TGeoTranslation(0.,0.,0.)

TGeoTranslation	trSub

TGeoTranslation	trBox

TGeoTranslation	trCrystal

TGeoRotation	rotQuar = new TGeoRotation()

TGeoRotation	rotSub

TGeoRotation	rotBox

TGeoRotation	rotCrystal

TGeoShape *	QuarterShape

TGeoShape *	SubunitShape

TGeoShape *	BoxShape

TGeoShape *	CrystalShape = new TGeoBBox(name,0.5kWidth,0.5kHeight,kLength,0)

TGeoVolume *	QuarterVol

TGeoVolume *	SubunitVol

TGeoVolume *	SubunitVol1

TGeoVolume *	SubunitVol2

TGeoVolume *	SubunitVol3

TGeoVolume *	SubunitVol4

TGeoVolume *	SubunitVol5

TGeoVolume *	SubunitVol6

TGeoVolume *	SubunitVol7

TGeoVolume *	SubunitVol8

TGeoVolume *	SubunitVol9

TGeoVolume *	BoxVol

TGeoVolume *	BoxVol1

TGeoVolume *	BoxVol2

TGeoVolume *	BoxVol3

TGeoVolume *	BoxVol4

TGeoVolume *	BoxVol5

TGeoVolume *	BoxVol6

TGeoVolume *	CrystalVol = new TGeoVolume(name,CrystalShape,gGeoMan->GetMedium(medium))

TString	name = "CrystalShape"

TString	name1

TString	name2

TString	name3

TString	name4

TString	name5

TString	name6

TString	name7

TString	name8

TString	name9

cout << "-----------------------------------------------> Quarter VOLUME<<endl;name="QuarterShape";QuarterShape=newTGeoBBox(name,sizeOfQuar0.5,sizeOfQuar0.5,10.5,0);name="QuarterNewVol";TStringmedium="air";QuarterVol=newTGeoVolumeAssembly(name);name="SubunitShape";SubunitShape=newTGeoBBox(name,2kWidth+2kAlveoleThickness+kCrystalCrystal+2kCrystalAlveole,2kHeight+2kAlveoleThickness+kCrystalCrystal+2kCrystalAlveole,kLength,0);TStringmedium="air";name="SubunitVol";name1="SubunitVol1";name2="SubunitVol2";name3="SubunitVol3";name4="SubunitVol4";name5="SubunitVol5";name6="SubunitVol6";name7="SubunitVol7";name8="SubunitVol8";name9="SubunitVol9";SubunitVol=newTGeoVolumeAssembly(name);SubunitVol1=newTGeoVolumeAssembly(name1);SubunitVol2=newTGeoVolumeAssembly(name2);SubunitVol3=newTGeoVolumeAssembly(name3);SubunitVol4=newTGeoVolumeAssembly(name4);SubunitVol5=newTGeoVolumeAssembly(name5);SubunitVol6=newTGeoVolumeAssembly(name6);SubunitVol7=newTGeoVolumeAssembly(name7);SubunitVol8=newTGeoVolumeAssembly(name8);SubunitVol9=newTGeoVolumeAssembly(name9);name="BoxShape";BoxShape=newTGeoBBox(name,kWidth+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,kHeight+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,kLength,0);TStringmedium="air";name="BoxVol";BoxVol=newTGeoVolumeAssembly(name);name1="BoxVol1";name2="BoxVol2";name3="BoxVol3";name4="BoxVol4";name5="BoxVol5";name6="BoxVol6";BoxVol1=newTGeoVolumeAssembly(name1);BoxVol2=newTGeoVolumeAssembly(name2);BoxVol3=newTGeoVolumeAssembly(name3);BoxVol4=newTGeoVolumeAssembly(name4);BoxVol5=newTGeoVolumeAssembly(name5);BoxVol6=newTGeoVolumeAssembly(name6);for(Int_tb=0;b <kNumOfBoxes;b++){cout <<""<<endl;cout <<"---------------->BOXnumber:" <<b<<endl;if(b==0){trBox=newTGeoTranslation(kSizeOfCrys+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,kSizeOfCrys+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,0.);rotBox=newTGeoRotation();rotBox.RotateX(0.);rotBox.RotateY(0.);}if(b==1){trBox=newTGeoTranslation(-kSizeOfCrys-kAlveoleThickness-0.5kCrystalCrystal-kCrystalAlveole,-kSizeOfCrys-kAlveoleThickness-0.5kCrystalCrystal-kCrystalAlveole,0.);rotBox=newTGeoRotation();rotBox.RotateX(0.);rotBox.RotateY(0.);}if(b==2){trBox=newTGeoTranslation(kSizeOfCrys+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,-kSizeOfCrys-kAlveoleThickness-0.5kCrystalCrystal-kCrystalAlveole,0.);rotBox=newTGeoRotation();rotBox.RotateX(0.);rotBox.RotateY(0.);}if(b==3){trBox=newTGeoTranslation(-kSizeOfCrys-kAlveoleThickness-0.5kCrystalCrystal-kCrystalAlveole,kSizeOfCrys+kAlveoleThickness+0.5kCrystalCrystal+kCrystalAlveole,0.);rotBox=newTGeoRotation();rotBox.RotateX(0.);rotBox.RotateY(0.);}TGeoCombiTrans*trrotBox=newTGeoCombiTrans(trBox,rotBox);name="BoxVol";name+=b;trrotBox-> SetName(name);trrotBox-> RegisterYourself();SubunitVol-> AddNode(BoxVol,b,trrotBox);if(b==0){name+=b;trrotBox-> SetName(name);trrotBox-> RegisterYourself();SubunitVol1-> AddNode(BoxVol1,b,trrotBox);SubunitVol4-> AddNode(BoxVol2,b,trrotBox);SubunitVol6-> AddNode(BoxVol6,b,trrotBox);SubunitVol7-> AddNode(BoxVol4,b,trrotBox);SubunitVol8->	AddNode (BoxVol, b, trrotBox)

TString	medium ="PWO"

Int_t	jj =0

TGeoCombiTrans *	trrotQuar = new TGeoCombiTrans(trQuar,rotQuar)

Function Documentation

SubunitVol9 AddNode	(	BoxVol5	,
		b	,
		trrotBox
	)

top AddNode	(	QuarterVol	,
		1	,
		trrotQuar
	)

fi Close ( )

gGeoMan CloseGeometry ( )

top Draw ( )

Definition at line 64 of file NHitsPerEvent.C.

            {
 
         std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;
 
 //      gHitPoints->Clear();
         ExtractData(mvdHitArray, gHitPoints);
         ExtractData(mvdProjectedArray, gProjections);
         //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());
 
         std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
         gHitPoints->Draw("apl");
         gProjections->Draw("pl");
 }

gGeoManager Export ( outfile )

for ( )

Definition at line 294 of file createRootGeoFileBwEnd_Str24.4.C.

                                        {
     cout << "               " << endl;
     cout << "         ----->  CRYSTAL number: "<< k << endl;
      if (k==0){     
         trCrystal = new TGeoTranslation(kSizeOfCrys/2. + 0.5*kCrystalCrystal,kSizeOfCrys/2. + 0.5*kCrystalCrystal,0.);
         rotCrystal = new TGeoRotation(); // "right-upper" Crystal
      }else if (k==1){
         trCrystal= new TGeoTranslation(-kSizeOfCrys/2. - 0.5*kCrystalCrystal,-kSizeOfCrys/2. - 0.5*kCrystalCrystal,0.); 
         rotCrystal = new TGeoRotation();  // "left-lower" Crystal
         rotCrystal.RotateZ(180.);      
      }else if (k==2){
         trCrystal= new TGeoTranslation(-kSizeOfCrys/2. - 0.5*kCrystalCrystal,kSizeOfCrys/2. + 0.5*kCrystalCrystal,0.); 
         rotCrystal = new TGeoRotation();  // "left-upper" Crystal
         rotCrystal.RotateZ(90.);       
      }else if (k==3){
         trCrystal= new TGeoTranslation(kSizeOfCrys/2. + 0.5*kCrystalCrystal,-kSizeOfCrys/2. - 0.5*kCrystalCrystal,0.); 
         rotCrystal = new TGeoRotation(); // "right-lower" Crystal
         rotCrystal.RotateZ(270.);      
      }
 
     TGeoCombiTrans* trrotCrystal= new TGeoCombiTrans(trCrystal,rotCrystal); 
     
     name = "CrystalVol";
     name+=k;
     trrotCrystal->SetName(name);
     trrotCrystal->RegisterYourself();
 
     BoxVol->AddNode(CrystalVol,k,trrotCrystal);
 
     if (k==0) {
       name+=k;
       trrotCrystal->SetName(name);
       trrotCrystal->RegisterYourself();
       
       BoxVol4->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol5->AddNode(CrystalVol,k,trrotCrystal);
     }
     if (k==1) {
       name+=k;
       trrotCrystal->SetName(name);
       trrotCrystal->RegisterYourself();
       
       BoxVol1->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol2->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol3->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol6->AddNode(CrystalVol,k,trrotCrystal);
     }
     if (k==2) {
       name+=k;
       trrotCrystal->SetName(name);
       trrotCrystal->RegisterYourself();
       
       BoxVol3->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol4->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol6->AddNode(CrystalVol,k,trrotCrystal);
     }
     if (k==3) {
       name+=k;
       trrotCrystal->SetName(name);
       trrotCrystal->RegisterYourself();
       
       BoxVol1->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol3->AddNode(CrystalVol,k,trrotCrystal);
       BoxVol5->AddNode(CrystalVol,k,trrotCrystal);
     }
   }

for ( Int_t q = 1; q<=3; q++ )

Definition at line 438 of file createRootGeoFileBwEnd_Str24.4.C.

References reflection.

                             {
     TGeoCombiTrans reflection; 
     TGeoTranslation ttt;
 
     if (q==1){                // second copy: -X -Y
       reflection.ReflectY(1); // right-lower Quarter   
       reflection.ReflectX(1);
       ttt = new TGeoTranslation(0.,0.,0.);//-0.5*sizeOfQuar,-0.5*sizeOfQuar,0.);
     }
     if (q==2){                // third copy: +X -Y
       reflection.ReflectY(1); // left-lower Quarter
       reflection.ReflectX(1); 
       reflection.ReflectX(1);
       ttt = new TGeoTranslation(0.,0.,0.);//0.5*sizeOfQuar,-0.5*sizeOfQuar,0.);
     }
     if (q==3){                // fourth copy: +X +Y
       reflection.ReflectX(1); // left-upper Quarter      
       reflection.ReflectX(1); 
       ttt = new TGeoTranslation(0.,0.,0.);//0.5*sizeOfQuar,0.5*sizeOfQuar,0.);
     }
 
     top->AddNode(QuarterVol,q+1,new TGeoCombiTrans(ttt,reflection)); 
   }

if ( b = =1 )

Definition at line 245 of file createRootGeoFileBwEnd_Str24.4.C.

References b.

              {
       name+=b;
       trrotBox->SetName(name);
       trrotBox->RegisterYourself();
 
       SubunitVol1->AddNode(BoxVol,b,trrotBox);
       SubunitVol2->AddNode(BoxVol,b,trrotBox);
       SubunitVol3->AddNode(BoxVol2,b,trrotBox);
       SubunitVol4->AddNode(BoxVol,b,trrotBox);
       SubunitVol5->AddNode(BoxVol,b,trrotBox);
       SubunitVol6->AddNode(BoxVol,b,trrotBox);
     }