Functions
	basiclibs ()

gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPassive")

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

geoFace	readMedia ()

geoFace	print ()

geoBuild	createMedium (FairMediumAir)

geoBuild	createMedium (FairMediumVacuum)

geoBuild	createMedium (FairMediumFusedSil)

geoBuild	createMedium (FairMediumDIRCAir)

geoBuild	createMedium (FairMediumMirror)

cout<<"-I- overall top"<< endl;TGeoVolume vTop;TGeoBBox lTop=new TGeoBBox(200, 200, 200);vTop=new TGeoVolume("top", lTop, gGeoManager->GetMedium("air"));gGeoManager->	SetTopVolume (vTop)

cout<<"-I- calculating local mother"<< endl;TGeoVolume vLocalMother;TGeoPgon lLocalMother=new TGeoPgon(fPhi, 360., fEdges, 2);lLocalMother->	DefineSection (0,-fCoatThickness, fWindowHeightHalf-fAbsorberThickness-fCoatThickness, fDiskRMax+fMirrorHeight+fMcpHeight)

lLocalMother	DefineSection (1, fDiskThickness+fCoatThickness, fWindowHeightHalf-fAbsorberThickness-fCoatThickness, fDiskRMax+fMirrorHeight+fMcpHeight)

vTop	AddNode (vLocalMother, 0, new TGeoCombiTrans(0, 0, fDiskDistanceZ, new TGeoRotation(0)))

cout<<"-I- calculating radiator disk"<< endl;TGeoPgon *lDiskGlass=new TGeoPgon("DG", fPhi, 360., fEdges, 2);lDiskGlass->	DefineSection (0, 0., fDiskRMin, fDiskRMax)

lDiskGlass	DefineSection (1, fDiskThickness, fDiskRMin, fDiskRMax)

	if (fWindowIsBox)

trDW	RegisterYourself ()

vDisk	SetLineColor (colYellow)

vLocalMother	AddNode (vDisk, 0, new TGeoCombiTrans(0., 0., 0., new TGeoRotation(0)))

	if (fAbsorberUse)

