10 Oct
October 2009 posts
2009.10.02 Jason vs. CVS EEMC geometry: sampling fraction and shower shapes
Tests with Jason geometry file (ecalgeo.g23)
Monte-Carlo setup:
- One photon per event
- EEMC only geometry with LOW_EM option
- Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Color coding:
- Photon with Jason geometry (single particle MC)
- Photon with CVS (cAir fix) geometry (single particle MC)
- Eta-meson [pp2006 data] (single photons from eta-meson decay)
Sampling fraction
Figure 1: Sampling fraction vs. thrown energy (upper plot)
and vs. azimuthal angle (lower left) and rapidity (lower right)
Shower shapes
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)
Figure 2: Shower shapes
Shower shapes sorted by pre-shower energy
Pre-shower bins:
- Ep1 = 0, Ep2 = 0 (no energy in both EEMC pre-shower layers)
- Ep1 = 0, Ep2 > 0
- 0 < Ep1 < 4 MeV
- 4 < Ep1 < 10 MeV
- Ep1 > 10 MeV
- All pre-shower bins combined
Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch
Figure 3: Shower shapes (left) and their ratio (right)
Figure 4: Shower shape ratios
2009.10.05 Fix to the Jason geometry file
Why volume numbers has changed in Jason geometry file?
The number of nested volumes (nv),
is the total number of parent volumes for the sensitive volume
(sensitive volume is indicated by the HITS in the tree structure below).
For the Jason and CVS files this nv number seems to be the same
(see block tree structures below).
Then why volume ids id in g2t tables has changed?
The answer I found (which seems trivial to me know)
is that in the original (CVS) file ECAL
block has been instantiated (positioned) twice.
The second appearance is the prototype (East) version of the Endcap
(Original ecalgeo.g from CVS)
if (emcg_OnOff==1 | emcg_OnOff==3) then Position ECAL in CAVE z=+center endif if (emcg_OnOff==2 | emcg_OnOff==3) then Position ECAL in CAVE z=-center ThetaZ=180 endif
In Jason version the second appearance has been removed
(what seems natural and it should not have any effect)
(ecalgeo.g Jason edits, g23):
IF (emcg_OnOff>0) THEN Create ECAL ..... IF (emcg_OnOff==2 ) THEN Prin1 ('East Endcap has been removed from the geometry' ) ENDIF EndIF! emcg_OnOff
Unfortunately, this affects the way GEANT counts nested volumes
(effectively the total number was reduced by 1, from 8 to 7)
and this is the reason why the volume numbering scheme
in g2t tables has been affected.
I propose to put back these East Endcap line back,
since in this case it will not require any additional
changes to the EEMC decoder and g2t tables.
Block tree of the geometry file
blue - added volumes in Jason file
red - G10 volume removed in Jason file
HITS - sensitive volumes
---- Jason file ----
ECAL
EAGA
|EMSS
| -EFLP
| |ECVO
| | |EMOD
| | | |ESEC
| | | | |ERAD
| | | | | -ELED
| | | | |EMGT
| | | | | |EPER
| | | | | | |ETAR
| | | | | | | -EALP
| | | | | | | -ESCI -> HITS
| |ESHM
| | |ESPL
| | | |EXSG
| | | | -EXPS
| | | | -EHMS -> HITS
| | | | -EBLS
| | | | -EFLS
| | |ERSM
| -ESSP
| -ERCM
| -EPSB
|ECGH
| -ECHC
---- CVS file ----
ECAL
EAGA
|EMSS
| -EFLP
| |ECVO
| | |EMOD
| | | |ESEC
| | | | |ERAD
| | | | | -ELED
| | | | |EMGT
| | | | | |EPER
| | | | | | |ETAR
| | | | | | | -EALP
| | | | | | | -ESCI -> HITS
| |ESHM
| | |ESPL
| | | |EXSG
| | | | -EHMS -> HITS
| | | -EXGT
| | -ERSM
| -ESSP
| -ERCM
| -EPSB
|ECGH
| -ECHC
Block definitions
Jason geometry file
Create ECAL Block ECAL is one EMC EndCap wheel Create and Position EAGA AlphaZ=halfi EndBlock Block EAGA IS HALF OF WHEEL AIR VOLUME FORTHE ENDCAP MODULE Create AND Position EMSS konly='MANY' Create AND Position ECGH alphaz=90 kOnly='ONLY' EndBlock Block EMSS is the steel support of the endcap module Create AND Position EFLP z=zslice-center+zwidth/2 Create AND Position ECVO z=zslice-center+zwidth/2 Create AND Position ESHM z=zslice-center+zwidth/2 kOnly='MANY' Create AND Position ECVO z=zslice-center+zwidth/2 Create AND Position ESSP z=zslice-center+zwidth/2 Create ERCM Create EPSB EndBlock Block ECVO is one of endcap volume with megatiles and radiators Create AND Position EMOD alphaz=d3 ncopy=i_sector EndBlock Block ESHM is the shower maxsection Create and Position ESPL z=currentk Only='MANY' Create ERSM EndBlock Block ECGH is air gap between endcap half wheels Create ECHC EndBlock Block ECHC is steel endcap half cover EndBlock Block ESSP is stainless steelback plate EndBlock Block EPSB IS A PROJECTILE STAINLESS STEEL BAR EndBlock Block ERCM is stainless steel tie rod in calorimeter sections EndBlock Block ERSM is stainless steel tie rod in shower max EndBlock Block EMOD (fsect,lsect) IS ONE MODULEOF THE EM ENDCAP Create AND Position ESEC z=section-curr+secwid/2 EndBlock Block ESEC is a single em section Create AND Position ERAD z=length+(cell)/2+esec_deltaz Create AND Position EMGT z=length+(gap+cell)/2+esec_deltaz Create AND Position ERAD z=length+cell/2+esec_deltaz EndBlock Block EMGT is a 30 degree megatile Create AND Position EPER alphaz=myPhi EndBlock Block EPER is a 5 degree slice of a 30 degree megatile (subsector) Create and Position ETAR x=(rbot+rtop)/2ort=yzx EndBlock Block ETAR is a single calorimeter cell, containing scintillator, fiber router, etc... Create AND Position EALP y=(-megatile+emcs_alincell)/2 Create AND Position ESCI y=(-megatile+g10)/2+emcs_alincell _ EndBlock Block ESCI is the active scintillator (polystyrene) layer EndBlock Block ERAD is the lead radiator with stainless steel cladding Create AND Position ELED EndBlock Block ELED is a lead absorber plate EndBlock Block EFLP is the aluminum (aluminium) front plate of the endcap EndBlock Block EALP is the thin aluminium plate in calorimeter cell EndBlock Block ESPL is the logical volume containing an SMD plane Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' EndBlock Block EXSG Is another logical volume... this one acutally creates the planes Create and Position EXPS kONLY='MANY' Create and Position EHMS x=xc y=yc alphaz=-45 kOnly='ONLY' Create and Position EBLS x=xc y=yc z=(+esmd_apex/2+esmd_back_layer/2) alphaz=-45 kOnly='ONLY' Create and Position EHMS x=xc y=yc alphaz=-45 ort=x-y-z kOnly='ONLY' Create and Position EFLS x=xc y=yc z=(-esmd_apex/2-esmd_front_layer/2) alphaz=-45 ort=x-y-z kOnly='ONLY' EndBlock Block EHMS defines the triangular SMD strips Endblock! EHMS Block EFLS is the layer of material on the front of the SMD planes EndBlock! EFLS Block EBLS is the layer of material on the back of the SMD planes EndBlock! EFLS Block EXPS is the plastic spacer in the shower maximum section EndBlock
CVS geometry file
Create ECAL Block ECAL is one EMC EndCap wheel Create and Position EAGA AlphaZ=halfi EndBlock Block EAGA is half of wheel air volume forthe EndCap module Create and Position EMSS konly='MANY' Create and Position ECGH AlphaZ=90 konly='ONLY' EndBlock Block EMSS is steel support of the EndCap module Create and Position EFLP z=zslice-center+slcwid/2 Create and Position ECVO z=zslice-center+slcwid/2 Create and Position ESHM z=zslice-center+slcwid/2 Create and Position ECVO z=zslice-center+slcwid/2 Create and Position ESSP z=zslice-center+slcwid/2 Create ERCM Create EPSB EndBlock Block ECVO is one of EndCap Volume with megatiles and radiators Create and Position EMOD AlphaZ=d3 Ncopy=J_section EndBlock Block ESHM is the SHower Maxsection Create and Position ESPL z=current Create ERSM Endblock Block ECGH is air Gap between endcap Half wheels Create ECHC EndBlock Block ECHC is steel EndCap Half Cover EndBlock Block ESSP is Stainless Steelback Plate endblock Block EPSB is Projectile Stainless steel Bar endblock Block ERCM is stainless steel tie Rod in CaloriMeter sections endblock Block ERSM is stainless steel tie Rod in Shower Max endblock Block EMOD is one moduleof the EM EndCap Create and Position ESEC z=section-curr+secwid/2 endblock Block ESEC is a single EM section Create and Position ERAD z=len + (cell)/2 Create and Position EMGT z=len +(gap+cell)/2 Create and Position ERAD z=len + cell/2 Endblock Block EMGT is a megatile EM section Create and Position EPER AlphaZ=(emcs_Nslices/2-isec+0.5)*dphi Endblock Block EPER is a EM subsection period (super layer) Create and Position ETAR x=(RBot+RTop)/2ORT=YZX EndBlock Block ETAR is one CELL of scintillator, fiber and plastic Create and Position EALP y=(-mgt+emcs_AlinCell)/2 Create and Position ESCI y=(-mgt+G10)/2+emcs_AlinCell _ EndBlock Block ESCI is the active scintillator (polystyren) layer endblock Block ERAD is radiator Create and PositionELED endblock Block ELED is lead absorber Plate endblock Block EFLP is First Aluminium plate endblock Block EALP is ALuminiumPlate in calorimeter cell endblock Block ESPL is one of the Shower maxPLanes Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 Endblock Block EXSG is the Shower maxGap for scintillator strips Create EHMS endblock Block EHMS is sHower Max Strip Endblock Block EXGT is the G10 layer in the Shower Max EndBlock
Original (ecalgeo.g) file from CVS
Original (ecalgeo.g) file from CVS
****************************************************************************** Module ECALGEO is the EM EndCap Calorimeter GEOmetry Created 26 jan 1996 Author Rashid Mehdiyev * * Version 1.1, W.J. Llope * - changed sensitive medium names... * * Version 2.0, R.R. Mehdiyev 16.04.97 * - Support walls included * - intercell and intermodule gaps width updated * - G10 layers inserted * Version 2.1, R.R. Mehdiyev 23.04.97 * - Shower Max Detector geometry added * - Variable eta grid step size introduced * Version 2.2, R.R. Mehdiyev 03.12.97 * - Eta grid corrected * - Several changes in volume's dimensions * - Material changes in SMD * * Version 3.0, O. Rogachevsky 28.11.99 * - New proposal for calorimeter SN 0401 * * Version 4.1, O.Akio 3 Jan 01 * - Include forward pion detectors * Version 5.0, O. Rogachevsky 20.11.01 * - FPD is eliminated in this version * - More closed to proposal description * of calorimeter and SMD structure * ****************************************************************************** +CDE,AGECOM,GCONST,GCUNIT. * Content EAGA,EALP,ECAL,ECHC,ECVO,ECGH,EFLP,EHMS, ELED,EMGT,EMOD,EPER,EPSB,ERAD,ERCM,ERSM, ESHM,ESEC,ESCI,ESGH,ESPL,ESSP,EMSS, ETAR,EXGT,EXSG * Structure EMCG { Version, int Onoff, int fillMode} Structure EMCS { Type,ZOrig,ZEnd,EtaMin,EtaMax, PhiMin,PhiMax,Offset, Nsupsec,Nsector,Nsection,Nslices, Front,AlinCell,Frplast,Bkplast,PbPlate,LamPlate, BckPlate,Hub,Rmshift,SMShift,GapPlt,GapCel, GapSMD,SMDcentr,TieRod(2),Bckfrnt,GapHalf,Cover} * Structure EETR { Type,Etagr,Phigr,Neta,EtaBin(13)} * Structure ESEC { Isect, FPlmat, Cell, Scint, Nlayer } * Structure EMXG {Version,Sapex,Sbase,Rin,Rout,F4} * Structure EXSE {Jsect,Zshift,Sectype(6)} * Integer I_section,J_section,Ie,is,isec,i_str,Nstr,Type,ii,jj, cut,fsect,lsect,ihalf,filled * Real center,Plate,Cell,G10,diff,halfi, tan_low,tan_upp,Tanf,RBot,Rtop,Deta,etax,sq2,sq3, dup,dd,d2,d3,rshift,dphi,radiator,orgkeep,endkeep * Real maxcnt,msecwd,mxgten,curr,Secwid,Section, curcl,EtaTop,EtaBot,slcwid,zslice,Gap,mgt, xleft,xright,yleft,yright,current, rth,len,p,xc,yc,xx,yy,rbotrad, Rdel,dxy,ddn,ddup Integer N Parameter (N=12) * Tanf(etax) = tan(2*atan(exp(-etax))) * * ---------------------------------------------------------------------------- * * FillMode =1 only 2-5 sectors (in the first half) filled with scintillators * FillMode =2 all sectors filled (still only one half of one side) * FillMode =3 both halves (ie all 12 sectors are filled) Fill EMCG ! EM EndCAp Calorimeter basic data Version = 5.0 ! Geometry version OnOff = 3 ! Configurations 0-no, 1-west 2-east 3-both FillMode = 3 ! sectors fill mode Fill EMCS ! EM Endcap Calorimeter geometry Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrig = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.400 ! space for SMD detector SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half * Rmshift = 2.121 ! radial shift of module * -------------------------------------------------------------------------- Fill EETR ! Eta and Phi grid values Type = 1 ! =1 endcap, =2 fpd EtaGr = 1.0536 ! eta_top/eta_bot tower granularity PhiGr = 0.0981747 ! Phi granularity (radians) NEta = 12 ! Eta granularity EtaBin = {2.0,1.9008,1.8065,1.7168,1.6317,1.5507,1.4738, 1.4007,1.3312,1.2651,1.2023,1.1427,1.086}! Eta rapidities *--------------------------------------------------------------------------- Fill ESEC ! First EM section ISect = 1 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness * Fill ESEC ! First EM section ISect = 2 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness * Fill ESEC ! Second EM section ISect = 3 ! Section number Nlayer = 4 ! Number of Sci layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness * Fill ESEC ! Third EM section ISect = 4 ! Section Nlayer = 18 ! Number of layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness * Fill ESEC ! 4th EM section ISect = 5 ! Section Nlayer = 1 ! Number of layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness *---------------------------------------------------------------------------- Fill EMXG ! EM Endcap SMD basic data Version = 1 ! Geometry version Sapex = 0.7 ! Scintillator strip apex Sbase = 1.0 ! Scintillator strip base Rin = 77.41 ! inner radius of SMD plane Rout = 213.922 ! outer radius of SMD plane F4 = .15 ! F4 thickness *---------------------------------------------------------------------------- Fill EXSE ! First SMD section JSect = 1 ! Section number Zshift = -1.215 ! Section width sectype = {4,1,0,2,1,0} ! 1-V,2-U,3-cutV,4-cutU * Fill EXSE ! Second SMD section JSect = 2 ! Section number Zshift = 0. ! Section width sectype = {0,2,1,0,2,3} ! 1-V,2-U,3-cutV,4-cutU * Fill EXSE ! Third SMD section JSect = 3 ! Section number Zshift = 1.215 ! Section width sectype = {1,0,2,1,0,2} ! 1-V,2-U,3-cutV,4-cutU *---------------------------------------------------------------------------- * Use EMCG * sq3 = sqrt(3.) sq2 = sqrt(2.) prin1 emcg_version ('ECALGEO version ', F4.2) * Endcap USE EMCS type=1 USE EETR type=1 orgkeep = emcs_ZOrig endkeep = emcs_ZEnd if(emcg_OnOff>0) then diff = 0.0 center = (emcs_ZOrig+emcs_ZEnd)/2 Tan_Upp = tanf(emcs_EtaMin) Tan_Low = tanf(emcs_EtaMax) rth = sqrt(1. + Tan_Low*Tan_Low) rshift = emcs_Hub * rth dup=emcs_Rmshift*Tan_Upp dd=emcs_Rmshift*rth d2=rshift + dd radiator = emcs_Pbplate + 2*emcs_LamPlate * d3=emcs_Rmshift-2*emcs_smshift dphi = (emcs_PhiMax-emcs_PhiMin)/emcs_Nsector Create ECAL if (emcg_OnOff==1 | emcg_OnOff==3) then Position ECAL in CAVE z=+center endif if (emcg_OnOff==2 | emcg_OnOff==3) then Position ECAL in CAVE z=-center ThetaZ=180 endif if(section > emcs_Zend) then prin0 section,emcs_Zend (' ECALGEO error: sum of sections exceeds maximum ',2F12.4) endif prin1 section (' EndCap calorimeter total depth ',F12.4) endif prin1 ('ECALGEO finished') * * ---------------------------------------------------------------------------- Block ECAL is one EMC EndCap wheel Material Air Medium standard Attribute ECAL seen=1 colo=7 ! lightblue shape CONE dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup do ihalf=1,2 filled=1 halfi = -105 + (ihalf-1)*180 if (ihalf=2 & emcg_FillMode<3) filled = 0 Create and Position EAGA AlphaZ=halfi enddo * EndBlock * ---------------------------------------------------------------------------- Block EAGA is half of wheel air volume for the EndCap module Attribute EAGA seen=1 colo=1 serial=filled ! black Material Air shape CONS dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax if (filled=1) then Create and Position EMSS konly='MANY' curr = orgkeep ; curcl = endkeep Create and position ECGH AlphaZ=90 konly='ONLY' endif EndBlock * ---------------------------------------------------------------------------- Block EMSS is steel support of the EndCap module Attribute EMSS seen=1 colo=1 ! black Material Iron shape CONS dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax zslice = emcs_ZOrig prin1 zslice (' Front Al plane starts at: ',F12.4) slcwid = emcs_Front Create and Position EFLP z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' First calorimeter starts at: ',F12.4) fsect = 1; lsect = 3 slcwid = emcs_SMDcentr - emcs_GapSMD/2 - zslice * Create and Position ECVO z=zslice-center+slcwid/2 slcwid = emcs_GapSMD zslice = emcs_SMDcentr - emcs_GapSMD/2 prin1 section,zslice (' 1st calorimeter ends, SMD starts at: ',2F10.5) Create and Position ESHM z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' SMD ends at: ',F10.5) * slcwid = 0 fsect = 4; lsect = 5 do I_section =fsect,lsect USE ESEC Isect=I_section Slcwid = slcwid + esec_cell*esec_Nlayer enddo slcwid = emcs_bckfrnt - zslice * Create and Position ECVO z = zslice-center+slcwid/2 zslice = emcs_bckfrnt prin1 section,zslice (' 2nd calorimeter ends, Back plate starts at: ',2F10.5) slcwid = emcs_BckPlate * Create and Position ESSP z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' BackPlate ends at: ',F10.5) slcwid = emcs_Zend-emcs_ZOrig Create ERCM do i_str = 1,2 do is = 1,5 xx = emcs_phimin + is*30 yy = xx*degrad xc = cos(yy)*emcs_TieRod(i_str) yc = sin(yy)*emcs_TieRod(i_str) Position ERCM z=0 x=xc y=yc enddo enddo rth = orgkeep*Tan_Upp+dup + 2.5/2 xc = (endkeep - orgkeep)*Tan_Upp len = .5*(endkeep + orgkeep)*Tan_Upp + dup + 2.5/2 yc = emcs_Zend-emcs_ZOrig p = atan(xc/yc)/degrad Create EPSB do is = 1,6 xx = -75 + (is-1)*30 yy = xx*degrad xc = cos(yy)*len yc = sin(yy)*len Position EPSB x=xc y=yc AlphaZ=xx enddo EndBlock * ---------------------------------------------------------------------------- Block ECVO is one of EndCap Volume with megatiles and radiators Material Air Attribute ECVO seen=1 colo=3 ! green shape CONS dz=slcwid/2, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup do J_section = 1,6 if (1 < J_section < 6 | emcg_FillMode > 1)then filled = 1 else filled = 0 endif d3 = 75 - (J_section-1)*30 Create and Position EMOD AlphaZ=d3 Ncopy=J_section enddo * EndBlock * ---------------------------------------------------------------------------- Block ESHM is the SHower Max section * Material Air Attribute ESHM seen=1 colo=4 ! blue Shape CONS dz=SlcWid/2, rmn1=zslice*Tan_Low-dd, rmn2=(zslice+slcwid)*Tan_Low-dd, rmx1=(zslice)*Tan_Upp+dup, rmx2=(zslice+slcwid)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax USE EMXG Version=1 maxcnt = emcs_SMDcentr prin1 zslice,section,center (' Z start for SMD,section: ',3F12.4) * do J_section = 1,3 USE EXSE Jsect=J_section * current = exse_Zshift secwid = emxg_Sapex + 2.*emxg_F4 section = maxcnt + exse_zshift prin1 j_section,current,section,secwid (' layer, Z, width : ',i3,3F12.4) rbot=section*Tan_Low rtop=section*Tan_Upp prin1 j_section,rbot,rtop (' layer, rbot,rtop : ',i3,2F12.4) Create and position ESPL z=current * end do Create ERSM do i_str = 1,2 do is = 1,5 xx = emcs_phimin + (is)*30 yy = xx*degrad xc = cos(yy)*emcs_TieRod(i_str) yc = sin(yy)*emcs_TieRod(i_str) Position ERSM z=0 x=xc y=yc enddo enddo Endblock * ---------------------------------------------------------------------------- Block ECGH is air Gap between endcap Half wheels Material Air Medium standard Attribute ECGH seen=0 colo=7 ! lightblue shape TRD1 dz=(emcs_Zend-emcs_ZOrig)/2, dy =(emcs_gaphalf+emcs_cover)/2, dx1=orgkeep*Tan_Upp+dup, dx2=endkeep*Tan_Upp+dup rth = emcs_GapHalf + emcs_cover xx=curr*Tan_Low-d2 xleft = sqrt(xx*xx - rth*rth) yy=curr*Tan_Upp+dup xright = sqrt(yy*yy - rth*rth) secwid = yy - xx xx=curcl*Tan_Low-d2 yleft = sqrt(xx*xx - rth*rth) yy=curcl*Tan_Upp+dup yright = sqrt(yy*yy - rth*rth) slcwid = yy - xx xx=(xleft+xright)/2 yy=(yleft + yright)/2 xc = yy - xx len = (xx+yy)/2 yc = curcl - curr p = atan(xc/yc)/degrad rth = -(emcs_GapHalf + emcs_cover)/2 Create ECHC Position ECHC x=len y=rth Position ECHC x=-len y=rth AlphaZ=180 EndBlock * ---------------------------------------------------------------------------- Block ECHC is steel EndCap Half Cover Attribute ECHC seen=1 colo=1 ! black Material Iron shape TRAP dz=(curcl-curr)/2, thet=p, bl1=secwid/2, tl1=secwid/2, bl2=slcwid/2, tl2=slcwid/2, h1=emcs_cover/2 h2=emcs_cover/2, phi=0 alp1=0 alp2=0 EndBlock * ---------------------------------------------------------------------------- Block ESSP is Stainless Steel back Plate * Material Iron Attribute ESSP seen=1 colo=6 fill=1 shape CONS dz=emcs_BckPlate/2, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax endblock * ---------------------------------------------------------------------------- Block EPSB is Projectile Stainless steel Bar * Material Iron Attribute EPSB seen=1 colo=6 fill=1 shape TRAP dz=(emcs_Zend-emcs_ZOrig)/2, thet=p, bl1=2.5/2, tl1=2.5/2, bl2=2.5/2, tl2=2.5/2, h1=2.0/2 h2=2.0/2, phi=0 alp1=0 alp2=0 endblock * ---------------------------------------------------------------------------- Block ERCM is stainless steel tie Rod in CaloriMeter sections * Material Iron Attribute ERSM seen=1 colo=6 fill=1 shape TUBE dz=slcwid/2, rmin=0, rmax=1.0425 ! nobody knows exactly endblock * ---------------------------------------------------------------------------- Block ERSM is stainless steel tie Rod in Shower Max * Material Iron Attribute ERSM seen=1 colo=6 fill=1 shape TUBE dz=slcwid/2, rmin=0, rmax=1.0425 endblock * ---------------------------------------------------------------------------- Block EMOD is one module of the EM EndCap Attribute EMOD seen=1 colo=3 serial=filled ! green Material Air Shape CONS dz=slcwid/2, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup * * Running parameter 'section' contains the position of the current section * It should not be modified in daughters, use 'current' variable instead. * SecWid is used in all 'CONS' daughters to define dimensions. * * section = zslice curr = zslice + slcwid/2 Do I_section =fsect,lsect USE ESEC Isect=I_section * Secwid = esec_cell*esec_Nlayer if (I_section = 3 | I_section = 5) then ! no last radiator Secwid = Secwid - radiator else if (I_section = 4) then ! add one more radiator Secwid = Secwid - esec_cell + radiator endif Create and position ESEC z=section-curr+secwid/2 section = section + secwid * enddo endblock * ---------------------------------------------------------------------------- Block ESEC is a single EM section Attribute ESEC seen=1 colo=1 serial=filled Material Air Medium standard * Shape CONS dz=secwid/2, rmn1=(section-diff)*Tan_Low-dd, rmn2=(section+secwid-diff)*Tan_Low-dd, rmx1=(section-diff)*Tan_Upp+dup, rmx2=(section+secwid-diff)*Tan_Upp+dup * len = -secwid/2 current = section mgt = esec_scint + emcs_AlinCell _ + emcs_FrPlast + emcs_BkPlast gap = esec_cell - radiator - mgt prin2 I_section,section (' ESEC:I_section,section',i3,F12.4) Do is = 1,esec_Nlayer * define actual cell thickness: Cell = esec_cell plate = radiator * if (is=nint(esec_Nlayer) & (I_section = 3 | I_section = 5)) then Cell = mgt + gap Plate=0 else if (I_section = 4 & is = 1) then ! radiator only Cell = radiator endif * prin2 I_section,is,len,cell,current (' ESEC:I_section,is,len,cell,current ',2i3,3F12.4) if (I_section = 4 & is = 1) then ! radiator only cell = radiator + .14 Create and Position ERAD z=len + (cell)/2 len = len + cell current = current + cell else cell = mgt if(filled = 1) then Create and Position EMGT z=len +(gap+cell)/2 xx = current + (gap+cell)/2 prin2 I_section,is,xx (' MEGA I_section,is ',2i3,F10.4) endif len = len + cell + gap current = current + cell + gap if (Plate>0) then cell = radiator Create and Position ERAD z=len + cell/2 len = len + cell current = current + cell end if end if end do Endblock * ---------------------------------------------------------------------------- Block EMGT is a megatile EM section Attribute EMGT seen=1 colo=1 Material Air Medium standard * Shape CONS dz=mgt/2, rmn1=(current-diff)*Tan_Low-dd, rmn2=(current+mgt-diff)*Tan_Low-dd, rmx1=(current-diff)*Tan_Upp+dup, rmx2=(current+mgt-diff)*Tan_Upp+dup if (I_section=1 | I_section=2 | I_section=5) then Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) else Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) end if * Do isec=1,nint(emcs_Nslices) Create and Position EPER AlphaZ=(emcs_Nslices/2-isec+0.5)*dphi End Do Endblock *--------------------------------------------------------------------------- Block EPER is a EM subsection period (super layer) * Material POLYSTYREN Attribute EPER seen=1 colo=1 Shape CONS dz=mgt/2, phi1=emcs_PhiMin/emcs_Nsector, phi2=+emcs_PhiMax/emcs_Nsector, rmn1=(current-diff)*Tan_Low-dd, rmn2=(current+mgt-diff)*Tan_Low-dd, rmx1=(current-diff)*Tan_Upp+dup, rmx2=(current+mgt-diff)*Tan_Upp+dup * curcl = current+mgt/2 Do ie = 1,nint(eetr_NEta) EtaBot = eetr_EtaBin(ie) EtaTop = eetr_EtaBin(ie+1) RBot=(curcl-diff)*Tanf(EtaBot) * if(Plate > 0) then ! Ordinary Sci layer RTop=min((curcl-diff)*Tanf(EtaTop), _ ((current-diff)*Tan_Upp+dup)) else ! last Sci layer in section RTop=min((curcl-diff)*Tanf(EtaTop), _ ((current-diff)*Tan_Upp+dup)) endif check RBot<RTop * xx=tan(pi*emcs_PhiMax/180.0/emcs_Nsector) yy=cos(pi*emcs_PhiMax/180.0/emcs_Nsector) Create and Position ETAR x=(RBot+RTop)/2 ORT=YZX prin2 ie,EtaTop,EtaBot,rbot,rtop (' EPER : ie,EtaTop,EtaBot,rbot,rtop ',i3,4F12.4) enddo * EndBlock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ETAR is one CELL of scintillator, fiber and plastic * Attribute ETAR seen=1 colo=4 ! blue * local z goes along the radius, y is the thickness Shape TRD1 dy=mgt/2 dz=(RTop-RBot)/2, dx1=RBot*xx-emcs_GapCel/yy, dx2=RTop*xx-emcs_GapCel/yy * Create and Position EALP y=(-mgt+emcs_AlinCell)/2 G10 = esec_scint Create and Position ESCI y=(-mgt+G10)/2+emcs_AlinCell _ +emcs_FrPlast EndBlock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ESCI is the active scintillator (polystyren) layer * Material POLYSTYREN Material Cpolystyren Isvol=1 Attribute ESCI seen=1 colo=7 fill=0 ! lightblue * local z goes along the radius, y is the thickness Shape TRD1 dy=esec_scint/2, dz=(RTop-RBot)/2-emcs_GapCel Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) Call GSTPAR (ag_imed,'CUTNEU',0.001) Call GSTPAR (ag_imed,'CUTHAD',0.001) Call GSTPAR (ag_imed,'CUTMUO',0.001) * define Birks law parameters Call GSTPAR (ag_imed,'BIRK1',1.) Call GSTPAR (ag_imed,'BIRK2',0.013) Call GSTPAR (ag_imed,'BIRK3',9.6E-6) * HITS ESCI Birk:0:(0,10) * xx:16:H(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: endblock * ---------------------------------------------------------------------------- Block ERAD is radiator * Material Iron Attribute ERAD seen=1 colo=6 fill=1 ! violet Shape CONS dz=radiator/2, rmn1=(current)*Tan_Low-dd, rmn2=(current+cell)*Tan_Low-dd, rmx1=(current)*Tan_Upp+dup, rmx2=(current+radiator)*Tan_Upp+dup Create and Position ELED endblock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ELED is lead absorber Plate * Material Lead Attribute ELED seen=1 colo=4 fill=1 Shape TUBS dz=emcs_Pbplate/2, rmin=(current)*Tan_Low, rmax=(current+emcs_Pbplate)*Tan_Upp, Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) Call GSTPAR (ag_imed,'CUTNEU',0.001) Call GSTPAR (ag_imed,'CUTHAD',0.001) Call GSTPAR (ag_imed,'CUTMUO',0.001) endblock * ---------------------------------------------------------------------------- Block EFLP is First Aluminium plate * Material Aluminium Attribute EFLP seen=1 colo=3 fill=1 ! green shape CONS dz=emcs_Front/2, Rmn1=68.813 Rmn2=68.813, Rmx1=(zslice-diff)*Tan_Upp+dup, Rmx2=(zslice + slcwid-diff)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax endblock * ---------------------------------------------------------------------------- Block EALP is ALuminium Plate in calorimeter cell * Material Aluminium Material StrAluminium isvol=0 Attribute EALP seen=1 colo=1 Shape TRD1 dy=emcs_AlinCell/2 dz=(RTop-RBot)/2 Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) Call GSTPAR (ag_imed,'LOSS',1.) Call GSTPAR (ag_imed,'STRA',1.) endblock * ---------------------------------------------------------------------------- Block ESPL is one of the Shower max PLanes * Material Air Attribute ESPL seen=1 colo=3 ! blue Shape TUBS dz=SecWid/2, rmin=section*Tan_Low-1.526, rmax=(section-secwid/2)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax USE EMXG Version=1 msecwd = (emxg_Sapex+emxg_F4)/2 do isec=1,6 cut=1 d3 = 75 - (isec-1)*30 if (exse_sectype(isec) = 0 | (emcg_FillMode=1 & (isec=6 | isec=1))) then cut = 0 Create and position EXSG AlphaZ=d3 Ncopy=isec else if(exse_sectype(isec) = 1) then ! V Create and position EXSG AlphaZ=d3 Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 else if(exse_sectype(isec) = 2) then ! U Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 else if(exse_sectype(isec) = 3) then ! cut V cut=2 Create and position EXSG AlphaZ=d3 Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 else if(exse_sectype(isec) = 4) then ! cut U cut=2 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 endif enddo Endblock * ---------------------------------------------------------------------------- Block EXSG is the Shower max Gap for scintillator strips * Attribute EXSG seen=1 colo=7 serial=cut ! black Material Air Shape TUBS dz=SecWid/2, rmin=section*Tan_Low-1.526, rmax=(section-secwid/2)*Tan_Upp+dup, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec * Rbot = emxg_Rin Rtop = emxg_Rout if(cut > 0) then if(cut = 1) then Rdel = 3.938 Nstr = 288 else Rdel = -.475 Nstr = 285 endif rth = .53*rdel ! .53 --- tentatavily ddn = sq3*1.713 + Rdel ddup = .5*1.846 + 1.713 prin2 Rbot,Rtop,Nstr (' EXSG: Rbot,Rtop,Nstr',2F12.4,I5) mgt = emxg_Sbase + .01 do i_str = 1,nstr p = .5*(i_str-1)*mgt + 41.3655 * if (p <= (.5*rbot*sq3 + rth)) then dxy = 1.9375*sq2 xleft = .5*sq2*p*(sq3 + 1.) - dxy yleft = .5*sq2*p*(sq3 - 1.) - dxy yright = .5*sq2*(sqrt( rbot*rbot - p*p) - p) xright = sq2*p + yright else if ((.5*rbot*sq3 + rth) < p <= (.5*Rtop + 1.5)) then prin2 i_str,p (' EXSG: 2 - -i_str,p:',i3,F12.4) dxy = 1.9375*sq2 xleft = .5*sq2*p*(sq3 + 1.) - dxy yleft = .5*sq2*p*(sq3 - 1.) - dxy dxy = rdel*sq2/sq3 yright = .5*sq2*p*(1.- 1./sq3) xright = sq2*p - yright - dxy yright = -yright - dxy else if (p > (.5*rtop +1.5)) then prin2 i_str,p (' EXSG: 3 - - i_str,p:',i3,F12.4) yleft = (sqrt(rtop*rtop - p*p) - p)/sq2 xleft = sq2*p + yleft dxy = rdel*sq2/sq3 yright = .5*sq2*p*(1.- 1./sq3) xright = sq2*p - yright - dxy yright = -yright - dxy dxy = 0. if ((.5*sq3*160.- ddn) < p <= (.5*sq3*160.+ ddup) ) then prin2 i_str,p (' EXSG: 4 - - i_str,p:',i3,F12.4) xc = .5*(sq3*160.+1.846) yc = xc - .5*sq3*1.713 if (p > yc) then dxy = .5*sq2*(2/sq3*rdel + .5*sq3*1.846 +_ sqrt(1.713*1.713 - (p-xc)*(p-xc))) else dxy = sq2/sq3*(p - .5*sq3* 160. + ddn) endif else if ((.5*sq3*195.- ddn) < p <= (.5*sq3*195. + ddup) ) then prin2 i_str,p (' EXSG: 5 - - i_str,p:',i3,F12.4) xc = .5*(sq3*195.+1.846) yc = xc - .5*sq3*1.713 if (p > yc) then dxy = .5*sq2*(2/sq3*rdel + .5*sq3*1.846 +_ sqrt(1.713*1.713 - (p-xc)*(p-xc))) else dxy = sq2/sq3*(p - .5*sq3*195. + ddn) endif endif xright = xright + dxy yright = yright + dxy endif dxy = section*Tan_Upp - Rtop xc = .5*(xright+xleft) + dxy yc = .5*(yright+yleft) xx = .5*sq2*(xleft+yleft) yy = .5*sq2*(xright+yright) len = xx-yy prin2 i_str,p,yy,xx,len,xc,yc (' EXSG: i_str,x,y1,y2,len,xc,yc:',i3,6F12.4) * Create EHMS if (mod(i_str,2) != 0 ) then Position EHMS x=xc y=yc AlphaZ=-45 else Position EHMS x=xc y=yc AlphaZ=-45 ORT=X-Y-Z endif end do endif * dcut exsg z 0 0 10 0.1 0.1 * dcut exsg y 0 10 -50 0.7 0.7 endblock * ---------------------------------------------------------------------------- Block EHMS is sHower Max Strip * Material POLYSTYREN Material Cpolystyren Isvol=1 Attribute EHMS seen=1 colo=2 serial=cut ! red Shape TRD1 dx1=0 dx2=emxg_Sbase/2 dy=len/2 dz=emxg_Sapex/2 Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) * define Birks law parameters Call GSTPAR (ag_imed,'BIRK1',1.) Call GSTPAR (ag_imed,'BIRK2',0.0130) Call GSTPAR (ag_imed,'BIRK3',9.6E-6) * HITS EHMS Birk:0:(0,10) * xx:16:SH(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: Eloss:0:(0,10) * Endblock * ---------------------------------------------------------------------------- Block EXGT is the G10 layer in the Shower Max * * G10 is about 60% SiO2 and 40% epoxy Component Si A=28.08 Z=14 W=0.6*1*28./60. Component O A=16 Z=8 W=0.6*2*16./60. Component C A=12 Z=6 W=0.4*8*12./174. Component H A=1 Z=1 W=0.4*14*1./174. Component O A=16 Z=8 W=0.4*4*16./174. Mixture g10 Dens=1.7 Attribute EXGT seen=1 colo=7 Shape TUBS dz=emxg_F4/2, rmin=(section-diff)*Tan_Low-1.526, rmax=(section+msecwd-diff)*Tan_Upp, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) EndBlock * ---------------------------------------------------------------------------- * ECAL nice views: dcut ecvo x 1 10 -5 .5 .1 * draw emdi 105 0 160 2 13 .2 .1 * draw emdi 120 180 150 1 14 .12 .12 * --------------------------------------------------------------------------- end
ecalgeo.g geometry file (Jason edits, g23)
ecalgeo.g geometry file (Jason Webb edits, g23)
c***************************************************************************** Module ECALGEO is the EM EndCap Calorimeter GEOmetry c-- Created 26 jan 1996 Author Rashid Mehdiyev c-- c Version 1.1, W.J. Llope c - changed sensitive medium names... c c Version 2.0, R.R. Mehdiyev 16.04.97 c - Support walls included c - intercell and intermodule gaps width updated c - G10 layers inserted c Version 2.1, R.R. Mehdiyev 23.04.97 c - Shower Max Detector geometry added c - Variable eta grid step size introduced c Version 2.2, R.R. Mehdiyev 03.12.97 c - Eta grid corrected c - Several changes in volumes dimensions c - Material changes in SMD c c Version 3.0, O. Rogachevsky 28.11.99 c - New proposal for calorimeter SN 0401 c c Version 4.1, O.Akio 3 Jan 01 c - Include forward pion detectors c c Version 5.0, O. Rogachevsky 20.11.01 c - FPD is eliminated in this version c - More closed to proposal description c of calorimeter and SMD structure c c***************************************************************************** +CDE,AGECOM,GCONST,GCUNIT. * Content EAGA,EALP,ECAL,ECHC,ECVO,ECGH,EFLP,EHMS, ELED,EMGT,EMOD,EPER,EPSB,ERAD,ERCM,ERSM, ESHM,ESEC,ESCI,ESGH,ESPL,ESSP,EMSS,ETAR, EXGT,EXSG,EXPS,EFLS,EBLS Structure EMCG { Version, int Onoff, int fillMode} Structure EMCS { Version,Type,zorg,zend,EtaMin,EtaMax, PhiMin,PhiMax,Offset, Nsupsec,Nsector,Nsection,Nslices, Front,AlinCell,Frplast,Bkplast,PbPlate,LamPlate, BckPlate,Hub,Rmshift,SMShift,GapPlt,GapCel, GapSMD,SMDcentr,TieRod(2),Bckfrnt,GapHalf,Cover, Rtie,slop} Structure EETR { Type,Etagr,Phigr,Neta,EtaBin(13)} Structure ESEC { Isect, FPlmat, Cell, Scint, Nlayer, deltaz, Jiggle(18) } Structure EMXG {Version,Sapex,Sbase,Rin,Rout,F4} Structure EXSE {Jsect,Zshift,Sectype(6)} Structure ESMD {Version, front_layer, back_layer, spacer_layer, base, apex } Integer I_section,J_section,Ie,is,isec,istrip,Nstr,Type,ii,jj, cut,fsect,lsect,ihalf,filled,i,j,k,i_sector Real center,Plate,Cell,G10,halfi, tan_low,tan_upp,Tanf,RBot,Rtop,Deta,etax,sq2,sq3, dup,dd,d2,d3,rshift,dphi,radiator Real maxcnt,msecwd,mxgten,curr,Secwid,Section, curcl,EtaTop,EtaBot,zwidth,zslice,Gap,megatile, xleft,xright,yleft,yright,current, rth,length,p,xc,yc,xx,yy,rdel,dxy,ddn,ddup Real myPhi Integer N Parameter (N=12) Tanf(etax) = tan(2*atan(exp(-etax))) c-------------------------------------------------------------------------------- c Data c c FillMode =1 only 2-5 sectors (in the first half) filled with scintillators c FillMode =2 all sectors filled (still only one half of one side) c FillMode =3 both halves (ie all 12 sectors are filled) c c OnOff =0 Do not build geometry c OnOff =1 Build West Endcap c OnOff =2 Build East Endcap (disabled) c OnOff =3 Build Both Endcaps (east disabled) c c Note: Fill EMCG ! EM EndCAp Calorimeter basic data Version = 5.0 ! Geometry version OnOff = 3 ! Configurations 0-no, 1-west 2-east 3-both FillMode = 3 ! sectors fill mode c-- Fill EMCS ! EM Endcap Calorimeter geometry Version = 1 ! Versioning Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrg = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.400 ! space for SMD detector << version 2 -- 3.600 >> SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half Rtie = 1.0425 ! Radius of tie rod Slop = 0.1400 ! Added to cell containing radiator 6 (formerly hardcoded in geom) c-- Fill EMCS ! EM Endcap Calorimeter geometry Version = 2 ! Versioning Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrg = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.600 ! space for SMD detector (* from master_geom_bmp.xls *) SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half Rtie = 0.75 ! Radius of tie rod Slop = 0.0000 ! Added to cell containing radiator 6 (formerly hardcoded in geom) c-- c--------------------------------------------------------------------------- c-- c-- Supporting documentation: c-- http://drupal.star.bnl.gov/STAR/system/files/SMD_module_stack.pdf c-- Fill ESMD ! shower maximum detector information Version = 1 ! versioning information front_layer = 0.161 ! thickness of front layer back_layer = 0.210 ! thickness of back layer base = 1.0 ! base of the SMD strip apex = 0.7 ! apex of the SMD strip spacer_layer = 1.2 ! spacer layer c-- Fill EETR ! Eta and Phi grid values Type = 1 ! =1 endcap, =2 fpd EtaGr = 1.0536 ! eta_top/eta_bot tower granularity PhiGr = 0.0981747 ! Phi granularity (radians) NEta = 12 ! Eta granularity EtaBin = {2.0,1.9008,1.8065,1.7168,1.6317,1.5507,1.4738, 1.4007,1.3312,1.2651,1.2023,1.1427,1.086}! Eta rapidities c-- c--------------------------------------------------------------------------- c-- Fill ESEC ! Preshower 1 / Radiator 1 ISect = 1 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.475 ! Sci layer thickness (4.75mm Bicron) deltaz = -0.014 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- c-- Note: Jiggle allows one to shift each megatile by Jiggle(i) degrees, where c-- i indicates the layer within the section of the calorimeter. This feature c-- has only been crudely tested... i.e. it compiles and creates a reasonable c-- set of pictures, but I have not verified that every scintillator shows up... c-- There could be volume conflicts and this would need to be checked. --JW c-- Fill ESEC ! Preshower 2 / Radiator 2 ISect = 2 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.475 ! Sci layer thickness (4.75mm Bicron) deltaz = -0.0182 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Megatiles 3-6 / Radiators 3-5 ISect = 3 ! Section number Nlayer = 4 ! Number of Sci layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness deltaz = -0.0145 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Megatiles 7-23 / Radiators 6-23 ISect = 4 ! Section Nlayer = 18 ! Number of layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness deltaz = +0.0336 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Postshower ISect = 5 ! Section Nlayer = 1 ! Number of layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness (5.0mm Kurarary) deltaz = +0.036 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- c---------------------------------------------------------------------------- c-- Fill EMXG ! EM Endcap SMD basic data Version = 1 ! Geometry version Sapex = 0.7 ! Scintillator strip apex Sbase = 1.0 ! Scintillator strip base Rin = 77.41 ! inner radius of SMD plane Rout = 213.922 ! outer radius of SMD plane F4 = .15 ! F4 thickness c-- c---------------------------------------------------------------------------- c-- Fill EXSE ! First SMD section JSect = 1 ! Section number Zshift = -1.215 ! Section width sectype = {4,1,0,2,1,0} ! 1-V,2-U,3-cutV,4-cutU c-- Fill EXSE ! Second SMD section JSect = 2 ! Section number Zshift = 0. ! Section width sectype = {0,2,1,0,2,3} ! 1-V,2-U,3-cutV,4-cutU c-- Fill EXSE ! Third SMD section JSect = 3 ! Section number Zshift = 1.215 ! Section width sectype = {1,0,2,1,0,2} ! 1-V,2-U,3-cutV,4-cutU c-- c---------------------------------------------------------------------------- c-- Materials c-- c-- PVC used in the SMD spacer layers c-- Component H A=1 Z=1 W=3.0*1.0/62.453 Component C A=12 Z=6 W=2.0*12.0/62.453 Component Cl A=35.453 Z=17 W=1.0*35.453/62.453 Mixture PVC_Spacer Dens=1.390*(1.20/1.00) c-- c-- Lead alloy used in the radiators c-- Component Sn A=118.710 Z=50 W=0.014 Component Ca A=40.0780 Z=20 W=0.00075 Component Al A=26.9815 Z=13 W=0.0003 Component Pb A=207.190 Z=82 W=0.98495 Mixture PbAlloy DENS=11.35 c-- c-- Stainless Steel used in various places c-- Component Cr A=51.9960 Z=24 W=0.19 Component Ni A=58.6934 Z=28 W=0.09 Component Fe A=55.8450 Z=26 W=0.72 Mixture Steel DENS=8.03 c-- c-- Aluminized mylar. According to information which I dug up on a google c-- search, this is typically mylar coated with a thin (1000 angstrom) layer c-- of aluminium on each side. c-- c-- http://www.eljentechnology.com/datasheets/EJ590-B10HH%20data%20sheet.pdf c-- Component Mylar A=12.875 Z=6.4580 w=0.999 Component Al A=26.980 Z=13.000 w=0.001 Mixture AlMylar dens=1.390 c-- c-- G10 Epoxy used in various places c-- Component Si A=28.08 Z=14 W=0.6*1*28./60. Component O A=16 Z=8 W=0.6*2*16./60. Component C A=12 Z=6 W=0.4*8*12./174. Component H A=1 Z=1 W=0.4*14*1./174. Component O A=16 Z=8 W=0.4*4*16./174. Mixture G10 Dens=1.7 c-- c-- Fibreglass cloth used in SMD stackup. I googled this one too... a self- c-- described expert quotes typical densities and percent by volume c-- http://en.allexperts.com/q/Composite-Materials-2430/fiberglass-1.htm c-- c-- glass fiber: 2.6 g/cm3 (17.6%) resin: 1.3 g/cm3 (82.4%) c-- c-- Fiberglass density = 1.529 g/cm3 c-- c-- I will assume that G10 epoxy is close enough to the typical resins c-- used, at least in terms of chemical composition. Then c-- Component G10 A=18.017 Z=9.013 W=1.3*0.824/(1.3*0.824+2.6*0.176) Component Si A=28.08 Z=14 W=2.6*0.176/(1.3*0.824+2.6*0.176)*28.08/60.08 Component O A=16 Z=8 W=2.6*0.176/(1.3*0.824+2.6*0.176)*32.00/60.08 Mixture Fiberglass dens=1.53 c-- c-- c---------------------------------------------------------------------------- c-- Select versions of various geometry data c-- Use EMCG Use EMCS Version=2 Use EETR c-- c---------------------------------------------------------------------------- c-- Calculate frequently used quantities c-- sq3 = sqrt(3.) ! 1/tan(30deg) = sq3 sq2 = sqrt(2.) c-- c-- center = (emcs_zorg+emcs_zend)/2 ! center of the calorimeter tan_upp = tanf(emcs_etamin) ! think this is angle pointing to top of calo tan_low = tanf(emcs_etamax) ! think this is angle pointing to bot of calo rth = sqrt(1. + tan_low*tan_low) ! ?? rshift = emcs_hub * rth ! ?? dup = emcs_rmshift*tan_upp ! dd = emcs_rmshift*rth ! d2 = rshift + dd ! radiator = emcs_pbplate + 2*emcs_lamplate ! thickness of radiator assembly dphi = (emcs_phimax-emcs_phimin)/emcs_nsector ! single endcap sector c-- c---------------------------------------------------------------------------- c---------------------------------------------------------------------------- c-- BEGIN Prin1 emcg_version ('ecalgeo version: ',F4.2) c-- IF (emcg_OnOff>0) THEN c-- c-- Build the EEMC geometry for one half wheel c-- Create ECAL c-- c-- Position the two halves. Bottom half installed in 2003, top c-- half in 2004... so we allow logic to allow for the time c-- evolution of the calorimeter c-- c-- c-- West Endcap c-- IF (emcg_OnOff==1 | emcg_OnOff==3) THEN Position ECAL in CAVE z=+center ENDIF IF (section > emcs_zend) THEN Prin1 section, emcs_zend (' ECALGEO error: sum of sections exceeds maximum ',2F12.4) ENDIF IF (emcg_OnOff==2 ) THEN Prin1 ('East Endcap has been removed from the geometry' ) ENDIF c-- EndIF! emcg_OnOff c-- Prin1 ('ECALGEO finished') c-- c-- END c---------------------------------------------------------------------------- c----------------------------------------------------------------- Block ECAL -- c-- Block ECAL is one EMC EndCap wheel c-- c-- The EEMC is built from two 180 degree half-wheels tilted at an angle c-- with respect to zero in the STAR reference frame. This block is serves c-- as a logical volume which creates the two half wheels. c-- c-- Creates: c-- + EAGA c-- Material Air Attribute ECAL seen=0 colo=7 ! lightblue c-- Shape CONE dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2, rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup, rmx2=emcs_zend*tan_upp+dup c-- c-- DO ihalf=1,2 c-- filled = 1 halfi = -105 + (ihalf-1)*180 if (ihalf=2 & emcg_FillMode<3) filled = 0 c-- Create and Position EAGA AlphaZ=halfi c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block EAGA -- c-- Block EAGA IS HALF OF WHEEL AIR VOLUME FOR THE ENDCAP MODULE c-- c-- The eemc is divided into two halves. one half installed for 2003 run, c-- second half added for 2004 and beyond. the eaga block represents one c-- of these half-wheels. it is an air volume which will be filled in c-- with additional detector components. c-- c-- Creates: c-- + EMSS -- steel support block c-- + ECGH -- air gap between the two halves c-- C-- Material AIR Attribute EAGA seen=0 colo=1 serial=FILLED ! BLACK C-- Shape CONS dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2 rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup rmx2=emcs_zend*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- c-- IF ( FILLED .EQ. 1 ) THEN c-- Create AND Position EMSS konly='MANY' c-- curr = emcs_zorg curcl = emcs_zend c-- Create AND Position ECGH alphaz=90 kOnly='ONLY' c-- ENDIF c-- EndBlock c----------------------------------------------------------------- Block EMSS -- c-- Block EMSS is the steel support of the endcap module c-- c-- Creates: c-- + EFLP -- ALUMINIUM FRONT PLATE c-- + ECVO -- VOLUMES TO CONTAIN RADIATORS AND MEGATILES c-- + ESHM -- SHOWER MAX DETECTOR VOLUME c-- + ESSP -- STAINLESS STEEL BACKPLATE c-- + ERCM -- STAINLESS STEEL TIE-RODS PENETRATING ECVO c-- c-- Material Steel c-- Attribute EMSS seen=1 colo=1 ! BLACK Shape CONS dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2 rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup rmx2=emcs_zend*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- c-- Aluminium front plate C-- zslice = emcs_zorg zwidth = emcs_front c-- Prin1 zslice+zwidth/2 (' Front Al plate centered at: ', F12.4 ) c-- Create AND Position EFLP z=zslice-center+zwidth/2 zslice = zslice + zwidth C-- Prin1 zslice (' FIRST CALORIMETER STARTS AT: ',F12.4) c-- c-- Preshower 1, preshower 2, and calorimeter tiles up to c-- megatile number six. c-- fsect = 1 ! first section lsect = 3 ! last section c-- zwidth = emcs_smdcentr - emcs_gapsmd/2 - zslice ! width of current slice c-- Prin1 zslice+zwidth/2 ('Sections 1-3 positioned at: ', F12.4 ) c-- Create AND Position ECVO z=zslice-center+zwidth/2 c-- zwidth = emcs_gapsmd zslice = emcs_smdcentr - emcs_gapsmd/2 c-- Prin1 section, zslice (' 1st calorimeter ends, smd starts at: ',2f10.5) Prin1 zwidth (' smd width = ',f10.5 ) c-- Prin1 zslice+zwidth/2 ('SMD section centered at: ', F12.4 ) c-- Do not kill neighbors Create AND Position ESHM z=zslice-center+zwidth/2 kOnly='MANY' zslice = zslice + zwidth c-- Prin1 zslice (' SMD ends at: ',f10.5) c-- c-- fsect = 4 ! first section lsect = 5 ! last section c-- c-- Calculate the width of the last two calorimeter sections c-- zwidth = 0 DO i_section = fsect,lsect c-- USE ESEC isect=i_section zwidth = zwidth + esec_cell*esec_nlayer c-- ENDDO c-- c-- ============================================================= c-- c-- Total width will be between the back plate and the current c-- position... this effectively turns the geometry into an c-- accordian... whatever was defined earlier will compress c-- / expand this section. so correcting the smd gap will c-- result in some small, sub-mm shifts of radiators and c-- megatiles... one would like to actually place these c-- into their absolute positions. c-- c-- ============================================================== c-- zwidth = emcs_bckfrnt - zslice c-- Prin1 zslice+zwidth/2 ('Sections 4-5 positioned at: ', F12.4 ) c-- Create AND Position ECVO z=zslice-center+zwidth/2 c-- zslice = emcs_bckfrnt c-- Prin1 section,zslice (' 2nd calorimeter ends, back plate starts at: ',2f10.5) c-- zwidth = emcs_bckplate c-- Create AND Position ESSP z=zslice-center+zwidth/2 c-- zslice = zslice + zwidth c-- Prin1 zslice ('EEMC Al backplate ends at: ',F12.4 ) c-- c-- Done with the calorimeter stackup. now go back and cut through the c-- calorimeter stack with the tie rods c-- c-- slice width will be full calorimeter depth zwidth = emcs_zend-emcs_zorg c-- Create ERCM c-- DO i = 1,2 ! two tie rods along DO j = 1,5 ! each gap between sectors (5 gaps) xx = emcs_phimin + j*30 yy = xx*degrad xc = cos(yy)*emcs_tierod(i) yc = sin(yy)*emcs_tierod(i) Position ERCM z=0 x=xc y=yc ENDDO ENDDO c-- c-- Now add in projective steel bars which form part of the support c-- structure of the eemc c-- rth = emcs_zorg*tan_upp+dup + 2.5/2 xc = (emcs_zend - emcs_zorg)*tan_upp length = .5*(emcs_zend + emcs_zorg)*tan_upp + dup + 2.5/2 yc = emcs_zend-emcs_zorg p = atan(xc/yc)/degrad c-- Create EPSB DO i = 1,6 c-- xx = -75 + (i-1)*30 yy = xx*degrad xc = cos(yy)*length yc = sin(yy)*length c-- Position EPSB X=XC Y=YC ALPHAZ=XX c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ECVO -- c-- Block ECVO is one of endcap volume with megatiles and radiators c-- c-- CreateS: c-- + EMOD -- Responsible for creating esec which, in a glorious example c-- of spaghetti code, turns around and creates esec, which is c-- responsible for creating the radiators before and after the c-- smd layers. C-- Material AIR Attribute ECVO seen=1 colo=3 ! GREEN Shape CONS dz=zwidth/2, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup c-- c-- Loop over the SIX SECTORS in the current half-wheel. determine c-- whether the sector is filled or not, and create the "module". c-- By "module", we really mean endcap sector. (Lots of code in the c-- EEMC borrows from the barrel, and so barrel modlues get mapped c-- to EEMC sectors). c-- DO i_sector = 1,6 c-- IF (1 < I_SECTOR < 6 | EMCG_FILLMODE > 1) THEN filled = 1 ELSE filled = 0 ENDIF c-- d3 = 75 - (i_sector-1)*30 Create AND Position EMOD alphaz=d3 ncopy=i_sector c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ESHM -- c-- Block ESHM is the shower max section c-- c-- CreateS: c-- + ESPL -- SHOWER MAXIMUM DETECTOR PLANES c-- + ERSM -- TIE RODS W/IN THE SHOWER MAXIMUM DETECTOR c-- Material AIR Attribute ESHM seen=1 colo=4 ! BLUE c-- Shape CONS dz=zwidth/2, rmn1=(zslice*tan_low)-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=(zslice)*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- USE EMXG c-- maxcnt = emcs_smdcentr Prin1 zslice, section, center (' === z start for smd,section: ',3f12.4) c-- c-- Loop over the three possible locations for the smd planes and c-- create them. note that code w/in espl will decide which of c-- 5 types of smd planes are created... u, v, cutu,cutv or spacer. c-- DO j_section = 1,3 c-- USE EXSE jsect=j_section c-- current = exse_zshift secwid = emxg_sapex + 2.*emxg_f4 section = maxcnt + exse_zshift c-- Prin1 j_section,current,section,secwid (' layer, z, width : ',i3,3f12.4) c-- rbot=section*tan_low rtop=section*tan_upp c-- Prin1 j_section,rbot,rtop (' layer, rbot,rtop : ',i3,2f12.4) c-- Prin1 j_section, center+current (' smd layer=',I1,' z=',F12.4 ) c-- Do not kill neighbors Create and Position ESPL z=current kOnly='MANY' c-- ENDDO c-- c-- Add in the tie rods which penetrate the SMD layers c-- Create ERSM c-- DO i = 1,2 DO j = 1,5 xx = emcs_phimin + j*30 yy = xx*degrad xc = cos(yy)*emcs_tierod(i) yc = sin(yy)*emcs_tierod(i) Position ERSM Z=0 X=XC Y=YC END DO END DO C-- EndBlock c----------------------------------------------------------------- Block ECGH -- c-- Block ECGH is air gap between endcap half wheels c-- c-- Creates: c-- + ECHC -- THE STAINLESS STEEL COVER FOR 1/2 OF THE EEMC. c-- Material AIR Medium standard Attribute ECGH seen=0 colo=7 ! LIGHTBLUE Shape TRD1 dz=(emcs_zend-emcs_zorg)/2, dy =(emcs_gaphalf+emcs_cover)/2, dx1=emcs_zorg*tan_upp+dup, dx2=emcs_zend*tan_upp+dup c-- c-- rth = emcs_gaphalf + emcs_cover xx=curr*tan_low-d2 xleft = sqrt(xx*xx - rth*rth) yy=curr*tan_upp+dup xright = sqrt(yy*yy - rth*rth) secwid = yy - xx xx=curcl*tan_low-d2 yleft = sqrt(xx*xx - rth*rth) yy=curcl*tan_upp+dup yright = sqrt(yy*yy - rth*rth) zwidth = yy - xx xx=(xleft+xright)/2 yy=(yleft + yright)/2 xc = yy - xx length = (xx+yy)/2 yc = curcl - curr p = atan(xc/yc)/degrad rth = -(emcs_gaphalf + emcs_cover)/2 c-- Create ECHC c-- Position ECHC X=+LENGTH Y=RTH Position ECHC X=-LENGTH Y=RTH ALPHAZ=180 c-- EndBlock c----------------------------------------------------------------- Block ECHC -- c-- Block ECHC is steel endcap half cover c-- Material steel Attribute ECHC seen=1 colo=1 ! BLACK c-- Shape TRAP dz=(curcl-curr)/2, thet=p, bl1=secwid/2, tl1=secwid/2, bl2=zwidth/2, tl2=zwidth/2, h1=emcs_cover/2, h2=emcs_cover/2, phi=0, alp1=0, alp2=0 c-- EndBlock c----------------------------------------------------------------- Block ESSP -- c-- Block ESSP is stainless steel back plate c-- Material steel Attribute ESSP seen=1 colo=6 fill=1 Shape CONS dz=emcs_bckplate/2, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin, phi2=emcs_phimax c-- EndBlock c----------------------------------------------------------------- Block EPSB -- c-- Block EPSB IS A PROJECTILE STAINLESS STEEL BAR C-- Material Steel Attribute EPSB seen=1 colo=6 FILL=1 Shape TRAP dz=(emcs_zend-emcs_zorg)/2, thet=p, bl1=2.5/2, tl1=2.5/2, bl2=2.5/2, tl2=2.5/2, h1=2.0/2, h2=2.0/2, phi=0, alp1=0, alp2=0 c-- c-- EndBlock c----------------------------------------------------------------- Block ERCM -- c-- Block ERCM is stainless steel tie rod in calorimeter sections c-- Material Steel Attribute ERSM seen=1 colo=6 FILL=1 c-- Shape TUBE dz=zwidth/2, rmin=0, rmax=emcs_rtie c-- c-- Looks like the tie rods are meant to engage the 1.525 cm diameter holes c-- piercing the ears of the smd spacer... 1.5 cm may be a better approximation c-- here. c-- c-- http://drupal.star.bnl.gov/star/system/files/smd_spacer_drawings.pdf c-- EndBlock c----------------------------------------------------------------- Block ERSM -- c-- Block ERSM is stainless steel tie rod in shower max c-- Material Steel Attribute ERSM seen=1 colo=6 FILL=1 c-- Shape TUBE dz=zwidth/2, rmin=0, rmax=emcs_rtie c-- c-- see comments above c-- EndBlock c----------------------------------------------------------------- Block EMOD -- c-- Block EMOD (fsect,lsect) IS ONE MODULE OF THE EM ENDCAP c-- c-- Arguements: (do be defined prior to the creation of this block) c-- c-- fsect -- first section to create c-- lsect -- last section to create c-- Attribute EMOD seen=1 colo=3 serial=FILLED ! GREEN Material Air Shape CONS dz=zwidth/2, phi1=emcs_phimin/emcs_nsupsec, phi2=emcs_phimax/emcs_nsupsec, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup c-- c-- Running parameter 'section' contains the position of the current section c-- it should not be modified in daughters, use 'current' variable instead. c-- secwid is used in all 'cons' daughters to define dimensions. c-- section = zslice curr = zslice + zwidth/2 c-- c-- DO i_section = fsect, lsect USE ESEC isect=i_section c-- secwid = esec_cell*esec_nlayer c-- c-- Section 3 precedes the smd. section 5 is the post shower. in c-- both cases these sections end with a scintillator layer and no c-- radiator. c-- IF (I_SECTION = 3 | I_SECTION = 5) THEN secwid = secwid - radiator ELSE IF (I_SECTION = 4) THEN ! add one more radiator secwid = secwid - esec_cell + radiator ENDIF c-- Prin1 i_section, section-curr+secwid/2 ('+ ECVO isection=',I1,' zcenter=', F12.4) c-- Create AND Position ESEC z=section-curr+secwid/2 c-- section = section + secwid c-- ENDDO! Loop over sections c-- EndBlock c----------------------------------------------------------------- Block ESEC -- c-- Block ESEC is a single em section Material AIR Medium standard Attribute ESEC seen=1 colo=1 serial=filled lsty=2 c-- Shape CONS dz=secwid/2, rmn1=(section)*tan_low-dd, rmn2=(section+secwid)*tan_low-dd, rmx1=(section)*tan_upp+dup, rmx2=(section+secwid)*tan_upp+dup c-- length = -secwid/2 current = section c-- megatile = esec_scint+emcs_alincell+emcs_frplast+emcs_bkplast c-- gap = esec_cell - radiator - megatile Prin2 i_section,section (' ESEC:i_section,section',i3,f12.4) c-- c-- Loop over all layers in this section c-- DO is = 1,esec_nlayer c-- c-- Define actual cell thickness: cell = esec_cell plate = radiator c-- IF (is=nint(esec_nlayer) & (i_section = 3 | i_section = 5)) THEN c-- cell = megatile + gap plate=0 c-- ELSE IF (i_section = 4 & is = 1) THEN ! RADIATOR ONLY c-- cell = radiator c-- ENDIF c-- Prin2 i_section,is,length,cell,current (' esec:i_section,is,length,cell,current ',2i3,3f12.4) C-- C-- This handles the special case in the section after the smd. c-- this section begins with a lead radiator. the previous section c-- ended with a plastic scintillator c-- IF (i_section = 4 & is = 1) THEN ! radiator only c-- c$$$ cell = radiator + .14 cell = radiator + emcs_slop ! ^^^^ probably the fiber router layer... but is this needed here? c-- Prin1 is, current + cell/2+esec_deltaz ( ' + ESEC radiator ilayer=',I2,' z=',F12.4 ) Create AND Position ERAD z=length+(cell)/2+esec_deltaz c-- length = length + cell current = current + cell c-- c-- All other cases are standard radiator followed by scintillator c-- ELSE c-- cell = megatile IF (FILLED = 1) THEN c-- Create AND Position EMGT z=length+(gap+cell)/2+esec_deltaz c-- xx = current + (gap+cell)/2+esec_deltaz prin2 i_section,is,xx (' mega i_section,is ',2i3,f10.4) Prin1 is, xx (' + ESEC megatile ilayer=',I2,' z=',F12.4) c-- ENDIF c-- length = length + cell + gap current = current + cell + gap c-- IF (PLATE>0) THEN c-- cell = radiator Prin1 is, current + cell/2+esec_deltaz ( ' + ESEC radiator ilayer=',I2,' z=',F12.4 ) Create AND Position ERAD z=length+cell/2+esec_deltaz c-- length = length + cell current = current + cell c-- ENDIF c-- ENDIF c-- ENDDO c-- c-- EndBlock c----------------------------------------------------------------- Block EMGT -- c-- Block EMGT is a 30 degree megatile c-- Material Air Medium Standard Attribute EMGT seen=1 colo=1 lsty=2 c-- Shape CONS dz=megatile/2, rmn1=(current)*tan_low-dd, rmn2=(current+megatile)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+megatile)*tan_upp+dup c-- c-- DO isec=1,nint(emcs_nslices) c-- myPhi = (emcs_nslices/2-isec+0.5)*dphi + esec_jiggle(is) c-- Create AND Position EPER alphaz=myPhi c-- END DO c-- EndBlock c----------------------------------------------------------------- Block EPER -- c-- Block EPER is a 5 degree slice of a 30 degree megatile (subsector) c-- c-- Creates: c-- + ETAR -- The pseudo-rapidity divivisions in the megatiles c-- Material Polystyren Attribute EPER seen=1 colo=1 lsty=1 c-- c-- c-- Shape CONS dz=megatile/2, phi1=emcs_phimin/emcs_nsector, phi2=emcs_phimax/emcs_nsector, rmn1=(current)*tan_low-dd, rmn2=(current+megatile)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+megatile)*tan_upp+dup c-- curcl = current+megatile/2 DO ie = 1, nint(eetr_neta) c-- etabot = eetr_etabin(ie) etatop = eetr_etabin(ie+1) rbot=(curcl)*tanf(etabot) rtop=min((curcl)*tanf(etatop), ((current)*tan_upp+dup)) c-- check rbot<rtop c-- xx=tan(pi*emcs_phimax/180.0/emcs_nsector) yy=cos(pi*emcs_phimax/180.0/emcs_nsector) Create and Position ETAR x=(rbot+rtop)/2 ort=yzx prin2 ie,etatop,etabot,rbot,rtop (' EPER : ie,etatop,etabot,rbot,rtop ',i3,4f12.4) c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ETAR -- c-- c-- ETAR is a single cell of scintillator, including fiber router, plastic, c-- etc... c-- c-- local z is radially outward in star c-- local y is the thickness of the layer c-- Block ETAR is a single calorimeter cell, containing scintillator, fiber router, etc... c-- Material POLYSTYREN Attribute ETAR seen=1 colo=4 lsty=1 ! BLUE c-- Shape TRD1 dy=megatile/2 dz=(rtop-rbot)/2, dx1=rbot*xx-emcs_gapcel/yy, dx2=rtop*xx-emcs_gapcel/yy c-- Create AND Position EALP y=(-megatile+emcs_alincell)/2 g10 = esec_scint Create AND Position ESCI y=(-megatile+g10)/2+emcs_alincell _ +emcs_frplast c-- EndBlock c----------------------------------------------------------------- Block ESCI -- c-- Block ESCI is the active scintillator (polystyrene) layer c-- c-- Obtain the definition of polystyrene on this line, next line clones Material Polystyren Material Ecal_scint isvol=1 Medium Ecal_active isvol=1 c-- Attribute ESCI seen=1 colo=7 fill=0 lsty=1 ! LIGHTBLUE c-- local z goes along the radius, y is the thickness Shape TRD1 dy=esec_scint/2, dz=(rtop-rbot)/2-emcs_gapcel c-- c-- Call ecal_set_cuts( ag_imed, 'detector' ) c-- c-- HITS ESCI BIRK:0:(0,10) c-- c-- EndBlock c----------------------------------------------------------------- Block ERAD -- c-- Block ERAD is the lead radiator with stainless steel cladding c-- c-- Creates: c-- + ELED -- the business end of the calorimeter... c-- Material STEEL c-- Attribute ERAD seen=1 colo=6 fill=1 lsty=1 ! VIOLET Shape CONS dz=radiator/2, rmn1=(current)*tan_low-dd, rmn2=(current+cell)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+radiator)*tan_upp+dup c-- Create AND Position ELED c-- EndBlock c------------------------------------------------------------------------- c----------------------------------------------------------------- Block ELED -- c-- Block ELED is a lead absorber plate c-- c-- Material PbAlloy Medium Ecal_lead Attribute ELED seen=1 colo=4 fill=1 lsty=1 c-- Shape TUBS dz=emcs_pbplate/2, rmin=(current)*tan_low, rmax=(current+emcs_pbplate)*tan_upp, c-- Call ecal_set_cuts( ag_imed, 'radiator' ) c-- EndBlock c-- c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EFLP -- c-- Block EFLP is the aluminum (aluminium) front plate of the endcap c-- Material ALUMINIUM Attribute EFLP seen=1 colo=3 fill=1 lsty=1 ! GREEN Shape CONS dz=emcs_front/2, rmn1=68.813 rmn2=68.813, rmx1=(zslice)*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- EndBlock c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EALP -- c-- Block EALP is the thin aluminium plate in calorimeter cell c-- c-- Material Aluminium Attribute EALP seen=1 colo=1 lsty=1 c-- c-- Shape TRD1 dy=emcs_alincell/2 dz=(rtop-rbot)/2 c-- c-- Thin aluminium plate in each calorimeter cell. The energy-loss c-- fluctuations are restricted in this thin material. c-- CALL GsTPar (AG_IMED,'CUTGAM',0.00001) CALL GsTPar (AG_IMED,'CUTELE',0.00001) CALL GsTPar (AG_IMED,'LOSS',1.) CALL GsTPar (AG_IMED,'STRA',1.) c-- EndBlock c----------------------------------------------------------------- Block ESPL -- c-- Block ESPL is the logical volume containing an SMD plane c-- Material Air Attribute ESPL seen=1 colo=4 lsty=4 Shape TUBS dz=emcs_gapsmd/3/2, rmin=section*tan_low-1.526, rmax=(section-secwid/2)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- USE EMXG version=1 msecwd = (emxg_sapex+emxg_f4)/2 c-- ^^^^^^ what is this used for? --jw c-- looks like the g10 layer which we are retiring c-- c-- loop over the six sectors in an endcap half wheel c-- DO isec=1,6 cut=1 d3 = 75 - (isec-1)*30 c-- IF (exse_sectype(isec)=0|(emcg_fillmode=1&(isec=6|isec=1))) THEN cut = 0 c -- come back and build spacers -- ElseIF (exse_sectype(isec) = 1) then ! v c-- Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 2) then ! u c-- Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 3) then ! cut v c-- cut=2 Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 4) then ! cut u c-- cut=2 Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' c-- EndIF c-- EndDO! loop over six sectors in eemc half wheel c-- c-- repeat the loop and add in the spacer layers c-- DO isec=1,6 d3=75 - (isec-1)*30 IF (exse_sectype(isec)=0|(emcg_fillmode=1&(isec=6|isec=1))) then cut = 0 c-- Do not kill neighbors Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ c potential side effect... may screw up the mapping c of the smd strips into the tables? c EndIF EndDO c-- EndBlock c----------------------------------------------------------------- Block EXSG -- c-- Block EXSG Is another logical volume... this one acutally creates the planes c-- c-- Creates: c-- + EHMS -- shower max strips c-- + EFLS -- front cover for SMD planes c-- + EBLS -- back cover for SMD planes c-- Attribute EXSG seen=1 colo=7 serial=cut lsty=3 ! MEH Material Air c$$$ Medium TMED_EXSG stemax=0.01 Shape TUBS dz=emcs_gapsmd/3/2, rmin=section*tan_low-1.526, rmax=(section-secwid/2)*tan_upp+dup, phi1=emcs_phimin/emcs_nsupsec-5, phi2=emcs_phimax/emcs_nsupsec+5 c-- rbot = emxg_rin rtop = emxg_rout c-- c-- Code to handle smd spacers c-- IF ( cut .eq. 0 ) THEN Create and Position EXPS kONLY='MANY' ENDIF c-- c-- Code to handle smd planes c-- IF (cut > 0) THEN c-- c-- setup which plane we are utilizing c-- IF (cut = 1) THEN nstr = 288 ELSE nstr = 285 ENDIF c-- c-- loop over all smd strips and place them w/in this smd plane c-- DO istrip = 1,nstr c-- Call ecal_get_strip( section, cut, istrip, xc, yc, length ) c-- IF (mod(istrip,2) != 0 ) THEN Create and Position EHMS x=xc y=yc alphaz=-45 kOnly='ONLY' Create and Position EBLS x=xc y=yc z=(+esmd_apex/2+esmd_back_layer/2) alphaz=-45 kOnly='ONLY' ELSE Create and Position EHMS x=xc y=yc alphaz=-45 ort=x-y-z kOnly='ONLY' Create and Position EFLS x=xc y=yc z=(-esmd_apex/2-esmd_front_layer/2) alphaz=-45 ort=x-y-z kOnly='ONLY' ENDIF c-- Prin1 istrip, xc, yc, length ( 'SMD Plane: strip=',I3,' xc=',F5.1,' yc=,'F5.1,' length=',F5.1 ) c-- ENDDO c-- ENDIF c-- c-- * dcut exsg z 0 0 10 0.1 0.1 * dcut exsg y 0 10 -50 0.7 0.7 c-- EndBlock c-- c-- c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EHMS -- c-- Block EHMS defines the triangular SMD strips c-- Material Ecal_scint Medium Ecal_active isvol=1 Attribute EHMS seen=1 colo=2 serial=cut lsty=1 ! red c-- Shape TRD1 dx1=0 dx2=emxg_Sbase/2 dy=length/2 dz=emxg_Sapex/2 c-- HITS EHMS Birk:0:(0,10) c-- Endblock! EHMS c----------------------------------------------------------------------------- c--- c-- Several thin layers of material are applied to the front and back of the c-- SMD planes to provide structural support. We combine these layers into c-- a single effective volume, which is affixed to the base of the SMD c-- strips. As with the SMD strips, z along the depth, y is length c-- c-- http://drupal.star.bnl.gov/STAR/system/files/SMD_module_stack.pdf c-- c-- 1.19 mm G10 c-- 0.25 mm Fiberglass and epoxy c-- 0.17 mm Aluminized mylar c-- c-- Weight in mixture by mass = (depth)*(Area) c-- c-- Weighted density is given by sum (density)_i * (depth)_i / sum (depth)_i c-- c----------------------------------------------------------------- Block EFLS -- c-- Block EFLS is the layer of material on the front of the SMD planes c-- c-- Component G10 A=18.017 Z=9.013 w=1.19*1.700/(1.19*1.700+0.25*1.530+0.17*1.390) Component Fiberglass A=19.103 Z=9.549 w=0.25*1.530/(1.19*1.700+0.25*1.530+0.17*1.390) Component AlMylar A=12.889 Z=6.465 w=0.17*1.390/(1.19*1.700+0.25*1.530+0.17*1.390) Mixture EFLS dens=(1.19*1.7+0.25*1.53+0.17*1.39)/(1.19+0.25+0.17) Attribute EFLS seen=1 colo=22 lsty=1 Shape TRD1 dz=esmd_front_layer/2 dy=length/2 dx1=esmd_base/2 dx2=esmd_base/2 c-- EndBlock! EFLS c-- c-- see link above for documentation c-- c-- 0.10 mm aluminized mylar c-- 0.25 mm fiberglass and epoxy c-- 1.50 mm WLS fiber router layer (polystyrene) c-- 0.25 mm aluminum c-- c----------------------------------------------------------------- Block EBLS -- c-- Block EBLS is the layer of material on the back of the SMD planes c-- Component AlMylar A=12.889 Z=6.465 w=0.10*1.390/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Fiberglass A=19.103 Z=9.549 w=0.25*1.530/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Polystyren A=11.154 Z=5.615 w=1.50*1.032/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Al A=28.08 Z=14.00 w=0.25*2.699/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Mixture EBLS dens=(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699)/(0.10+0.25+1.50+0.25) c-- Attribute EFLS seen=1 colo=22 lsty=1 Shape TRD1 dz=esmd_back_layer/2 dy=length/2 dx1=esmd_base/2 dx2=esmd_base/2 c-- EndBlock! EFLS c----------------------------------------------------------------- Block EXPS -- c-- Block EXPS is the plastic spacer in the shower maximum section c-- c-- Simple implementation of the spacer in the shwoer maximum detector. c-- This implmentation neglects the ears and the source tube. c-- c-- n.b. There may be a side effect in the way this gets created... c-- it could overwrite SMD strips which extend into this plane. c-- Probably need to go with a different approach here. c-- c-- Scanned Drawings: c-- + http://drupal.star.bnl.gov/STAR/system/files/SMD_spacer_drawings.pdf c-- c-- thickness is 1.2 cm, as given by detail B and C... but I do not want c-- to do alot of complicated recoding of the geometry. So I am limiting c-- it to be the same width as a normal SMD volume. c-- Material PVC_Spacer Attribute EXPS seen=1 colo=6 lsty=1 lwid=2 c-- c-- Spacer layers are extended by +/- 5 degrees into the adjacent sectors. c-- The kONLY='Many' option at creation time should mean that conflicts c-- in volume will be resolved in favor of the SMD strips. c-- Shape TUBS dz=esmd_apex/2, rmin=(section)*Tan_Low-1.526, rmax=(section+msecwd)*Tan_Upp, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec c-- EndBlock c-- END c----------------------------------------------------------------- End Module -- c------------------------------------------------------------------------------ c-- Helper subroutines and functions c------------------------------------------------------------------------------ c-- c-- Subroutine ecal_set_cuts(id, medium) c-- c-- id -- integer ID idetifying the current tracking medium c-- medium -- character switch selecting the type of cuts to be c-- used in this tracking volumne c-- c------------------------------------------------------------------------------ Subroutine ecal_set_cuts(id,medium) c-- Implicit NONE Integer id Character medium*(*) c-- Integer radiator, megatile, detector Save radiator, megatile, detector c-- IF ( medium == 'print' ) THEN c-- Write (*,400) radiator Write (*,401) megatile Write (*,402) detector c-- Call GpTMed( +radiator ) Call GpTMed( -megatile ) Call GpTMed( -detector ) c-- Return c-- ENDIF c-- 400 Format('radiator cuts set for ag_imed=',I3) 401 Format('megatile cuts set for ag_imed=',I3) 402 Format('detector cuts set for ag_imed=',I3) c-- c-- c-- Setup common cuts for neutrons, hadrons and muons c-- Call GsTPar (id,'CUTNEU',0.001) Call GsTPar (id,'CUTHAD',0.001) Call GsTPar (id,'CUTMUO',0.001) c-- IF ( medium == 'radiator' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) radiator = id C-- c-- ELSEIF ( medium == 'megatile' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) megatile = id c-- c-- ELSEIF ( medium == 'detector' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) c-- Call GsTPar (id,'BIRK1',1.) Call GsTPar (id,'BIRK2',0.0130) Call GsTPar (id,'BIRK3',9.6E-6) detector = id c-- c-- ELSE Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) Write(*,300) 300 Format('Warning: unknown medium[',A20,'] in ecal_set_cuts') c-- c-- ENDIF c-- Return End c----------------------------------------------------------------------- c----------------------------------------------------------------------- c-- c-- Subroutine ecal_get_strip( section, cut, istrip, xcenter, ycenter, length ) c-- in in in out out out Implicit NONE c-- Real section Integer cut ! 0=no plane 1=normal plane 2=cut plane Integer istrip ! strip index Real xcenter ! output Real ycenter ! output Real length ! output c-- Integer nstrips Real rdel ! shift in radius (?) Real rth Real ddn, ddup Real megatile, p c-- Real xleft, yleft, xright, yright Real dxy, xx, yy Real sqrt2, sqrt3 c-- c-- SMD data copied from data structures above c-- Real base, apex Data base, apex / 1.0, 0.7/ !cm c-- Real Rbot, Rtop Data Rbot, Rtop / 77.41, 213.922 / c-- Real EtaMin, EtaMax Data EtaMin, EtaMax / 1.086, 2.000 / c-- Real tan_theta_min, tan_theta_max c-- Real tanf, eta tanf(eta) = tan(2*atan(exp(-eta))) c-- tan_theta_min = tanf( EtaMax ) tan_theta_max = tanf( EtaMin ) c-- IF (cut = 1) THEN rdel = 3.938 nstrips = 288 ELSE rdel = -.475 nstrips = 285 ENDIF c-- xcenter=0. ycenter=0. length=0. c-- IF ( cut = 0 ) THEN RETURN ENDIF c-- sqrt2 = sqrt(2.0) sqrt3 = sqrt(3.0) c-- rth = .53*rdel ! .53 --- tentatavily jcw-- wtf? ddn = sqrt(3.0)*1.713 + rdel ddup = .5*1.846 + 1.713 megatile = base + .01 c-- p = .5*(istrip-1)*megatile + 41.3655 IF (p <= (.5*rbot*sqrt3 + rth)) THEN dxy = 1.9375*sqrt2 xleft = .5*sqrt2*p*(sqrt3 + 1.) - dxy yleft = .5*sqrt2*p*(sqrt3 - 1.) - dxy yright = .5*sqrt2*(sqrt( rbot*rbot - p*p) - p) xright = sqrt2*p + yright ELSEIF ((.5*rbot*sqrt3 + rth) < p <= (.5*rtop + 1.5)) THEN dxy = 1.9375*sqrt2 xleft = .5*sqrt2*p*(sqrt3 + 1.) - dxy yleft = .5*sqrt2*p*(sqrt3 - 1.) - dxy dxy = rdel*sqrt2/sqrt3 yright = .5*sqrt2*p*(1.- 1./sqrt3) xright = sqrt2*p - yright - dxy yright = -yright - dxy ELSEIF (p > (.5*rtop +1.5)) THEN yleft = (sqrt(rtop*rtop - p*p) - p)/sqrt2 xleft = sqrt2*p + yleft dxy = rdel*sqrt2/sqrt3 yright = .5*sqrt2*p*(1.- 1./sqrt3) xright = sqrt2*p - yright - dxy yright = -yright - dxy dxy = 0. c-- IF ((.5*sqrt3*160.- ddn) < p <= (.5*sqrt3*160.+ ddup) ) THEN xcenter = .5*(sqrt3*160.+1.846) ycenter = xcenter - .5*sqrt3*1.713 IF (p > ycenter) THEN dxy = .5*sqrt2*(2/sqrt3*rdel + .5*sqrt3*1.846 +_ sqrt(1.713*1.713 - (p-xcenter)*(p-xcenter))) ELSE dxy = sqrt2/sqrt3*(p - .5*sqrt3* 160. + ddn) ENDIF ELSEIF ((.5*sqrt3*195.- ddn) < p <= (.5*sqrt3*195. + ddup) ) THEN xcenter = .5*(sqrt3*195.+1.846) ycenter = xcenter - .5*sqrt3*1.713 IF (p > ycenter) THEN dxy = .5*sqrt2*(2/sqrt3*rdel + .5*sqrt3*1.846 +_ sqrt(1.713*1.713 - (p-xcenter)*(p-xcenter))) ELSE dxy = sqrt2/sqrt3*(p - .5*sqrt3*195. + ddn) ENDIF ENDIF xright = xright + dxy yright = yright + dxy ENDIF dxy = section*tan_theta_max - rtop xcenter = .5*(xright+xleft) + dxy ycenter = .5*(yright+yleft) xx = .5*sqrt2*(xleft+yleft) yy = .5*sqrt2*(xright+yright) length = xx-yy c-- c-- Return c-- End! Subroutine smd_strip c-- * ---------------------------------------------------------------------------- * ECAL nice views: dcut ecvo x 1 10 -5 .5 .1 * draw emdi 105 0 160 2 13 .2 .1 * draw emdi 120 180 150 1 14 .12 .12 * --------------------------------------------------------------------------- c-- examples of HITS * HITS EHMS Birk:0:(0,10) * xx:16:SH(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: Eloss:0:(0,10) *
2009.10.12 Jason EEMC geometry: effects of adding layers
Effect of added layers in Jason geometry file (ecalgeo.g23)
Monte-Carlo setup:
- One photon per event
- EEMC only geometry with LOW_EM option
- Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Cuts for shower shapes:
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
All shapes are normalized to 1 at peak (central strip)
Added layer definition from Jason file:
- EXPS is the plastic spacer in the shower maximum section
- EBLS is the layer of material on the back of the SMD planes
- EFLS is the layer of material on the front of the SMD planes
Some comments:
- Figs. 1-2 show that I can reproduce
sampling fraction and shower shapes
which I see with geometry file from CVS
if I disable all three added layers in Jason geometry file
(this assumes/shows that G10 layer have tiny effect).
This a good starting point, since it indicate that
all other (cosmetic) code modifications
are most probably done correctly and has no
effect on simulated detector response. - Fig. 3 shows effect of each added layer
(plastic spacers and layers in front/back of SMD)
on the sampling fraction and 2x1/3x3 energy profile:- Each layer contributes more or less equally to the sampling fraction.
- Energy profile (E2x1 / E3x3) does not affected by the added layers
- Fig. 4 shows effect of each added layer on the shower shapes:
- Back SMD layer does not contribute much (as expected).
- Front and spacers introduce equal amount of "shape narrowing".
- Figs. 5-6 show pre-shower sorted shower shapes
and comparison with eta-meson shapes.
No layers and G10 removed
Figure 1: Sampling fraction vs. thrown energy
Figure 2: Shower shapes
Adding new laters (spacer, front, back)
Figure 3: Sampling fraction vs. thrown energy (left), 2x1/3x3 energy ratio (right)
See legend for details
Figure 4: Shower shapes. See legend for details
Shower shapes sorted by pre-shower energy
Pre-shower bins:
- Ep1 = 0, Ep2 = 0 (no energy in both EEMC pre-shower layers)
- Ep1 = 0, Ep2 > 0
- 0 < Ep1 < 4 MeV
- 4 < Ep1 < 10 MeV
- Ep1 > 10 MeV
- All pre-shower bins combined
Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch
Figure 5: Shower shapes (left) and their ratio (right)
Figure 6: Shower shape ratios
2009.10.13 Jason EEMC geometry: position correlations
Effect of added layers in Jason geometry file (ecalgeo.g23)
Monte-Carlo setup:
- One photon per event
- EEMC only geometry with LOW_EM option
- Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Added layer definition from Jason file:
- EXPS is the plastic spacer in the shower maximum section
- EBLS is the layer of material on the back (routing layers) of the SMD planes
- EFLS is the layer of material on the front (G10, etc) of the SMD planes
Geometry configurations and notations (shown in the center of the plot):
- j-noLayers: Jason geometry: no EXPS, EBLS, EFLS
- j-back: Jason geometry, EBLS only
- j-front: Jason geometry, EFLS only
- j-spacer: Jason geometry, EXPS only
- j-all: Jason geometry, all new layers included
- geom-cvs geometry file from CVS after cAir bug fixed
cross section of 1st SMD plane labeled with "SUV" ordering
Note: u-v ordering scheme can be found here (Fig. 9-11)
Figure 1: Average number of SMD u-strip fired vs. thrown photon's (x,y)
Figure 2:Average number of SMD v-strip fired vs. thrown photon's (x,y)
Figure 3:Average SMD u-energy vs. thrown photon's (x,y)
Figure 4:Average SMD v-energy vs. thrown photon's (x,y)
2009.10.16 Jason geometry file: Full STAR simulations
Monte-Carlo setup:
- One photon per event
- EEMC only and Full STAR geometry configurations with LOW_EM option
(using Victor's geometry fix) - Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- eemc-cvs: EEMC only with geometry file from CVS (cAir-fixed)
- full-cvs: Full STAR with geometry file from CVS (cAir-fixed)
- eemc-j: EEMC only with Jason geometry file
- full-j: Full STAR with Jason geometry file
Figure 1: Sampling fraction
Figure 2: Total energy distribution
Figure 3: Shower shapes (left) and shape ratios (right) for 0 < pre-shower1 < 4MeV
Pre-shower sorted shapes (for completeness)
Figure 4: Shower shapes (all pre-shower bins)
Figure 5: Shower shapes ratio (all pre-shower bins)
2009.10.20 Sampling fraction problem: full STAr vs. EEMC stand alone geometry
For the previous study click here
Monte-Carlo setup:
- One photon per event
- EEMC only and Full STAR geometry configurations with LOW_EM option
(using Victor's geometry fix) - Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- eemc-cvs: EEMC only with geometry file from CVS (cAir-fixed)
- full-cvs: Full STAR with geometry file from CVS (cAir-fixed)
- eemc-j: EEMC only with Jason geometry file
- full-j: Full STAR with Jason geometry file
Figure 1: Average energy in SMD-u plane vs. position of the thrown photon
SMD v (left) and u (right) sampling fraction (E_smd/E_thrown) vs. E_thrown
Figure 2: Sampling fraction (E_tower^total/E_thrown) vs. position of the thrown photon
Sampling fraction (E_tower^total/E_thrown) vs. E_thrown
Figure 3: Number of towers above threshold vs. position of the thrown photon
Number of towers above threshold vs. E_thrown
Other EEMC layers: pre-shower, postshower
Figure 4: (left) Pre-shower1 and (right) Pre-shower2 sampling fraction vs. E_thrown
Figure 5: (left) High tower sampling fraction and (right) residual energy, [E_tot-E_3x3]/E_thrown, vs. E_thrown
2009.10.26 Jason vs. CVS EEMC: removed SMD layers
Monte-Carlo setup:
- One photon per event
- Disabled new SMD layers (EXPS EBLS EFLS) in Jason geometry
- EEMC only and Full STAR geometry configurations with LOW_EM option
Note: LOW_EM option seems not to work for EEMC only configuration (double checking)
(using Victor's geometry fix) - Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- eemc-cvs: EEMC only with geometry file from CVS (cAir-fixed)
- full-cvs: Full STAR with geometry file from CVS (cAir-fixed)
- eemc-j-noL: EEMC only with Jason geometry file (disabled 3-new SMD layers)
- full-j-noL: Full STAR with Jason geometry file (disabled 3-new SMD layers)
Figure 1: number of post-shower tiles
Figure 2: number of pre-1-shower tiles
Figure 3: number of pre-2-shower tiles
Figure 4: number of towers
2D
Figure 5: Average pre-shower1 energy
Figure 6: Average pre-shower2 energy
Figure 7: Average number of SMD-u strips
Figure 8: Average number of SMD-v strips
Figure 9: Average post-shower energy
Sampling fraction
Figure 10: Sampling fraction 1x1 vs. thrown energy
Figure 11: Sampling fraction 2x1 vs. thrown energy
Figure 12: Sampling fraction 3x3 vs. thrown energy
Figure 13: Sampling fraction (total energy) vs. thrown energy
Figure 14: Sampling fraction 1x1
Figure 15: Sampling fraction 2x1
Figure 16: Sampling fraction 3x3
SMD shower shapes
Figure 17: SMD shower shape (v-plane)
2009.10.27 Jason EEMC geometry: effect of removing new SMD layers
Monte-Carlo setup:
- One photon per event
- Disabled/Enabled new SMD layers (EXPS EBLS EFLS) in Jason geometry
- EEMC only and Full STAR geometry configurations with LOW_EM option
(using Victor's geometry fix) - Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- eemc-j: EEMC only with Jason geometry file
- full-j: Full STAR with Jason geometry file
- eemc-j-noL: EEMC only with Jason geometry file (disabled 3-new SMD layers)
- full-j-noL: Full STAR with Jason geometry file (disabled 3-new SMD layers)
Effect of removing SMD layers on SMD strips
Figure 1: Average number of SMD-u strips
Figure 2: Average number of SMD-v strips
Effect of removing SMD layers on sampling fraction
Figure 3: distribution of 1x1 sampling fraction
Figure 4: distribution of 2x1 sampling fraction
Figure 5: distribution of 3x3 sampling fraction
Figure 6: 1x1 sampling fraction vs. thrown energy
Figure 7: 2x1 sampling fraction vs. thrown energy
Figure 8: 3x3 sampling fraction vs. thrown energy
2009.10.27: Jason EEMC geometry: comparison without LOW_EM option
Monte-Carlo setup:
- One photon per event
- Disabled new SMD layers (EXPS EBLS EFLS) in Jason geometry
- EEMC only and Full STAR geometry configurations without LOW_EM option
(using Victor's geometry fix) - Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- eemc-cvs: EEMC only with geometry file from CVS (cAir-fixed)
- full-cvs: Full STAR with geometry file from CVS (cAir-fixed)
- eemc-j-noL: EEMC only with Jason geometry file (disabled 3-new SMD layers)
- full-j-noL: Full STAR with Jason geometry file (disabled 3-new SMD layers)
Figure 1: Sampling fraction 1x1
Figure 2: Sampling fraction 2x1
Figure 3: Sampling fraction 3x3
Figure 4: Sampling fraction total energy
Figure 5: Sampling fraction pre1-shower
Figure 6: Sampling fraction pre2-shower
Figure 7: Sampling fraction smd-u
Figure 8: Sampling fraction smd-v
Figure 9: Sampling fraction post-shower
Sampling fraction vs. thrown energy
Figure 10: Sampling fraction 1x1 vs. thrown energy
Figure 11: Sampling fraction 2x1 vs. thrown energy
Figure 12: Sampling fraction 3x3 vs. thrown energy
Figure 13: Sampling fraction (tatal energy) vs. thrown energy
2009.10.30: Jason EEMC geometry: Jason with ELED block from CVS file
FYI: Alice blog on ELED block study
Monte-Carlo setup:
- One photon per event
- Disabled SMD layers (EXPS EBLS EFLS) in Jason geometry
- Put ELED block from CVS file into Jason geometry
geometry configurations without LOW_EM option
(using Victor's geometry fix)- Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
- Using A2Emaker to get reconstructed Tower/SMD energy
(no EEMC SlowSimulator in chain) - Vertex z=0
- ~50K/per particle type
- Non-zero energy: 3 sigma above pedestal
Geometry configurations and notations (shown in the center of the plot):
- full-cvs: Full STAR with geometry file from CVS (cAir-fixed)
- full-j: EEMC only with Jason geometry file (disabled 3-new SMD layers, ELED block replaced with that from CVS)
Figure 1: Sampling fraction 1x1 (up-left), 2x1 (up-right), 3x3 (low-left), total energy (low-right)
Figure 2: Sampling fraction pre1 (up-left), pre2 (up-right), SMD-u (low-left), post (low-right)
Figure 3: Shower shapes (left) and shower shape ratio (right)
