gen_crt_frags Namespace Reference

Functions
def	get_mod_id
def	get_mod_id_num
def	beam
def	beamVol
def	strip
def	module
def	minosSideTagger
def	minosNorthTagger
def	minosSouthTagger
def	DCTagger
def	cernTopTagger
def	cernLatRimTagger
def	cernLongRimTagger
def	detectorEnclosure

Variables
	testmode = False
	printModIds = False
float	WVWIDTH = 1031.8
	Parameters ##################### warm vessel (cm), increased as of 7/19/2019 to reflect full extent of WV profile in detector hall. More...
float	WVHEIGHT = 627.4
float	WVLENGTH = 2268.8
float	WVFOOTELEVATION = 10.16
float	ISLANDWIDTH = 118.0
float	IVLENGTH = 1996.0
float	TOPLEDGERISERTOFLOOR = 991.9
float	LEDGERISERHEIGTH = 46.99
float	TOPCRTBEAMTOFLOOR = 970.8
float	CRTBEAMSPACING = 92.71
float	BOTTOMCRTROLLERHEIGHT = 3.02
float	SIDECRTWVOFFSET = 4.13
float	SIDECRTPOSTWIDTH = 4.13
float	SIDECRTPOSTSPACING = 4.13
int	SIDECRTSTACKDEPTH = 3
float	SIDECRTSHELFTHICK = 0.56
float	overlap = 0.0065
float	gap = 0.005
float	crosstwomodule = 10.7
float	separatetwomodule = 0.6
float	offsetontopmodule = 48.26
float	verticalmodule = 485.14
float	CRTBEAMLENGTH = 1132.84
float	CRTBEAMHEIGHT = 25.3
float	CRTBEAMWIDTH = 25.4
float	CRTBEAMFLANGETHICK = 1.434
float	CRTBEAMWEBTHICK = 0.89407
int	CRTBEAMMASSDENS = 73
float	YM = 4.1
float	ZM = 800.0
float	XM = 1.0
float	MINOSSTRAIGHTSNOUT = 37.5
float	MINOSBENDSNOUT = 26.5
float	MINOSSNOUTLENGTH = 0.5
float	XC = 23.0
float	ZC = 184.0
float	YCTOP = 1.0
float	YCBOT = 1.5
float	XD = 5.0
float	ZD = 322.5
float	YD = 1.0
int	NXM = 20
int	NXC = 8
int	NXD = 32
float	PADM = 0.05
float	PADC = 0.1
float	PADD = 0.05
float	PADModule = 0.1
float	PADStrip = 0.01
float	PADTagger = 0.001
int	mModL = ZM+2
int	mModW = YM*NXM+2
int	mModH = XM+2
int	cModW = XC*NXC+2
int	cModH = YCTOP+YCBOT+2
int	dModH = YD*2
float	LAYERSPACE = 8.27
int	NMODSTACK = 9
int	NMODSTACKSOUTHY = 10
float	SLIDERSPACE = 18.5
float	STACKOVERLAP = 50.0
float	SIDECRTROLLOFFSET = 44.29
float	SIDECRTNORTHWALLTOFLOOR = 152.2
float	SIDECRTSOUTHWALLLATOFFSET = 1.1
float	DCSPACER = 32.6
tuple	LONGOFF5 = (3*ISLANDWIDTH+481.8)
tuple	LONGOFF2 = (ISLANDWIDTH+181.8)
float	CERNMODSPACE = 0.2
int	NTOPX = 6
int	NTOPZ = 14
int	NSLOPELAT = 14
int	NSLOPEFRONT = 6
float	SLOPEINCLINATION = 90.0
tuple	CERNROOFL = NTOPZ*cModW+(NTOPZ-1)
float	SHELLY = 1.1
float	MINOSSOUTHY = -0.5
	MINOSLATFIXY = MINOSSOUTHY
float	MINOSLATROLLY = MINOSLATFIXY-0.5
int	MINOSLATSOUTHACTIVEOVERHANG = 2
float	MINOSLATSOUTHZ = -0.5
float	MINOSLATNORTHZ = 0.5
float	MINOSLATCENTZ = 0.5
tuple	posDCInDetEncl = (0,-480.135, 0)
float	CERNRIMSWVOFFSET = 34.0
float	CERNTOPY = SHELLY*0.5
float	CERNRIMSY = CERNTOPY-0.5
float	CERNRIMSZ = -0.5
float	CERNRIMNY = CERNRIMSY+29.0
float	CERNRIMNZ = CERNROOFL+CERNRIMSZ+0.6
float	CERNTOPZ = CERNRIMSZ+0.5
float	CERNRIMLATX = 0.5
	CERNRIMLATY = CERNRIMSY
	CERNRIMLATZ = CERNTOPZ