cout<<"-I- calculating detectors (mother for mirror/mcp)" << endl;TGeoBBox lDetector[fDetectorTypes];TGeoVolume vDetector[fDetectorTypes];for(Int_t i=0;i< fDetectorTypes;i++){lDetector[i]=new TGeoBBox(fDetectorWidth/2., fDetectorHeight/2., fDiskThickness/2.);name="detector";name+=(i);vDetector[i]=new TGeoVolume(name.Data(), lDetector[i], gGeoManager-> GetMedium("DIRCair"));vDetector[i]->	SetLineColor (fDetectorColor[i])

cout<<"-I- calculating mirrors" << endl;TGeoBBox lMirror[fDetectorTypes];TGeoVolume vMirror[fDetectorTypes];for(Int_t i=0;i< fDetectorTypes;i++){lMirror[i]=new TGeoBBox(fDetectorWidth/2., fMirrorHeight/2., fDiskThickness/2.);name="mirror";name+=(i);medium="Mirror";vMirror[i]=new TGeoVolume(name.Data(), lMirror[i], gGeoManager-> GetMedium(medium.Data()));vMirror[i]->	SetLineColor (fMirrorColor[i])

vDetector[i]	AddNode (vMirror[i], 0, new TGeoCombiTrans(0.,-fDetectorHeight/2.+fMirrorHeight/2., 0., new TGeoRotation(0)))

cout<<"-I- calculating mcps" << endl;TGeoBBox lMcp=new TGeoBBox(fDetectorWidth/2., fMcpHeight/2., fDiskThickness/2.);TGeoVolume vMcp[fDetectorTypes];for(Int_t i=0;i< fDetectorTypes;i++){name="mcp";name+=(i);vMcp[i]=new TGeoVolume(name.Data(), lMcp, gGeoManager->GetMedium("FusedSil"));vMcp[i]->	SetLineColor (colGray)

vDetector[i]	AddNode (vMcp[i], 0, new TGeoCombiTrans(0.,-fDetectorHeight/2.+fMirrorHeight+fMcpHeight/2., 0., new TGeoRotation(0)))

gGeoManager	CloseGeometry ()

vTop	Write ()

fi	Close ()

cout<< "Done."<< endl;vTop->	Raytrace ()

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

FairGeoInterface *	geoFace = geoLoad->getGeoInterface()

FairGeoMedia *	geoMedia = geoFace->getMedia()

FairGeoBuilder *	geoBuild = geoLoad->getGeoBuilder()

TString	fGeoFile = "../../geometry/dsk.root"

Double_t const	cm = 1.

Double_t const	mm = cm/10.

Int_t const	colBlack = 1

Int_t const	colYellow = 5

Int_t const	colRed = 2

Int_t const	colBlue = 4

Int_t const	colGreen = 3

Int_t const	colGray = 15

Int_t const	fEdges = 8

Int_t const	fDetectorTypes = 2

Int_t const	fDetectorColor [fDetectorTypes] = {colRed,colBlue}

Int_t const	fMirrorColor [fDetectorTypes] = {colRed,colBlue}

Int_t const	fDetectorsPerEdge = 120

Double_t const	fAlpha = 360. /fEdges

Double_t const	fAlphaRad = fAlpha * TMath::DegToRad()

Double_t const	fPhi = fAlpha /2.

Double_t const	fCoatThickness = 1. *cm

Double_t const	fDiskThickness = 2. *cm

Double_t const	fDiskRMin = 3. *cm

Double_t const	fDiskDistanceZ = 196. *cm

Double_t const	fDiskRMax = fDiskDistanceZ * TMath::Tan(22.*TMath::DegToRad())

Double_t const	fWindowHeightHalf = fDiskDistanceZ * TMath::Tan( 5.*TMath::DegToRad())

Double_t const	fWindowWidthHalf = fDiskDistanceZ * TMath::Tan(10.*TMath::DegToRad())

Bool_t const	fWindowIsBox = kTRUE

Double_t const	fAbsorberThickness = 10 *mm

Bool_t const	fAbsorberUse = kTRUE

Double_t const	fMirrorHeight = 0.5 *mm

Double_t const	fMcpHeight = 1. *cm

Double_t const	fEdgeWidthHalf = fDiskRMax * TMath::Tan(fPhi * TMath::DegToRad())

Double_t const	fEdgeWidth = 2. * fEdgeWidthHalf

Double_t const	fDetectorWidth = fEdgeWidth/fDetectorsPerEdge

Double_t const	fDetectorHeight = fMirrorHeight + fMcpHeight

TVector3	edgeStartPos

TVector3	corner

TGeoRotation	nullRotation

TGeoRotation	rotation

TString	name

TString	medium ="air"

FairGeoMedium *	FairMediumAir = geoMedia->getMedium("air")

FairGeoMedium *	FairMediumVacuum = geoMedia->getMedium("vacuum")

FairGeoMedium *	FairMediumFusedSil = geoMedia->getMedium("FusedSil")

FairGeoMedium *	FairMediumDIRCAir = geoMedia->getMedium("DIRCair")

FairGeoMedium *	FairMediumMirror = geoMedia->getMedium("Mirror")

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

	vLocalMother = new TGeoVolume("mother", lLocalMother, gGeoManager->GetMedium("DIRCair"))

	else

TGeoTranslation *	trDW = new TGeoTranslation("trDW",0., 0., fDiskThickness/2.)

TGeoCompositeShape *	lDisk = new TGeoCompositeShape("DG - DW:trDW")

TGeoVolume *	vDisk = new TGeoVolume("radiator", lDisk, gGeoManager->GetMedium("FusedSil"))

cout<<"-I- placing detectors" << endl;edgeStartPos.SetXYZ(fEdgeWidthHalf-fDetectorWidth/2., fDiskRMax+fDetectorHeight/2., fDiskThickness/2.);for(Int_t iEdge=0;iEdge< fEdges;iEdge++){for(Int_t iDet=0;iDet< fDetectorsPerEdge;iDet++){corner=edgeStartPos;corner.SetX(corner.X()-iDet fDetectorWidth);corner.RotateZ(iEdge fAlphaRad);rotation=nullRotation;rotation.RotateZ(iEdge *fAlpha);vLocalMother->	AddNode (vDetector[(fDetectorsPerEdge iEdge+iDet)%fDetectorTypes], fDetectorsPerEdge iEdge+iDet, new TGeoCombiTrans(corner.X(), corner.Y(), corner.Z(), new TGeoRotation(rotation)))

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

Function Documentation

vTop AddNode	(	vLocalMother	,
		0	,
		new	TGeoCombiTrans0, 0, fDiskDistanceZ, new TGeoRotation(0)
	)

vLocalMother AddNode	(	vDisk	,
		0	,
		new	TGeoCombiTrans0., 0., 0., new TGeoRotation(0)
	)

vDetector [i] AddNode	(	vMirror	[i],
		0	,
		new	TGeoCombiTrans0.,-fDetectorHeight/2.+fMirrorHeight/2., 0., new TGeoRotation(0)
	)

vDetector [i] AddNode	(	vMcp	[i],
		0	,
		new	TGeoCombiTrans0.,-fDetectorHeight/2.+fMirrorHeight+fMcpHeight/2., 0., new TGeoRotation(0)
	)

basiclibs ( )

create root-geo-file for disk DIRC

fi Close ( )

gGeoManager CloseGeometry ( )

geoBuild createMedium ( FairMediumAir )

geoBuild createMedium ( FairMediumVacuum )

geoBuild createMedium ( FairMediumFusedSil )

geoBuild createMedium ( FairMediumDIRCAir )

geoBuild createMedium ( FairMediumMirror )

cout<<"-I- calculating local mother"<<endl; TGeoVolume* vLocalMother; TGeoPgon* lLocalMother = new TGeoPgon(fPhi, 360., fEdges, 2); lLocalMother-> DefineSection	(	0	,
		-	fCoatThickness,
		fWindowHeightHalf-fAbsorberThickness-	fCoatThickness,
		fDiskRMax+fMirrorHeight+	fMcpHeight
	)

lLocalMother DefineSection	(	1	,
		fDiskThickness+	fCoatThickness,
		fWindowHeightHalf-fAbsorberThickness-	fCoatThickness,
		fDiskRMax+fMirrorHeight+	fMcpHeight
	)

cout<<"-I- calculating radiator disk"<<endl; TGeoPgon* lDiskGlass = new TGeoPgon("DG",fPhi,360.,fEdges, 2); lDiskGlass-> DefineSection	(	0	,
		0.	,
		fDiskRMin	,
		fDiskRMax
	)

lDiskGlass DefineSection	(	1	,
		fDiskThickness	,
		fDiskRMin	,
		fDiskRMax
	)

if ( fWindowIsBox )

Definition at line 134 of file createRootGeoFile.C.

Referenced by PndFilteredPrimaryGenerator::AddFilter(), Check_particle_path(), PndSoftTriggerTask::Init(), PndSoftTriggerTask::ReadConfiguration(), PndFsmCmpDet::readParameters(), PndFsmCombiDet::readParameters(), and PndMvdMQFileSamplerBursts::Run().

                     {
     TGeoBBox* lDiskWindow = new TGeoBBox("DW",fWindowWidthHalf,fWindowHeightHalf,fDiskThickness);
   } else {

if ( fAbsorberUse )

Definition at line 149 of file createRootGeoFile.C.

                     {
   cout<<"-I- calculating absorber"<<endl;
   if (fWindowIsBox) {
   // in the window needs to be an absorber, (black colour) so the photons wont do
   // reflections back into the disk.
   // z = 0 ..... fDiskThickness
   TGeoBBox* lAbsorberOuter = new TGeoBBox("AO",fWindowWidthHalf,fWindowHeightHalf,fDiskThickness/2.);
   TGeoBBox* lAbsorberInner = new TGeoBBox("AI",fWindowWidthHalf-fAbsorberThickness,
                                                fWindowHeightHalf-fAbsorberThickness,fDiskThickness);
   } else {
   TGeoEltu* lAbsorberOuter = new TGeoEltu("AO",fWindowWidthHalf,fWindowHeightHalf,fDiskThickness/2.);
   TGeoEltu* lAbsorberInner = new TGeoEltu("AI",fWindowWidthHalf-fAbsorberThickness,
                                                fWindowHeightHalf-fAbsorberThickness,fDiskThickness);
   }
   TGeoCompositeShape* lAbsorber = new TGeoCompositeShape("AO - AI");
   TGeoVolume* vAbsorber = new TGeoVolume("absorber", lAbsorber, gGeoManager->GetMedium("vacuum"));
   vAbsorber->SetLineColor(colBlack);
   // again we need to shift them by the half fDiskThickness so place the center where it belongs
   vLocalMother->AddNode(vAbsorber,0,new TGeoCombiTrans(0., 0.,fDiskThickness/2., new TGeoRotation(0)));
   }

gSystem Load ( "libGeoBase" )

gSystem Load ( "libParBase" )

gSystem Load ( "libBase" )

gSystem Load ( "libPassive" )

rtdb print ( )

Referenced by DecayTreeFitter::Fitter::fillStream(), DecayTreeFitter::Constraint::filter(), and DecayTreeFitter::Fitter::fit().

cout<< "Done." << endl; vTop-> Raytrace ( )

geoFace readMedia ( )

trc4b RegisterYourself ( )

Referenced by createRootGeoFileBwEnd_2017(), hypGeGeoCOSYsetup2014(), hypGeGeoCOSYsetup2014Actives(), hypGeGeoCOSYsetup2014Passives(), and hypGeGeoCOSYTarget().

vDisk SetLineColor ( colYellow )

Referenced by create1StationGem_firstStation_v2(), create1StationGem_lastStation_v2(), create1StationsGem_firstStation(), create1StationsGem_lastStation(), create3StationsGem(), create3StationsGem_Tube(), create3StationsGem_v2(), create3StationsGem_v2_1(), create4StationsGem(), createRootGeoFileBwEnd_2017(), and PndGeoHypGeModifiedBeamPipe::PndGeoHypGeModifiedBeamPipe().

cout<<"-I- calculating detectors (mother for mirror/mcp)"<<endl; TGeoBBox* lDetector[fDetectorTypes]; TGeoVolume* vDetector[fDetectorTypes]; for (Int_t i=0; i<fDetectorTypes; i++) { lDetector[i] = new TGeoBBox(fDetectorWidth/2., fDetectorHeight/2., fDiskThickness/2.); name="detector"; name+=(i); vDetector[i] = new TGeoVolume(name.Data(),lDetector[i],gGeoManager->GetMedium("DIRCair")); vDetector[i]-> SetLineColor ( fDetectorColor [i] )

cout<<"-I- calculating mirrors"<<endl; TGeoBBox* lMirror[fDetectorTypes]; TGeoVolume* vMirror[fDetectorTypes]; for (Int_t i=0; i<fDetectorTypes; i++) { lMirror[i] = new TGeoBBox(fDetectorWidth/2., fMirrorHeight/2., fDiskThickness/2.); name="mirror"; name+=(i); medium="Mirror"; vMirror[i] = new TGeoVolume(name.Data(),lMirror[i],gGeoManager->GetMedium(medium.Data())); vMirror[i]-> SetLineColor ( fMirrorColor [i] )

cout<<"-I- calculating mcps"<<endl; TGeoBBox* lMcp = new TGeoBBox(fDetectorWidth/2.,fMcpHeight/2.,fDiskThickness/2.); TGeoVolume* vMcp[fDetectorTypes]; for (Int_t i=0; i<fDetectorTypes; i++) { name="mcp"; name+=(i); vMcp[i] = new TGeoVolume(name.Data(),lMcp,gGeoManager->GetMedium("FusedSil")); vMcp[i]-> SetLineColor ( colGray )

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

cout<<"-I- overall top"<<endl; TGeoVolume* vTop; TGeoBBox* lTop = new TGeoBBox(200,200,200); vTop = new TGeoVolume("top", lTop, gGeoManager->GetMedium("air")); gGeoManager-> SetTopVolume ( vTop )