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
Hide whitespace changes
Inline Side-by-side
Showing
with 2503 additions and 198 deletions
compact/ecal_backward_hybrid.xml 0 → 100644
View file @ 1008a184
<lccdd>
<define>
<comment>
Transition area.
The idea behind this parametrization is that:
one glass module with its wrap is always
a size of 4 crystal modules with its wraps.
Then the transition area (where glass meets crystals) has no gaps
+----------------+----------------+
| +----+ +----+ | +----------+ |
| | | | | | | | |
| +----+ +----+ | | | |
| +----+ +----+ | | | |
| | | | | | | | |
| +----+ +----+ | +----------+ |
+----------------+----------------+
crystal glass
This implies that:
GlassModule_wrap = 2*CrystalModule_wrap
GlassModule_sx = 2*CrystalModule_sx
GlassModule_sy = 2*CrystalModule_sy
</comment>
<constant name="CrystalModule_width" value="20.00*mm"/>
<constant name="CrystalModule_length" value="200.00*mm"/>
<constant name="CrystalModule_wrap" value="0.50*mm"/>
<constant name="CrystalModule_z0" value="0.0*cm"/>
<constant name="GlassModule_width" value="2*CrystalModule_width"/>
<constant name="GlassModule_length" value="55.00*cm"/>
<constant name="GlassModule_wrap" value="2*CrystalModule_wrap"/>
<constant name="GlassModule_z0" value="0.0*cm"/>
<constant name="EcalEndcapNIonCutout_dphi" value="30*degree"/>
<constant name="EcalEndcapN_thickness" value="GlassModule_length"/>
<constant name="EcalEndcapN_z0" value="-EcalEndcapN_zmin - EcalEndcapN_thickness/2"/>
<constant name="EcalEndcapNCrystal_rmax" value="40*cm"/>
<constant name="CrystalModule_distance" value="CrystalModule_width + CrystalModule_wrap"/>
<constant name="GlassModule_distance" value="GlassModule_width + GlassModule_wrap"/>
</define>
<display>
<vis name="HybridEcalOuterVis" alpha="0.5" r= "0.3" g="0.3" b="0.3" showDaughters="true" visible="true"/>
</display>
<detectors>
<documentation level="10">
#### Backwards Endcap EM Calorimeter
Backwards Endcap EM Calorimeter, placements generated by script
</documentation>
<detector id="ECalEndcapN_ID" name="EcalEndcapN" type="HybridCalorimeter" readout="EcalEndcapNHits">
<position x="0" y="0" z="EcalEndcapN_z0"/>
<rotation x="0" y="0" z="0"/>
</detector>
</detectors>
<readouts>
<comment>Effectively no segmentation, the segmentation is used to provide cell dimension info</comment>
<readout name="EcalEndcapNHits">
<segmentation type="MultiSegmentation" key="sector">
<segmentation name="CrystalSeg" key_value="1" type="CartesianGridXY" grid_size_x="CrystalModule_distance" grid_size_y="CrystalModule_distance"/>
<segmentation name="GlassSeg" key_value="2" type="CartesianGridXY" grid_size_x="GlassModule_distance" grid_size_y="GlassModule_distance"/>
</segmentation>
<id>system:8,sector:4,module:20,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
compact/ecal_barrel.xml
View file @ 1008a184
<lccdd>
<display>
<vis name="EcalBarrelEnvelope_vis" alpha="0.9" r="0.99" g="0.5" b="0" showDaughters="true" visible="false" />
<vis name="EcalBarrelStave_vis" alpha="0.9" r="0.99" g="0.5" b="0" showDaughters="true" visible="false" />
</display>
<display>
</display>
<define>
<comment>
---------------------------------------
......@@ -23,7 +21,7 @@
<constant name="EcalBarrel_ModRepeat" value="CaloSides"/>
<constant name="EcalBarrel_ModLength" value="0.5*m"/>
<constant name="EcalBarrel_ModWidth" value="0.5*m"/>
<constant name="EcalBarrel_AvailThickness" value="EcalBarrel_TotalThickness-EcalBarrel_Support_thickness"/>
<constant name="EcalBarrel_AvailThickness" value="EcalBarrelRegion_thickness-EcalBarrel_Support_thickness"/>
<constant name="EcalBarrel_ReadoutLayerThickness"
value="EcalBarrel_SiliconThickness
+ EcalBarrel_ElectronicsThickness
......@@ -47,38 +45,37 @@
<detectors>
<comment>
---------------------
Barrel EM Calorimeter
---------------------
<documentation>
### Barrel EM Calorimeter
A layered EM calorimeter with tungsten and silicon (AstroPix)
</comment>
</documentation>
<detector
id="ECalBarrel_ID"
name="EcalBarrel"
type="athena_EcalBarrel"
readout="EcalBarrelHits"
calorimeterType="EM_BARREL"
vis="EcalBarrelEnvelope_vis"
vis="EcalBarrelEnvelopeVis"
offset="EcalBarrel_offset">
<dimensions numsides="EcalBarrel_ModRepeat"
rmin="EcalBarrel_rmin"
z="EcalBarrel_length"/>
<staves vis="EcalBarrelStave_vis">
<support inside="true" material="Steel235" vis="AnlOrange"
<staves vis="EcalBarrelStaveVis">
<support inside="true" material="Steel235" vis="EcalBarrelSupportVis"
thickness="EcalBarrel_Support_thickness"
n_beams="3" grid_size="25.0*cm" >
</support>
</staves>
<layer repeat="min(EcalBarrel_ReadoutLayerNumber, EcalBarrelLayers_nMax)" vis="AnlBlue">
<slice material="Silicon" thickness="EcalBarrel_SiliconThickness" sensitive="yes" limits="cal_limits" vis="AnlGray"/>
<slice material="Silicon" thickness="EcalBarrel_ElectronicsThickness" vis="AnlGold"/>
<slice material="Copper" thickness="EcalBarrel_CopperThickness" vis="AnlGray"/>
<slice material="Kapton" thickness="EcalBarrel_KaptonThickness" vis="AnlGold"/>
<slice material="Epoxy" thickness="EcalBarrel_EpoxyThickness" vis="AnlGray"/>
<slice material="CarbonFiber" thickness="EcalBarrel_CarbonThickness" vis="AnlGold"/>
<slice material="TungstenDens24" thickness="EcalBarrel_TungstenThickness" vis="AnlGray"/>
<slice material="Air" thickness="EcalBarrel_LayerSpacing - EcalBarrel_TungstenThickness" vis="AnlGold"/>
<layer repeat="min(EcalBarrel_ReadoutLayerNumber, EcalBarrelLayers_nMax)" vis="EcalBarrelLayerVis">
<slice material="Silicon" thickness="EcalBarrel_SiliconThickness" sensitive="yes" limits="cal_limits" vis="EcalBarrelSliceVis"/>
<slice material="Silicon" thickness="EcalBarrel_ElectronicsThickness" vis="EcalBarrelSliceVis"/>
<slice material="Copper" thickness="EcalBarrel_CopperThickness" vis="EcalBarrelSliceVis"/>
<slice material="Kapton" thickness="EcalBarrel_KaptonThickness" vis="EcalBarrelSliceVis"/>
<slice material="Epoxy" thickness="EcalBarrel_EpoxyThickness" vis="EcalBarrelSliceVis"/>
<slice material="CarbonFiber" thickness="EcalBarrel_CarbonThickness" vis="EcalBarrelSliceVis"/>
<slice material="TungstenDens24" thickness="EcalBarrel_TungstenThickness" vis="EcalBarrelSliceVis"/>
<slice material="Air" thickness="EcalBarrel_LayerSpacing - EcalBarrel_TungstenThickness" vis="EcalBarrelSliceVis"/>
</layer>
</detector>
</detectors>
......
compact/ecal_barrel_interlayers.xml 0 → 100644
View file @ 1008a184
<lccdd>
<display>
</display>
<define>
<comment>
---------------------------------------
EM Calorimeter Parameters with AstroPix
---------------------------------------
</comment>
<comment>
To change the number of imaging layers from, e.g., 9 to 6,
change EcalBarrelImagingLayers_nMax variable
</comment>
<constant name="EcalBarrelImagingLayers_nMax" value="6"/>
<constant name="EcalBarrel_Support_thickness" value="5*cm"/>
<constant name="EcalBarrel_SiliconThickness" value="500*um"/>
<constant name="EcalBarrel_ElectronicsThickness" value="150*um"/>
<constant name="EcalBarrel_CopperThickness" value="100*um"/>
<constant name="EcalBarrel_KaptonThickness" value="200*um"/>
<constant name="EcalBarrel_EpoxyThickness" value="100*um"/>
<constant name="EcalBarrel_CarbonThickness" value="0.5*mm"/>
<constant name="EcalBarrel_CarbonSpacerWidth" value="4*mm"/>
<constant name="EcalBarrel_LayerSpacing" value="10.0*mm"/>
<constant name="EcalBarrel_FiberRadius" value="0.5*mm"/>
<constant name="EcalBarrel_FiberXSpacing" value="1.34*mm"/>
<constant name="EcalBarrel_FiberZSpacing" value="1.22*mm"/>
<constant name="EcalBarrel_SpaceBetween" value="0.1*mm"/>
<comment>
For Pb/SiFi (GlueX): X0 ~ 1.45 cm
For W/SiFi (sPHENIX): X0 ~ 0.7 cm (but different fiber orientation)
</comment>
<constant name="EcalBarrel_RadiatorThickness" value="EcalBarrel_FiberZSpacing*13"/>
<constant name="EcalBarrel_ModRepeat" value="CaloSides"/>
<constant name="EcalBarrel_ModLength" value="0.5*m"/>
<constant name="EcalBarrel_ModWidth" value="0.5*m"/>
<constant name="EcalBarrel_AvailThickness" value="EcalBarrelRegion_thickness-EcalBarrel_Support_thickness"/>
<constant name="EcalBarrel_ImagingLayerThickness"
value="EcalBarrel_SiliconThickness
+ EcalBarrel_ElectronicsThickness
+ EcalBarrel_CopperThickness
+ EcalBarrel_KaptonThickness
+ EcalBarrel_EpoxyThickness
+ EcalBarrel_CarbonThickness"/>
<constant name="EcalBarrel_ScFiLayerThickness"
value="EcalBarrel_RadiatorThickness
+ EcalBarrel_CarbonThickness
+ EcalBarrel_LayerSpacing"/>
<constant name="EcalBarrelImagingLayers_num"
value="min(EcalBarrelImagingLayers_nMax,
floor((EcalBarrel_AvailThickness-EcalBarrel_ImagingLayerThickness)/
(EcalBarrel_ImagingLayerThickness + EcalBarrel_ScFiLayerThickness +
EcalBarrel_SpaceBetween)))"/>
<comment>
Adjusting size of the ScFi back chunk to match number of imiaging layers
and 20 radiation lengths in total
</comment>
<constant name="EcalBarrel_FiberChunkLayers_num" value = "20-EcalBarrelImagingLayers_num+1"/>
<constant name="EcalBarrel_ImagingPartThickness"
value="(EcalBarrelImagingLayers_num-1)*(EcalBarrel_ImagingLayerThickness + EcalBarrel_ScFiLayerThickness + EcalBarrel_SpaceBetween)
+ EcalBarrel_ImagingLayerThickness + EcalBarrel_SpaceBetween"/>
<constant name="EcalBarrel_ScFiPartThickness_max"
value="max(0, EcalBarrel_AvailThickness-EcalBarrel_ImagingPartThickness)"/>
<constant name="EcalBarrel_ScFiPartThickness"
value="min(EcalBarrel_ScFiPartThickness_max,
EcalBarrel_FiberZSpacing*13*EcalBarrel_FiberChunkLayers_num)"/>
<constant name="EcalBarrel_SensitiveLayers_rmax"
value="EcalBarrel_rmin + EcalBarrel_ImagingPartThickness + EcalBarrel_ScFiPartThickness"/>
</define>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<detectors>
<comment>
---------------------------------------
Imaging Layers of Barrel EM Calorimeter
---------------------------------------
Silicon (Astropix) readout layers for imaging 3d showers
</comment>
<detector
id="ECalBarrel_ID"
name="EcalBarrelImaging"
type="athena_EcalBarrelInterlayers"
readout="EcalBarrelHits"
calorimeterType="EM_BARREL"
vis="EcalBarrelEnvelopeVis"
offset="EcalBarrel_offset">
<dimensions numsides="EcalBarrel_ModRepeat"
rmin="EcalBarrel_rmin"
z="EcalBarrel_length"/>
<layer repeat="EcalBarrelImagingLayers_num" vis="EcalBarrelLayerVis"
space_between="EcalBarrel_ScFiLayerThickness + EcalBarrel_SpaceBetween"
space_before="0.*cm">
<slice material="Silicon" thickness="EcalBarrel_SiliconThickness" sensitive="yes" limits="cal_limits" vis="EcalBarrelSliceVis"/>
<slice material="Silicon" thickness="EcalBarrel_ElectronicsThickness" vis="EcalBarrelSliceVis"/>
<slice material="Copper" thickness="EcalBarrel_CopperThickness" vis="EcalBarrelSliceVis"/>
<slice material="Kapton" thickness="EcalBarrel_KaptonThickness" vis="EcalBarrelSliceVis"/>
<slice material="Epoxy" thickness="EcalBarrel_EpoxyThickness" vis="EcalBarrelSliceVis"/>
<slice material="CarbonFiber" thickness="EcalBarrel_CarbonThickness" vis="EcalBarrelSliceVis"/>
</layer>
</detector>
<comment>
---------------------------------------
Pb/ScFi Layers of Barrel EM Calorimeter
---------------------------------------
</comment>
<detector
id="ECalBarrel2_ID"
name="EcalBarrelScFi"
type="athena_EcalBarrelInterlayers"
readout="EcalBarrelScFiHits"
calorimeterType="EM_BARREL"
vis="EcalBarrelEnvelopeVis"
offset="EcalBarrel_offset">
<dimensions numsides="EcalBarrel_ModRepeat"
rmin="EcalBarrel_rmin"
z="EcalBarrel_length"/>
<staves vis="EcalBarrelStaveVis">
<support material="Steel235" vis="EcalBarrelSupportVis" n_beams="3" grid_size="25.0*cm"
thickness="EcalBarrel_Support_thickness" beam_thickness="EcalBarrel_Support_thickness/4" >
</support>
</staves>
<layer repeat="EcalBarrelImagingLayers_num-1" vis="EcalBarrelLayerVis"
space_between="EcalBarrel_ImagingLayerThickness + EcalBarrel_SpaceBetween"
space_before="EcalBarrel_ImagingLayerThickness + EcalBarrel_SpaceBetween/2.">
<slice material="Lead" thickness="EcalBarrel_RadiatorThickness" vis="EcalBarrelFiberLayerVis">
<fiber material="PlasticScint"
sensitive="yes"
radius="EcalBarrel_FiberRadius"
spacing_x="EcalBarrel_FiberXSpacing"
spacing_z="EcalBarrel_FiberZSpacing"
vis="EcalBarrelFiberLayerVis">
</fiber>
</slice>
<slice material="CarbonFiber" thickness="EcalBarrel_CarbonThickness" vis="EcalBarrelSliceVis"/>
<slice material="Air" thickness="EcalBarrel_LayerSpacing " vis="EcalBarrelSliceVis"/>
</layer>
<layer repeat="EcalBarrel_FiberChunkLayers_num" vis="EcalBarrelLayerVis"
space_before="EcalBarrel_ImagingLayerThickness + EcalBarrel_SpaceBetween">
<slice material="Lead"
thickness="EcalBarrel_ScFiPartThickness/EcalBarrel_FiberChunkLayers_num"
vis="EcalBarrelFiberLayerVis">
<fiber material="PlasticScint"
sensitive="yes"
radius="EcalBarrel_FiberRadius"
spacing_x="EcalBarrel_FiberXSpacing"
spacing_z="EcalBarrel_FiberZSpacing"
vis="EcalBarrelFiberLayerVis">
</fiber>
</slice>
</layer>
</detector>
</detectors>
<readouts>
<readout name="EcalBarrelHits">
<segmentation type="CartesianGridXY" grid_size_x="0.5 * mm" grid_size_y="0.5 * mm"/>
<id>system:8,module:8,layer:8,slice:8,x:32:-16,y:-16</id>
</readout>
<readout name="EcalBarrelScFiHits">
<segmentation type="NoSegmentation"/>
<id>system:8,module:8,layer:8,slice:8,grid:16,fiber:16</id>
</readout>
</readouts>
</lccdd>
compact/ecal_forward_homogenous.xml 0 → 100644
View file @ 1008a184
<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>
compact/ecal_forward_scfi.xml 0 → 100644
View file @ 1008a184
<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>
compact/fake_dirc.xml 0 → 100644
View file @ 1008a184
<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>
compact/gem_tracker_endcap.xml
View file @ 1008a184
<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>
......
compact/hcal.xml
View file @ 1008a184
<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>
......
compact/helmholtz.xml
View file @ 1008a184
......@@ -12,7 +12,7 @@
<constant name="SolenoidBarrelInnerThermalShieldThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelInnerFirstMLILayerThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumVesselThickness" value=" 20.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumThickness" value=" 0.0 * mm"/> <!-- modification -->
<constant name="SolenoidBarrelInnerHeliumThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelCoilFormerThickness" value=" 50.0 * mm"/>
<constant name="SolenoidBarrelCoilThickness" value="210.0 * mm"/> <!-- modification -->
<constant name="SolenoidBarrelBlockThickness" value="210.0 * mm"/> <!-- addition -->
......@@ -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"
......@@ -329,9 +329,20 @@
</detectors>
<fields>
<field type="solenoid" name="GlobalSolenoid" inner_field="3.0 * tesla" outer_field="0.5*tesla"
<!--field type="solenoid" name="GlobalSolenoid" inner_field="3.0 * tesla" outer_field="0.5*tesla"
zmin="-SolenoidBarrelCoil_zmax" zmax="SolenoidBarrelCoil_zmax"
inner_radius="SolenoidalFieldRadius" outer_radius="SolenoidalFieldRadius+20*cm"/>
inner_radius="SolenoidalFieldRadius" outer_radius="SolenoidalFieldRadius+20*cm"/-->
<field type="FieldMapBrBz" name="GlobalSolenoid" field_type="magnetic"
field_map="fieldmaps/EIC_v.0.1.0_Magnetic_Field_Map_2021_05_28_radial_coords__cm___T_.401301.line.Bmap"
url="https://eicweb.phy.anl.gov/EIC/detectors/athena/uploads/2bb404a58b09e9b9fbd33befded5a04b/EIC_v.0.1.0_Magnetic_Field_Map_2021_05_28_radial_coords__cm___T_.401301.line.Bmap"
scale="-1.0">
<dimensions>
<transverse step="2.0*cm" rmin="0*cm" rmax="1000*cm" />
<longitudinal step="2.0*cm" zmin="-800*cm" zmax="800*cm" />
<translation x="0.0*cm" y="0.0*cm" z="0.0*cm" />
<rotation x="0" y="0" z="0" />
</dimensions>
</field>
</fields>
</lccdd>
compact/materials.xml
View file @ 1008a184
......@@ -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>
compact/mrich.xml 0 → 100644
View file @ 1008a184
<lccdd>
<comment> MRICH </comment>
<define>
<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"/>
<constant name="MRICHFresnelLens_thickness" value="0.06*inch"/>
<constant name="MRICHAerogelLensGap_thickness" value="2*mm"/>
<constant name="MRICHLensMirrorGap_thickness" value="2*mm"/>
<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="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"/>
<constant name="MRICHMirror_width1" value="MRICHAerogel_width "/>
<constant name="MRICHMirror_width2" value="MRICHGlassWindow_width"/>
<constant name="MRICHFresnelLensEffectiveDiameter" value="6.0*inch"/>
<constant name="MRICHFresnelLensGroove_pitch" value="inch/125"/>
<constant name="MRICHCarbonFrame_thickness" value="1.0*mm"/>
<constant name="MRICHCarbonFrame_width" value="MRICHAerogel_width+2.0*MRICHFoam_thickness + 2.0*MRICHCarbonFrame_thickness"/>
<constant name="MRICHModules_nx" value="floor((MRICH_rmax-MRICH_rmin)/MRICHCarbonFrame_width)"/>
<constant name="MRICHModules_ny" value="floor((MRICH_rmax-MRICH_rmin)/MRICHCarbonFrame_width)"/>
<constant name="MRICHCarbonFrame_length"
value="MRICHAerogel_thickness
+ 2.0*MRICHCarbonFrame_thickness
+ 2.0*MRICHFoam_thickness
+ MRICHAerogelLensGap_thickness
+ MRICHFresnelLens_thickness
+ MRICHLensPhotoDet_length
+ MRICHGlassPhotoDet_thickness
+ 2.0*MRICHPhotoDetMCPlate_thickness
+ MRICHPhotoDetAnode_thickness
+ MRICHPhotoDetPCB_thickness
+ MRICHPhotoDetCopper_thickness
+ MRICHPhotoDetKapton_thickness
+ MRICHRearExtraSpace_thickness "/>
</define>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<detectors>
<detector id="MRICH_ID" name="MRICH" type="athena_MRICH"
readout="MRICHHits"
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>
<readouts>
<readout name="MRICHHits">
<segmentation type="CartesianGridXY" grid_size_x="3*mm" grid_size_y="3*mm" />
<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>
compact/optical_materials.xml 0 → 100644
View file @ 1008a184
<lccdd>
<properties>
<matrix name="RINDEX__Vacuum" coldim="2" values="
1.0*eV 1.0
5.1*eV 1.0
"/>
<matrix name="RINDEX__Air" coldim="2" values="
1.0*eV 1.00029
5.1*eV 1.00029
"/>
<matrix name="RINDEX__Quartz" coldim="2" values="
1.0*eV 1.46
5.1*eV 1.46
"/>
<matrix name="RINDEX__N2" coldim="2" values="
1.0*eV 1.00033
4.0*eV 1.00033
5.1*eV 1.00033
"/>
<matrix name="RINDEX__Pyrex" coldim="2" values="
1.0*eV 1.5
4.0*eV 1.5
5.1*eV 1.5
"/>
<matrix name="ABSLENGTH__Pyrex" coldim="2" values="
1.0*eV 10.0*cm
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
5.1*eV 0.9
"/>
<matrix name="RINDEX__Aerogel" coldim="2" values="
1.0*eV 1.030
4.0*eV 1.030
5.1*eV 1.030
"/>
<matrix name="ABSLENGTH__Aerogel" coldim="2" values="
1.0*eV 4.0*cm
4.0*eV 4.0*cm
5.1*eV 4.0*cm
"/>
<matrix name="RINDEX__Acrylic" coldim="2" values="
1240*eV/1100 1.49
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">
<D type="density" unit="g/cm3" value="0.0012"/>
<fraction n="0.754" ref="N"/>
<fraction n="0.234" ref="O"/>
<fraction n="0.012" ref="Ar"/>
<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"/>
<property name="RINDEX" ref="RINDEX__N2"/>
</material>
<material name="PyrexGlassOptical">
<D type="density" value="2.23" 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"/>
<property name="RINDEX" ref="RINDEX__Pyrex"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Pyrex"/>
</material>
<material name="AerogelOptical">
<D value="0.2" unit="g / cm3"/>
<fraction n="0.625" ref="SiliconOxide"/>
<fraction n="0.374" ref="SiliconOxide"/>
<fraction n="0.1" ref="C"/>
<property name="RINDEX" ref="RINDEX__Aerogel"/>
<property name="ABSLENGTH" ref="ABSLENGTH__Aerogel"/>
</material>
<material name="AcrylicOptical">
<D type="density" value="1.18" unit="g/cm3"/>
<composite n="5" ref="C"/>
<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>
<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"/>
<property name="EFFICIENCY" coldim="2" values="2.034*eV 0.8 4.136*eV 1.0"/>
<property name="RINDEX" coldim="2" values="1.034*eV 1.5 4.136*eV 1.5"/>
</opticalsurface>
-->
<opticalsurface finish="polished" model="unified" name="MRICH_MirrorOpticalSurface" type="dielectric_metal" value="0">
</opticalsurface>
<opticalsurface finish="polished" model="unified" name="MRICH_LensOpticalSurface" type="dielectric_dielectric" value="0">
<property name="REFLECTIVITY" coldim="2" values="1240*eV/1100 0.08 1240*eV/400 0.08"/>
<!--
<property name="RINDEX" coldim="2" values="2.034*eV 1.56 4.136*eV 1.56"/>
<property name="SPECULARLOBECONSTANT" coldim="2" values="2.034*eV 0.3 4.136*eV 0.3 "/>
<property name="SPECULARSPIKECONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
<property name="BACKSCATTERCONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
-->
</opticalsurface>
<opticalsurface finish="polished" model="unified" name="MRICH_PhotoSensorOpticalSurface" type="dielectric_dielectric" value="0">
<!--
<property name="RINDEX" coldim="2" values="2.034*eV 1.35 4.136*eV 1.40"/>
<property name="SPECULARLOBECONSTANT" coldim="2" values="2.034*eV 0.3 4.136*eV 0.3 "/>
<property name="SPECULARSPIKECONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
<property name="BACKSCATTERCONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
-->
</opticalsurface>
<opticalsurface finish="polished" model="unified" name="MRICH_AerogelOpticalSurface" type="dielectric_dielectric" value="0">
<!--
<property name="RINDEX" coldim="2" values="2.034*eV 1.010 4.136*eV 1.010"/>
<property name="SPECULARLOBECONSTANT" coldim="2" values="2.034*eV 0.3 4.136*eV 0.3 "/>
<property name="SPECULARSPIKECONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
<property name="BACKSCATTERCONSTANT" coldim="2" values="2.034*eV 0.2 4.136*eV 0.2 "/>
-->
</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>
compact/pfrich.xml 0 → 100644
View file @ 1008a184
<?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>
compact/pid_config_acadia.xml 0 → 100644
View file @ 1008a184
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="mrich.xml"/>
<include ref="drich.xml"/>
</lccdd>
compact/pid_config_canyonlands.xml 0 → 100644
View file @ 1008a184
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="drich.xml"/>
<include ref="pfrich.xml"/>
</lccdd>
compact/pid_config_deathvalley.xml 0 → 100644
View file @ 1008a184
<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>
compact/solenoid.xml
View file @ 1008a184
......@@ -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"
......@@ -324,9 +323,20 @@
</detectors>
<fields>
<field type="solenoid" name="GlobalSolenoid" inner_field="3.0 * tesla" outer_field="0.5*tesla"
<!--field type="solenoid" name="GlobalSolenoid" inner_field="3.0 * tesla" outer_field="0.5*tesla"
zmin="-SolenoidBarrelCoil_zmax" zmax="SolenoidBarrelCoil_zmax"
inner_radius="SolenoidalFieldRadius" outer_radius="SolenoidalFieldRadius+20*cm"/>
inner_radius="SolenoidalFieldRadius" outer_radius="SolenoidalFieldRadius+20*cm"/-->
<field type="FieldMapBrBz" name="GlobalSolenoid" field_type="magnetic"
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="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>
</field>
</fields>
</lccdd>
compact/subsystem_views/README.md 0 → 100644
View file @ 1008a184
# 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.
compact/subsystem_views/README.mds 0 → 100644
View file @ 1008a184
# 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.
compact/subsystem_views/calorimeters.xml 0 → 100644
View file @ 1008a184
<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>