Functions
	for (Int_t i=0;i< 16;i++)

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

gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPndData")

gSystem	Load ("libPassive")

geoFace	setMediaFile ("../../geometry/media_pnd.geo")

geoFace	readMedia ()

geoFace	print ()

gGeoMan	SetTopVolume (top)

gGeoMan	CloseGeometry ()

top	Write ()

fi	Close ()

top	Draw ()

Variables
const Double_t	kSpaceInSub = 0.024

const Double_t	kSpaceSubGlue = 0.036

const Double_t	kAlveoleThickness = 0.018

const Int_t	kNumOfQuarters = 1

const Int_t	kNumOfBoxes = 4

const Int_t	kNumOfCrystals = 4

Double_t	tr = 2.51875

Double_t	b1 = 2.6000-0.5*tr

Double_t	b2 = 2.4375-0.5*tr

Double_t	dz = 10.000

Double_t	vert [20] = -0.5*tr

Double_t	vert0 [20]

Double_t	vertBox [20]

Double_t	vertSub [20]

Double_t	vertQuar [20]

Double_t	spacesQuar = 0.0

TString	outfile = "../../geometry/emc_module3new.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 = 96.192

TGeoVolume *	top = new TGeoVolumeAssembly("Emc3")

TGeoTranslation	trQuar

TGeoRotation	rotQuar

TGeoTranslation	trSub

TGeoTranslation	trBox

TGeoTranslation	trCrystal

TGeoRotation	rotSub

TGeoRotation	rotBox

TGeoRotation	rotCrystal

TGeoShape *	QuarterShape

TGeoShape *	SubunitShape

TGeoShape *	BoxShape

TGeoShape *	CrystalShape

TGeoVolume *	QuarterVol [4]

TGeoVolume *	SubunitVol [4]

TGeoVolume *	BoxVol [4]

TGeoVolume *	CrystalVol

TString	name

Function Documentation

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");
 }

for ( )

Definition at line 54 of file createRootGeoFileFwEnd_fixed.C.

References Double_t, i, and kSpaceInSub.

                              {
     Double_t spaces = 0.5*kSpaceInBox + kSpaceInSub;
     if (i==0 || i==1) 
       {
         vertBox[i] = (-vert[3] - spaces);   // lower-left corner back-side
       }
     else if (i==8 || i==9) 
       {
         vertBox[i] = (-vert[11] - spaces);  // lower-left corner front-side 
       }
     else if (i==7 || i==15)
       { 
         vertBox[i] = (-vert[i-1] - spaces); // lower-right points (front & back sides) 
       }
     else if (i==2 || i==10)                // upper-left points (front & back sides) 
       {
         vertBox[i] = (-vert[i+1] - spaces);
       }
     else 
       {
         vertBox[i] = (vert[i] + spaces);    // rest of those points
       }
   }