Functions
gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPndData")

gSystem	Load ("libPassive")

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

geoFace	readMedia ()

geoFace	print ()

gGeoMan	SetTopVolume (top)

SuperModuleVolume	SetLineColor (2)

TyvekVolume	SetLineColor (3)

SuperModuleVolume	AddNode (TyvekVolume, 0)

	for (int ix=1;ix<=nClsSuperModuleX;ix++)

LeadVolume	SetLineColor (4)

LayerHolderVolume	AddNode (LeadVolume, 0, iCombTrans)

SciVolume	SetLineColor (5)

LayerHolderVolume	AddNode (SciVolume, 0, iCombTrans)

FiberHoleVolume	AddNode (FiberVolume, 0, iCombTrans)

fb	Close ()

MechFsc	Print ()

MechFsc	PrintNodes ()

rotMechFsc	RotateX (-90.)

rotMechFsc	RotateY (180.)

top	AddNode (MFscVol, 0, new TGeoCombiTrans(11.0, 4.8,-38.9, new TGeoRotation(rotMechFsc)))

gGeoMan	CloseGeometry ()

fi	cd ()

top	Write ()

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

const Int_t	nClsSuperModuleX = 2

const Int_t	nClsSuperModuleY = 2

Double_t	cellxsize = TMath::Max(Absorb_xsize,Sci_xsize)

Double_t	cellysize = TMath::Max(Absorb_ysize,Sci_ysize)

Double_t	suprModXsize = 2.cellxsize + 2.crystalspace + Air_thickness

Double_t	suprModYsize = 2.cellysize + 2.crystalspace + Air_thickness

Double_t	xsize = suprModXsize*numcellsx/nClsSuperModuleX

Double_t	ysize = suprModYsize*numcellsy/nClsSuperModuleY

Double_t	cell_thickness = 2.0*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 = "new_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 = "FscSuperModuleShape"

TGeoShape *	SuperModuleShape = new TGeoBBox(name,suprModXsize/2.,suprModYsize/2.,thickness/2.)

TString	medium = "carbon"

TGeoVolume *	SuperModuleVolume = new TGeoVolume(name, SuperModuleShape, gGeoMan->GetMedium(medium))

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

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

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

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

TGeoCombiTrans *	iCombTrans = NULL

TGeoRotation	rot

Int_t	iCopy = 0

Double_t	x

Double_t	y

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

TGeoVolume *	LayerHolderVolume = new TGeoVolumeAssembly(name)

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

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

TString	FileName = "fsc_cad2root_mech_support.root"

TFile *	fb = new TFile(FileName)

TGeoVolume *	MechFsc =(TGeoVolume *)fb->Get("topNode")

TGeoNode *	n =MechFsc->GetNode(0)

TGeoVolume *	MFscVol = n->GetVolume()

TGeoRotation	rotMechFsc

Function Documentation

SuperModuleVolume AddNode	(	TyvekVolume	,
		0
	)

LayerHolderVolume AddNode	(	LeadVolume	,
		0	,
		iCombTrans
	)

LayerHolderVolume AddNode	(	SciVolume	,
		0	,
		iCombTrans
	)

FiberHoleVolume AddNode	(	FiberVolume	,
		0	,
		iCombTrans
	)

top AddNode	(	MFscVol	,
		0	,
		new	TGeoCombiTrans11.0, 4.8,-38.9, new TGeoRotation(rotMechFsc)
	)

fi cd ( )

fb Close ( )

gGeoMan CloseGeometry ( )

top Draw ( "ogl" )

for ( int ix = 1; ix <= nClsSuperModuleX; ix++ )

!NO OVERLAPS

Definition at line 131 of file modifyRootFscGeometryFile.C.

References Double_t, iCopy, and nClsSuperModuleY.

    {
        for (int iy = 1; iy <= nClsSuperModuleY; iy++)
        {
            x = -cellxsize + 2.*(cellxsize/(Double_t)nClsSuperModuleX)*(ix - 0.5);
            y = -cellxsize + 2.*(cellxsize/(Double_t)nClsSuperModuleY)*(iy - 0.5);
            cout<<"===Inside supermodule"<<endl;
            cout <<"x="<<x<<" y="<<y<<endl;
            cout <<"ix="<<ix<<" iy="<<iy<<endl;
 
            iCombTrans = new TGeoCombiTrans(x,y,0, new TGeoRotation(rot));
            name = "FscModuleVolume";
            name+=iCopy;
            iCombTrans->SetName(name);
            iCombTrans->RegisterYourself();
            TyvekVolume->AddNode(ModuleVolume,iCopy,iCombTrans);
            iCopy++;
        }
    }