Skip to content
Snippets Groups Projects

Compare revisions

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

Source

Select target project
No results found

Target

Select target project
  • EIC/detectors/athena
  • zwzhao/athena
  • FernandoTA/athena
  • palspeic/athena
4 results
Show changes
Showing
with 2590 additions and 536 deletions
<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="EcalBarrel_ReadoutLayerNumber" 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>
......
<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>
<lccdd>
<comment>
Simplified homogenous implementation of the forward ScFi calorimeter
</comment>
<define>
<constant name="EcalEndcapP_FiberRadius" value="0.235*cm"/>
<constant name="EcalEndcapP_FiberOffset" value="0.5*mm"/>
<constant name="EcalEndcapP_FiberSpaceX" value="0.265*mm"/>
<constant name="EcalEndcapP_FiberSpaceY" value="0.425*mm"/>
</define>
<limits>
</limits>
<regions>
</regions>
<!-- Common Generic visualization attributes -->
<comment>Common Generic visualization attributes</comment>
<display>
</display>
<detectors>
<comment>
------------------------------------------
Forward (Positive Z) Endcap EM Calorimeter
------------------------------------------
An EM calorimeter with ScFi modules
</comment>
<detector id="ECalEndcapP_ID"
name="EcalEndcapP"
type="ScFiCalorimeter"
vis="EcalEndcapVis"
readout="EcalEndcapPHits">
<position x="0" y="0" z="EcalEndcapP_zmin + EcalEndcapP_length/2."/>
<dimensions rmin="EcalEndcapP_rmin" rmax="EcalEndcapP_rmax" length="EcalEndcapP_length"/>
<module sizex="25*mm" sizey="25*mm" sizez="170*mm" material="AvgTungstenScFi" vis="EcalEndcapBlockVis">
</module>
</detector>
</detectors>
<!-- Definition of the readout segmentation/definition -->
<readouts>
<readout name="EcalEndcapPHits">
<segmentation type="NoSegmentation"/>
<id>system:8,ring:8,module:20,fiber_x:8,fiber_y:8</id>
</readout>
</readouts>
<plugins>
</plugins>
</lccdd>
<lccdd>
<define>
<constant name="EcalEndcapP_FiberRadius" value="0.235*cm"/>
<constant name="EcalEndcapP_FiberOffset" value="0.5*mm"/>
<constant name="EcalEndcapP_FiberSpaceX" value="0.265*mm"/>
<constant name="EcalEndcapP_FiberSpaceY" value="0.425*mm"/>
</define>
<limits>
</limits>
<regions>
</regions>
<!-- Common Generic visualization attributes -->
<comment>Common Generic visualization attributes</comment>
<display>
</display>
<detectors>
<comment>
------------------------------------------
Forward (Positive Z) Endcap EM Calorimeter
------------------------------------------
An EM calorimeter with ScFi modules
</comment>
<detector id="ECalEndcapP_ID"
name="EcalEndcapP"
type="ScFiCalorimeter"
vis="EcalEndcapVis"
readout="EcalEndcapPHits">
<position x="0" y="0" z="EcalEndcapP_zmin + EcalEndcapP_length/2."/>
<dimensions rmin="EcalEndcapP_rmin" rmax="EcalEndcapP_rmax" length="EcalEndcapP_length"/>
<module sizex="25*mm" sizey="25*mm" sizez="170*mm" material="TungstenDens24" vis="EcalEndcapBlockVis">
<fiber material="Polystyrene"
radius="EcalEndcapP_FiberRadius"
offset="EcalEndcapP_FiberOffset"
spacex="EcalEndcapP_FiberSpaceX"
spacey="EcalEndcapP_FiberSpaceY"/>
</module>
</detector>
</detectors>
<!-- Definition of the readout segmentation/definition -->
<readouts>
<readout name="EcalEndcapPHits">
<segmentation type="NoSegmentation"/>
<id>system:8,ring:8,module:20,fiber_x:8,fiber_y:8</id>
</readout>
</readouts>
<plugins>
</plugins>
</lccdd>
<lccdd>
<define>
<constant name="DIRCModule_count" value="CaloSides" />
<constant name="DIRCBar_thickness" value="17*mm" />
<constant name="DIRCBar_length" value="DIRC_length-8*cm" />
<constant name="DIRCSkinFront_thickness" value="2.5*mm" />
<constant name="DIRCSkinBack_thickness" value="2.5*mm" />
<constant name="DIRCFrame_width" value="20*mm"/>
<constant name="DIRCFrame_length" value="DIRC_length"/>
<constant name="DIRCFrame_thickness" value="40*mm"/>
<constant name="DIRCModule_rmax" value="DIRC_rmin + DIRCBar_thickness + DIRCSkinFront_thickness + DIRCSkinBack_thickness"/>
<constant name="DIRCModule_thickness" value="DIRCModule_rmax - DIRC_rmin" />
<constant name="DIRCFrame_rmax" value="DIRC_rmin + DIRCFrame_thickness" />
<constant name="DIRCModule_halfangle" value="180*degree/DIRCModule_count" />
<constant name="DIRCModule_width" value="2*DIRC_rmin * tan(DIRCModule_halfangle)"/>
<constant name="DIRCModule_rmin" value="DIRC_rmin + 1*cm"/>
</define>
<materials>
</materials>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<detectors>
<detector id="BarrelDIRC_ID" name="FakeDIRC" type="athena_FakeDIRC" readout="DIRCBarHits" vis="DIRCVis">
<dimensions rmin="DIRC_rmin" rmax="DIRC_rmax" length="DIRC_length" />
<position x="0" y="0" z="DIRC_offset" />
<comment> Fake DIRC modules </comment>
<module name="DIRCModule" vis="DIRCModuleVis">
<module_component name="FrontSkin"
material="CarbonFiber"
sensitive="false"
width="DIRCModule_width"
thickness="DIRCSkinFront_thickness"
vis="DIRCSupportVis"
length="DIRCBar_length" />
<module_component name="QuartzBar"
material="Quartz"
sensitive="true"
width="DIRCModule_width"
thickness="DIRCBar_thickness"
vis="DIRCBarVis"
length="DIRCBar_length" />
<module_component name="BackSkin"
material="CarbonFiber"
sensitive="false"
width="DIRCModule_width"
thickness="DIRCSkinBack_thickness"
vis="DIRCSupportVis"
length="DIRCBar_length" />
<frame material="StainlessSteel"
width="DIRCFrame_width"
length="DIRCFrame_length"
vis="DIRCFrameVis"
thickness="DIRCFrame_thickness" />
</module>
<comment> Fake DIRC layers </comment>
<layer module="DIRCModule" id="1" vis="DIRCLayerVis">
<barrel_envelope
inner_r="DIRC_rmin"
outer_r="DIRC_rmax"
z_length="DIRC_length" />
<rphi_layout
phi_tilt="0"
nphi="DIRCModule_count"
phi0="0"
rc="0.5*(DIRCModule_rmin+DIRCModule_rmax)"
dr="0" />
<z_layout
dr="0.0*mm"
z0="0.0*mm"
nz="1" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="DIRCBarHits">
<segmentation type="CartesianGridXY" grid_size_x="3.0*mm" grid_size_y="3.0*mm" />
<id>system:8,layer:4,module:8,section:4,x:32:-16,y:-16</id>
</readout>
</readouts>
<plugins>
</plugins>
<fields>
</fields>
</lccdd>
<lccdd>
<comment>
//////////////////////////////////////////////////
// Far Forward Ion Zero Degree Calorimeter - Ecal
// PWO crystal bars
//////////////////////////////////////////////////
</comment>
<detectors>
<detector id="ffi_ZDC_ECAL_ID" name="ffi_ZDC_ECAL" type="ffi_ZDC_ECAL" readout="ffi_ZDC_ECAL_Hits" vis="ffi_ZDC_ECAL_Vis">
<position x="ffi_ZDC_ECAL_x_pos" y="ffi_ZDC_ECAL_y_pos" z="ffi_ZDC_ECAL_z_pos"/>
<rotation x="ffi_ZDC_ECAL_rotateX_angle" y="ffi_ZDC_ECAL_rotateY_angle" z="ffi_ZDC_ECAL_rotateZ_angle"/>
<dimensions x="ffi_ZDC_ECAL_width" z="ffi_ZDC_ECAL_thickness"/>
<module name="ffi_ZDC_ECAL_module"
vis="ffi_ZDC_ECAL_module_Vis"
material="PbWO4"
thickness="ffi_ZDC_ECAL_module_thickness"
width="ffi_ZDC_ECAL_module_width"
gap="ffi_ZDC_ECAL_module_gap"
ntower="ffi_ZDC_ECAL_module_ntower"/>
</detector>
</detectors>
<readouts>
<readout name="ffi_ZDC_ECAL_Hits">
<segmentation type="CartesianGridXY" grid_size_x="5.0*cm" grid_size_y="5.0*cm" />
<id>system:8,module:14,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<comment>
//////////////////////////////////////////////////
// Far Forward Ion Zero Degree Calorimeter - Hcal
// Si+W sampling
//////////////////////////////////////////////////
</comment>
<define>
<constant name="ffi_ZDC_HCAL_NRepeat" value="40"/>
<constant name="ffiZDCHcalTungstenThickness" value="3.5 * mm"/>
<constant name="ffiZDCHcalSiliconThickness" value="500.0 * um"/>
<constant name="ffiZDCHcalAirThickness" value="1.0 * mm"/>
</define>
<detectors>
<detector id="ffi_ZDC_HCAL_ID" name="ffi_ZDC_HCAL" type="ffi_ZDC_HCAL" readout="ffi_ZDC_HCAL_Hits" vis="ffi_ZDC_HCAL_Vis">
<position x="ffi_ZDC_HCAL_x_pos" y="ffi_ZDC_HCAL_y_pos" z="ffi_ZDC_HCAL_z_pos"/>
<rotation x="ffi_ZDC_HCAL_rotateX_angle" y="ffi_ZDC_HCAL_rotateY_angle" z="ffi_ZDC_HCAL_rotateZ_angle"/>
<dimensions x="ffi_ZDC_HCAL_width" z="ffi_ZDC_HCAL_length"/>
<layer repeat="ffi_ZDC_HCAL_NRepeat">
<slice name="Tungsten_slice" material="TungstenDens24" thickness="ffiZDCHcalTungstenThickness" vis="AnlBlue"/>
<slice name="Silicon_slice" material="Silicon" thickness="ffiZDCHcalSiliconThickness" vis="AnlRed" sensitive="true"/>
<slice name="Air_slice" material="Air" thickness="ffiZDCHcalAirThickness" vis="AnlGold"/>
</layer>
</detector>
</detectors>
<readouts>
<readout name="ffi_ZDC_HCAL_Hits">
<segmentation type="CartesianGridXY" grid_size_x="1.0*cm" grid_size_y="1.0*cm" />
<id>system:8,layer:12,slice:12,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<comment>Far Forward Hadron Detectors</comment>
<define>
</define>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<include ref="far_forward_ZDC_Ecal.xml" />
<include ref="far_forward_ZDC_Hcal.xml" />
<detectors>
</detectors>
<readouts>
</readouts>
<plugins>
</plugins>
<fields>
</fields>
</lccdd>
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<comment>
----------------------------------------
Forward Off-Momentum Tracker Parameters
----------------------------------------
</comment>
<constant name="ForwardOffMTracker_zpos" value="B1APF_CenterPosition + B1APF_Length/2.0 + 10.0*cm"/>
<constant name="ForwardOffMTracker_xpos" value="B1APF_XPosition+40.0*cm"/>
</define>
<detectors>
<detector
id="ForwardOffMTracker_ID"
name="ForwardOffMTracker"
readout="ForwardOffMTrackerHits"
type="ref_RectangularTracker"
insideTrackingVolume="true"
reflect="false"
vis="AnlRed">
<position x="ForwardOffMTracker_xpos" y="0" z="ForwardOffMTracker_zpos"/>
<layer id="1"
vis="AnlViolet"
x="60.0*cm"
y="20.0*cm" >
<slice material="Air" thickness="1.0*cm" vis="AnlViolet" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlViolet" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlViolet" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlViolet" sensitive="true" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="ForwardOffMTrackerHits">
<segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
<id>system:8,barrel:3,layer:4,slice:4,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<comment>
---------------------------------
Far Forward B0 Tracker Parameters
---------------------------------
</comment>
<constant name="ForwardRomanPotA_zpos" value="B1APF_CenterPosition + B1APF_Length/2.0 + 20.0*cm"/>
<constant name="ForwardRomanPotA_xpos" value="B1APF_XPosition"/>
<constant name="ForwardRomanPotB_zpos" value="B2PF_CenterPosition - B2PF_Length/2.0 - 20.0*cm"/>
<constant name="ForwardRomanPotB_xpos" value="B2PF_XPosition "/>
</define>
<detectors>
<detector
id="ForwardRomanPotA_ID"
name="ForwardRomanPotA"
readout="ForwardRomanPotHits"
type="ref_RectangularTracker"
insideTrackingVolume="true"
reflect="false"
vis="AnlRed">
<position x="ForwardRomanPotA_xpos" y="0" z="ForwardRomanPotA_zpos"/>
<layer id="1"
vis="AnlBlue"
x="20.0*cm"
y="10.0*cm" >
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
</layer>
</detector>
<detector
id="ForwardRomanPotB_ID"
name="ForwardRomanPotB"
readout="ForwardRomanPotHits"
type="ref_RectangularTracker"
insideTrackingVolume="true"
reflect="false"
vis="AnlRed">
<position x="ForwardRomanPotB_xpos" y="0" z="ForwardRomanPotB_zpos"/>
<layer id="1"
vis="AnlBlue"
x="20.0*cm"
y="10.0*cm" >
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
<slice material="Air" thickness="1.0*cm" vis="AnlBlue" sensitive="true" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="ForwardRomanPotHits">
<segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
<id>system:8,barrel:3,layer:4,slice:4,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<define>
</define>
<limits>
</limits>
<regions>
</regions>
<display>
</display>
<!-- Define detector -->
<detectors>
<detector
id="HCalBarrel_ID"
name="HcalBarrel"
type="DD4hep_PolyhedraBarrelCalorimeter2"
readout="HcalBarrelHits"
calorimeterType="HAD_BARREL"
gap="0.*cm"
material="Steel235">
<dimensions
numsides="CaloSides"
rmin="HcalBarrel_rmin"
z="HcalBarrelLength"/>
<staves vis="HcalBarrelVis"/>
<layer repeat="HcalBarrelLayers">
<slice material="Steel235" thickness="HcalSteelThickness"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="RPCGasDefault" thickness="HcalRPCGasThickness" sensitive="yes" limits="cal_limits"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="G10" thickness="HcalG10Thickness"/>
<slice material="Air" thickness="HcalAirThickness"/>
</layer>
</detector>
<detector
id="HCalEndcapP_ID"
name="HcalEndcapP"
type="refdet_PolyhedraEndcapCalorimeter2"
readout="HcalEndcapHits"
vis="HcalEndcapVis"
calorimeterType="HAD_ENDCAP"
reflect="false">
<dimensions
numsides="CaloSides"
zmin="HcalEndcap_zmin"
rmin="HcalEndcapP_rmin"
rmax="HcalBarrel_rmax"/>
<layer repeat="HcalEndcapLayers" vis="HcalEndcapVis" >
<slice material="Steel235" thickness="HcalSteelThickness"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="RPCGasDefault" thickness="HcalRPCGasThickness" sensitive="yes" limits="cal_limits"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="G10" thickness="HcalG10Thickness"/>
<slice material="Air" thickness="HcalAirThickness"/>
</layer>
</detector>
<detector
id="HCalEndcapN_ID"
name="HcalEndcapN"
type="refdet_PolyhedraEndcapCalorimeter2"
readout="HcalEndcapHits"
vis="HcalEndcapVis"
calorimeterType="HAD_ENDCAP" reflect="true">
<dimensions
numsides="CaloSides"
zmin="HcalEndcap_zmin"
rmin="HcalEndcapN_rmin"
rmax="HcalBarrel_rmax"/>
<layer repeat="HcalEndcapLayers" vis="HcalEndcapVis" >
<slice material="Steel235" thickness="HcalSteelThickness"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="RPCGasDefault" thickness="HcalRPCGasThickness" sensitive="yes" limits="cal_limits"/>
<slice material="PyrexGlass" thickness="HcalPyrexThickness"/>
<slice material="G10" thickness="HcalG10Thickness"/>
<slice material="Air" thickness="HcalAirThickness"/>
</layer>
</detector>
</detectors>
<!-- Definition of the readout segmentation/definition -->
<readouts>
</readouts>
<plugins>
</plugins>
<fields>
</fields>
</lccdd>
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<constant name="ForwardRICH_zmin" value="BarrelTracking_length/2.0 + ForwardTracking_length "/>
<constant name="ForwardRICH_rmin" value="ForwardPID_rmin1"/>
<constant name="ForwardRICH_rmax0" value="TrackerBarrel_rmax"/>
<constant name="ForwardRICH_rmax1" value="EcalBarrel_rmin"/>
<constant name="ForwardRICH_rmax2" value="Solenoid_rmin-6*cm"/>
<constant name="ForwardRICHDepth" value="0.9*m"/>
</define>
<detectors>
<detector
id="ForwardRICH_ID"
name="ForwardRICH"
type="refdet_ForwardRICH"
readout="ForwardRICHHits"
vis="RICHVis">
<dimensions
z0="ForwardRICH_zmin"
snout_length="ForwardRICH_length - ForwardRICHDepth"
length="ForwardRICH_length"
rmin="ForwardRICH_rmin"
rmax0="ForwardRICH_rmax0"
rmax1="ForwardRICH_rmax1"
rmax2="ForwardRICH_rmax2"/>
<radiator material="N2cherenkov" />
<tank
length="ForwardRICHDepth"
gas="N2cherenkov"
vis="GreenVis"
rmin="ForwardRICH_rmin"
rmax1="ForwardRICH_rmax1"
rmax2="ForwardRICH_rmax2" />
<comment> What are the following MCP-PMT parameters?</comment>
<mcppmt
z0="0.0*cm"
rmin="ForwardRICH_rmax1"
rmax="ForwardRICH_rmax2"
rtol="1.0*cm"
vis="BlueVis"
module_size="10*cm"
module_gap="0.2*cm"
thickness="1.0*cm"
material="Quartz" />
<mirror
z0="ForwardRICH_length - 40*cm"
thickness="1*mm"
material="PyrexGlass"
vis="GrayVis">
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="0*degree" />
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="60*degree" />
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="120*degree" />
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="180*degree" />
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="240*degree" />
<slice focus="10*cm" curve="300*cm" rmin="ForwardRICH_rmin" rmax="ForwardRICH_rmax1" phiw="59*degree" rotz="300*degree" />
</mirror>
</detector>
</detectors>
<readouts>
<readout name="ForwardRICHHits">
<segmentation type="CartesianGridXY" grid_size_x="3*mm" grid_size_y="3*mm" />
<id>system:8,layer:4,piece:4,module:14,x:32:-16,y:-16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<comment> Endcap Tracker </comment>
<comment> Endcap Tracker for Acadia </comment>
<define>
<constant name="GEMTrackerEndcap_zmin" value="BarrelTracking_zmax"/>
<constant name="GEMTrackerEndcap_NLayers" value="1"/>
<constant name="GEMTrackerEndcapFoilOpeningAngle" value="30.1*degree"/>
<constant name="GEMTrackerEndcapP_thickness" value="min(10.0*cm,ForwardTrackingRegion_length)" />
<constant name="GEMTrackerEndcapN_thickness" value="min(10.0*cm,BackwardTrackingRegion_length)" />
<constant name="GEMTrackerEndcapP_zmin" value="ForwardTrackingRegion_zmin"/>
<constant name="GEMTrackerEndcapN_zmin" value="BackwardTrackingRegion_zmin"/>
<constant name="GEMTrackerEndcapP_rmin" value="19*cm" />
<constant name="GEMTrackerEndcapN_rmin" value="10*cm" />
<constant name="GEMTrackerEndcapP_rmax" value="200*cm" />
<constant name="GEMTrackerEndcapN_rmax" value="BackwardPIDRegion_rmax" />
<constant name="GEMTrackerEndcapOpeningAngle" value="30.4*degree"/>
<constant name="GEMTrackerEndcap_dz" value="15*mm" />
<comment>
____X2____
\ /
\ / Y (=z for the trd)
\_X1_/
</comment>
<constant name="GEMTrackerEndcapFoilX1" value="43.0 * mm"/>
<constant name="GEMTrackerEndcapLayer_thickness" value="ForwardTrackingGEMLength/GEMTrackerEndcap_NLayers"/>
<constant name="GEMTrackerEndcapFoil_rmin" value="GEMTrackerEndcapFoilX1/(tan(GEMTrackerEndcapFoilOpeningAngle/2.0)*2.0)"/>
<constant name="GEMTrackerEndcapFoilY" value="BarrelTracking_rmax- GEMTrackerEndcapFoil_rmin- 1.0*cm"/>
<constant name="GEMTrackerEndcapFoilX2" value="2*GEMTrackerEndcapFoilY*tan(GEMTrackerEndcapFoilOpeningAngle/2.0) + GEMTrackerEndcapFoilX1"/>
<constant name="GEMTrackerEndcapFoil_rmax" value="(GEMTrackerEndcapFoilY+GEMTrackerEndcapFoil_rmin)/(cos(GEMTrackerEndcapFoilOpeningAngle/2.0))"/>
<constant name="GEMTrackerEndcapP_x1" value="2.0*GEMTrackerEndcapP_rmin*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapP_x2" value="2.0*GEMTrackerEndcapP_rmax*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapP_y" value="GEMTrackerEndcapP_rmax*cos(GEMTrackerEndcapOpeningAngle/2.0) - GEMTrackerEndcapP_rmin"/>
<constant name="GEMTrackerEndcapN_x1" value="2.0*GEMTrackerEndcapN_rmin*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapN_x2" value="2.0*GEMTrackerEndcapN_rmax*sin(GEMTrackerEndcapOpeningAngle/2.0)"/>
<constant name="GEMTrackerEndcapN_y" value="GEMTrackerEndcapN_rmax*cos(GEMTrackerEndcapOpeningAngle/2.0) - GEMTrackerEndcapN_rmin"/>
<comment>
Needs implimentation:
Carbon fiber support ring and spokes
Initial preliminary details can be found here: https://www.dropbox.com/s/grt9p8feqlsh52f/GEM-RD-GEMWheel.pdf?dl=0
This design is still very much preliminary
// Frame parameters accoring to Kondo's sbsCrossSection.pdf file),
<constant name="GEMTrackerEndcapFrame_thickness" value="18.0 * mm"/>
//frame thickness is from 1*3mm (drift gap) + 5*2mm frames
</comment>
<constant name="GEMTrackerEndcapFrame_thickness" value="13.0 * mm"/>
<constant name="GEMTrackerEndcapFrameBotEdge_width" value="10.0 * mm"/>
<constant name="GEMTrackerEndcapFrameTopEdge_width" value="10.0 * mm"/>
<constant name="GEMTrackerEndcapFrameTopEdge_width" value="15.0 * mm"/>
<constant name="GEMTrackerEndcapFrameSideEdge_width" value=" 8.0 * mm"/>
<comment> Account for the fact that GEMs have holes, so they are not solid volumes </comment>
<constant name="GEMAreaFactor" value="0.80" />
<comment>
From https://github.com/eic/EicToyModel/blob/master/source/eicroot/GemGeoParData.cxx
From https://github.com/eic/EicToyModel/blob/master/source/eicroot/GEMGeoParData.cxx
// FIXME: put aluminum layer later as well),
mEntranceWindowMaterial ( "GemKapton"),
mEntranceWindowMaterial ( "GEMKapton"),
mEntranceWindowThickness ( 50.0 * etm::um),
// Use evaristo.pdf p.10 for the foil parameters:
......@@ -39,16 +55,16 @@
// - readout foils : 30um kapton + 3um copper total)
mDriftFoilKaptonThickness ( 50.0 * etm::um),
mDriftFoilCopperThickness ( 3.0 * etm::um),
mGemFoilAreaFraction ( 0.80),
mGemFoilKaptonThickness ( 30.0 * etm::um),
mGemFoilCopperThickness ( 3.0 * etm::um),
mGEMFoilAreaFraction ( 0.80),
mGEMFoilKaptonThickness ( 30.0 * etm::um),
mGEMFoilCopperThickness ( 3.0 * etm::um),
mReadoutG10Thickness ( 0.0 * etm::mm),
mReadoutKaptonThickness ( 30.0 * etm::um),
mReadoutCopperThickness ( 3.0 * etm::um),
// 3mm thick Nomex honeycomb for SBS GEMs),
mReadoutSupportMaterial ( "GemNomex"),
mReadoutSupportMaterial ( "GEMNomex"),
mReadoutSupportThickness ( 3.0 * etm::mm),
mEntranceRegionLength ( 3.0 * etm::mm),
......@@ -71,136 +87,138 @@
<detectors>
<detector
id="GEMTrackerEndcap_ID"
name="GEMTrackerEndcap"
type="refdet_GEMTrackerEndcap"
id="TrackerEndcapP_3_ID"
name="GEMTrackerEndcapP"
type="athena_GEMTrackerEndcap"
readout="GEMTrackerEndcapHits"
vis="AnlRed"
vis="TrackerVis"
reflect="false">
<module name="GEMModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<module name="GEMModule1" vis="TrackerLayerVis">
<trd x1="GEMTrackerEndcapP_x1/2.0" x2="GEMTrackerEndcapP_x2/2.0" z="GEMTrackerEndcapP_y/2"/>
<comment>
Frame disabled as it gives issue with the ACTS geometry conversion
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapFoilX1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
x2="GEMTrackerEndcapFoilX2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
z="GEMTrackerEndcapFoilY/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<trd x1="GEMTrackerEndcapP_x1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
x2="GEMTrackerEndcapP_x2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
z="GEMTrackerEndcapP_y/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</comment>
<comment> Going from HV side to readout side</comment>
<module_component thickness="0.127 * mm" material="Mylar"/>
<module_component thickness="50.0*um" material="Kapton" name="entrance_window"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="entrance region" />
<module_component thickness="50.0*um" material="Kapton"/>
<module_component thickness=" 3.0*um" material="Copper"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="drift region"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region I"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region II"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction region"/>
<module_component thickness="30.0*um" material="Kapton" name="readout" sensitive="true"/>
<module_component thickness=" 3.0*um" material="Copper" name="readout_Cu"/>
<module_component thickness="127.0*um" material="Mylar"/>
<module_component thickness="200.0*um" material="Epoxy" sensitive="true" vis="AnlProcess_Blue"/>
<module_component thickness="25.0*um" material="Mylar" name="entrance_window" vis="TrackerLayerVis"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="entrance_region" />
<module_component thickness="50.0*um" material="Kapton" name="cathode_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="cathode_cu"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" sensitive="true" name="drift_region" vis="AnlProcess_Blue"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_1"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_1_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_I"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_2"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_2_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_II"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_3"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_3_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction_region"/>
<module_component thickness="50.0*um" material="Kapton" name="readout_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="readout_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="exit_region"/>
<module_component thickness="25.0*um" material="Mylar" name="exit_window" vis="TrackerLayerVis"/>
</module>
<module name="GEMSupportModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapFoilX1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
x2="GEMTrackerEndcapFoilX2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapFoilOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapFoilOpeningAngle)"
z="GEMTrackerEndcapFoilY/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</module>
<module name="GEMSupportModule2" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapFoilY/2"/>
<module_component thickness="1.0*mm" material="Mylar"/>
<module name="GEMSupportSpokes" vis="TrackerSupportVis">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapP_y/2"/>
<module_component thickness="1.0*mm" material="Mylar" vis="TrackerSupportVis"/>
</module>
<layer id="1" >
<ring vis="AnlRed" module="GEMModule1"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm" />
<!--
<ring vis="AnlRed" module="GEMSupportModule1"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm"/>
-->
<ring vis="AnlBlue" phi0="15.0*degree" module="GEMSupportModule2"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="0 * mm" />
<envelope vis="TrackerVis"
rmin = "GEMTrackerEndcapP_rmin"
rmax = "GEMTrackerEndcapP_rmax"
length="GEMTrackerEndcapP_thickness"
zstart="GEMTrackerEndcapP_zmin" />
<layer_material surface="representing" binning="binPhi,binR" bins0="40*12" bins1="1024"/>
<ring vis="TrackerVis" module="GEMModule1"
r="GEMTrackerEndcapP_rmin+GEMTrackerEndcapP_y/2.0"
zstart="0"
nmodules="12" dz="GEMTrackerEndcap_dz" />
<ring vis="AnlBlue" phi0="GEMTrackerEndcapOpeningAngle/2." module="GEMSupportSpokes"
r="GEMTrackerEndcapP_rmin+GEMTrackerEndcapP_y/2.0"
zstart="0"
nmodules="12" dz="0*mm" />
</layer>
</detector>
<detector
id="GEMTrackerEndcapN_ID"
id="TrackerEndcapN_3_ID"
name="GEMTrackerEndcapN"
type="refdet_GEMTrackerEndcap"
type="athena_GEMTrackerEndcap"
readout="GEMTrackerEndcapHits"
vis="AnlRed"
vis="TrackerVis"
reflect="true">
<module name="GEMModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX1/2.0" x2="GEMTrackerEndcapFoilX2/2.0" z="GEMTrackerEndcapFoilY/2"/>
<module name="GEMModule1" vis="TrackerLayerVis">
<trd x1="GEMTrackerEndcapN_x1/2.0" x2="GEMTrackerEndcapN_x2/2.0" z="GEMTrackerEndcapN_y/2"/>
<comment>
Frame disabled as it gives issue with the ACTS geometry conversion
<frame material="G10" thickness="GEMTrackerEndcapFrame_thickness" vis="TrackerSupportVis">
<position x="0" y="0" z="0" />
<trd x1="GEMTrackerEndcapN_x1/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
+ GEMTrackerEndcapFrameBotEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
x2="GEMTrackerEndcapN_x2/2.0 - GEMTrackerEndcapFrameSideEdge_width/cos(GEMTrackerEndcapOpeningAngle)
- GEMTrackerEndcapFrameTopEdge_width*tan(GEMTrackerEndcapOpeningAngle)"
z="GEMTrackerEndcapN_y/2 - (GEMTrackerEndcapFrameBotEdge_width +GEMTrackerEndcapFrameTopEdge_width)/2.0 ">
<position x="0" y="0" z="0" />
</trd>
</frame>
</comment>
<comment> Going from HV side to readout side</comment>
<module_component thickness="0.127 * mm" material="Mylar"/>
<module_component thickness="50.0*um" material="Kapton" name="entrance_window"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="entrance region" />
<module_component thickness="50.0*um" material="Kapton"/>
<module_component thickness=" 3.0*um" material="Copper"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" name="drift region"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region I"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer region II"/>
<module_component thickness="30.0*um" material="Kapton" name="gem_foil"/>
<module_component thickness=" 3.0*um" material="Copper" name="gem_foil_Cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction region"/>
<module_component thickness="30.0*um" material="Kapton" name="readout" sensitive="true"/>
<module_component thickness=" 3.0*um" material="Copper" name="readout_Cu"/>
<module_component thickness="127.0*um" material="Mylar"/>
<module_component thickness="200.0*um" material="Epoxy" sensitive="true" vis="AnlProcess_Blue"/>
<module_component thickness="25.0*um" material="Mylar" name="entrance_window" vis="TrackerLayerVis"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="entrance_region" />
<module_component thickness="50.0*um" material="Kapton" name="cathode_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="cathode_cu"/>
<module_component thickness=" 3.0*mm" material="Ar10CO2" sensitive="true" name="drift_region" vis="AnlProcess_Blue"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_1"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_1_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_I"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_2"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_2_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="transfer_region_II"/>
<module_component thickness="GEMAreaFactor*50.0*um" material="Kapton" name="gem_foil_3"/>
<module_component thickness="GEMAreaFactor*10.0*um" material="Copper" name="gem_foil_3_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="induction_region"/>
<module_component thickness="50.0*um" material="Kapton" name="readout_foil"/>
<module_component thickness=" 5.0*um" material="Copper" name="readout_cu"/>
<module_component thickness=" 2.0*mm" material="Ar10CO2" name="exit_region"/>
<module_component thickness="25.0*um" material="Mylar" name="exit_window" vis="TrackerLayerVis"/>
</module>
<module name="GEMSupportModule1" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFoilX2/2.0" x2="GEMTrackerEndcapFoilX1/2.0" z="GEMTrackerEndcapFrameBotEdge_width"/>
<module_component thickness="GEMTrackerEndcapFrame_thickness" material="Mylar"/>
</module>
<module name="GEMSupportModule2" vis="AnlProcess_Blue">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapFoilY/2"/>
<module_component thickness="4.0*mm" material="Mylar"/>
<module name="GEMSupportSpokes" vis="TrackerSupportVis">
<trd x1="GEMTrackerEndcapFrameSideEdge_width" x2="GEMTrackerEndcapFrameSideEdge_width" z="GEMTrackerEndcapN_y/2"/>
<module_component thickness="1.0*mm" material="Mylar" vis="TrackerSupportVis"/>
</module>
<layer id="1" >
<ring vis="AnlRed"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="10 * mm" module="GEMModule1" />
<ring vis="AnlBlue" phi0="15.0*degree"
r="GEMTrackerEndcapFoil_rmin+GEMTrackerEndcapFoilY/2.0"
zstart="GEMTrackerEndcap_zmin + 0.5*GEMTrackerEndcapLayer_thickness"
nmodules="12" dz="0 * mm" module="GEMSupportModule2" />
<envelope vis="TrackerVis"
rmin = "GEMTrackerEndcapN_rmin"
rmax = "GEMTrackerEndcapN_rmax"
length="GEMTrackerEndcapN_thickness"
zstart="GEMTrackerEndcapN_zmin" />
<layer_material surface="representing" binning="binPhi,binR" bins0="40*12" bins1="1024"/>
<ring vis="TrackerVis" module="GEMModule1"
r="GEMTrackerEndcapN_rmin+GEMTrackerEndcapN_y/2.0"
zstart="0.1"
nmodules="12" dz="GEMTrackerEndcap_dz" />
<ring vis="AnlBlue" phi0="GEMTrackerEndcapOpeningAngle/2." module="GEMSupportSpokes"
r="GEMTrackerEndcapN_rmin+GEMTrackerEndcapN_y/2.0"
zstart="0.1"
nmodules="12" dz="0*mm" />
</layer>
</detector>
</detectors>
<readouts>
<readout name="GEMTrackerEndcapHits">
<segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
<id>system:8,barrel:2,layer:4,module:12,sensor:2,x:32:-16,y:-16</id>
<segmentation type="CartesianGridXZ" grid_size_x="50*um" grid_size_z="250*um" />
<id>system:8,layer:4,module:12,sensor:2,x:32:-16,z:-16</id>
</readout>
</readouts>
......
<lccdd>
<define>
<constant name="HcalBarrel_offset" value="(HcalEndcapP_zmin-HcalEndcapN_zmin)/2.0" />
<constant name="HcalBarrel_length" value="HcalBarrelForward_zmax + HcalBarrelBackward_zmax"/>
<constant name="HcalBarrel_offset" value="(HcalBarrelForward_zmax - HcalBarrelBackward_zmax)/2"/>
<constant name="HcalBarrelAvailThickness" value="HcalBarrel_thickness"/>
<documentation>
#### Material Thickness
</documentation>
<constant name="HcalBarrelSteelThickness" value="2.0 * cm"/>
<constant name="HcalBarrelPolystyreneThickness" value="0.5 * cm"/>
<constant name="HcalEndcapNSteelThickness" value="HcalBarrelSteelThickness"/>
<constant name="HcalEndcapNPolystyreneThickness" value="HcalBarrelPolystyreneThickness"/>
<constant name="HcalEndcapPSteelThickness" value="HcalBarrelSteelThickness"/>
<constant name="HcalEndcapPPolystyreneThickness" value="0.3 * cm"/>
<documentation>
- Hcal Barrel Layers and computed Thickness
</documentation>
<constant name="HcalBarrelSingleLayerThickness"
value="HcalBarrelSteelThickness + HcalBarrelPolystyreneThickness"/>
<constant name="HcalBarrelLayer_NRepeat" value="floor(HcalBarrelAvailThickness / HcalBarrelSingleLayerThickness)"/>
<constant name="HcalBarrelThickness" value="HcalBarrelLayer_NRepeat * HcalBarrelSingleLayerThickness"/>
<documentation>
- Hcal Endcap N Layers and computed Thickness
</documentation>
<constant name="HcalEndcapNSingleLayerThickness"
value="HcalEndcapNSteelThickness + HcalEndcapNPolystyreneThickness"/>
<constant name="HcalEndcapNLayer_NRepeat" value="floor(HcalEndcapN_length / HcalEndcapNSingleLayerThickness)"/>
<constant name="HcalEndcapN_thickness" value="HcalEndcapNLayer_NRepeat * HcalEndcapNSingleLayerThickness"/>
<documentation>
- Hcal Endcap P Layers and computed Thickness
</documentation>
<constant name="HcalEndcapPSingleLayerThickness"
value="HcalEndcapPSteelThickness + HcalEndcapPPolystyreneThickness"/>
<constant name="HcalEndcapPLayer_NRepeat" value="floor(HcalEndcapP_length / HcalEndcapPSingleLayerThickness)"/>
<constant name="HcalEndcapP_thickness" value="HcalEndcapPLayer_NRepeat * HcalEndcapPSingleLayerThickness"/>
</define>
<limits>
......@@ -11,11 +48,6 @@
</regions>
<display>
<vis name="HcalVis" alpha="0.9" r="0.75" g="0.0" b="0.0" showDaughters="false"/>
<vis name="HcalBarrelVis" alpha="0.9" r="0.0" g="0.3" b="0.8" showDaughters="true" visible="true" />
<vis name="HcalEndcapVis" alpha="0.9" r="0.0" g="0.3" b="0.8" showDaughters="true" visible="true" />
<vis name="HcalSensorVis" alpha="0.9" r="1.0" g="0.0" b="0.2" showDaughters="true" visible="true" />
<vis name="HcalAbsorberVis" alpha="0.9" r="0.4" g="0.4" b="0.6" showDaughters="true" visible="true" />
</display>
<!-- Define detector -->
......@@ -26,6 +58,7 @@
name="HcalBarrel"
type="athena_EcalBarrel"
readout="HcalBarrelHits"
vis="HcalBarrelVis"
calorimeterType="HAD_BARREL"
offset="HcalBarrel_offset"
gap="0.*cm"
......@@ -33,19 +66,19 @@
<dimensions
numsides="CaloSides"
rmin="HcalBarrel_rmin"
z="HcalBarrelLength"/>
<staves vis="HcalBarrelVis"/>
<layer repeat="HcalBarrelLayer_NRepeat">
<slice material="Steel235" thickness="2.0*cm" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="0.5*cm" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
z="HcalBarrel_length"/>
<staves vis="HcalBarrelStaveVis"/>
<layer repeat="HcalBarrelLayer_NRepeat" vis="HcalBarrelLayerVis">
<slice material="Steel235" thickness="HcalBarrelSteelThickness" vis="HcalAbsorberVis"/>
<slice material="Polystyrene" thickness="HcalBarrelPolystyreneThickness" vis="HcalSensorVis" sensitive="yes" limits="cal_limits"/>
</layer>
</detector>
<detector
id="HCalEndcapP_ID"
name="HcalEndcapP"
type="refdet_PolyhedraEndcapCalorimeter2"
readout="HcalEndcapHits"
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="HcalEndcapHits"
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,10 +138,14 @@
<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="HcalEndcapHits">
<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="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>
</readouts>
<plugins>
......
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<comment>
--------------------------
Solenoid Magnet Parameters
--------------------------
</comment>
<comment>Solenoid Barrel Parameters</comment>
<constant name="SolenoidBarrelInnerVacuumVesselThickness" value=" 20.0 * mm"/>
<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="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 -->
<constant name="SolenoidBarrelCoilOverbindThickness" value="200.0 * mm"/>
<constant name="SolenoidBarrelOuterHeliumThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelOuterHeliumVesselThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelOuterFirstMLILayerThickness" value=" 15.0 * mm"/>
<constant name="SolenoidBarrelOuterThermalShieldThickness" value=" 10.0 * mm"/>
<constant name="SolenoidBarrelOuterSecondMLILayerThickness" value=" 20.0 * mm"/>
<constant name="SolenoidBarrelOuterVacuumVesselThickness" value=" 25.0 * mm"/>
<constant name="SolenoidBarrelInnerVacuumVessel_zmax" value="1920.0 * mm"/>
<constant name="SolenoidBarrelInnerSecondMLILayer_zmax" value="1895.0 * mm"/>
<constant name="SolenoidBarrelInnerThermalShield_zmax" value="1875.0 * mm"/>
<constant name="SolenoidBarrelInnerFirstMLILayer_zmax" value="1855.0 * mm"/>
<constant name="SolenoidBarrelInnerHeliumVessel_zmax" value="1840.0 * mm"/>
<constant name="SolenoidBarrelInnerHelium_zmax" value="1820.0 * mm"/>
<constant name="SolenoidBarrelCoilFormer_zmax" value="1800.0 * mm"/>
<constant name="SolenoidBarrelCoil_zmax" value="1800.0 * mm"/>
<constant name="SolenoidBarrelBlock_zmax" value=" 545.0 * mm"/> <!-- addition -->
<constant name="SolenoidBarrelCoilOverbind_zmax" value="1800.0 * mm"/>
<constant name="SolenoidBarrelOuterHelium_zmax" value="1820.0 * mm"/>
<constant name="SolenoidBarrelOuterHeliumVessel_zmax" value="1840.0 * mm"/>
<constant name="SolenoidBarrelOuterFirstMLILayer_zmax" value="1855.0 * mm"/>
<constant name="SolenoidBarrelOuterThermalShield_zmax" value="1875.0 * mm"/>
<constant name="SolenoidBarrelOuterSecondMLILayer_zmax" value="1895.0 * mm"/>
<constant name="SolenoidBarrelOuterVacuumVessel_zmax" value="1920.0 * mm"/>
<constant name="SolenoidBarrelInnerVacuumVessel_rmin" value="SolenoidBarrel_rmin"/>
<constant name="SolenoidBarrelInnerVacuumVessel_rmax" value="SolenoidBarrelInnerVacuumVessel_rmin + SolenoidBarrelInnerVacuumVesselThickness"/>
<constant name="SolenoidBarrelInnerSecondMLILayer_rmin" value="SolenoidBarrelInnerVacuumVessel_rmax"/>
<constant name="SolenoidBarrelInnerSecondMLILayer_rmax" value="SolenoidBarrelInnerSecondMLILayer_rmin + SolenoidBarrelInnerSecondMLILayerThickness"/>
<constant name="SolenoidBarrelInnerThermalShield_rmin" value="SolenoidBarrelInnerSecondMLILayer_rmax"/>
<constant name="SolenoidBarrelInnerThermalShield_rmax" value="SolenoidBarrelInnerThermalShield_rmin + SolenoidBarrelInnerThermalShieldThickness"/>
<constant name="SolenoidBarrelInnerFirstMLILayer_rmin" value="SolenoidBarrelInnerThermalShield_rmax"/>
<constant name="SolenoidBarrelInnerFirstMLILayer_rmax" value="SolenoidBarrelInnerFirstMLILayer_rmin + SolenoidBarrelInnerFirstMLILayerThickness"/>
<constant name="SolenoidBarrelInnerHeliumVessel_rmin" value="SolenoidBarrelInnerFirstMLILayer_rmax"/>
<constant name="SolenoidBarrelInnerHeliumVessel_rmax" value="SolenoidBarrelInnerHeliumVessel_rmin + SolenoidBarrelInnerHeliumVesselThickness"/>
<constant name="SolenoidBarrelInnerHelium_rmin" value="SolenoidBarrelInnerHeliumVessel_rmax"/>
<constant name="SolenoidBarrelInnerHelium_rmax" value="SolenoidBarrelInnerHelium_rmin + SolenoidBarrelInnerHeliumThickness"/>
<constant name="SolenoidBarrelCoilFormer_rmin" value="SolenoidBarrelInnerHelium_rmax"/>
<constant name="SolenoidBarrelCoilFormer_rmax" value="SolenoidBarrelCoilFormer_rmin + SolenoidBarrelCoilFormerThickness"/>
<constant name="SolenoidBarrelCoil_rmin" value="SolenoidBarrelCoilFormer_rmax"/>
<constant name="SolenoidBarrelCoil_rmax" value="SolenoidBarrelCoil_rmin + SolenoidBarrelCoilThickness"/>
<constant name="SolenoidBarrelBlock_rmin" value="SolenoidBarrelCoilFormer_rmax"/> <!-- addition -->
<constant name="SolenoidBarrelBlock_rmax" value="SolenoidBarrelBlock_rmin + SolenoidBarrelBlockThickness"/> <!-- addition -->
<constant name="SolenoidBarrelCoilOverbind_rmin" value="max(SolenoidBarrelCoil_rmax,SolenoidBarrelBlock_rmax)"/> <!-- modification -->
<constant name="SolenoidBarrelCoilOverbind_rmax" value="SolenoidBarrelCoilOverbind_rmin + SolenoidBarrelCoilOverbindThickness"/>
<constant name="SolenoidBarrelOuterHelium_rmin" value="SolenoidBarrelCoilOverbind_rmax"/>
<constant name="SolenoidBarrelOuterHelium_rmax" value="SolenoidBarrelOuterHelium_rmin + SolenoidBarrelOuterHeliumThickness"/>
<constant name="SolenoidBarrelOuterHeliumVessel_rmin" value="SolenoidBarrelOuterHelium_rmax"/>
<constant name="SolenoidBarrelOuterHeliumVessel_rmax" value="SolenoidBarrelOuterHeliumVessel_rmin + SolenoidBarrelOuterHeliumVesselThickness"/>
<constant name="SolenoidBarrelOuterFirstMLILayer_rmin" value="SolenoidBarrelOuterHeliumVessel_rmax"/>
<constant name="SolenoidBarrelOuterFirstMLILayer_rmax" value="SolenoidBarrelOuterFirstMLILayer_rmin + SolenoidBarrelOuterFirstMLILayerThickness"/>
<constant name="SolenoidBarrelOuterThermalShield_rmin" value="SolenoidBarrelOuterFirstMLILayer_rmax"/>
<constant name="SolenoidBarrelOuterThermalShield_rmax" value="SolenoidBarrelOuterThermalShield_rmin + SolenoidBarrelOuterThermalShieldThickness"/>
<constant name="SolenoidBarrelOuterSecondMLILayer_rmin" value="SolenoidBarrelOuterThermalShield_rmax"/>
<constant name="SolenoidBarrelOuterSecondMLILayer_rmax" value="SolenoidBarrelOuterSecondMLILayer_rmin + SolenoidBarrelOuterSecondMLILayerThickness"/>
<constant name="SolenoidBarrelOuterVacuumVessel_rmin" value="SolenoidBarrelOuterSecondMLILayer_rmax"/>
<constant name="SolenoidBarrelOuterVacuumVessel_rmax" value="SolenoidBarrelOuterVacuumVessel_rmin + SolenoidBarrelOuterVacuumVesselThickness"/>
<constant name="SolenoidalFieldRadius"
value="(SolenoidBarrelCoil_rmin + SolenoidBarrelCoil_rmax) / 2.0"/>
<comment> Solenoid Endcap </comment>
<comment>At the center, the end caps are outside the barrels</comment>
<constant name="SolenoidEndcapCoilEndSupport_rmin" value="SolenoidBarrelCoilFormer_rmin"/>
<constant name="SolenoidEndcapCoilEndSupport_rmax" value="SolenoidBarrelCoilOverbind_rmax"/>
<constant name="SolenoidEndcapHelium_rmin" value="SolenoidBarrelInnerHelium_rmin"/>
<constant name="SolenoidEndcapHelium_rmax" value="SolenoidBarrelOuterHelium_rmax"/>
<constant name="SolenoidEndcapHeliumVessel_rmin" value="SolenoidBarrelInnerHeliumVessel_rmin"/>
<constant name="SolenoidEndcapHeliumVessel_rmax" value="SolenoidBarrelOuterHeliumVessel_rmax"/>
<constant name="SolenoidEndcapFirstMLILayer_rmin" value="SolenoidBarrelInnerFirstMLILayer_rmin"/>
<constant name="SolenoidEndcapFirstMLILayer_rmax" value="SolenoidBarrelOuterFirstMLILayer_rmax"/>
<comment>Below here, the end caps are inside the barrels</comment>
<constant name="SolenoidEndcapThermalShield_rmin" value="SolenoidBarrelInnerThermalShield_rmax"/>
<constant name="SolenoidEndcapThermalShield_rmax" value="SolenoidBarrelOuterThermalShield_rmin"/>
<constant name="SolenoidEndcapSecondMLILayer_rmin" value="SolenoidBarrelInnerSecondMLILayer_rmax"/>
<constant name="SolenoidEndcapSecondMLILayer_rmax" value="SolenoidBarrelOuterSecondMLILayer_rmin"/>
<constant name="SolenoidEndcapVacuumVessel_rmin" value="SolenoidBarrelInnerVacuumVessel_rmax"/>
<constant name="SolenoidEndcapVacuumVessel_rmax" value="SolenoidBarrelOuterVacuumVessel_rmin"/>
<constant name="SolenoidEndcapCoilEndSupport_zmin" value="1800.0 * mm"/>
<constant name="SolenoidEndcapHelium_zmin" value="1820.0 * mm"/>
<constant name="SolenoidEndcapHeliumVessel_zmin" value="1840.0 * mm"/>
<constant name="SolenoidEndcapFirstMLILayer_zmin" value="1855.0 * mm"/>
<constant name="SolenoidEndcapThermalShield_zmin" value="1865.0 * mm"/>
<constant name="SolenoidEndcapSecondMLILayer_zmin" value="1875.0 * mm"/>
<constant name="SolenoidEndcapVacuumVessel_zmin" value="1895.0 * mm"/>
<constant name="SolenoidEndcapCoilEndSupportThickness" value="20.0 * mm"/>
<constant name="SolenoidEndcapHeliumThickness" value="20.0 * mm"/>
<constant name="SolenoidEndcapHeliumVesselThickness" value="10.0 * mm"/>
<constant name="SolenoidEndcapFirstMLILayerThickness" value="10.0 * mm"/>
<constant name="SolenoidEndcapThermalShieldThickness" value="10.0 * mm"/>
<constant name="SolenoidEndcapSecondMLILayerThickness" value="20.0 * mm"/>
<constant name="SolenoidEndcapVacuumVesselThickness" value="25.0 * mm"/>
</define>
<detectors>
<comment>
---------------
Solenoid Magnet
---------------
The solenoid superconducting magnet and surrounding cryostat and
yoke endcaps.
</comment>
<detector
id="Solenoid_ID"
name="SolenoidCoilBarrel"
type="athena_SolenoidCoil"
insideTrackingVolume="false">
<position x="0" y="0" z="Solenoid_offset"/>
<layer id="1"
inner_r="SolenoidBarrelInnerVacuumVessel_rmin"
outer_z="SolenoidBarrelInnerVacuumVessel_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelInnerVacuumVesselThickness"/>
</layer>
<layer id="2"
inner_r="SolenoidBarrelInnerSecondMLILayer_rmin"
outer_z="SolenoidBarrelInnerSecondMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelInnerSecondMLILayerThickness"/>
</layer>
<layer id="3"
inner_r="SolenoidBarrelInnerThermalShield_rmin"
outer_z="SolenoidBarrelInnerThermalShield_zmax"
vis="SolenoidCryostatVis">
<slice material="Copper" thickness="SolenoidBarrelInnerThermalShieldThickness"/>
</layer>
<layer id="4"
inner_r="SolenoidBarrelInnerFirstMLILayer_rmin"
outer_z="SolenoidBarrelInnerFirstMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelInnerFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_r="SolenoidBarrelInnerHeliumVessel_rmin"
outer_z="SolenoidBarrelInnerHeliumVessel_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelInnerHeliumVesselThickness"/>
</layer>
<layer id="6"
inner_r="SolenoidBarrelInnerHelium_rmin"
outer_z="SolenoidBarrelInnerHelium_zmax"
vis="SolenoidCryostatVis">
<slice material="Helium" thickness="SolenoidBarrelInnerHeliumThickness"/>
</layer>
<layer id="7"
inner_r="SolenoidBarrelCoilFormer_rmin"
outer_z="SolenoidBarrelCoilFormer_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelCoilFormerThickness"/>
</layer>
<layer id="8"
inner_r="SolenoidBarrelCoil_rmin"
outer_z="max(SolenoidBarrelCoil_zmax,SolenoidBarrelBlock_zmax)"
vis="SolenoidCryostatVis">
<slice material="Copper" thickness="SolenoidBarrelCoilThickness" outer_z="SolenoidBarrelCoil_zmax" inner_z="SolenoidBarrelBlock_zmax"/>
<slice material="Steel235" thickness="SolenoidBarrelBlockThickness" outer_z="SolenoidBarrelBlock_zmax"/>
</layer>
<layer id="9"
inner_r="SolenoidBarrelCoilOverbind_rmin"
outer_z="SolenoidBarrelCoilOverbind_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelCoilOverbindThickness"/>
</layer>
<layer id="10"
inner_r="SolenoidBarrelOuterHelium_rmin"
outer_z="SolenoidBarrelOuterHelium_zmax"
vis="SolenoidCryostatVis">
<slice material="Helium" thickness="SolenoidBarrelOuterHeliumThickness"/>
</layer>
<layer id="11"
inner_r="SolenoidBarrelOuterHeliumVessel_rmin"
outer_z="SolenoidBarrelOuterHeliumVessel_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelOuterHeliumVesselThickness"/>
</layer>
<layer id="12"
inner_r="SolenoidBarrelOuterFirstMLILayer_rmin"
outer_z="SolenoidBarrelOuterFirstMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelOuterFirstMLILayerThickness"/>
</layer>
<layer id="13"
inner_r="SolenoidBarrelOuterThermalShield_rmin"
outer_z="SolenoidBarrelOuterThermalShield_zmax"
vis="SolenoidCryostatVis">
<slice material="Copper" thickness="SolenoidBarrelOuterThermalShieldThickness"/>
</layer>
<layer id="14"
inner_r="SolenoidBarrelOuterSecondMLILayer_rmin"
outer_z="SolenoidBarrelOuterSecondMLILayer_zmax"
vis="SolenoidCryostatVis">
<slice material="Vacuum" thickness="SolenoidBarrelOuterSecondMLILayerThickness"/>
</layer>
<layer id="15"
inner_r="SolenoidBarrelOuterVacuumVessel_rmin"
outer_z="SolenoidBarrelOuterVacuumVessel_zmax"
vis="SolenoidCryostatVis">
<slice material="Steel235" thickness="SolenoidBarrelOuterVacuumVesselThickness"/>
</layer>
</detector>
<detector
id="Solenoid_ID"
name="SolenoidCoilEndcapP"
type="ref_SolenoidEndcap"
insideTrackingVolume="false"
reflect="false"
vis="SolenoidCryostatVis">
<position x="0" y="0" z="Solenoid_offset"/>
<layer id="1"
inner_z="SolenoidEndcapCoilEndSupport_zmin"
inner_r="SolenoidEndcapCoilEndSupport_rmin"
outer_r="SolenoidEndcapCoilEndSupport_rmax">
<slice material="Steel235" thickness="SolenoidEndcapCoilEndSupportThickness"/>
</layer>
<layer id="2"
inner_z="SolenoidEndcapHelium_zmin"
inner_r="SolenoidEndcapHelium_rmin"
outer_r="SolenoidEndcapHelium_rmax">
<slice material="Helium" thickness="SolenoidEndcapHeliumThickness"/>
</layer>
<layer id="3"
inner_z="SolenoidEndcapHeliumVessel_zmin"
inner_r="SolenoidEndcapHeliumVessel_rmin"
outer_r="SolenoidEndcapHeliumVessel_rmax">
<slice material="Steel235" thickness="SolenoidEndcapHeliumVesselThickness"/>
</layer>
<layer id="4"
inner_z="SolenoidEndcapFirstMLILayer_zmin"
inner_r="SolenoidEndcapFirstMLILayer_rmin"
outer_r="SolenoidEndcapFirstMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_z="SolenoidEndcapThermalShield_zmin"
inner_r="SolenoidEndcapThermalShield_rmin"
outer_r="SolenoidEndcapThermalShield_rmax">
<slice material="Copper" thickness="SolenoidEndcapThermalShieldThickness"/>
</layer>
<layer id="6"
inner_z="SolenoidEndcapSecondMLILayer_zmin"
inner_r="SolenoidEndcapSecondMLILayer_rmin"
outer_r="SolenoidEndcapSecondMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapSecondMLILayerThickness"/>
</layer>
<layer id="7"
inner_z="SolenoidEndcapVacuumVessel_zmin"
inner_r="SolenoidEndcapVacuumVessel_rmin"
outer_r="SolenoidEndcapVacuumVessel_rmax">
<slice material="Steel235" thickness="SolenoidEndcapVacuumVesselThickness"/>
</layer>
</detector>
<detector
id="Solenoid_ID"
name="SolenoidCoilEndcapN"
type="ref_SolenoidEndcap"
insideTrackingVolume="false"
reflect="true"
vis="SolenoidCryostatVis">
<position x="0" y="0" z="Solenoid_offset"/>
<layer id="1"
inner_z="SolenoidEndcapCoilEndSupport_zmin"
inner_r="SolenoidEndcapCoilEndSupport_rmin"
outer_r="SolenoidEndcapCoilEndSupport_rmax">
<slice material="Steel235" thickness="SolenoidEndcapCoilEndSupportThickness"/>
</layer>
<layer id="2"
inner_z="SolenoidEndcapHelium_zmin"
inner_r="SolenoidEndcapHelium_rmin"
outer_r="SolenoidEndcapHelium_rmax">
<slice material="Helium" thickness="SolenoidEndcapHeliumThickness"/>
</layer>
<layer id="3"
inner_z="SolenoidEndcapHeliumVessel_zmin"
inner_r="SolenoidEndcapHeliumVessel_rmin"
outer_r="SolenoidEndcapHeliumVessel_rmax">
<slice material="Steel235" thickness="SolenoidEndcapHeliumVesselThickness"/>
</layer>
<layer id="4"
inner_z="SolenoidEndcapFirstMLILayer_zmin"
inner_r="SolenoidEndcapFirstMLILayer_rmin"
outer_r="SolenoidEndcapFirstMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapFirstMLILayerThickness"/>
</layer>
<layer id="5"
inner_z="SolenoidEndcapThermalShield_zmin"
inner_r="SolenoidEndcapThermalShield_rmin"
outer_r="SolenoidEndcapThermalShield_rmax">
<slice material="Copper" thickness="SolenoidEndcapThermalShieldThickness"/>
</layer>
<layer id="6"
inner_z="SolenoidEndcapSecondMLILayer_zmin"
inner_r="SolenoidEndcapSecondMLILayer_rmin"
outer_r="SolenoidEndcapSecondMLILayer_rmax">
<slice material="Vacuum" thickness="SolenoidEndcapSecondMLILayerThickness"/>
</layer>
<layer id="7"
inner_z="SolenoidEndcapVacuumVessel_zmin"
inner_r="SolenoidEndcapVacuumVessel_rmin"
outer_r="SolenoidEndcapVacuumVessel_rmax">
<slice material="Steel235" thickness="SolenoidEndcapVacuumVesselThickness"/>
</layer>
</detector>
</detectors>
<fields>
<!--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"/-->
<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>
......@@ -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"/>
......@@ -147,6 +199,12 @@
<composite n="12" ref="H"/>
<composite n="3" ref="O"/>
</material>
<material name="TungstenPowder">
<D value="11.25" unit="g / cm3"/>
<fraction n="0.954" ref="W"/>
<fraction n="0.040" ref="Ni"/>
<fraction n="0.006" ref="Fe"/>
</material>
<material name="TungstenDens23">
<D value="17.7" unit="g / cm3"/>
<fraction n="0.925" ref="W"/>
......@@ -165,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"/>
......@@ -221,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"/>
......@@ -248,5 +318,68 @@
<D type="density" value="2.21" unit="g/cm3"/>
<composite n="1" ref="C"/>
</material>
<material name="PolyvinylAcetate">
<D type="density" value="1.19" unit="g/cm3"/>
<composite n="4" ref="C"/>
<composite n="6" ref="H"/>
<composite n="2" ref="O"/>
</material>
<material name="Plexiglass">
<D type="density" value="1.19" unit="g/cm3"/>
<composite n="5" ref="C"/>
<composite n="8" ref="H"/>
<composite n="2" ref="O"/>
</material>
<material name="StainlessSteel">
<D type="density" value="8.3" unit="g / cm3"/>
<fraction n="0.74" ref="Fe"/>
<fraction n="0.18" ref="Cr"/>
<fraction n="0.08" ref="Ni"/>
</material>
<material name="ZDC_Tungsten">
<D type="density" value="18.73" unit="g / cm3"/>
<fraction n="0.943" ref="W"/>
<fraction n="0.038" ref="Ni"/>
<fraction n="0.019" ref="Cu"/>
</material>
<material name="ZDC_EpoxyGlue">
<D type="density" value="1.38" unit="g/cm3"/>
<composite n="8" ref="H"/>
<composite n="10" ref="C"/>
<composite n="4" ref="O"/>
</material>
<material name="ZDC_Scintillator">
<D type="density" value="1.032" unit="g/cm3"/>
<composite n="8" ref="C"/>
<composite n="8" ref="H"/>
</material>
<material name="MMGAS_ResistivePaste">
<D type="density" unit="g/cm3" value="0.77906"/>
<composite n="1" ref="C"/>
</material>
<material name="MMGAS_InoxForMesh" state="solid">
<D type="density" unit="g/cm3" value="2.8548"/>
<fraction n="0.19" ref="Cr"/>
<fraction n="0.68" ref="Fe"/>
<fraction n="0.02" ref="Mn"/>
<fraction n="0.10" ref="Ni"/>
<fraction n="0.01" ref="Si"/>
</material>
<material name="Ar90IsoButane" state="gas">
<D type="density" unit="g/cm3" value="0.00170335"/>
<fraction n="0.90" ref="Ar"/>
<fraction n="0.0826586" ref="C"/>
<fraction n="0.0173414" ref="H"/>
</material>
<material name="ZDC_WSciFiMixture">
<D type="density" value="10.17" unit="g/cm3"/>
<composite n="0.950" ref="W"/>
<composite n="0.009" ref="H"/>
<composite n="0.041" ref="C"/>
</material>
<material name="WSciFi_UCLA_Abs">
<D type="density" value="12.4" unit="g / cm3"/>
<fraction n="0.96" ref="W"/>
<fraction n="0.04" ref="Polystyrene"/>
</material>
</materials>
<lccdd>
<comment> MRICH </comment>
<define>
<constant name="MRICH_rmin" value="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>
<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>
<?xml version="1.0" encoding="UTF-8"?>
<lccdd>
<define>
<!-- vessel geometry -->
<constant name="PFRICH_Length" value="BackwardRICHRegion_length"/> <!-- vessel z-length -->
<constant name="PFRICH_zmin" value="-BackwardRICHRegion_zmin"/> <!-- vessel front -->
<constant name="PFRICH_zmax" value="PFRICH_zmin - PFRICH_Length"/> <!-- vessel back -->
<constant name="PFRICH_rmin0" value="BackwardPIDRegion_tan * BackwardRICHRegion_zmin"/> <!-- bore radius at vessel frontplane -->
<constant name="PFRICH_rmin1" value="BackwardPIDRegion_tan * (BackwardRICHRegion_zmin + BackwardRICHRegion_length)"/> <!-- bore radius at vessel backplane -->
<constant name="PFRICH_rmax" value="BackwardPIDRegion_rmax - 2*cm"/> <!-- vessel backplane radius -->
<constant name="PFRICH_wall_thickness" value="0.5*cm"/> <!-- thickness of radial walls -->
<constant name="PFRICH_window_thickness" value="0.1*cm"/> <!-- thickness of entrance and exit walls -->
<!-- additional parameters -->
<constant name="PFRICH_aerogel_thickness" value="3.0*cm"/> <!-- aerogel thickness -->
<constant name="PFRICH_sensor_active_size_default" value="24.0*mm"/> <!-- sensor side length (effective area) -->
<constant name="PFRICH_sensor_full_size_default" value="25.8*mm"/> <!-- sensor side length (full size, with enclosure) -->
<!-- lores values are used for global vizualizations _only_ -->
<constant name="PFRICH_sensor_active_size_lores" value="PFRICH_sensor_active_size_default * 4"/>
<constant name="PFRICH_sensor_full_size_lores" value="PFRICH_sensor_full_size_default * 4"/>
<!-- actual values used by the parametrization -->
<constant name="PFRICH_sensor_active_size" value="PFRICH_sensor_active_size_default"/>
<constant name="PFRICH_sensor_full_size" value="PFRICH_sensor_full_size_default"/>
<constant name="PFRICH_sensor_thickness" value="1.5*mm"/> <!-- sensor thickness -->
<constant name="PFRICH_sensor_dist" value="40*cm"/> <!-- distance between aerogel exit plane and sensor entrance plane -->
<constant name="PFRICH_num_px" value="8"/> <!-- number of pixels along one side of the sensor -->
<!-- debugging switches -->
<comment>
- `PFRICH_debug_optics`: 1 = all components become vacuum; test opticalphotons from IP
2 = all components become vacuum, except for `gasvol`, test charged particles from IP
0 = off
</comment>
<constant name="PFRICH_debug_optics" value="0"/>
</define>
<detectors>
<!-- /detectors/detector -->
<documentation level="10">
### pfRICH: Proximity Focusing RICH
</documentation>
<detector
id="BackwardRICH_ID"
name="PFRICH"
type="athena_PFRICH"
readout="PFRICHHits"
gas="C4F10_PFRICH"
material="Aluminum"
vis_vessel="DRICH_vessel_vis"
vis_gas="DRICH_gas_vis"
debug_optics="PFRICH_debug_optics"
>
<!-- /detectors/detector/dimensions -->
<documentation level="10">
#### Vessel
- dimensions:
- `zmin`: z-position of vessel front plane
- `length`: overall z-length of the full vessel
- `rmin0` and `rmin1`: bore radius at front plane and back plane, respectively
- `rmax0` and `rmax1`: outer radius of vessel, at front plane and back plane, respectively
- `wall_thickness`: thickness of radial walls
- `window_thickness`: thickness of entrance and exit disks
</documentation>
<dimensions
zmin="PFRICH_zmin"
zmax="PFRICH_zmax"
length="PFRICH_Length"
rmin0="PFRICH_rmin0"
rmin1="PFRICH_rmin1"
rmax0="PFRICH_rmax"
rmax1="PFRICH_rmax"
wall_thickness="PFRICH_wall_thickness"
window_thickness="PFRICH_window_thickness"
/>
<!-- /detectors/detector/radiator -->
<documentation level="10">
#### Radiator
- radiator is defined in a wedge of azimuthal space, composed of aerogel and a
filter; the filter is applied to the back of the aerogel, so that it separates
the aerogel and gas radiators
- dimensions:
- `frontplane`: front of the aerogel, w.r.t. front plane of the vessel envelope
- `rmin` and `rmax`: inner and outer radius (at the front plane; radial bounds are conical)
- `phiw`: azimuthal width of wedge
- `thickness`: radiator thickness, defined separately for aerogel and filter
- `pitch`: controls the angle of the radiator (0=vertical)
</documentation>
<radiator
frontplane="-PFRICH_window_thickness"
rmin="PFRICH_rmin0 + PFRICH_wall_thickness + 0.2*cm"
rmax="(PFRICH_rmax/PFRICH_zmax)*PFRICH_zmin + 8.0*cm"
phiw="60*degree"
pitch="0*degree"
>
<aerogel
material="Aerogel_DRICH"
vis="DRICH_aerogel_vis"
thickness="PFRICH_aerogel_thickness"
/>
<filter
material="Acrylic_DRICH"
vis="DRICH_filter_vis"
thickness="0.3*mm"
/>
</radiator>
<!-- /detectors/detector/sensors -->
<documentation level="10">
#### Sensors
</documentation>
<sensors>
<!-- /detectors/detector/sensors/module -->
<documentation level="10">
##### Sensor module
- based on Hamamatsu S13361-3050AE-08
- effective area: 24.0 x 24.0 mm
- enclosure size: 25.8 x 25.8 mm
- pixel size: 3x3 mm (8x8 matrix)
- thickness: 1.5 mm
- dimensions:
- `side`: side length of the square module
- `thickness`: thickness of the sensor module
- `gap`: provides room between the squares, to help prevent them from overlapping
- note: the values of `side` and `gap` will determine how many sensors there are, since the
sensor placement algorithm will try to place as many as it can in the specified region
</documentation>
<module
material="Silicon"
surface="SensorSurface_DRICH"
vis="PFRICH_sensor_vis"
side="PFRICH_sensor_active_size"
thickness="PFRICH_sensor_thickness"
gap="0.5*(PFRICH_sensor_full_size-PFRICH_sensor_active_size) + 0.5*mm"
/>
<!-- /detectors/detector/sensors/plane -->
<documentation level="10">
##### Sensor plane
- sensors will be placed on a plane
- plane dimensions:
- `sensordist`: distance between sensor plane active surface (e.g., photocathode) and aerogel backplane
- `rmin`: minimum radial position of a sensor's centroid
- `rmax`: maximum radial position of a sensor's centroid
</documentation>
<plane
sensordist="PFRICH_sensor_dist"
rmin="PFRICH_rmin1 + 2*cm"
rmax="PFRICH_rmax - 4*cm"
/>
<services>
<comment>
Material should be equivallent with 3x0.5cm Al, spread over the entire available distance.
This was confirmed by Silvia.
</comment>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
<component name="air" thickness="40*mm" vis="InvisibleNoDaughters" material="Air"/>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
<component name="air" thickness="40*mm" vis="InvisibleNoDaughters" material="Air"/>
<component name="aluminum" thickness="5*mm" vis="RICHServiceVis" material="Aluminum"/>
</services>
</sensors>
</detector>
</detectors>
<documentation level="10">
#### Readout
- segmentation: square matrix of pixels
- `grid_size_x,y`: size of each sensor, but note we must divide sensor size
by 1 less than the number of pixels, to account for fenceposting
- `offset_x,y`: specified such that the `x` and `y` indicators are unsigned
- indicators and `cellID` bits:
| indicator | offset | length |
|-----------|--------|--------|
| pfRICH ID | 0 | 8 |
| sensor | 8 | 12 |
| x pixel | 20 | 16 |
| y pixel | 36 | 16 |
</documentation>
<readouts>
<readout name="PFRICHHits">
<segmentation
type="CartesianGridXY"
grid_size_x="PFRICH_sensor_active_size/(PFRICH_num_px-1)"
grid_size_y="PFRICH_sensor_active_size/(PFRICH_num_px-1)"
offset_x="-PFRICH_sensor_active_size/2.0"
offset_y="-PFRICH_sensor_active_size/2.0"
/>
<id>system:8,module:12,x:20:16,y:16</id>
</readout>
</readouts>
</lccdd>
<lccdd>
<include ref="fake_dirc.xml"/>
<!--include ref="dirc.xml"/-->
<include ref="mrich.xml"/>
<include ref="drich.xml"/>
</lccdd>