float	SHELLZ = 1.01
float	SHELLWVOFFSET = SHELLZ*0.5
tuple	minosCutModLengthNorth = (256.54, 309.9, 309.9, 508.19, 508.19, 508.19)
tuple	minosCutModLengthSoutheast = (497.84, 497.84, 497.84, 497.84, 497.84, 497.84, 497.84, 497.84, 497.84)
tuple	minosCutModLengthSouthwest = (497.84, 497.84, 497.84, 497.84, 497.84, 497.84, 325.12, 325.12, 325.12)
float	MINOSNORTHY = -0.5
tuple	gdml = ET.Element('gdml')
tuple	solids = ET.SubElement(gdml, 'solids')
tuple	structure = ET.SubElement(gdml, 'structure')
dictionary	solids_store = {}
tuple	modToFeb = dict()
int	feb_id = 0
int	mod_id = -1
int	beam_id = 0
int	nModM = 0
int	nModC = 0
int	nModD = 0
tuple	m = ET.SubElement(materials, 'element', name='aluminum', formula='Al', Z='13')
	now if test mode generate materials, CRT shell, world, gdml header else just generate CRT shell for use with master geometry generator script More...
tuple	ws = ET.SubElement(solids, 'box', name='World', lunit="cm", x='1500', y='1500', z='3000')
tuple	w = ET.SubElement(structure, 'volume', name='volWorld')
tuple	pv = ET.SubElement(w, 'physvol')
string	posname = 'pos'
tuple	setup = ET.SubElement(gdml, 'setup', name='Default', version='1.0')

Detailed Description

Generate the GDML for the SBND CRT geometry.

We define a global XML DOM which the function for each volume will append to.
Volume-building functions are called hierarchically, so if you call module(),
it will construct the module and all the parts that make it up, so you end
up with complete GDML for one module.

Each physical volume has a corresponding unique logical volume, as required
by LArG4 to keep track of energy depositions. The solids, however, can safely
be referenced many times, and so are stored only once (using a hash keyed on
the the linear dimensions).

The output of this code is a file "crt.gdml" which contains the GDML snippets
to paste into the full SBND geometry.

Created by A. Mastbaum <mastbaum@uchicago.edu>, 2016/10/27
Downloaded by C. Hilgenberg <Chris.Hilgenberg@colostate.edu> 2017/10/25
Modified by C. Hilgenberg 2017/10/26

Function Documentation

def gen_crt_frags.beam

(

)

build one wide flange beam for top CRT support

Definition at line 212 of file gen_crt_frags.py.

def beam():
    '''build one wide flange beam for top CRT support'''

    global beam_id
    beam_id += 1

    xx = str(CRTBEAMWIDTH)
    yy = str(CRTBEAMHEIGHT)
    zz = str(CRTBEAMLENGTH)
    xxsub = str(0.5*(CRTBEAMWIDTH-CRTBEAMWEBTHICK))
    yysub = str(CRTBEAMHEIGHT-2*CRTBEAMFLANGETHICK)
    zzsub = zz
    xsubpos = 0.5*(CRTBEAMWIDTH-0.5*(CRTBEAMWIDTH-CRTBEAMWEBTHICK))
    area = CRTBEAMWIDTH*CRTBEAMHEIGHT-(CRTBEAMWIDTH-CRTBEAMWEBTHICK)*(CRTBEAMHEIGHT-2*CRTBEAMFLANGETHICK)
    if beam_id==1: print('modeled - true beam areas (cm^2): '+str(area-92.90304))

    sname = 'TopCRTSupportBeam'
    vname = 'vol'+sname+'_'+str(beam_id)

    if not sname in solids_store:
        snameext = sname+'_external'
        snameint = sname+'_internal'
        snamesub = sname+'_firstsubtraction'
        sexternal= ET.SubElement(solids, 'box', name=snameext, lunit="cm", x=xx, y=yy, z=zz)
        sinternal = ET.SubElement(solids, 'box', name=snameint, lunit="cm", x=xxsub, y=yysub, z=zzsub)
        ssub = ET.SubElement(solids, 'subtraction', name=snamesub)
        #s = ET.SubElement(solids, 'subtraction', name=sname)
        ET.SubElement(ssub, 'first', ref=snameext)
        ET.SubElement(ssub, 'second', ref=snameint)
        ET.SubElement(ssub, 'position', name='beamsubpos1', unit='cm', x=str(xsubpos), y='0', z='0')
        s = ET.SubElement(solids, 'subtraction', name=sname)
        ET.SubElement(s, 'first', ref=snamesub)
        ET.SubElement(s, 'second', ref=snameint)
        ET.SubElement(s, 'position', name='beamsubpos2', unit='cm', x=str(-1*xsubpos), y='0', z='0')
        solids_store[sname] = s

    else:
        s = solids_store[sname]

    v = ET.SubElement(structure, 'volume', name=vname) #Logical volume
    ET.SubElement(v, 'materialref', ref='STEEL_A992')
    ET.SubElement(v, 'solidref', ref=sname)

    return (s,v)