Functions
gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPndData")

gSystem	Load ("libPassive")

geoFace	setMediaFile ("/home/freez/fairsoft_mar15/pandaroot/geometry/media_pnd.geo")

geoFace	readMedia ()

geoFace	print ()

gGeoMan	SetTopVolume (top)

ModuleVolume	SetLineColor (2)

TyvekVolume	SetLineColor (3)

ModuleVolume	AddNode (TyvekVolume, 0)

LeadVolume	SetLineColor (4)

LayerHolderVolume	AddNode (LeadVolume, 0, iCombTrans)

SciVolume	SetLineColor (5)

LayerHolderVolume	AddNode (SciVolume, 0, iCombTrans)

	for (int iz=1;iz<=EMCnumlayers;iz++)

FiberHoleVolume	AddNode (FiberVolume, 0, iCombTrans)

gGeoMan	CloseGeometry ()

top	Write ()

fi	Close ()

top	Draw ("ogl")

Variables
const Double_t	zpos = 760.0

const Int_t	numcellsx = 54

const Int_t	numcellsy = 28

const Double_t	tyvekbundle = 0.025

const Double_t	paperbundle = 0.025

Double_t	crystalspace = tyvekbundle+paperbundle

const Double_t	Absorb_xsize = 5.5

const Double_t	Absorb_ysize = 5.5

const Double_t	Air_thickness = 0.0

const Double_t	Pb_thickness = 0.0275

const Double_t	Sci_thickness = 0.15

const Double_t	Sci_xsize = 5.5

const Double_t	Sci_ysize = 5.5

const Int_t	EMCnumlayers = 380

const Int_t	numfibers = 6

const Double_t	fiber_r = 0.06

const Double_t	fiber_hole_r = 0.07

Double_t	cellxsize = TMath::Max(Absorb_xsize,Sci_xsize)+2.*crystalspace

Double_t	cellysize = TMath::Max(Absorb_ysize,Sci_ysize)+2.*crystalspace

Double_t	xsize = cellxsize*numcellsx

Double_t	ysize = cellysize*numcellsy

Double_t	cell_thickness = 2*Air_thickness+Pb_thickness+Sci_thickness

Double_t	thickness = cell_thicknessEMCnumlayers + 2.crystalspace

Double_t	holder_thickness = cell_thickness*EMCnumlayers

Double_t	hole_xsize = 6.*cellxsize

Double_t	hole_ysize = 6.*cellysize

Double_t	hole_xpos = 2.0*cellxsize

Double_t	hole_ypos = 1.0*cellxsize

TString	outfile = "emc_module5_fsc.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 *	MediumAir = Media->getMedium("air")

FairGeoMedium *	MediumTyvek = Media->getMedium("tyvek")

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

FairGeoMedium *	MediumLead = Media->getMedium("lead")

FairGeoMedium *	MediumScint = Media->getMedium("FscScint")

FairGeoMedium *	MediumFiber = Media->getMedium("FscFiber")

Int_t	nmed =geobuild->createMedium(MediumAir)

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

TGeoVolume *	top = new TGeoVolumeAssembly("Emc5")

TString	name = "FscModuleShape"

TGeoShape *	ModuleShape = new TGeoBBox(name,cellxsize/2.,cellysize/2.,thickness/2.)

TString	medium = "carbon"

TGeoVolume *	ModuleVolume = new TGeoVolume(name, ModuleShape, gGeoMan->GetMedium(medium))

TGeoShape *	TyvekShape = new TGeoBBox(name,cellxsize/2. - paperbundle,cellysize/2. - paperbundle,thickness/2. - paperbundle)

TGeoVolume *	TyvekVolume = new TGeoVolume(name, TyvekShape, gGeoMan->GetMedium(medium))

TGeoShape *	LayerHolderShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, cell_thickness/2.)

TGeoVolume *	LayerHolderVolume = new TGeoVolumeAssembly(name)

TGeoCombiTrans *	iCombTrans = NULL

TGeoRotation	rot

TGeoShape *	LeadShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, Pb_thickness/2.)

TGeoVolume *	LeadVolume = new TGeoVolume(name, LeadShape, gGeoMan->GetMedium(medium))

TGeoShape *	SciShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, Sci_thickness/2.)

TGeoVolume *	SciVolume = new TGeoVolume(name, SciShape, gGeoMan->GetMedium(medium))

Double_t	z

Int_t	iCopy = 0

Double_t	x

Double_t	y

TGeoShape *	FiberHoleShape = new TGeoTube(name, 0., fiber_hole_r, (holder_thickness)/2.)

TGeoVolume *	FiberHoleVolume = new TGeoVolume(name, FiberHoleShape, gGeoMan->GetMedium(medium))

TGeoShape *	FiberShape = new TGeoTube(name, 0., fiber_r, (holder_thickness)/2.)

TGeoVolume *	FiberVolume = new TGeoVolume(name, FiberShape, gGeoMan->GetMedium(medium))

Int_t	CopyCode =0

Function Documentation

ModuleVolume AddNode	(	TyvekVolume	,
		0
	)

LayerHolderVolume AddNode	(	LeadVolume	,
		0	,
		iCombTrans
	)

LayerHolderVolume AddNode	(	SciVolume	,
		0	,
		iCombTrans
	)

FiberHoleVolume AddNode	(	FiberVolume	,
		0	,
		iCombTrans
	)

fi Close ( )

gGeoMan CloseGeometry ( )

top Draw ( "ogl" )

for ( int iz = 1; iz <= EMCnumlayers; iz++ )

!NO OVERLAPS

Definition at line 149 of file createRootFscGeometryFile.C.

    {
            z = -holder_thickness/2. + 0.5 * cell_thickness *(2.*iz - 1.);
            iCombTrans = new TGeoCombiTrans(0,0,z, new TGeoRotation(rot));
            name = "FscLayerHolder";
            name += iz-1;
            TyvekVolume->AddNode(LayerHolderVolume,iz-1,iCombTrans);
    }