Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • EIC/detectors/athena
  • zwzhao/athena
  • FernandoTA/athena
  • palspeic/athena
4 results
Show changes
Showing
with 2340 additions and 214 deletions
<lccdd>
<comment>
Simplified homogenous implementation of the forward ScFi calorimeter
</comment>
<define>
<constant name="EcalEndcapP_FiberRadius" value="0.235*cm"/>
<constant name="EcalEndcapP_FiberOffset" value="0.5*mm"/>
<constant name="EcalEndcapP_FiberSpaceX" value="0.265*mm"/>
<constant name="EcalEndcapP_FiberSpaceY" value="0.425*mm"/>
</define>
<limits>
</limits>
<regions>
</regions>
<!-- Common Generic visualization attributes -->
<comment>Common Generic visualization attributes</comment>
<display>
</display>
<detectors>
<comment>
------------------------------------------
Forward (Positive Z) Endcap EM Calorimeter
------------------------------------------
An EM calorimeter with ScFi modules
</comment>
<detector id="ECalEndcapP_ID"
name="EcalEndcapP"
type="ScFiCalorimeter"
vis="EcalEndcapVis"
readout="EcalEndcapPHits">
<position x="0" y="0" z="EcalEndcapP_zmin + EcalEndcapP_length/2."/>
<dimensions rmin="EcalEndcapP_rmin" rmax="EcalEndcapP_rmax" length="EcalEndcapP_length"/>
<module sizex="25*mm" sizey="25*mm" sizez="170*mm" material="AvgTungstenScFi" vis="EcalEndcapBlockVis">
</module>
</detector>
</detectors>
<!-- Definition of the readout segmentation/definition -->
<readouts>
<readout name="EcalEndcapPHits">
<segmentation type="NoSegmentation"/>
<id>system:8,ring:8,module:20,fiber_x:8,fiber_y:8</id>
</readout>
</readouts>
<plugins>
</plugins>
</lccdd>
<lccdd>
<define>
<constant name="EcalEndcapP_FiberRadius" value="0.235*cm"/>
<constant name="EcalEndcapP_FiberOffset" value="0.5*mm"/>
<constant name="EcalEndcapP_FiberSpaceX" value="0.265*mm"/>
<constant name="EcalEndcapP_FiberSpaceY" value="0.425*mm"/>
</define>
<limits>
</limits>
<regions>
</regions>
<!-- Common Generic visualization attributes -->
<comment>Common Generic visualization attributes</comment>
<display>
</display>
<detectors>
<comment>
------------------------------------------
Forward (Positive Z) Endcap EM Calorimeter
------------------------------------------
An EM calorimeter with ScFi modules
</comment>
<detector id="ECalEndcapP_ID"
name="EcalEndcapP"
type="ScFiCalorimeter"
vis="EcalEndcapVis"
readout="EcalEndcapPHits">
<position x="0" y="0" z="EcalEndcapP_zmin + EcalEndcapP_length/2."/>
<dimensions rmin="EcalEndcapP_rmin" rmax="EcalEndcapP_rmax" length="EcalEndcapP_length"/>
<module sizex="25*mm" sizey="25*mm" sizez="170*mm" material="TungstenDens24" vis="EcalEndcapBlockVis">
<fiber material="Polystyrene"
radius="EcalEndcapP_FiberRadius"
offset="EcalEndcapP_FiberOffset"
spacex="EcalEndcapP_FiberSpaceX"
spacey="EcalEndcapP_FiberSpaceY"/>
</module>
</detector>
</detectors>
<!-- Definition of the readout segmentation/definition -->
<readouts>
<readout name="EcalEndcapPHits">
<segmentation type="NoSegmentation"/>
<id>system:8,ring:8,module:20,fiber_x:8,fiber_y:8</id>
</readout>
</readouts>
<plugins>
</plugins>
</lccdd>
<lccdd>
<define>
<constant name="DIRCModule_count" value="CaloSides" />
<constant name="DIRCBar_thickness" value="17*mm" />
<constant name="DIRCBar_length" value="DIRC_length-8*cm" />
<constant name="DIRCSkinFront_thickness" value="2.5*mm" />
<constant name="DIRCSkinBack_thickness" value="2.5*mm" />
<constant name="DIRCFrame_width" value="20*mm"/>
<constant name="DIRCFrame_length" value="DIRC_length"/>
<constant name="DIRCFrame_thickness" value="40*mm"/>
<constant name="DIRCModule_rmax" value="DIRC_rmin + DIRCBar_thickness + DIRCSkinFront_thickness + DIRCSkinBack_thickness"/>
<constant name="DIRCModule_thickness" value="DIRCModule_rmax - DIRC_rmin" />
<constant name="DIRCFrame_rmax" value="DIRC_rmin + DIRCFrame_thickness" />
<constant name="DIRCModule_halfangle" value="180*degree/DIRCModule_count" />
<constant name="DIRCModule_width" value="2*DIRC_rmin * tan(DIRCModule_halfangle)"/>
<constant name="DIRCModule_rmin" value="DIRC_rmin + 1*cm"/>
</define>
<materials>
</materials>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<detectors>
<detector id="BarrelDIRC_ID" name="FakeDIRC" type="athena_FakeDIRC" readout="DIRCBarHits" vis="DIRCVis">
<dimensions rmin="DIRC_rmin" rmax="DIRC_rmax" length="DIRC_length" />
<position x="0" y="0" z="DIRC_offset" />
<comment> Fake DIRC modules </comment>
<module name="DIRCModule" vis="DIRCModuleVis">
<module_component name="FrontSkin"
material="CarbonFiber"
sensitive="false"
width="DIRCModule_width"
thickness="DIRCSkinFront_thickness"
vis="DIRCSupportVis"
length="DIRCBar_length" />
<module_component name="QuartzBar"
material="Quartz"
sensitive="true"
width="DIRCModule_width"
thickness="DIRCBar_thickness"
vis="DIRCBarVis"
length="DIRCBar_length" />
<module_component name="BackSkin"
material="CarbonFiber"
sensitive="false"
width="DIRCModule_width"
thickness="DIRCSkinBack_thickness"
vis="DIRCSupportVis"
length="DIRCBar_length" />
<frame material="StainlessSteel"
width="DIRCFrame_width"
length="DIRCFrame_length"
vis="DIRCFrameVis"
thickness="DIRCFrame_thickness" />
</module>
<comment> Fake DIRC layers </comment>
<layer module="DIRCModule" id="1" vis="DIRCLayerVis">
<barrel_envelope
inner_r="DIRC_rmin"
outer_r="DIRC_rmax"
z_length="DIRC_length" />
<rphi_layout
phi_tilt="0"
nphi="DIRCModule_count"
phi0="0"
rc="0.5*(DIRCModule_rmin+DIRCModule_rmax)"
dr="0" />
<z_layout
dr="0.0*mm"
z0="0.0*mm"
nz="1" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="DIRCBarHits">
<segmentation type="CartesianGridXY" grid_size_x="3.0*mm" grid_size_y="3.0*mm" />
<id>system:8,layer:4,module:8,section:4,x:32:-16,y:-16</id>
</readout>
</readouts>
<plugins>
</plugins>
<fields>
</fields>
</lccdd>
<lccdd>
<comment> Endcap Tracker </comment>
<comment> Endcap Tracker for Acadia </comment>
<define>
<constant name="GEMTrackerEndcap_zmin" value="BarrelTracking_zmax"/>
<constant name="GEMTrackerEndcap_NLayers" value="1"/>
<constant name="GEMTrackerEndcapFoilOpeningAngle" value="30.1*degree"/>
<constant name="GEMTrackerEndcapP_thickness" value="min(10.0*cm,ForwardTrackingRegion_length)" />
<constant name="GEMTrackerEndcapN_thickness" value="min(10.0*cm,BackwardTrackingRegion_length)" />
<constant name="GEMTrackerEndcapP_zmin" value="ForwardTrackingRegion_zmin"/>
<constant name="GEMTrackerEndcapN_zmin" value="BackwardTrackingRegion_zmin"/>
<constant name="GEMTrackerEndcapP_rmin" value="19*cm" />
<constant name="GEMTrackerEndcapN_rmin" value="10*cm" />
<constant name="GEMTrackerEndcapP_rmax" value="200*cm" />
<constant name="GEMTrackerEndcapN_rmax" value="BackwardPIDRegion_rmax" />
<constant name="GEMTrackerEndcapOpeningAngle" value="30.4*degree"/>
<constant name="GEMTrackerEndcap_dz" value="15*mm" />
<comment>
____X2____
\ /
\ / Y (=z for the trd)
\_X1_/
</comment>
<constant name="GEMTrackerEndcapFoilX1" value="43.0 * mm"/>
<constant name="GEMTrackerEndcapLayer_thickness" value="ForwardTrackingGEMLength/GEMTrackerEndcap_NLayers"/>
<constant name="GEMTrackerEndcapFoil_rmin" value="GEMTrackerEndcapFoilX1/(tan(GEMTrackerEndcapFoilOpeningAngle/2.0)*2.0)"/>
<constant name="GEMTrackerEndcapFoilY" value="BarrelTracking_rmax- GEMTrackerEndcapFoil_rmin- 1.0*cm"/>
<constant name="GEMTrackerEndcapFoilX2" value="2*GEMTrackerEndcapFoilY*tan(GEMTrackerEndcapFoilOpeningAngle/2.0) + GEMTrackerEndcapFoilX1"/>
<constant name="GEMTrackerEndcapFoil_rmax" value="(GEMTrackerEndcapFoilY+GEMTrackerEndcapFoil_rmin)/(cos(GEMTrackerEndcapFoilOpeningAngle/2.0))"/>
<constant name="GEMTrackerEndcapP_x1" value="2.0*GEMTrackerEndcapP_rmin*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapP_x2" value="2.0*GEMTrackerEndcapP_rmax*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapP_y" value="GEMTrackerEndcapP_rmax*cos(GEMTrackerEndcapOpeningAngle/2.0) - GEMTrackerEndcapP_rmin"/>
<constant name="GEMTrackerEndcapN_x1" value="2.0*GEMTrackerEndcapN_rmin*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapN_x2" value="2.0*GEMTrackerEndcapN_rmax*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapN_y" value="GEMTrackerEndcapN_rmax*cos(GEMTrackerEndcapOpeningAngle/2.0) - GEMTrackerEndcapN_rmin"/>
<comment>
Needs implimentation:
Carbon fiber support ring and spokes
Initial preliminary details can be found here: https://www.dropbox.com/s/grt9p8feqlsh52f/GEM-RD-GEMWheel.pdf?dl=0
This design is still very much preliminary
// Frame parameters accoring to Kondo's sbsCrossSection.pdf file),
<constant name="GEMTrackerEndcapFrame_thickness" value="18.0 * mm"/>
//frame thickness is from 1*3mm (drift gap) + 5*2mm frames
</comment>
<constant name="GEMTrackerEndcapFrame_thickness" value="13.0 * mm"/>
<constant name="GEMTrackerEndcapFrameBotEdge_width" value="10.0 * mm"/>
<constant name="GEMTrackerEndcapFrameTopEdge_width" value="10.0 * mm"/>
<constant name="GEMTrackerEndcapFrameTopEdge_width" value="15.0 * mm"/>
<constant name="GEMTrackerEndcapFrameSideEdge_width" value=" 8.0 * mm"/>
<comment> Account for the fact that GEMs have holes, so they are not solid volumes </comment>
<constant name="GEMAreaFactor" value="0.80" />
<comment>
From https://github.com/eic/EicToyModel/blob/master/source/eicroot/GemGeoParData.cxx
From https://github.com/eic/EicToyModel/blob/master/source/eicroot/GEMGeoParData.cxx
// FIXME: put aluminum layer later as well),
mEntranceWindowMaterial ( "GemKapton"),
mEntranceWindowMaterial ( "GEMKapton"),
mEntranceWindowThickness ( 50.0 * etm::um),
// Use evaristo.pdf p.10 for the foil parameters:
......@@ -39,16 +55,16 @@
// - readout foils : 30um kapton + 3um copper total)
mDriftFoilKaptonThickness ( 50.0 * etm::um),
mDriftFoilCopperThickness ( 3.0 * etm::um),
mGemFoilAreaFraction ( 0.80),
mGemFoilKaptonThickness ( 30.0 * etm::um),
mGemFoilCopperThickness ( 3.0 * etm::um),
mGEMFoilAreaFraction ( 0.80),
mGEMFoilKaptonThickness ( 30.0 * etm::um),
mGEMFoilCopperThickness ( 3.0 * etm::um),
mReadoutG10Thickness ( 0.0 * etm::mm),
mReadoutKaptonThickness ( 30.0 * etm::um),
mReadoutCopperThickness ( 3.0 * etm::um),
// 3mm thick Nomex honeycomb for SBS GEMs),
mReadoutSupportMaterial ( "GemNomex"),
mReadoutSupportMaterial ( "GEMNomex"),
mReadoutSupportThickness ( 3.0 * etm::mm),
mEntranceRegionLength ( 3.0 * etm::mm),
......@@ -71,136 +87,138 @@
<detectors>
<detector
id="GEMTrackerEndcap_ID"
name="GEMTrackerEndcap"
type="refdet_GEMTrackerEndcap"
id="TrackerEndcapP_3_ID"
name="GEMTrackerEndcapP"
type="athena_GEMTrackerEndcap"
readout="GEMTrackerEndcapHits"
vis="AnlRed"
vis="TrackerVis"
reflect="false">
<module name="GEMModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<module name="GEMModule1" vis="TrackerLayerVis">
<trd x1="GEMTrackerEndcapP_x1/2.0" x2="GEMTrackerEndcapP_x2/2.0" z="GEMTrackerEndcapP_y/2"/>
<comment>
Frame disabled as it gives issue with the ACTS geometry conversion
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapFoilX1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
x2="GEMTrackerEndcapFoilX2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
z="GEMTrackerEndcapFoilY/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<trd x1="GEMTrackerEndcapP_x1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
x2="GEMTrackerEndcapP_x2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
z="GEMTrackerEndcapP_y/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</comment>
<comment> Going from HV side to readout side</comment>
<module_component thickness="0.127 * mm" material="Mylar"/>
<module_component thickness="50.0*um" material="Kapton" name="entrance_window"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="entrance region" />
<module_component thickness="50.0*um" material="Kapton"/>
<module_component thickness=" 3.0*um" material="Copper"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="drift region"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region I"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region II"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction region"/>
<module_component thickness="30.0*um" material="Kapton" name="readout" sensitive="true"/>
<module_component thickness=" 3.0*um" material="Copper" name="readout_Cu"/>
<module_component thickness="127.0*um" material="Mylar"/>
<module_component thickness="200.0*um" material="Epoxy" sensitive="true" vis="AnlProcess_Blue"/>
<module_component thickness="25.0*um" material="Mylar" name="entrance_window" vis="TrackerLayerVis"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="entrance_region" />
<module_component thickness="50.0*um" material="Kapton" name="cathode_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="cathode_cu"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" sensitive="true" name="drift_region" vis="AnlProcess_Blue"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_1"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_1_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_I"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_2"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_2_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_II"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_3"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_3_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction_region"/>
<module_component thickness="50.0*um" material="Kapton" name="readout_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="readout_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="exit_region"/>
<module_component thickness="25.0*um" material="Mylar" name="exit_window" vis="TrackerLayerVis"/>
</module>
<module name="GEMSupportModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapFoilX1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
x2="GEMTrackerEndcapFoilX2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
z="GEMTrackerEndcapFoilY/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</module>
<module name="GEMSupportModule2" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapFoilY/2"/>
<module_component thickness="1.0*mm" material="Mylar"/>
<module name="GEMSupportSpokes" vis="TrackerSupportVis">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapP_y/2"/>
<module_component thickness="1.0*mm" material="Mylar" vis="TrackerSupportVis"/>
</module>
<layer id="1" >
<ring vis="AnlRed" module="GEMModule1"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm" />
<!--
<ring vis="AnlRed" module="GEMSupportModule1"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm"/>
-->
<ring vis="AnlBlue" phi0="15.0*degree" module="GEMSupportModule2"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="0 * mm" />
<envelope vis="TrackerVis"
rmin = "GEMTrackerEndcapP_rmin"
rmax = "GEMTrackerEndcapP_rmax"
length="GEMTrackerEndcapP_thickness"
zstart="GEMTrackerEndcapP_zmin" />
<layer_material surface="representing" binning="binPhi,binR" bins0="40*12" bins1="1024"/>
<ring vis="TrackerVis" module="GEMModule1"
r="GEMTrackerEndcapP_rmin+GEMTrackerEndcapP_y/2.0"
zstart="0"
nmodules="12" dz="GEMTrackerEndcap_dz" />
<ring vis="AnlBlue" phi0="GEMTrackerEndcapOpeningAngle/2." module="GEMSupportSpokes"
r="GEMTrackerEndcapP_rmin+GEMTrackerEndcapP_y/2.0"
zstart="0"
nmodules="12" dz="0*mm" />
</layer>
</detector>
<detector
id="GEMTrackerEndcapN_ID"
id="TrackerEndcapN_3_ID"
name="GEMTrackerEndcapN"
type="refdet_GEMTrackerEndcap"
type="athena_GEMTrackerEndcap"
readout="GEMTrackerEndcapHits"
vis="AnlRed"
vis="TrackerVis"
reflect="true">
<module name="GEMModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<module name="GEMModule1" vis="TrackerLayerVis">
<trd x1="GEMTrackerEndcapN_x1/2.0" x2="GEMTrackerEndcapN_x2/2.0" z="GEMTrackerEndcapN_y/2"/>
<comment>
Frame disabled as it gives issue with the ACTS geometry conversion
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapN_x1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
x2="GEMTrackerEndcapN_x2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
z="GEMTrackerEndcapN_y/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</comment>
<comment> Going from HV side to readout side</comment>
<module_component thickness="0.127 * mm" material="Mylar"/>
<module_component thickness="50.0*um" material="Kapton" name="entrance_window"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="entrance region" />
<module_component thickness="50.0*um" material="Kapton"/>
<module_component thickness=" 3.0*um" material="Copper"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="drift region"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region I"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region II"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction region"/>
<module_component thickness="30.0*um" material="Kapton" name="readout" sensitive="true"/>
<module_component thickness=" 3.0*um" material="Copper" name="readout_Cu"/>
<module_component thickness="127.0*um" material="Mylar"/>
<module_component thickness="200.0*um" material="Epoxy" sensitive="true" vis="AnlProcess_Blue"/>
<module_component thickness="25.0*um" material="Mylar" name="entrance_window" vis="TrackerLayerVis"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="entrance_region" />
<module_component thickness="50.0*um" material="Kapton" name="cathode_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="cathode_cu"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" sensitive="true" name="drift_region" vis="AnlProcess_Blue"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_1"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_1_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_I"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_2"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_2_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_II"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_3"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_3_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction_region"/>
<module_component thickness="50.0*um" material="Kapton" name="readout_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="readout_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="exit_region"/>
<module_component thickness="25.0*um" material="Mylar" name="exit_window" vis="TrackerLayerVis"/>
</module>
<module name="GEMSupportModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX2/2.0" x2="GEMTrackerEndcapFoilX1/2.0" z="GEMTrackerEndcapFrameBotEdge_width"/>
<module_component thickness="GEMTrackerEndcapFrame_thickness" material="Mylar"/>
</module>
<module name="GEMSupportModule2" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapFoilY/2"/>
<module_component thickness="4.0*mm" material="Mylar"/>
<module name="GEMSupportSpokes" vis="TrackerSupportVis">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapN_y/2"/>
<module_component thickness="1.0*mm" material="Mylar" vis="TrackerSupportVis"/>
</module>
<layer id="1" >
<ring vis="AnlRed"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm" module="GEMModule1" />
<ring vis="AnlBlue" phi0="15.0*degree"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="0 * mm" module="GEMSupportModule2" />
<envelope vis="TrackerVis"
rmin = "GEMTrackerEndcapN_rmin"
rmax = "GEMTrackerEndcapN_rmax"
length="GEMTrackerEndcapN_thickness"
zstart="GEMTrackerEndcapN_zmin" />
<layer_material surface="representing" binning="binPhi,binR" bins0="40*12" bins1="1024"/>
<ring vis="TrackerVis" module="GEMModule1"
r="GEMTrackerEndcapN_rmin+GEMTrackerEndcapN_y/2.0"
zstart="0.1"
nmodules="12" dz="GEMTrackerEndcap_dz" />
<ring vis="AnlBlue" phi0="GEMTrackerEndcapOpeningAngle/2." module="GEMSupportSpokes"
r="GEMTrackerEndcapN_rmin+GEMTrackerEndcapN_y/2.0"
zstart="0.1"
nmodules="12" dz="0*mm" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="GEMTrackerEndcapHits">
<segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
<id>system:8,barrel:2,layer:4,module:12,sensor:2,x:32:-16,y:-16</id>
<segmentation type="CartesianGridXZ" grid_size_x="50*um" grid_size_z="250*um" />
<id>system:8,layer:4,module:12,sensor:2,x:32:-16,z:-16</id>
</readout>
</readouts>
......
<lccdd>
<define>
<constant name="HcalBarrel_offset" value="(HcalEndcapP_zmin-HcalEndcapN_zmin)/2.0" />
<constant name="HcalBarrel_length" value="HcalBarrelForward_zmax + HcalBarrelBackward_zmax"/>
<constant name="HcalBarrel_offset" value="(HcalBarrelForward_zmax - HcalBarrelBackward_zmax)/2"/>
<constant name="HcalBarrelAvailThickness" value="HcalBarrel_thickness"/>
<documentation>
#### Material Thickness
</documentation>
<constant name="HcalBarrelSteelThickness" value="2.0 * cm"/>
<constant name="HcalBarrelPolystyreneThickness" value="0.5 * cm"/>
<constant name="HcalEndcapNSteelThickness" value="HcalBarrelSteelThickness"/>
<constant name="HcalEndcapNPolystyreneThickness" value="HcalBarrelPolystyreneThickness"/>
<constant name="HcalEndcapPSteelThickness" value="HcalBarrelSteelThickness"/>
<constant name="HcalEndcapPPolystyreneThickness" value="0.3 * cm"/>
<documentation>
- Hcal Barrel Layers and computed Thickness
</documentation>
<constant name="HcalBarrelSingleLayerThickness"
value="HcalBarrelSteelThickness + HcalBarrelPolystyreneThickness"/>
<constant name="HcalBarrelLayer_NRepeat" value="floor(HcalBarrelAvailThickness / HcalBarrelSingleLayerThickness)"/>
<constant name="HcalBarrelThickness" value="HcalBarrelLayer_NRepeat * HcalBarrelSingleLayerThickness"/>
<documentation>
- Hcal Endcap N Layers and computed Thickness
</documentation>
<constant name="HcalEndcapNSingleLayerThickness"
value="HcalEndcapNSteelThickness + HcalEndcapNPolystyreneThickness"/>
<constant name="HcalEndcapNLayer_NRepeat" value="floor(HcalEndcapN_length / HcalEndcapNSingleLayerThickness)"/>
<constant name="HcalEndcapN_thickness" value="HcalEndcapNLayer_NRepeat * HcalEndcapNSingleLayerThickness"/>
<documentation>
- Hcal Endcap P Layers and computed Thickness
</documentation>
<constant name="HcalEndcapPSingleLayerThickness"
value="HcalEndcapPSteelThickness + HcalEndcapPPolystyreneThickness"/>
<constant name="HcalEndcapPLayer_NRepeat" value="floor(HcalEndcapP_length / HcalEndcapPSingleLayerThickness)"/>
<constant name="HcalEndcapP_thickness" value="HcalEndcapPLayer_NRepeat * HcalEndcapPSingleLayerThickness"/>
</define>
<limits>
......@@ -11,11 +48,6 @@
</regions>
<display>
<vis name="HcalVis" alpha="0.9" r="0.75" g="0.0" b="0.0" showDaughters="false"/>
<vis name="HcalBarrelVis" alpha="0.9" r="0.0" g="0.3" b="0.8" showDaughters="true" visible="true" />
<vis name="HcalEndcapVis" alpha="0.9" r="0.0" g="0.3" b="0.8" showDaughters="true" visible="true" />
<vis name="HcalSensorVis" alpha="0.9" r="1.0" g="0.0" b="0.2" showDaughters="true" visible="true" />
<vis name="HcalAbsorberVis" alpha="0.9" r="0.4" g="0.4" b="0.6" showDaughters="true" visible="true" />
</display>
<!-- Define detector -->
......@@ -26,6 +58,7 @@
name="HcalBarrel"
type="athena_EcalBarrel"
readout="HcalBarrelHits"
vis="HcalBarrelVis"
calorimeterType="HAD_BARREL"
offset="HcalBarrel_offset"
gap="0.*cm"
......@@ -33,19 +66,19 @@
<dimensions
numsides="CaloSides"
rmin="HcalBarrel_rmin"
z="HcalBarrelLength"/>
<staves vis="HcalBarrelVis"/>
<layer repeat="HcalBarrelLayer_NRepeat">
<slice material="Steel235" thickness="2.0*cm" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="0.5*cm" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
z="HcalBarrel_length"/>
<staves vis="HcalBarrelStaveVis"/>
<layer repeat="HcalBarrelLayer_NRepeat" vis="HcalBarrelLayerVis">
<slice material="Steel235" thickness="HcalBarrelSteelThickness" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="HcalBarrelPolystyreneThickness" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
</layer>
</detector>
<detector
id="HCalEndcapP_ID"
name="HcalEndcapP"
type="refdet_PolyhedraEndcapCalorimeter2"
readout="HcalHadronEndcapHits"
type="athena_PolyhedraEndcapCalorimeter2"
readout="HcalEndcapPHits"
vis="HcalEndcapVis"
calorimeterType="HAD_ENDCAP"
reflect="false">
......@@ -55,17 +88,35 @@
zmin="HcalEndcapP_zmin"
rmin="HcalEndcapP_rmin"
rmax="HcalBarrel_rmax"/>
<layer repeat="51" vis="HcalEndcapVis" >
<slice material="Steel235" thickness="2.0" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="0.3*cm" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
<layer repeat="HcalEndcapPLayer_NRepeat" vis="HcalEndcapLayerVis" >
<slice material="Steel235" thickness="HcalEndcapPSteelThickness" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="HcalEndcapPPolystyreneThickness" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
</layer>
</detector>
<detector
id="PassiveSteelRingEndcapP_ID"
name="PassiveSteelRingEndcapP"
type="athena_PolyhedraEndcapCalorimeter2"
vis="PassiveSteelVis"
calorimeterType="HAD_ENDCAP"
reflect="false">
<position x="0" y="0" z="0"/>
<dimensions
numsides="CaloSides"
zmin="EcalEndcapP_zmin"
rmin="EcalEndcapP_rmax"
rmax="HcalBarrel_rmax"/>
<layer repeat="1">
<slice material="Steel235" thickness="EcalEndcapP_length" vis="HcalAbsorberVis"/>
</layer>
</detector>
<detector
id="HCalEndcapN_ID"
name="HcalEndcapN"
type="refdet_PolyhedraEndcapCalorimeter2"
readout="HcalElectronEndcapHits"
type="athena_PolyhedraEndcapCalorimeter2"
readout="HcalEndcapNHits"
vis="HcalEndcapVis"
calorimeterType="HAD_ENDCAP" reflect="true">
<position x="0" y="0" z="0"/>
......@@ -74,9 +125,9 @@
zmin="HcalEndcapN_zmin"
rmin="HcalEndcapN_rmin"
rmax="HcalBarrel_rmax"/>
<layer repeat="HcalEndcapNLayer_NRepeat" vis="HcalEndcapVis" >
<slice material="Steel235" thickness="HcalSteelThickness" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="0.5*cm" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
<layer repeat="HcalEndcapNLayer_NRepeat" vis="HcalEndcapLayerVis" >
<slice material="Steel235" thickness="HcalEndcapNSteelThickness" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="HcalEndcapNPolystyreneThickness" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
</layer>
</detector>
......@@ -87,11 +138,11 @@
<segmentation type="CartesianGridXY" grid_size_x="100* mm" grid_size_y="100 * mm"/>
<id>system:8,barrel:3,module:4,layer:8,slice:5,x:32:-16,y:-16</id>
</readout>
<readout name="HcalHadronEndcapHits">
<readout name="HcalEndcapPHits">
<segmentation type="CartesianGridXY" grid_size_x="100 * mm" grid_size_y="100 * mm"/>
<id>system:8,barrel:3,module:4,layer:8,slice:5,x:32:-16,y:-16</id>
</readout>
<readout name="HcalElectronEndcapHits">
<readout name="HcalEndcapNHits">
<segmentation type="CartesianGridXY" grid_size_x="100 * mm" grid_size_y ="100 * mm"/>
<id>system:8,barrel:3,module:4,layer:8,slice:5,x:32:-16,y:-16</id>
</readout>
......
......@@ -122,7 +122,7 @@
<detector
id="Solenoid_ID"
name="SolenoidCoilBarrel"
type="refdet_SolenoidCoil"
type="athena_SolenoidCoil"
insideTrackingVolume="false">
<position x="0" y="0" z="Solenoid_offset"/>
<layer id="1"
......
......@@ -24,7 +24,7 @@
<composite n="15" ref="C"/>
<composite n="7" ref="O"/>
</material>
<material name="Quartz">
<material name="Quartz"> <!-- different density than `SiliconDioxide` -->
<D type="density" value="2.2" unit="g/cm3"/>
<composite n="1" ref="Si"/>
<composite n="2" ref="O"/>
......@@ -45,11 +45,16 @@
<fraction n="0.998" ref="Fe"/>
<fraction n=".002" ref="C"/>
</material>
<material name="SiliconOxide">
<material name="SiliconOxide"> <!-- different density than `SiliconDioxide` -->
<D type="density" value="2.65" unit="g/cm3"/>
<composite n="1" ref="Si"/>
<composite n="2" ref="O"/>
</material>
<material name="SiliconDioxide"> <!-- density from `G4_SILICON_DIOXIDE` (NIST DB) -->
<D type="density" value="2.32" unit="g/cm3"/>
<composite n="1" ref="Si"/>
<composite n="2" ref="O"/>
</material>
<material name="BoronOxide">
<D type="density" value="2.46" unit="g/cm3"/>
<composite n="2" ref="B"/>
......@@ -92,6 +97,13 @@
<fraction n="0.040" ref="SodiumOxide"/>
<fraction n="0.023" ref="AluminumOxide"/>
</material>
<material name="PyrexGlass25">
<D type="density" value="0.5575" unit="g/cm3"/>
<fraction n="0.806" ref="SiliconOxide"/>
<fraction n="0.130" ref="BoronOxide"/>
<fraction n="0.040" ref="SodiumOxide"/>
<fraction n="0.023" ref="AluminumOxide"/>
</material>
<material name="PbGlass">
<D type="density" value="3.86" unit="g / cm3"/>
<fraction n="0.5080" ref="LeadOxide"/>
......@@ -100,6 +112,46 @@
<fraction n="0.0278" ref="SodiumOxide"/>
<fraction n="0.0050" ref="ArsenicOxide"/>
</material>
<material name="SciGlass">
<D type="density" value="4.22" unit="g / cm3"/>
<fraction n="0.3875" ref="Ba"/>
<fraction n="0.2146" ref="Gd"/>
<fraction n="0.1369" ref="Si"/>
<fraction n="0.2610" ref="O"/>
</material>
<documentation level="3">
#### Material for TOF
a level 3 doc
</documentation>
<material name="CarbonFoam">
<D type="density" value="0.09" unit="g/cm3"/>
<fraction n="1.00" ref="C"/>
</material>
<material name="CFRPMix">
<D type="density" value="0.03" unit="g/cm3"/>
<fraction n="0.074" ref="O"/>
<fraction n="0.870" ref="C"/>
<fraction n="0.023" ref="H"/>
<fraction n="0.033" ref="Cl"/>
</material>
<material name="CFRPMix2">
<D type="density" value="1.60" unit="g/cm3"/>
<fraction n="0.074" ref="O"/>
<fraction n="0.870" ref="C"/>
<fraction n="0.023" ref="H"/>
<fraction n="0.033" ref="Cl"/>
</material>
<material name="NOVEC7200">
<D type="density" value="1.44" unit="g/cm3"/>
<fraction n="0.061" ref="O"/>
<fraction n="0.019" ref="H"/>
<fraction n="0.272" ref="C"/>
<fraction n="0.648" ref="F"/>
</material>
<documentation level="3">
#### Carbon fiber
a level 3 doc
</documentation>
<material name="CarbonFiber">
<D type="density" value="1.5" unit="g/cm3"/>
<fraction n="0.65" ref="C"/>
......@@ -171,6 +223,12 @@
<fraction n="0.044" ref="Ni"/>
<fraction n="0.006" ref="Fe"/>
</material>
<material name="AvgTungstenScFi">
<D value="10.15" unit="g / cm3"/>
<fraction n="0.948" ref="W"/>
<fraction n="0.052*19./40" ref="C"/>
<fraction n="0.052*21./40" ref="H"/>
</material>
<material name="CarbonFiber_25percent">
<D type="density" value="0.375" unit="g / cm3"/>
<fraction n="1.0" ref="CarbonFiber"/>
......@@ -227,6 +285,12 @@
<fraction n="0.041960452" ref="H"/>
<fraction n="0.33301847" ref="O"/>
</material>
<material name="MylarMLI">
<D type="density" value="0.08" unit="g/cm3" />
<fraction n="0.62502108" ref="C"/>
<fraction n="0.041960452" ref="H"/>
<fraction n="0.33301847" ref="O"/>
</material>
<material name="Nomex">
<D type="density" value="0.029" unit="g/cm3" />
<composite n="14" ref="C"/>
......@@ -254,5 +318,68 @@
<D type="density" value="2.21" unit="g/cm3"/>
<composite n="1" ref="C"/>
</material>
<material name="PolyvinylAcetate">
<D type="density" value="1.19" unit="g/cm3"/>
<composite n="4" ref="C"/>
<composite n="6" ref="H"/>
<composite n="2" ref="O"/>
</material>
<material name="Plexiglass">
<D type="density" value="1.19" unit="g/cm3"/>
<composite n="5" ref="C"/>
<composite n="8" ref="H"/>
<composite n="2" ref="O"/>
</material>
<material name="StainlessSteel">
<D type="density" value="8.3" unit="g / cm3"/>
<fraction n="0.74" ref="Fe"/>
<fraction n="0.18" ref="Cr"/>
<fraction n="0.08" ref="Ni"/>
</material>
<material name="ZDC_Tungsten">
<D type="density" value="18.73" unit="g / cm3"/>
<fraction n="0.943" ref="W"/>
<fraction n="0.038" ref="Ni"/>
<fraction n="0.019" ref="Cu"/>
</material>
<material name="ZDC_EpoxyGlue">
<D type="density" value="1.38" unit="g/cm3"/>
<composite n="8" ref="H"/>
<composite n="10" ref="C"/>
<composite n="4" ref="O"/>
</material>
<material name="ZDC_Scintillator">
<D type="density" value="1.032" unit="g/cm3"/>
<composite n="8" ref="C"/>
<composite n="8" ref="H"/>
</material>
<material name="MMGAS_ResistivePaste">
<D type="density" unit="g/cm3" value="0.77906"/>
<composite n="1" ref="C"/>
</material>
<material name="MMGAS_InoxForMesh" state="solid">
<D type="density" unit="g/cm3" value="2.8548"/>
<fraction n="0.19" ref="Cr"/>
<fraction n="0.68" ref="Fe"/>
<fraction n="0.02" ref="Mn"/>
<fraction n="0.10" ref="Ni"/>
<fraction n="0.01" ref="Si"/>
</material>
<material name="Ar90IsoButane" state="gas">
<D type="density" unit="g/cm3" value="0.00170335"/>
<fraction n="0.90" ref="Ar"/>
<fraction n="0.0826586" ref="C"/>
<fraction n="0.0173414" ref="H"/>
</material>
<material name="ZDC_WSciFiMixture">
<D type="density" value="10.17" unit="g/cm3"/>
<composite n="0.950" ref="W"/>
<composite n="0.009" ref="H"/>
<composite n="0.041" ref="C"/>
</material>
<material name="WSciFi_UCLA_Abs">
<D type="density" value="12.4" unit="g / cm3"/>
<fraction n="0.96" ref="W"/>
<fraction n="0.04" ref="Polystyrene"/>
</material>
</materials>
<lccdd>
<comment> MRICH </comment>
<define>
<constant name="MRICH_rmin" value="15*cm"/>
<constant name="MRICH_rmax" value="BarrelTracking_rmax"/>
<constant name="MRICH_length" value="BackwardCherenkov_length"/>
<constant name="MRICH_zmin" value="-BarrelTracking_length/2.0 - BackwardTracking_length"/>
<constant name="MRICH_rmin" value="10*cm"/>
<constant name="MRICH_rmax" value="BackwardPIDRegion_rmax"/>
<constant name="MRICH_length" value="BackwardPIDRegion_length"/>
<constant name="MRICH_zmin" value="BackwardPIDRegion_zmin"/>
<constant name="MRICHAerogel_thickness" value="30.0*mm"/>
<constant name="MRICHAerogel_width" value="126.5*mm"/>
<constant name="MRICHFoam_thickness" value="2*mm"/>
......@@ -14,8 +15,13 @@
<constant name="MRICHPhotoDet_thickness" value="1.5*mm"/>
<constant name="MRICHPhotoDet_length" value="48.5*mm"/>
<constant name="MRICHGlassWindow_width" value="103.5*mm"/>
<constant name="MRICHGlassPhotoDet_thickness" value="2.0*mm"/>
<constant name="MRICHRearExtraSpace_thickness" value="10.0*mm"/>
<constant name="MRICHGlassPhotoDet_thickness" value="5.0*mm"/>
<constant name="MRICHPhotoDetMCPlate_thickness" value="0.3*mm"/> <!-- FIXME: should be 1.2*mm with PyrexGlass25 -->
<constant name="MRICHPhotoDetAnode_thickness" value="3.8*mm"/>
<constant name="MRICHPhotoDetPCB_thickness" value="2.0*mm"/>
<constant name="MRICHPhotoDetCopper_thickness" value="0.1*mm"/>
<constant name="MRICHPhotoDetKapton_thickness" value="0.2*mm"/>
<constant name="MRICHRearExtraSpace_thickness" value="0.3*mm"/>
<constant name="MRICHLensPhotoDet_length" value="136.4*mm"/>
<constant name="MRICHMirror_thickness" value="2.0*mm"/>
<constant name="MRICHMirror_length" value="MRICHLensPhotoDet_length - MRICHLensMirrorGap_thickness"/>
......@@ -39,6 +45,11 @@
+ MRICHFresnelLens_thickness
+ MRICHLensPhotoDet_length
+ MRICHGlassPhotoDet_thickness
+ 2.0*MRICHPhotoDetMCPlate_thickness
+ MRICHPhotoDetAnode_thickness
+ MRICHPhotoDetPCB_thickness
+ MRICHPhotoDetCopper_thickness
+ MRICHPhotoDetKapton_thickness
+ MRICHRearExtraSpace_thickness "/>
</define>
......@@ -54,35 +65,94 @@
<detectors>
<detector id="MRICH_ID" name="MRICH" type="athena_MRICH"
readout="MRICHHits"
projective="false"
vis="InvisibleWithDaughters" material="Air">
<dimensions rmin="MRICH_rmin" rmax="MRICH_rmax" length="MRICHCarbonFrame_length" zmin="-1.0*abs(MRICH_zmin)"/>
<module name="MRICH_module1" vis="InvisibleWithDaughters"
width="MRICHCarbonFrame_width"
height="MRICHCarbonFrame_width"
length="MRICHCarbonFrame_length">
<frame vis="AnlGray" thickness="MRICHCarbonFrame_thickness" material="CarbonFiber"/>
<aerogel vis="AnlTeal" length="MRICHAerogel_thickness" width="MRICHAerogel_width" material="AerogelOptical">
<frame vis="AnlGold_1" thickness="MRICHFoam_thickness" material="PolystyreneFoam" />
</aerogel>
<lens vis="AnlViolet" thickness="MRICHFresnelLens_thickness"
pitch="MRICHFresnelLensGroove_pitch" focal_length="6.0*inch"
effective_diameter="MRICHFresnelLensEffectiveDiameter"
width="MRICHAerogel_width"
material="AcrylicOptical"/>
<mirror vis="AnlGray" x1="MRICHMirror_width1" x2="MRICHMirror_width2" length="MRICHMirror_length"
surface="MRICH_MirrorOpticalSurface" thickness="MRICHMirror_thickness" material="AluminumOxide"/>
<photodet width="MRICHGlassWindow_width" thickness="MRICHGlassPhotoDet_thickness" material="PyrexGlassOptical">
<sensor nx="2" ny="2" thickness="MRICHPhotoDet_thickness" width="MRICHPhotoDet_length" material="SiliconOxide"/>
</photodet>
</module>
<!--
<layer id="1">
<array nx="4" ny="4" module="MRICH_module1">
<position x="0" y="0" z="0"/>
</array>
</layer>
-->
reflect="true"
projective="true"
vis="InvisibleWithDaughters"
material="Air">
<dimensions rmin="MRICH_rmin" rmax="MRICH_rmax" length="abs(MRICH_length)" zmin="MRICH_zmin"/>
<module name="MRICH_module1"
vis="InvisibleWithDaughters"
width="MRICHCarbonFrame_width"
height="MRICHCarbonFrame_width"
length="MRICHCarbonFrame_length">
<frame vis="AnlGray" thickness="MRICHCarbonFrame_thickness" material="CarbonFiber"/>
<aerogel vis="MRICH_aerogel_vis"
length="MRICHAerogel_thickness"
width="MRICHAerogel_width"
material="AerogelOptical">
<frame vis="MRICH_frame_vis" thickness="MRICHFoam_thickness" material="PolystyreneFoam" />
</aerogel>
<lens vis="MRICH_lens_vis" thickness="MRICHFresnelLens_thickness"
pitch="MRICHFresnelLensGroove_pitch" focal_length="6.0*inch"
effective_diameter="MRICHFresnelLensEffectiveDiameter"
width="MRICHAerogel_width"
material="AcrylicOptical"/>
<mirror vis="MRICH_mirror_vis"
x1="MRICHMirror_width1" x2="MRICHMirror_width2" length="MRICHMirror_length"
surface="MRICH_MirrorOpticalSurface" thickness="MRICHMirror_thickness"
material="AluminumOxide"/>
<photodet width="MRICHGlassWindow_width" thickness="MRICHGlassPhotoDet_thickness" material="PyrexGlassOptical">
<sensor nx="2" ny="2" thickness="MRICHPhotoDet_thickness" width="MRICHPhotoDet_length" material="SiliconOxide"/>
<layer thickness="MRICHPhotoDetMCPlate_thickness" material="PyrexGlass"/> <!-- FIXME: should be PyrexGlass25 with 1.2*mm thickness -->
<layer thickness="MRICHPhotoDetMCPlate_thickness" material="PyrexGlass"/> <!-- FIXME: should be PyrexGlass25 with 1.2*mm thickness -->
<layer thickness="MRICHPhotoDetAnode_thickness" material="AluminumOxide"/>
<layer thickness="MRICHPhotoDetPCB_thickness" material="Fr4"/>
<layer thickness="MRICHPhotoDetCopper_thickness" material="Copper"/>
<layer thickness="MRICHPhotoDetKapton_thickness" material="Kapton"/>
</photodet>
</module>
<comment>
Modules are only listed here for one quadrant. The modules point back to the
z0 position (should be their actual position, but breaks parametrization).
</comment>
<positions scale="1.03" z0="157.5*cm">
<position x="-41.3250000000" y="178.6750000000"/>
<position x=" 96.0250000001" y="178.6750000000"/>
</positions>
<positions scale="1.036" z0="157.5*cm">
<position x="-41.3250000000" y="316.0250000001"/>
<position x=" 96.0250000001" y="316.0250000001"/>
<position x="233.3750000002" y="178.6750000000"/>
<position x="233.3750000002" y="316.0250000001"/>
</positions>
<positions scale="1.052" z0="157.5*cm">
<position x="-41.3250000000" y="453.3750000002"/>
<position x=" 96.0250000001" y="453.3750000002"/>
<position x="233.3750000002" y="453.3750000002"/>
<position x="370.7250000003" y="453.3750000002"/>
<position x="370.7250000003" y="316.0250000001"/>
<position x="370.7250000003" y="178.6750000000"/>
</positions>
<positions scale="1.078" z0="157.5*cm">
<position x="-41.3250000000" y="590.7250000003"/>
<position x=" 96.0250000001" y="590.7250000003"/>
<position x="233.3750000002" y="590.7250000003"/>
<position x="370.7250000003" y="590.7250000003"/>
<!--position x="508.0750000004" y="590.7250000003"/-->
<position x="508.0750000004" y="178.6750000000"/>
<position x="508.0750000004" y="316.0250000001"/>
<position x="508.0750000004" y="453.3750000002"/>
</positions>
<positions scale="1.09" z0="157.5*cm">
<position x="-41.3250000000" y="728.0750000004"/>
<position x=" 96.0250000001" y="728.0750000004"/>
<!--position x="233.3750000002" y="728.0750000004"/-->
<!--position x="370.7250000003" y="728.0750000004"/-->
<!--position x="508.0750000004" y="728.0750000004"/-->
<!--position x="645.4250000005" y="590.7250000003"/-->
<!--position x="645.4250000005" y="453.3750000002"/-->
<position x="645.4250000005" y="316.0250000001"/>
<position x="645.4250000005" y="178.6750000000"/>
</positions>
<positions scale="1.095" z0="157.5*cm">
<!--position x="-41.3250000000" y="865.4250000005"/-->
<!--position x=" 96.0250000001" y="865.4250000005"/-->
<!--position x="233.3750000002" y="865.4250000005"/-->
<!--position x="782.7750000006" y="178.6750000000"/-->
<!--position x="782.7750000006" y="316.0250000001"/-->
<!--position x="782.7750000006" y="453.3750000002"/-->
</positions>
<layer thickness="MRICHCarbonFrame_thickness" material="CarbonFiber"/>
</detector>
</detectors>
......@@ -92,5 +162,11 @@
<id>system:8,module:14,sensor:8,x:32:-16,y:-16</id>
</readout>
</readouts>
<!--Globals>
<Parameter Name="mrichInfo" Value="mrichmod/mrich_1_geoparams-0-0-4294967295-1527211159.xml"/>
</Globals-->
</lccdd>
......@@ -27,6 +27,317 @@
4.0*eV 10.0*cm
5.1*eV 10.0*cm
"/>
<matrix name="RINDEX__DIRCQuartz" coldim="2" values="
1.90745*eV 1.45653
1.93725*eV 1.45681
1.968*eV 1.4571
1.99974*eV 1.4574
2.03253*eV 1.45771
2.0664*eV 1.45804
2.10143*eV 1.45838
2.13766*eV 1.45873
2.17516*eV 1.45911
2.214*eV 1.4595
2.25426*eV 1.45991
2.296*eV 1.46034
2.33932*eV 1.4608
2.38431*eV 1.46128
2.43106*eV 1.46179
2.47968*eV 1.46233
2.53029*eV 1.4629
2.583*eV 1.4635
2.63796*eV 1.46415
2.69531*eV 1.46483
2.7552*eV 1.46557
2.81782*eV 1.46635
2.88335*eV 1.46719
2.952*eV 1.46809
3.024*eV 1.46907
3.0996*eV 1.47012
3.17908*eV 1.47125
3.26274*eV 1.47249
3.35092*eV 1.47383
3.44401*eV 1.47529
3.54241*eV 1.47689
3.64659*eV 1.47865
3.7571*eV 1.48059
3.87451*eV 1.48274
3.99949*eV 1.48513
4.13281*eV 1.48779
"/>
<matrix name="ABSLENGTH__DIRCQuartz" coldim="2" values="
1.90745*eV 2.33615e+06*mm
1.93725*eV 2.19567e+06*mm
1.968*eV 2.06162e+06*mm
1.99974*eV 1.93381e+06*mm
2.03253*eV 1.81203e+06*mm
2.0664*eV 1.6961e+06*mm
2.10143*eV 1.58582e+06*mm
2.13766*eV 1.48101e+06*mm
2.17516*eV 1.38148e+06*mm
2.214*eV 1.28706e+06*mm
2.25426*eV 1.19756e+06*mm
2.296*eV 1.11281e+06*mm
2.33932*eV 1.03264e+06*mm
2.38431*eV 956885*mm
2.43106*eV 885374*mm
2.47968*eV 817949*mm
2.53029*eV 754450*mm
2.583*eV 694722*mm
2.63796*eV 638613*mm
2.69531*eV 585973*mm
2.7552*eV 536656*mm
2.81782*eV 490520*mm
2.88335*eV 447425*mm
2.952*eV 407233*mm
3.024*eV 369813*mm
3.0996*eV 335032*mm
3.17908*eV 302764*mm
3.26274*eV 272886*mm
3.35092*eV 245275*mm
3.44401*eV 219814*mm
3.54241*eV 196390*mm
3.64659*eV 174889*mm
3.7571*eV 155204*mm
3.87451*eV 137229*mm
3.99949*eV 120863*mm
4.13281*eV 106006*mm
"/>
<matrix name="RINDEX__DIRCEpotek" coldim="2" values="
1.90745*eV 1.55403
1.93725*eV 1.55557
1.968*eV 1.55698
1.99974*eV 1.55827
2.03253*eV 1.55945
2.0664*eV 1.56056
2.10143*eV 1.5616
2.13766*eV 1.5626
2.17516*eV 1.56358
2.214*eV 1.56455
2.25426*eV 1.56553
2.296*eV 1.56654
2.33932*eV 1.5676
2.38431*eV 1.56872
2.43106*eV 1.56993
2.47968*eV 1.57125
2.53029*eV 1.57269
2.583*eV 1.57427
2.63796*eV 1.57601
2.69531*eV 1.57793
2.7552*eV 1.58005
2.81782*eV 1.58238
2.88335*eV 1.58495
2.952*eV 1.58777
3.024*eV 1.59086
3.0996*eV 1.59424
3.17908*eV 1.59793
3.26274*eV 1.60195
3.35092*eV 1.60631
3.44401*eV 1.61103
3.54241*eV 1.61614
3.64659*eV 1.62165
3.7571*eV 1.62758
3.87451*eV 1.63395
3.99949*eV 1.64077
4.13281*eV 1.64807
"/>
<matrix name="ABSLENGTH__DIRCEpotek" coldim="2" values="
1.90745*eV 2.33615e+06*mm
1.93725*eV 2.19567e+06*mm
1.968*eV 2.06162e+06*mm
1.99974*eV 1.93381e+06*mm
2.03253*eV 1.81203e+06*mm
2.0664*eV 1.6961e+06*mm
2.10143*eV 1.58582e+06*mm
2.13766*eV 1.48101e+06*mm
2.17516*eV 1.38148e+06*mm
2.214*eV 1.28706e+06*mm
2.25426*eV 1.19756e+06*mm
2.296*eV 1.11281e+06*mm
2.33932*eV 1.03264e+06*mm
2.38431*eV 956885*mm
2.43106*eV 885374*mm
2.47968*eV 817949*mm
2.53029*eV 754450*mm
2.583*eV 694722*mm
2.63796*eV 638613*mm
2.69531*eV 585973*mm
2.7552*eV 536656*mm
2.81782*eV 490520*mm
2.88335*eV 447425*mm
2.952*eV 407233*mm
3.024*eV 369813*mm
3.0996*eV 335032*mm
3.17908*eV 302764*mm
3.26274*eV 272886*mm
3.35092*eV 245275*mm
3.44401*eV 219814*mm
3.54241*eV 196390*mm
3.64659*eV 174889*mm
3.7571*eV 155204*mm
3.87451*eV 137229*mm
3.99949*eV 120863*mm
4.13281*eV 106006*mm
"/>
<matrix name="RINDEX__DIRCNlak33a" coldim="2" values="
1.0*eV 1.73816
1.2511*eV 1.73836
1.26386*eV 1.73858
1.27687*eV 1.73881
1.29016*eV 1.73904
1.30372*eV 1.73928
1.31758*eV 1.73952
1.33173*eV 1.73976
1.34619*eV 1.74001
1.36097*eV 1.74026
1.37607*eV 1.74052
1.39152*eV 1.74078
1.40731*eV 1.74105
1.42347*eV 1.74132
1.44*eV 1.7416
1.45692*eV 1.74189
1.47425*eV 1.74218
1.49199*eV 1.74249
1.51016*eV 1.74279
1.52878*eV 1.74311
1.54787*eV 1.74344
1.56744*eV 1.74378
1.58751*eV 1.74412
1.6081*eV 1.74448
1.62923*eV 1.74485
1.65092*eV 1.74522
1.6732*eV 1.74562
1.69609*eV 1.74602
1.71961*eV 1.74644
1.7438*eV 1.74687
1.76868*eV 1.74732
1.79427*eV 1.74779
1.82062*eV 1.74827
1.84775*eV 1.74878
1.87571*eV 1.7493
1.90452*eV 1.74985
1.93423*eV 1.75042
1.96488*eV 1.75101
1.99652*eV 1.75163
2.0292*eV 1.75228
2.06296*eV 1.75296
2.09787*eV 1.75368
2.13398*eV 1.75443
2.17135*eV 1.75521
2.21006*eV 1.75604
2.25017*eV 1.75692
2.29176*eV 1.75784
2.33492*eV 1.75882
2.37973*eV 1.75985
2.42631*eV 1.76095
2.47473*eV 1.76211
2.52514*eV 1.76335
2.57763*eV 1.76467
2.63236*eV 1.76608
2.68946*eV 1.76758
2.7491*eV 1.7692
2.81143*eV 1.77093
2.87666*eV 1.77279
2.94499*eV 1.7748
3.01665*eV 1.77698
3.09187*eV 1.77934
3.17095*eV 1.7819
3.25418*eV 1.7847
3.34189*eV 1.78775
3.43446*eV 1.79111
3.53231*eV 1.79481
3.6359*eV 1.79889
3.74575*eV 1.80343
3.86244*eV 1.8085
3.98663*eV 1.81419
4.11908*eV 1.82061
4.26062*eV 1.8279
4.41225*eV 1.83625
4.57506*eV 1.84589
4.75035*eV 1.85713
4.93961*eV 1.87039
"/>
<matrix name="ABSLENGTH__DIRCNlak33a" coldim="2" values="
1.0*eV 371813*mm
1.2511*eV 352095*mm
1.26386*eV 331021*mm
1.27687*eV 310814*mm
1.29016*eV 291458*mm
1.30372*eV 272937*mm
1.31758*eV 255238*mm
1.33173*eV 238342*mm
1.34619*eV 222234*mm
1.36097*eV 206897*mm
1.37607*eV 192313*mm
1.39152*eV 178463*mm
1.40731*eV 165331*mm
1.42347*eV 152896*mm
1.44*eV 141140*mm
1.45692*eV 130043*mm
1.47425*eV 119585*mm
1.49199*eV 109747*mm
1.51016*eV 100507*mm
1.52878*eV 91846.3*mm
1.54787*eV 83743.1*mm
1.56744*eV 76176.7*mm
1.58751*eV 69126.1*mm
1.6081*eV 62570.2*mm
1.62923*eV 56488.0*mm
1.65092*eV 50858.3*mm
1.6732*eV 45660.1*mm
1.69609*eV 40872.4*mm
1.71961*eV 36474.6*mm
1.7438*eV 32445.8*mm
1.76868*eV 28765.9*mm
1.79427*eV 25414.6*mm
1.82062*eV 22372.2*mm
1.84775*eV 19619.3*mm
1.87571*eV 17136.9*mm
1.90452*eV 14906.5*mm
1.93423*eV 12910.2*mm
1.96488*eV 11130.3*mm
1.99652*eV 9550.13*mm
2.0292*eV 8153.3*mm
2.06296*eV 6924.25*mm
2.09787*eV 5848.04*mm
2.13398*eV 4910.46*mm
2.17135*eV 4098.04*mm
2.21006*eV 3398.06*mm
2.25017*eV 2798.54*mm
2.29176*eV 2288.32*mm
2.33492*eV 1856.99*mm
2.37973*eV 1494.92*mm
2.42631*eV 1193.28*mm
2.47473*eV 943.973*mm
2.52514*eV 739.657*mm
2.57763*eV 573.715*mm
2.63236*eV 440.228*mm
2.68946*eV 333.94*mm
2.7491*eV 250.229*mm
2.81143*eV 185.064*mm
2.87666*eV 134.967*mm
2.94499*eV 96.9664*mm
3.01665*eV 68.5529*mm
3.09187*eV 47.6343*mm
3.17095*eV 32.4882*mm
3.25418*eV 21.7174*mm
3.34189*eV 14.2056*mm
3.43446*eV 9.07612*mm
3.53231*eV 5.65267*mm
3.6359*eV 3.4241*mm
3.74575*eV 2.01226*mm
3.86244*eV 1.14403*mm
3.98663*eV 0.62722*mm
4.11908*eV 0.330414*mm
4.26062*eV 0.166558*mm
4.41225*eV 0.0799649*mm
4.57506*eV 0.0363677*mm
4.75035*eV 0.0155708*mm
4.93961*eV 0.00623089*mm
"/>
<matrix name= "REFLECTIVITY_mirror" coldim="2" values="
1.0*eV 0.9
4.0*eV 0.9
......@@ -47,6 +358,494 @@
1240*eV/600 1.49
1240*eV/400 1.49
"/>
<matrix name="ABSLENGTH__Acrylic" coldim="2" values="
1240*eV/1100 2631*mm
1240*eV/1000 2631*mm
1240*eV/900 2631*mm
1240*eV/800 2631*mm
1240*eV/700 2500*mm
1240*eV/600 2272*mm
1240*eV/500 2000*mm
1240*eV/400 1315*mm
1240*eV/300 1613*mm
1240*eV/250 740*mm
1240*eV/225 125*mm
1240*eV/210 10*mm
1240*eV/200 0*mm
"/>
<!-- BEGIN dRICh and pfRICh material properties
- dumped from fun4all implementation
- see https://github.com/cisbani/dRICh/blob/main/share/source/g4dRIChOptics.hh
-->
<!-- dRICh gas -->
<matrix name="RINDEX__C2F6_DRICH" coldim="2" values="
1.7712*eV 1.00075
1.92389*eV 1.00075
2.10539*eV 1.00075
2.3247*eV 1.00075
2.59502*eV 1.00076
2.93647*eV 1.00076
3.38139*eV 1.00077
3.98521*eV 1.00078
4.85156*eV 1.0008
6.19921*eV 1.00084
"/>
<matrix name="ABSLENGTH__C2F6_DRICH" coldim="2" values="
1.7712*eV 10.0*m
1.92389*eV 10.0*m
2.10539*eV 10.0*m
2.3247*eV 10.0*m
2.59502*eV 10.0*m
2.93647*eV 10.0*m
3.38139*eV 10.0*m
3.98521*eV 10.0*m
4.85156*eV 10.0*m
6.19921*eV 10.0*m
"/>
<!-- pfRICh gas -->
<matrix name="RINDEX__C4F10_PFRICH" coldim="2" values="
1.7712*eV 1.0013
1.92389*eV 1.0013
2.10539*eV 1.00131
2.3247*eV 1.00131
2.59502*eV 1.00132
2.93647*eV 1.00133
3.38139*eV 1.00134
3.98521*eV 1.00136
4.85156*eV 1.0014
6.19921*eV 1.00149
"/>
<matrix name="ABSLENGTH__C4F10_PFRICH" coldim="2" values="
1.7712*eV 6.0*m
1.92389*eV 6.0*m
2.10539*eV 6.0*m
2.3247*eV 6.0*m
2.59502*eV 6.0*m
2.93647*eV 6.0*m
3.38139*eV 6.0*m
3.98521*eV 6.0*m
4.85156*eV 6.0*m
6.19921*eV 6.0*m
"/>
<!-- dRICh aerogel, for density=0.11 g/cm3 (OLD versions are for 0.1g/cm3) -->
<matrix name="RINDEX__Aerogel_DRICH_OLD" coldim="2" values="
1.87855*eV 1.01638
1.96673*eV 1.01642
2.0549*eV 1.01647
2.14308*eV 1.01652
2.23126*eV 1.01657
2.31943*eV 1.01662
2.40761*eV 1.01667
2.49579*eV 1.01673
2.58396*eV 1.01679
2.67214*eV 1.01685
2.76032*eV 1.01691
2.84849*eV 1.01698
2.93667*eV 1.01705
3.02485*eV 1.01712
3.11302*eV 1.01719
3.2012*eV 1.01727
3.28938*eV 1.01734
3.37755*eV 1.01742
3.46573*eV 1.01751
3.55391*eV 1.01759
3.64208*eV 1.01768
3.73026*eV 1.01777
3.81844*eV 1.01787
3.90661*eV 1.01796
3.99479*eV 1.01806
4.08297*eV 1.01816
4.17114*eV 1.01827
4.25932*eV 1.01838
4.3475*eV 1.01849
4.43567*eV 1.0186
4.52385*eV 1.01872
4.61203*eV 1.01884
4.7002*eV 1.01897
4.78838*eV 1.01909
4.87656*eV 1.01922
4.96473*eV 1.01936
5.05291*eV 1.0195
5.14109*eV 1.01964
5.22927*eV 1.01979
5.31744*eV 1.01994
5.40562*eV 1.02009
5.4938*eV 1.02025
5.58197*eV 1.02041
5.67015*eV 1.02057
5.75833*eV 1.02074
5.8465*eV 1.02092
5.93468*eV 1.0211
6.02286*eV 1.02128
6.11103*eV 1.02147
6.19921*eV 1.02167
"/>
<matrix name="RINDEX__Aerogel_DRICH" coldim="2" values="
1.87855*eV 1.01852
1.96673*eV 1.01856
2.05490*eV 1.01861
2.14308*eV 1.01866
2.23126*eV 1.01871
2.31943*eV 1.01876
2.40761*eV 1.01881
2.49579*eV 1.01887
2.58396*eV 1.01893
2.67214*eV 1.01899
2.76032*eV 1.01905
2.84849*eV 1.01912
2.93667*eV 1.01919
3.02485*eV 1.01926
3.11302*eV 1.01933
3.20120*eV 1.01941
3.28938*eV 1.01948
3.37755*eV 1.01956
3.46573*eV 1.01965
3.55391*eV 1.01973
3.64208*eV 1.01982
3.73026*eV 1.01991
3.81844*eV 1.02001
3.90661*eV 1.02010
3.99479*eV 1.02020
4.08297*eV 1.02030
4.17114*eV 1.02041
4.25932*eV 1.02052
4.34750*eV 1.02063
4.43567*eV 1.02074
4.52385*eV 1.02086
4.61203*eV 1.02098
4.70020*eV 1.02111
4.78838*eV 1.02123
4.87656*eV 1.02136
4.96473*eV 1.02150
5.05291*eV 1.02164
5.14109*eV 1.02178
5.22927*eV 1.02193
5.31744*eV 1.02208
5.40562*eV 1.02223
5.49380*eV 1.02239
5.58197*eV 1.02255
5.67015*eV 1.02271
5.75833*eV 1.02288
5.84650*eV 1.02306
5.93468*eV 1.02324
6.02286*eV 1.02342
6.11103*eV 1.02361
6.19921*eV 1.02381
"/>
<matrix name="ABSLENGTH__Aerogel_DRICH_OLD" coldim="2" values="
1.87855*eV 154*mm
1.96673*eV 156.17*mm
2.0549*eV 158.154*mm
2.14308*eV 159.974*mm
2.23126*eV 161.651*mm
2.31943*eV 163.2*mm
2.40761*eV 164.636*mm
2.49579*eV 165.97*mm
2.58396*eV 167.213*mm
2.67214*eV 168.374*mm
2.76032*eV 169.461*mm
2.84849*eV 170.481*mm
2.93667*eV 171.439*mm
3.02485*eV 172.342*mm
3.11302*eV 173.193*mm
3.2012*eV 173.998*mm
3.28938*eV 174.759*mm
3.37755*eV 175.481*mm
3.46573*eV 176.165*mm
3.55391*eV 176.817*mm
3.64208*eV 162.708*mm
3.73026*eV 140.953*mm
3.81844*eV 122.516*mm
3.90661*eV 106.833*mm
3.99479*eV 93.4428*mm
4.08297*eV 81.9694*mm
4.17114*eV 72.1072*mm
4.25932*eV 63.6011*mm
4.3475*eV 56.2434*mm
4.43567*eV 49.8599*mm
4.52385*eV 44.3054*mm
4.61203*eV 39.4601*mm
4.7002*eV 35.2211*mm
4.78838*eV 31.5049*mm
4.87656*eV 28.2374*mm
4.96473*eV 25.359*mm
5.05291*eV 22.8166*mm
5.14109*eV 20.5674*mm
5.22927*eV 18.5724*mm
5.31744*eV 16.7992*mm
5.40562*eV 15.2205*mm
5.4938*eV 13.8125*mm
5.58197*eV 12.5541*mm
5.67015*eV 11.4277*mm
5.75833*eV 10.4166*mm
5.8465*eV 9.50928*mm
5.93468*eV 8.69176*mm
6.02286*eV 7.95608*mm
6.11103*eV 7.29168*mm
6.19921*eV 6.69064*mm
"/>
<matrix name="ABSLENGTH__Aerogel_DRICH" coldim="2" values="
1.87855*eV 140.000*mm
1.96673*eV 141.973*mm
2.05490*eV 143.776*mm
2.14308*eV 145.431*mm
2.23126*eV 146.955*mm
2.31943*eV 148.364*mm
2.40761*eV 149.669*mm
2.49579*eV 150.882*mm
2.58396*eV 152.012*mm
2.67214*eV 153.067*mm
2.76032*eV 154.055*mm
2.84849*eV 154.982*mm
2.93667*eV 155.854*mm
3.02485*eV 156.674*mm
3.11302*eV 157.448*mm
3.20120*eV 158.180*mm
3.28938*eV 158.872*mm
3.37755*eV 159.528*mm
3.46573*eV 160.150*mm
3.55391*eV 160.742*mm
3.64208*eV 147.916*mm
3.73026*eV 128.139*mm
3.81844*eV 111.378*mm
3.90661*eV 97.121*mm
3.99479*eV 84.948*mm
4.08297*eV 74.518*mm
4.17114*eV 65.552*mm
4.25932*eV 57.819*mm
4.34750*eV 51.130*mm
4.43567*eV 45.327*mm
4.52385*eV 40.278*mm
4.61203*eV 35.873*mm
4.70020*eV 32.019*mm
4.78838*eV 28.641*mm
4.87656*eV 25.670*mm
4.96473*eV 23.054*mm
5.05291*eV 20.742*mm
5.14109*eV 18.698*mm
5.22927*eV 16.884*mm
5.31744*eV 15.272*mm
5.40562*eV 13.837*mm
5.49380*eV 12.557*mm
5.58197*eV 11.413*mm
5.67015*eV 10.389*mm
5.75833*eV 9.470*mm
5.84650*eV 8.645*mm
5.93468*eV 7.902*mm
6.02286*eV 7.233*mm
6.11103*eV 6.629*mm
6.19921*eV 6.082*mm
"/>
<matrix name="RAYLEIGH__Aerogel_DRICH_OLD" coldim="2" values="
1.87855*eV 309.218*mm
1.96673*eV 257.383*mm
2.0549*eV 215.968*mm
2.14308*eV 182.559*mm
2.23126*eV 155.367*mm
2.31943*eV 133.053*mm
2.40761*eV 114.606*mm
2.49579*eV 99.2482*mm
2.58396*eV 86.3795*mm
2.67214*eV 75.5291*mm
2.76032*eV 66.3313*mm
2.84849*eV 58.4918*mm
2.93667*eV 51.777*mm
3.02485*eV 45.998*mm
3.11302*eV 41.0045*mm
3.2012*eV 36.6696*mm
3.28938*eV 32.8931*mm
3.37755*eV 29.5904*mm
3.46573*eV 26.6917*mm
3.55391*eV 24.1402*mm
3.64208*eV 21.8856*mm
3.73026*eV 19.8884*mm
3.81844*eV 18.1144*mm
3.90661*eV 16.533*mm
3.99479*eV 15.1206*mm
4.08297*eV 13.856*mm
4.17114*eV 12.7208*mm
4.25932*eV 11.7005*mm
4.3475*eV 10.7791*mm
4.43567*eV 9.94752*mm
4.52385*eV 9.1938*mm
4.61203*eV 8.51136*mm
4.7002*eV 7.88964*mm
4.78838*eV 7.32468*mm
4.87656*eV 6.80988*mm
4.96473*eV 6.33864*mm
5.05291*eV 5.907*mm
5.14109*eV 5.51232*mm
5.22927*eV 5.14932*mm
5.31744*eV 4.81668*mm
5.40562*eV 4.51044*mm
5.4938*eV 4.22796*mm
5.58197*eV 3.9666*mm
5.67015*eV 3.72504*mm
5.75833*eV 3.50196*mm
5.8465*eV 3.29604*mm
5.93468*eV 3.10464*mm
6.02286*eV 2.92644*mm
6.11103*eV 2.76144*mm
6.19921*eV 2.607*mm
"/>
<matrix name="RAYLEIGH__Aerogel_DRICH" coldim="2" values="
1.87855*eV 281.107*mm
1.96673*eV 233.984*mm
2.05490*eV 196.334*mm
2.14308*eV 165.962*mm
2.23126*eV 141.242*mm
2.31943*eV 120.958*mm
2.40761*eV 104.188*mm
2.49579*eV 90.226*mm
2.58396*eV 78.527*mm
2.67214*eV 68.663*mm
2.76032*eV 60.301*mm
2.84849*eV 53.174*mm
2.93667*eV 47.070*mm
3.02485*eV 41.816*mm
3.11302*eV 37.277*mm
3.20120*eV 33.336*mm
3.28938*eV 29.903*mm
3.37755*eV 26.900*mm
3.46573*eV 24.265*mm
3.55391*eV 21.946*mm
3.64208*eV 19.896*mm
3.73026*eV 18.080*mm
3.81844*eV 16.468*mm
3.90661*eV 15.030*mm
3.99479*eV 13.746*mm
4.08297*eV 12.596*mm
4.17114*eV 11.564*mm
4.25932*eV 10.637*mm
4.34750*eV 9.799*mm
4.43567*eV 9.043*mm
4.52385*eV 8.358*mm
4.61203*eV 7.738*mm
4.70020*eV 7.172*mm
4.78838*eV 6.659*mm
4.87656*eV 6.191*mm
4.96473*eV 5.762*mm
5.05291*eV 5.370*mm
5.14109*eV 5.011*mm
5.22927*eV 4.681*mm
5.31744*eV 4.379*mm
5.40562*eV 4.100*mm
5.49380*eV 3.844*mm
5.58197*eV 3.606*mm
5.67015*eV 3.386*mm
5.75833*eV 3.184*mm
5.84650*eV 2.996*mm
5.93468*eV 2.822*mm
6.02286*eV 2.660*mm
6.11103*eV 2.510*mm
6.19921*eV 2.370*mm
"/>
<!-- dRICh acrylic -->
<matrix name="RINDEX__Acrylic_DRICH" coldim="2" values="
4.13281*eV 1.5017
4.22099*eV 1.5017
4.30916*eV 1.5017
4.39734*eV 1.5017
4.48552*eV 1.5017
4.57369*eV 1.5017
4.66187*eV 1.5017
4.75005*eV 1.5017
4.83822*eV 1.5017
4.9264*eV 1.5017
5.01458*eV 1.5017
5.10275*eV 1.5017
5.19093*eV 1.5017
5.27911*eV 1.5017
5.36728*eV 1.5017
5.45546*eV 1.5017
5.54364*eV 1.5017
5.63181*eV 1.5017
5.71999*eV 1.5017
5.80817*eV 1.5017
5.89634*eV 1.5017
5.98452*eV 1.5017
6.0727*eV 1.5017
6.16087*eV 1.5017
6.24905*eV 1.5017
6.33723*eV 1.5017
6.4254*eV 1.5017
6.51358*eV 1.5017
6.60176*eV 1.5017
6.68993*eV 1.5017
6.77811*eV 1.5017
6.86629*eV 1.5017
6.95446*eV 1.5017
7.04264*eV 1.5017
7.13082*eV 1.5017
7.21899*eV 1.5017
7.30717*eV 1.5017
7.39535*eV 1.5017
7.48353*eV 1.5017
7.5717*eV 1.5017
7.65988*eV 1.5017
7.74806*eV 1.5017
7.83623*eV 1.5017
7.92441*eV 1.5017
8.01259*eV 1.5017
8.10076*eV 1.5017
8.18894*eV 1.5017
8.27712*eV 1.5017
8.36529*eV 1.5017
8.45347*eV 1.5017
"/>
<matrix name="ABSLENGTH__Acrylic_DRICH" coldim="2" values="
4.13281*eV 82.0704*mm
4.22099*eV 36.9138*mm
4.30916*eV 13.3325*mm
4.39734*eV 5.03627*mm
4.48552*eV 2.3393*mm
4.57369*eV 1.36177*mm
4.66187*eV 0.933192*mm
4.75005*eV 0.708268*mm
4.83822*eV 0.573082*mm
4.9264*eV 0.483641*mm
5.01458*eV 0.420282*mm
5.10275*eV 0.373102*mm
5.19093*eV 0.33662*mm
5.27911*eV 0.307572*mm
5.36728*eV 0.283902*mm
5.45546*eV 0.264235*mm
5.54364*eV 0.247641*mm
5.63181*eV 0.233453*mm
5.71999*eV 0.221177*mm
5.80817*eV 0.210456*mm
5.89634*eV 0.201012*mm
5.98452*eV 0.192627*mm
6.0727*eV 0.185134*mm
6.16087*eV 0.178399*mm
6.24905*eV 0.172309*mm
6.33723*eV 0.166779*mm
6.4254*eV 0.166779*mm
6.51358*eV 0.166779*mm
6.60176*eV 0.166779*mm
6.68993*eV 0.166779*mm
6.77811*eV 0.166779*mm
6.86629*eV 0.166779*mm
6.95446*eV 0.166779*mm
7.04264*eV 0.166779*mm
7.13082*eV 0.166779*mm
7.21899*eV 0.166779*mm
7.30717*eV 0.166779*mm
7.39535*eV 0.166779*mm
7.48353*eV 0.166779*mm
7.5717*eV 0.166779*mm
7.65988*eV 0.166779*mm
7.74806*eV 0.166779*mm
7.83623*eV 0.166779*mm
7.92441*eV 0.166779*mm
8.01259*eV 0.166779*mm
8.10076*eV 0.166779*mm
8.18894*eV 0.166779*mm
8.27712*eV 0.166779*mm
8.36529*eV 0.166779*mm
8.45347*eV 0.166779*mm
"/>
<!-- END dRICh optical properties -->
</properties>
<materials>
<material name="AirOptical">
......@@ -57,6 +856,14 @@
<property name="RINDEX" ref="RINDEX__Air"/>
<property name="ABSLENGTH" coldim="2" values="1*eV 200*m 5*eV 200*m"/>
</material>
<material name="VacuumOptical">
<D type="density" unit="g/cm3" value="0.0000000001"/>
<fraction n="0.754" ref="N"/>
<fraction n="0.234" ref="O"/>
<fraction n="0.012" ref="Ar"/>
<property name="RINDEX" ref="RINDEX__Vacuum"/>
<property name="ABSLENGTH" coldim="2" values="1*eV 2000*m 5*eV 2000*m"/>
</material>
<material name="N2cherenkov">
<D type="density" value="0.00125" unit="g/cm3"/>
<composite n="1" ref="N"/>
......@@ -85,18 +892,88 @@
<composite n="2" ref="O"/>
<composite n="8" ref="H"/>
<property name="RINDEX" ref="RINDEX__Acrylic"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Acrylic"/>
</material>
<!-- BEGIN dRICh and pfRICh material definitions -->
<material name="C2F6_DRICH">
<D type="density" value="0.005734" unit="g/cm3"/>
<composite n="2" ref="C"/>
<composite n="6" ref="F"/>
<property name="RINDEX" ref="RINDEX__C2F6_DRICH"/>
<property name="ABSLENGTH" ref="ABSLENGTH__C2F6_DRICH"/>
</material>
<material name="C4F10_PFRICH">
<D type="density" value="0.009935" unit="g/cm3"/>
<composite n="4" ref="C"/>
<composite n="10" ref="F"/>
<property name="RINDEX" ref="RINDEX__C4F10_PFRICH"/>
<property name="ABSLENGTH" ref="ABSLENGTH__C4F10_PFRICH"/>
</material>
<material name="Aerogel_DRICH_OLD">
<D type="density" value="0.100" unit="g/cm3"/>
<comment> n_air = [dens(Si02)-dens(aerogel)] / [dens(Si02)-dens(Air) ] </comment>
<fraction n=" (2.32-0.1) / (2.32-0.0012)" ref="Air"/>
<fraction n="1 - (2.32-0.1) / (2.32-0.0012)" ref="SiliconDioxide"/>
<property name="RINDEX" ref="RINDEX__Aerogel_DRICH_OLD"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Aerogel_DRICH_OLD"/>
<property name="RAYLEIGH" ref="RAYLEIGH__Aerogel_DRICH_OLD"/>
</material>
<material name="Aerogel_DRICH">
<D type="density" value="0.110" unit="g/cm3"/>
<comment> n_air = [dens(Si02)-dens(aerogel)] / [dens(Si02)-dens(Air) ] </comment>
<fraction n=" (2.32-0.11) / (2.32-0.0012)" ref="Air"/>
<fraction n="1 - (2.32-0.11) / (2.32-0.0012)" ref="SiliconDioxide"/>
<property name="RINDEX" ref="RINDEX__Aerogel_DRICH"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Aerogel_DRICH"/>
<property name="RAYLEIGH" ref="RAYLEIGH__Aerogel_DRICH"/>
</material>
<material name="Acrylic_DRICH">
<D type="density" value="1.19" unit="g/cm3"/>
<comment> TO BE IMPROVED </comment>
<fraction n="0.99" ref="Plexiglass"/>
<fraction n="0.01" ref="PolyvinylAcetate"/>
<property name="RINDEX" ref="RINDEX__Acrylic_DRICH"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Acrylic_DRICH"/>
</material>
<!-- END dRICh material definitions -->
<material name="QuartzOptical">
<D type="density" value="2.2" unit="g/cm3"/>
<composite n="1" ref="Si"/>
<composite n="2" ref="O"/>
<property name="RINDEX" ref="RINDEX__DIRCQuartz"/>
<property name="ABSLENGTH" ref="ABSLENGTH__DIRCQuartz"/>
</material>
<material name="EpotekOptical">
<D type="density" value="1.2" unit="g/cm3"/>
<composite n="5" ref="H"/>
<composite n="3" ref="C"/>
<composite n="2" ref="O"/>
<property name="RINDEX" ref="RINDEX__DIRCEpotek"/>
<property name="ABSLENGTH" ref="ABSLENGTH__DIRCEpotek"/>
<!-- G4Material *Epotek = new G4Material("Epotek", density = 1.2 * g / cm3, ncomponents = 3);
Epotek->AddElement(C, natoms = 3);
Epotek->AddElement(H, natoms = 5);
Epotek->AddElement(O, natoms = 2);-->
</material>
<material name="Nlak33aOptical">
<D type="density" value="4.220" unit="g/cm3"/>
<composite n="1" ref="Si"/>
<composite n="2" ref="O"/>
<property name="RINDEX" ref="RINDEX__DIRCNlak33a"/>
<property name="ABSLENGTH" ref="ABSLENGTH__DIRCNlak33a"/>
<!--
Nlak33aMaterial = new G4Material("Nlak33a", density = 4.220 * g / cm3, ncomponents = 2);
Nlak33aMaterial->AddElement(Si, natoms = 1);
Nlak33aMaterial->AddElement(O, natoms = 2);
-->
</material>
</materials>
<surfaces>
<comment> For the values of "finish", model and type, see TGeoOpticalSurface.h !
</comment>
<opticalsurface finish="polished" model="glisur" name="MirrorOpticalSurface" type="dielectric_metal" value="0">
<property name="REFLECTIVITY" ref="REFLECTIVITY_mirror"/>
<property name="RINDEX" coldim="2" values="1.034*eV 1.5 4.136*eV 1.5"/>
<!--<property name="EFFICIENCY" ref="EFFICIENCY0x8b77240"/>-->
</opticalsurface>
<!--
<opticalsurface name="mirror2" finish="polished" model="glisur" type="dielectric_dielectric">
<property name="REFLECTIVITY" coldim="2" values="1.034*eV 0.8 4.136*eV 0.9"/>
......@@ -135,5 +1012,70 @@
-->
</opticalsurface>
<!-- BEGIN dRICh surface definitions -->
<opticalsurface name="MirrorSurface_DRICH" model="unified" finish="polishedfrontpainted" type="dielectric_dielectric">
<property name="REFLECTIVITY" coldim="2" values="
2.04358*eV 0.867812
2.0664*eV 0.865156
2.09046*eV 0.863906
2.14023*eV 0.86375
2.16601*eV 0.864062
2.20587*eV 0.864531
2.23327*eV 0.864375
2.26137*eV 0.865625
2.31972*eV 0.865313
2.35005*eV 0.865
2.38116*eV 0.864844
2.41313*eV 0.863828
2.44598*eV 0.863516
2.47968*eV 0.863125
2.53081*eV 0.862187
2.58354*eV 0.861719
2.6194*eV 0.861328
2.69589*eV 0.861016
2.73515*eV 0.861094
2.79685*eV 0.861602
2.86139*eV 0.862305
2.95271*eV 0.86375
3.04884*eV 0.865586
3.12665*eV 0.867383
3.2393*eV 0.870059
3.39218*eV 0.874199
3.52508*eV 0.878105
3.66893*eV 0.88252
3.82396*eV 0.887617
3.99949*eV 0.893721
4.13281*eV 0.898184
4.27679*eV 0.902744
4.48244*eV 0.907837
4.65057*eV 0.9102
4.89476*eV 0.909316
5.02774*eV 0.906174
5.16816*eV 0.900422
5.31437*eV 0.891521
5.63821*eV 0.859954
5.90401*eV 0.820831
6.19921*eV 0.762502
"/>
</opticalsurface>
<opticalsurface name="SensorSurface_DRICH" model="glisur" finish="polished" type="dielectric_metal">
<property name="EFFICIENCY" coldim="2" values="
1*eV 1
4*eV 1
7*eV 1
"/>
<property name="IMAGINARYRINDEX" coldim="2" values="
1*eV 1.69
4*eV 1.69
7*eV 1.69
"/>
<property name="REALRINDEX" coldim="2" values="
1*eV 1.92
4*eV 1.92
7*eV 1.92
"/>
</opticalsurface>
<!-- END dRICh surface definitions -->
</surfaces>
</lccdd>
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<!-- vessel geometry -->
<constant name="PFRICH_Length" value="BackwardRICHRegion_length"/> <!-- vessel z-length -->
<constant name="PFRICH_zmin" value="-BackwardRICHRegion_zmin"/> <!-- vessel front -->
<constant name="PFRICH_zmax" value="PFRICH_zmin - PFRICH_Length"/> <!-- vessel back -->
<constant name="PFRICH_rmin0" value="BackwardPIDRegion_tan * BackwardRICHRegion_zmin"/> <!-- bore radius at vessel frontplane -->
<constant name="PFRICH_rmin1" value="BackwardPIDRegion_tan * (BackwardRICHRegion_zmin + BackwardRICHRegion_length)"/> <!-- bore radius at vessel backplane -->
<constant name="PFRICH_rmax" value="BackwardPIDRegion_rmax - 2*cm"/> <!-- vessel backplane radius -->
<constant name="PFRICH_wall_thickness" value="0.5*cm"/> <!-- thickness of radial walls -->
<constant name="PFRICH_window_thickness" value="0.1*cm"/> <!-- thickness of entrance and exit walls -->
<!-- additional parameters -->
<constant name="PFRICH_aerogel_thickness" value="3.0*cm"/> <!-- aerogel thickness -->
<constant name="PFRICH_sensor_active_size_default" value="24.0*mm"/> <!-- sensor side length (effective area) -->
<constant name="PFRICH_sensor_full_size_default" value="25.8*mm"/> <!-- sensor side length (full size, with enclosure) -->
<!-- lores values are used for global vizualizations _only_ -->
<constant name="PFRICH_sensor_active_size_lores" value="PFRICH_sensor_active_size_default * 4"/>
<constant name="PFRICH_sensor_full_size_lores" value="PFRICH_sensor_full_size_default * 4"/>
<!-- actual values used by the parametrization -->
<constant name="PFRICH_sensor_active_size" value="PFRICH_sensor_active_size_default"/>
<constant name="PFRICH_sensor_full_size" value="PFRICH_sensor_full_size_default"/>
<constant name="PFRICH_sensor_thickness" value="1.5*mm"/> <!-- sensor thickness -->
<constant name="PFRICH_sensor_dist" value="40*cm"/> <!-- distance between aerogel exit plane and sensor entrance plane -->
<constant name="PFRICH_num_px" value="8"/> <!-- number of pixels along one side of the sensor -->
<!-- debugging switches -->
<comment>
- `PFRICH_debug_optics`: 1 = all components become vacuum; test opticalphotons from IP
2 = all components become vacuum, except for `gasvol`, test charged particles from IP
0 = off
</comment>
<constant name="PFRICH_debug_optics" value="0"/>
</define>
<detectors>
<!-- /detectors/detector -->
<documentation level="10">
### pfRICH: Proximity Focusing RICH
</documentation>
<detector
id="BackwardRICH_ID"
name="PFRICH"
type="athena_PFRICH"
readout="PFRICHHits"
gas="C4F10_PFRICH"
material="Aluminum"
vis_vessel="DRICH_vessel_vis"
vis_gas="DRICH_gas_vis"
debug_optics="PFRICH_debug_optics"
>
<!-- /detectors/detector/dimensions -->
<documentation level="10">
#### Vessel
- dimensions:
- `zmin`: z-position of vessel front plane
- `length`: overall z-length of the full vessel
- `rmin0` and `rmin1`: bore radius at front plane and back plane, respectively
- `rmax0` and `rmax1`: outer radius of vessel, at front plane and back plane, respectively
- `wall_thickness`: thickness of radial walls
- `window_thickness`: thickness of entrance and exit disks
</documentation>
<dimensions
zmin="PFRICH_zmin"
zmax="PFRICH_zmax"
length="PFRICH_Length"
rmin0="PFRICH_rmin0"
rmin1="PFRICH_rmin1"
rmax0="PFRICH_rmax"
rmax1="PFRICH_rmax"
wall_thickness="PFRICH_wall_thickness"
window_thickness="PFRICH_window_thickness"
/>
<!-- /detectors/detector/radiator -->
<documentation level="10">
#### Radiator
- radiator is defined in a wedge of azimuthal space, composed of aerogel and a
filter; the filter is applied to the back of the aerogel, so that it separates
the aerogel and gas radiators
- dimensions:
- `frontplane`: front of the aerogel, w.r.t. front plane of the vessel envelope
- `rmin` and `rmax`: inner and outer radius (at the front plane; radial bounds are conical)
- `phiw`: azimuthal width of wedge
- `thickness`: radiator thickness, defined separately for aerogel and filter
- `pitch`: controls the angle of the radiator (0=vertical)
</documentation>
<radiator
frontplane="-PFRICH_window_thickness"
rmin="PFRICH_rmin0 + PFRICH_wall_thickness + 0.2*cm"
rmax="(PFRICH_rmax/PFRICH_zmax)*PFRICH_zmin + 8.0*cm"
phiw="60*degree"
pitch="0*degree"
>
<aerogel
material="Aerogel_DRICH"
vis="DRICH_aerogel_vis"
thickness="PFRICH_aerogel_thickness"
/>
<filter
material="Acrylic_DRICH"
vis="DRICH_filter_vis"
thickness="0.3*mm"
/>
</radiator>
<!-- /detectors/detector/sensors -->
<documentation level="10">
#### Sensors
</documentation>
<sensors>
<!-- /detectors/detector/sensors/module -->
<documentation level="10">
##### Sensor module
- based on Hamamatsu S13361-3050AE-08
- effective area: 24.0 x 24.0 mm
- enclosure size: 25.8 x 25.8 mm
- pixel size: 3x3 mm (8x8 matrix)
- thickness: 1.5 mm
- dimensions:
- `side`: side length of the square module
- `thickness`: thickness of the sensor module
- `gap`: provides room between the squares, to help prevent them from overlapping
- note: the values of `side` and `gap` will determine how many sensors there are, since the
sensor placement algorithm will try to place as many as it can in the specified region
</documentation>
<module
material="Silicon"
surface="SensorSurface_DRICH"
vis="PFRICH_sensor_vis"
side="PFRICH_sensor_active_size"
thickness="PFRICH_sensor_thickness"
gap="0.5*(PFRICH_sensor_full_size-PFRICH_sensor_active_size) + 0.5*mm"
/>
<!-- /detectors/detector/sensors/plane -->
<documentation level="10">
##### Sensor plane
- sensors will be placed on a plane
- plane dimensions:
- `sensordist`: distance between sensor plane active surface (e.g., photocathode) and aerogel backplane
- `rmin`: minimum radial position of a sensor's centroid
- `rmax`: maximum radial position of a sensor's centroid
</documentation>
<plane
sensordist="PFRICH_sensor_dist"
rmin="PFRICH_rmin1 + 2*cm"
rmax="PFRICH_rmax - 4*cm"
/>
<services>
<comment>
Material should be equivallent with 3x0.5cm Al, spread over the entire available distance.
This was confirmed by Silvia.
</comment>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
<component name="air" thickness="40*mm" vis="InvisibleNoDaughters" material="Air"/>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
<component name="air" thickness="40*mm" vis="InvisibleNoDaughters" material="Air"/>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
</services>
</sensors>
</detector>
</detectors>
<documentation level="10">
#### Readout
- segmentation: square matrix of pixels
- `grid_size_x,y`: size of each sensor, but note we must divide sensor size
by 1 less than the number of pixels, to account for fenceposting
- `offset_x,y`: specified such that the `x` and `y` indicators are unsigned
- indicators and `cellID` bits:
| indicator | offset | length |
|-----------|--------|--------|
| pfRICH ID | 0 | 8 |
| sensor | 8 | 12 |
| x pixel | 20 | 16 |
| y pixel | 36 | 16 |
</documentation>
<readouts>
<readout name="PFRICHHits">
<segmentation
type="CartesianGridXY"
grid_size_x="PFRICH_sensor_active_size/(PFRICH_num_px-1)"
grid_size_y="PFRICH_sensor_active_size/(PFRICH_num_px-1)"
offset_x="-PFRICH_sensor_active_size/2.0"
offset_y="-PFRICH_sensor_active_size/2.0"
/>
<id>system:8,module:12,x:20:16,y:16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="mrich.xml"/>
<include ref="drich.xml"/>
</lccdd>
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="drich.xml"/>
<include ref="pfrich.xml"/>
</lccdd>
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="drich.xml"/>
<include ref="pfrich.xml"/>
<!--include ref="tof_endcap.xml"/-->
<include ref="tof_barrel.xml"/>
</lccdd>
......@@ -11,11 +11,11 @@
<constant name="SolenoidBarrelInnerSecondMLILayerThickness" value=" 20.0 * mm"/>
<constant name="SolenoidBarrelInnerThermalShieldThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelInnerFirstMLILayerThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumVesselThickness" value=" 20.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumVesselThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelCoilFormerThickness" value=" 50.0 * mm"/>
<constant name="SolenoidBarrelCoilThickness" value="200.0 * mm"/>
<constant name="SolenoidBarrelCoilOverbindThickness" value="200.0 * mm"/>
<constant name="SolenoidBarrelCoilFormerThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelCoilThickness" value="210.0 * mm"/>
<constant name="SolenoidBarrelCoilOverbindThickness" value=" 50.0 * mm"/>
<constant name="SolenoidBarrelOuterHeliumThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelOuterHeliumVesselThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelOuterFirstMLILayerThickness" value=" 15.0 * mm"/>
......@@ -108,17 +108,16 @@
</define>
<detectors>
<comment>
---------------
Solenoid Magnet
<documentation level="6">
### Solenoid Magnet
---------------
The solenoid superconducting magnet and surrounding cryostat and
yoke endcaps.
</comment>
</documentation>
<detector
id="Solenoid_ID"
name="SolenoidCoilBarrel"
type="refdet_SolenoidCoil"
type="athena_SolenoidCoil"
insideTrackingVolume="false">
<position x="0" y="0" z="Solenoid_offset"/>
<layer id="1"
......@@ -131,7 +130,7 @@
inner_r="SolenoidBarrelInnerSecondMLILayer_rmin"
outer_z="SolenoidBarrelInnerSecondMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelInnerSecondMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidBarrelInnerSecondMLILayerThickness"/>
</layer>
<layer id="3"
inner_r="SolenoidBarrelInnerThermalShield_rmin"
......@@ -143,7 +142,7 @@
inner_r="SolenoidBarrelInnerFirstMLILayer_rmin"
outer_z="SolenoidBarrelInnerFirstMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelInnerFirstMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidBarrelInnerFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_r="SolenoidBarrelInnerHeliumVessel_rmin"
......@@ -163,19 +162,19 @@
inner_r="SolenoidBarrelCoilFormer_rmin"
outer_z="SolenoidBarrelCoilFormer_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelCoilFormerThickness"/>
<slice material="Aluminum" thickness="SolenoidBarrelCoilFormerThickness"/>
</layer>
<layer id="8"
inner_r="SolenoidBarrelCoil_rmin"
outer_z="SolenoidBarrelCoil_zmax"
vis="SolenoidCryostatVis">
<slice material="Copper" thickness="SolenoidBarrelCoilThickness"/>
<slice material="Aluminum" thickness="SolenoidBarrelCoilThickness"/>
</layer>
<layer id="9"
inner_r="SolenoidBarrelCoilOverbind_rmin"
outer_z="SolenoidBarrelCoilOverbind_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelCoilOverbindThickness"/>
<slice material="Aluminum" thickness="SolenoidBarrelCoilOverbindThickness"/>
</layer>
<layer id="10"
......@@ -195,7 +194,7 @@
inner_r="SolenoidBarrelOuterFirstMLILayer_rmin"
outer_z="SolenoidBarrelOuterFirstMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelOuterFirstMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidBarrelOuterFirstMLILayerThickness"/>
</layer>
<layer id="13"
inner_r="SolenoidBarrelOuterThermalShield_rmin"
......@@ -207,7 +206,7 @@
inner_r="SolenoidBarrelOuterSecondMLILayer_rmin"
outer_z="SolenoidBarrelOuterSecondMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelOuterSecondMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidBarrelOuterSecondMLILayerThickness"/>
</layer>
<layer id="15"
inner_r="SolenoidBarrelOuterVacuumVessel_rmin"
......@@ -229,7 +228,7 @@
inner_z="SolenoidEndcapCoilEndSupport_zmin"
inner_r="SolenoidEndcapCoilEndSupport_rmin"
outer_r="SolenoidEndcapCoilEndSupport_rmax">
<slice material="Steel235" thickness="SolenoidEndcapCoilEndSupportThickness"/>
<slice material="Aluminum" thickness="SolenoidEndcapCoilEndSupportThickness"/>
</layer>
<layer id="2"
inner_z="SolenoidEndcapHelium_zmin"
......@@ -247,7 +246,7 @@
inner_z="SolenoidEndcapFirstMLILayer_zmin"
inner_r="SolenoidEndcapFirstMLILayer_rmin"
outer_r="SolenoidEndcapFirstMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapFirstMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidEndcapFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_z="SolenoidEndcapThermalShield_zmin"
......@@ -259,7 +258,7 @@
inner_z="SolenoidEndcapSecondMLILayer_zmin"
inner_r="SolenoidEndcapSecondMLILayer_rmin"
outer_r="SolenoidEndcapSecondMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapSecondMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidEndcapSecondMLILayerThickness"/>
</layer>
<layer id="7"
inner_z="SolenoidEndcapVacuumVessel_zmin"
......@@ -281,7 +280,7 @@
inner_z="SolenoidEndcapCoilEndSupport_zmin"
inner_r="SolenoidEndcapCoilEndSupport_rmin"
outer_r="SolenoidEndcapCoilEndSupport_rmax">
<slice material="Steel235" thickness="SolenoidEndcapCoilEndSupportThickness"/>
<slice material="Aluminum" thickness="SolenoidEndcapCoilEndSupportThickness"/>
</layer>
<layer id="2"
inner_z="SolenoidEndcapHelium_zmin"
......@@ -299,7 +298,7 @@
inner_z="SolenoidEndcapFirstMLILayer_zmin"
inner_r="SolenoidEndcapFirstMLILayer_rmin"
outer_r="SolenoidEndcapFirstMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapFirstMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidEndcapFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_z="SolenoidEndcapThermalShield_zmin"
......@@ -311,7 +310,7 @@
inner_z="SolenoidEndcapSecondMLILayer_zmin"
inner_r="SolenoidEndcapSecondMLILayer_rmin"
outer_r="SolenoidEndcapSecondMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapSecondMLILayerThickness"/>
<slice material="MylarMLI" thickness="SolenoidEndcapSecondMLILayerThickness"/>
</layer>
<layer id="7"
inner_z="SolenoidEndcapVacuumVessel_zmin"
......@@ -328,12 +327,12 @@
zmin="-SolenoidBarrelCoil_zmax" zmax="SolenoidBarrelCoil_zmax"
inner_radius="SolenoidalFieldRadius" outer_radius="SolenoidalFieldRadius+20*cm"/-->
<field type="FieldMapBrBz" name="GlobalSolenoid" field_type="magnetic"
field_map="fieldmaps/EIC_Magnetic_Field_Map_2021_05_07_radial_coords__cm___T_.120000.lines.Bmap"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena/uploads/9179a98c5a80df1f2b98968c747c1c48/EIC_Magnetic_Field_Map_2021_05_07_radial_coords__cm___T_.120000.lines.Bmap"
field_map="fieldmaps/EIC_v.2.0.3_Magnetic_Field_Map_2021_09_28_radial_coords_cm_T.Bmap.txt"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena/uploads/fd64db33e1c99547ab1840d8847fd1c1/EIC_v.2.0.3_Magnetic_Field_Map_2021_09_28_radial_coords_cm_T.Bmap.txt"
scale="-1.0">
<dimensions>
<transverse step="2.0*cm" rmin="0*cm" rmax="398*cm" />
<longitudinal step="2.0*cm" zmin="-600*cm" zmax="598*cm" />
<transverse step="2.0*cm" rmin="0*cm" rmax="998*cm" />
<longitudinal step="2.0*cm" zmin="-800*cm" zmax="798*cm" />
<translation x="0.0*cm" y="0.0*cm" z="0.0*cm" />
<rotation x="0" y="0" z="0" />
</dimensions>
......
# subsytem_views
These compact detector files are copies of the main `athena.xml` compact description file, modified to only show certain subsystems.
To add or update a subsystem:
1. Copy `athena.xml` to this directory (`compact/subsystem_views`) with the appropriate name.
2. Remove all the undesired detectors from this xml file. **Do not modify any of the included xml files** -- only modify the top level xml file here.
3. Check that that you have the desired color scheme (eg. `colors.xml`) and display attributes (eg. `display_detailed.xml`) included.
These compact detector files should not be used for any real simulations or studies.
# subsytem_views
These compact detector files are copies of the main `athena.xml` compact description file, modified to only show certain subsystems.
To add or update a subsystem:
1. Copy `athena.xml` to this directory (`compact/subsystem_views`) with the appropriate name.
2. Remove all the undesired detectors from this xml file. **Do not modify any of the included xml files** -- only modify the top level xml file here.
3. Check that that you have the desired color scheme (eg. `colors.xml`) and display attributes (eg. `display_detailed.xml`) included.
These compact detector files should not be used for any real simulations or studies.
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>
<documentation level="-1">
# Athena Detector
- https://eicweb.phy.anl.gov/EIC/detectors/athena.git
- https://eicweb.phy.anl.gov/EIC/detectors/ip6.git
</documentation>
<!-- Some information about detector -->
<info name="Athena Detector" title="Athena Detector"
author="Athena Collaboration"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena.git"
status="development"
version="v1 2021-03-16">
<comment> Athena </comment>
</info>
<define>
<documentation level="2">
## Main Constant Definitions
The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="ip6/definitions.xml" />
<include ref="ip6/far_forward/fields_275.xml" />
<include ref="compact/definitions.xml" />
</define>
<includes>
<gdmlFile ref="compact/elements.xml"/>
<gdmlFile ref="compact/materials.xml"/>
<file ref="compact/optical_materials.xml"/>
</includes>
<limits>
<limitset name="EICBeamlineLimits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>
<display>
<include ref="compact/colors.xml" />
<include ref="compact/display_detailed.xml" />
</display>
<documentation level="0">
## Detector Subsystems
### IP Subsystems
The interaction point subsystems are included before the central detector subsystems.
This is becuase the IP subsystems, for exmaple the beampipe, will define paramters
which are subsquently used in the central detector construction -- e.g. the vertex tracker
uses the beampipe OD to help define its placement.
The IP subsystems include the Far forward and backward regions. The list of subsystem includes:
- Interaction region beampipe
- B0 tracker
- Off-momentum tracker
- Far forward roman pots
- Zero Degree Calorimeter
- Beam line magnets.
- and more...
</documentation>
<comment>
Beamline elements
-----------------
</comment>
<include ref="ip6/central_beampipe.xml" />
<documentation level="5">
## Main magnet
Note: When changing magnet, also select dimensions in definitions.xml.
</documentation>
<include ref="compact/solenoid.xml"/>
<documentation level="10">
## Central calorimetry
</documentation>
<include ref="compact/ecal.xml"/>
<include ref="compact/hcal.xml"/>
<readouts>
</readouts>
</lccdd>
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>
<documentation level="-1">
# Athena Detector
- https://eicweb.phy.anl.gov/EIC/detectors/athena.git
- https://eicweb.phy.anl.gov/EIC/detectors/ip6.git
</documentation>
<!-- Some information about detector -->
<info name="Athena Detector" title="Athena Detector"
author="Athena Collaboration"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena.git"
status="development"
version="v1 2021-03-16">
<comment> Athena </comment>
</info>
<define>
<documentation level="2">
## Main Constant Definitions
The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="ip6/definitions.xml" />
<include ref="ip6/far_forward/fields_275.xml" />
<include ref="compact/definitions.xml" />
</define>
<includes>
<gdmlFile ref="compact/elements.xml"/>
<gdmlFile ref="compact/materials.xml"/>
<file ref="compact/optical_materials.xml"/>
</includes>
<limits>
<limitset name="EICBeamlineLimits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>
<display>
<include ref="compact/colors.xml" />
<include ref="compact/display_detailed.xml" />
</display>
<documentation level="0">
## Detector Subsystems
### IP Subsystems
The interaction point subsystems are included before the central detector subsystems.
This is becuase the IP subsystems, for exmaple the beampipe, will define paramters
which are subsquently used in the central detector construction -- e.g. the vertex tracker
uses the beampipe OD to help define its placement.
The IP subsystems include the Far forward and backward regions. The list of subsystem includes:
- Interaction region beampipe
- B0 tracker
- Off-momentum tracker
- Far forward roman pots
- Zero Degree Calorimeter
- Beam line magnets.
- and more...
</documentation>
<!-- <include ref="ip6/forward_ion_beamline.xml"/>-->
<!-- <include ref="ip6/beampipe.xml"/>-->
<detectors>
<!-- <detector id="VertexBarrelSubAssembly_ID"-->
<!-- name="VertexBarrelSubAssembly"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TrackerSubAssemblyVis">-->
<!-- <composite name="VertexBarrel" />-->
<!-- </detector>-->
<!-- <detector id="VertexEndcapSubAssembly_ID"-->
<!-- name="VertexEndcapSubAssembly"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TrackerSubAssemblyVis">-->
<!-- <composite name="VertexEndcapN" />-->
<!-- <composite name="VertexEndcapP" />-->
<!-- </detector>-->
<!-- <detector id="TrackerBarrelSubAssembly_Inner_ID"-->
<!-- name="TrackerBarrelSubAssembly_Inner"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TrackerSubAssemblyVis">-->
<!-- <composite name="TrackerBarrel_Inner"/>-->
<!-- </detector>-->
<!-- <detector id="TrackerEndcapSubAssembly_Inner_ID"-->
<!-- name="TrackerEndcapSubAssembly_Inner"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TrackerSubAssemblyVis">-->
<!-- <composite name="TrackerEndcapN_Inner"/>-->
<!-- <composite name="TrackerEndcapP_Inner"/>-->
<!-- </detector>-->
<!-- <detector id="TrackerBarrelSubAssembly_Outer_ID"-->
<!-- name="TrackerSubAssembly_Outer"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TrackerSubAssemblyVis">-->
<!-- <composite name="TrackerBarrel_Outer"/>-->
<!-- <composite name="TrackerEndcapP_Outer"/>-->
<!-- <composite name="TrackerEndcapN_Outer"/>-->
<!-- </detector>-->
<!-- <detector id="TOFSubAssembly_ID"-->
<!-- name="TOFSubAssembly"-->
<!-- type="DD4hep_SubdetectorAssembly"-->
<!-- vis="TOFSubAssemblyVis">-->
<!-- <composite name="BarrelTOF"/>-->
<!-- <composite name="ForwardTOF"/>-->
<!-- <composite name="BackwardTOF"/>-->
<!-- </detector>-->
</detectors>
<include ref="compact/fake_dirc.xml"/>
<readouts>
</readouts>
</lccdd>
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>
<documentation level="-1">
# Athena Detector
- https://eicweb.phy.anl.gov/EIC/detectors/athena.git
- https://eicweb.phy.anl.gov/EIC/detectors/ip6.git
</documentation>
<!-- Some information about detector -->
<info name="Athena Detector" title="Athena Detector"
author="Athena Collaboration"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena.git"
status="development"
version="v1 2021-03-16">
<comment> Athena </comment>
</info>
<define>
<documentation level="2">
## Main Constant Definitions
The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="ip6/definitions.xml" />
<include ref="ip6/far_forward/fields_275.xml" />
<include ref="compact/definitions.xml" />
</define>
<includes>
<gdmlFile ref="compact/elements.xml"/>
<gdmlFile ref="compact/materials.xml"/>
<file ref="compact/optical_materials.xml"/>
</includes>
<limits>
<limitset name="EICBeamlineLimits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>
<display>
<include ref="compact/colors.xml" />
<include ref="compact/display_detailed.xml" />
</display>
<documentation level="10">
## dRICh only
</documentation>
<include ref="ip6/central_beampipe.xml" />
<include ref="compact/drich.xml" />
</lccdd>
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>
<documentation level="-1">
# Athena Detector
- https://eicweb.phy.anl.gov/EIC/detectors/athena.git
- https://eicweb.phy.anl.gov/EIC/detectors/ip6.git
</documentation>
<!-- Some information about detector -->
<info name="Athena Detector" title="Athena Detector"
author="Athena Collaboration"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena.git"
status="development"
version="v1 2021-03-16">
<comment> Athena </comment>
</info>
<define>
<documentation level="2">
## Main Constant Definitions
The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="ip6/definitions.xml" />
<include ref="ip6/far_forward/fields_275.xml" />
<include ref="compact/definitions.xml" />
</define>
<includes>
<gdmlFile ref="compact/elements.xml"/>
<gdmlFile ref="compact/materials.xml"/>
<file ref="compact/optical_materials.xml"/>
</includes>
<limits>
<limitset name="EICBeamlineLimits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>
<display>
<include ref="compact/colors.xml" />
<!--include ref="compact/colors2.xml"/-->
<include ref="compact/display.xml"/>
</display>
<documentation level="0">
## Detector Subsystems
### IP Subsystems
The interaction point subsystems are included before the central detector subsystems.
This is becuase the IP subsystems, for exmaple the beampipe, will define paramters
which are subsquently used in the central detector construction -- e.g. the vertex tracker
uses the beampipe OD to help define its placement.
The IP subsystems include the Far forward and backward regions. The list of subsystem includes:
- Interaction region beampipe
- B0 tracker
- Off-momentum tracker
- Far forward roman pots
- Zero Degree Calorimeter
- Beam line magnets.
- and more...
</documentation>
<comment>
Beamline elements
-----------------
</comment>
<include ref="ip6/central_beampipe.xml" />
<documentation level="5">
## Main magnet
Note: When changing magnet, also select dimensions in definitions.xml.
</documentation>
<documentation level="10">
## Central tracking detectors
</documentation>
<include ref="compact/tracking_config_deathvalley.xml"/>
<documentation level="10">
### PID detectors
</documentation>
<include ref="compact/pid_config_deathvalley.xml"/>
<readouts>
</readouts>
</lccdd>