Skip to content

Commit

Permalink
CLD_o4_v05: Add LAr ECal Barrel and adapt main file of the detector d…
Browse files Browse the repository at this point in the history
…imensions.

Add materials needed for the LAr

There however are overlaps inside of LAr bath which experts on LAr should maybe take care of
  • Loading branch information
SwathiSasikumar committed Aug 25, 2023
1 parent 971d333 commit 1a708a9
Show file tree
Hide file tree
Showing 3 changed files with 197 additions and 15 deletions.
30 changes: 15 additions & 15 deletions FCCee/CLD/compact/CLD_o4_v05/CLD_o4_v05.xml
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">

<info name="CLD_o2_v05"
<info name="CLD_o4_v05"
title="FCCee detector model option 2 version 05 (updated Beampipe and VXD)"
author="Andre Sailer"
url="http://ilcsoft.desy.de"
Expand All @@ -17,7 +17,7 @@
</includes>

<define>
<constant name="world_side" value="6100*mm"/>
<constant name="world_side" value="7000*mm"/>
<constant name="world_x" value="world_side"/>
<constant name="world_y" value="world_side"/>
<constant name="world_z" value="world_side"/>
Expand Down Expand Up @@ -137,12 +137,12 @@
<constant name="OuterTracker_half_length" value="2300*mm"/>

<constant name="ECalBarrel_inner_radius" value="2150*mm"/>
<constant name="ECalBarrel_outer_radius" value="2352*mm"/>
<constant name="ECalBarrel_outer_radius" value="2820*mm"/>
<constant name="ECalBarrel_half_length" value="2210*mm"/>
<constant name="ECalBarrel_symmetry" value="12"/>

<constant name="ECalEndcap_inner_radius" value="340*mm"/>
<constant name="ECalEndcap_outer_radius" value="2455*mm"/>
<constant name="ECalEndcap_outer_radius" value="2820*mm"/>
<constant name="ECalEndcap_min_z" value="2307*mm"/>
<constant name="ECalEndcap_max_z" value="2509*mm"/>
<constant name="ECalEndcap_outer_symmetry" value="12"/>
Expand All @@ -156,38 +156,38 @@
<constant name="ECalPlug_outer_symmetry" value="12"/>
<constant name="ECalPlug_inner_symmetry" value="12"/>

<constant name="HCalBarrel_inner_radius" value="2400*mm"/>
<constant name="HCalBarrel_outer_radius" value="3566*mm"/>
<constant name="HCalBarrel_inner_radius" value="2850*mm"/>
<constant name="HCalBarrel_outer_radius" value="4016*mm"/>
<constant name="HCalBarrel_half_length" value="2210*mm"/>
<constant name="HCalBarrel_symmetry" value="12"/>

<constant name="HCalEndcap_inner_radius" value="340*mm"/>
<constant name="HCalEndcap_outer_radius" value="3566*mm"/>
<constant name="HCalEndcap_outer_radius" value="4016*mm"/>
<constant name="HCalEndcap_min_z" value="2539*mm"/>
<constant name="HCalEndcap_max_z" value="3705*mm"/>
<constant name="HCalEndcap_symmetry" value="12"/>
<constant name="HCalEndcap_zcutout" value="200*mm"/>
<constant name="HCalEndcap_rcutout" value="0*mm"/>

<constant name="HCalRing_inner_radius" value="2475*mm"/>
<constant name="HCalRing_inner_radius" value="2850*mm"/>
<constant name="HCalRing_outer_radius" value="HCalEndcap_outer_radius"/>
<constant name="HCalRing_min_z" value="2353.5*mm"/>
<constant name="HCalRing_max_z" value="HCalEndcap_min_z"/>
<constant name="HCalRing_symmetry" value="12"/>

<constant name="Solenoid_inner_radius" value="3719*mm"/>
<constant name="Solenoid_outer_radius" value="4272*mm"/>
<constant name="Solenoid_inner_radius" value="4169*mm"/>
<constant name="Solenoid_outer_radius" value="4722*mm"/>
<constant name="Solenoid_half_length" value="3705*mm"/>
<constant name="Solenoid_Coil_half_length" value="3476*mm"/>
<constant name="Solenoid_Coil_radius" value="3930*mm"/>
<constant name="Solenoid_Coil_radius" value="4380*mm"/>

<constant name="YokeBarrel_inner_radius" value="4479*mm"/>
<constant name="YokeBarrel_outer_radius" value="6000*mm"/>
<constant name="YokeBarrel_inner_radius" value="4929*mm"/>
<constant name="YokeBarrel_outer_radius" value="6450*mm"/>
<constant name="YokeBarrel_half_length" value="3755*mm"/>
<constant name="YokeBarrel_symmetry" value="12"/>

<constant name="YokeEndcap_inner_radius" value="400*mm"/>
<constant name="YokeEndcap_outer_radius" value="6000*mm"/>
<constant name="YokeEndcap_outer_radius" value="6450*mm"/>
<constant name="YokeEndcap_min_z" value="3755*mm"/>
<constant name="YokeEndcap_max_z" value="5300*mm"/>
<constant name="YokeEndcap_outer_symmetry" value="12"/>
Expand Down Expand Up @@ -298,7 +298,7 @@
<include ref="InnerTracker_o2_v07.xml"/>
<include ref="OuterTracker_o2_v07.xml"/>

<include ref="ECalBarrel_o2_v01_03.xml"/>
<include ref="LAr_ECalBarrel.xml"/>
<include ref="ECalEndcap_o2_v01_03.xml"/>

<include ref="HCalBarrel_o1_v01_01.xml"/>
Expand Down
144 changes: 144 additions & 0 deletions FCCee/CLD/compact/CLD_o4_v05/LAr_ECalBarrel.xml
Original file line number Diff line number Diff line change
@@ -0,0 +1,144 @@
<?xml version="1.0" encoding="UTF-8"?>
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">

<info name="FCCee_ECalBarrel"
title="Settings for FCCee Inclined ECal Barrel Calorimeter"
author="M.Aleksa,J.Faltova,A.Zaborowska, V. Volkl"
url="no"
status="development"
version="1.0">
<comment>
Settings for the inclined EM calorimeter.
The barrel is filled with liquid argon. Passive material includes lead in the middle and steal on the outside, glued together.
Passive plates are inclined by a certain angle from the radial direction.
In between of two passive plates there is a readout.
Space between the plate and readout is of trapezoidal shape and filled with liquid argon.
Definition of sizes, visualization settings, readout and longitudinal segmentation are specified.
</comment>
</info>

<define>
<!-- Inclination angle of the lead plates -->
<constant name="InclinationAngle" value="50*degree"/>
<!-- thickness of active volume between two absorber plates at barrel Rmin, measured perpendicular to the readout plate -->
<constant name="LArGapThickness" value="1.239749*2*mm"/>

<!-- Air margin, thicknesses of cryostat and LAr bath -->
<constant name="AirMarginThickness" value="54*mm"/> <!-- Space holder for air gap between cryostat vessels -->

<constant name="CryoBarrelFrontWarm" value="10*mm"/> <!-- Al solid corresponding to 0.11 X0 -->
<constant name="CryoBarrelFrontCold" value="3.8*mm"/> <!-- Al solid equivalent of 0.043 X0 sandwich CFRP -->
<constant name="CryoBarrelFront" value="CryoBarrelFrontWarm+CryoBarrelFrontCold"/>

<constant name="CryoBarrelBackCold" value="30*mm"/> <!-- Al solid corresponding to 0.34 X0 -->
<constant name="CryoBarrelBackWarm" value="2.7*mm"/> <!-- Al solid equivalent of 0.03 X0 sandwich CFRP -->
<constant name="SolenoidBarrel" value="70*mm"/> <!-- Al solenoid with thickness of 0.8 X0 -->
<constant name="CryoBarrelBack" value="CryoBarrelBackWarm+SolenoidBarrel+CryoBarrelBackCold"/>

<constant name="CryoBarrelSideWarm" value="30*mm"/>
<constant name="CryoBarrelSideCold" value="3.8*mm"/>
<constant name="CryoBarrelSide" value="CryoBarrelSideWarm+CryoBarrelSideCold"/>

<constant name="LArBathThicknessFront" value="5*mm"/>
<constant name="LArBathThicknessBack" value="40*mm"/>

<!-- air margin around calorimeter -->
<constant name="BarCryoECal_rmin" value="ECalBarrel_inner_radius+AirMarginThickness"/>
<constant name="BarCryoECal_rmax" value="ECalBarrel_outer_radius-AirMarginThickness"/>
<constant name="BarCryoECal_dz" value="ECalBarrel_half_length"/>
<!-- calorimeter active volume -->
<constant name="EMBarrel_rmin" value="BarCryoECal_rmin+CryoBarrelFront+LArBathThicknessFront"/>
<constant name="EMBarrel_rmax" value="BarCryoECal_rmax-CryoBarrelBack-LArBathThicknessBack"/>
<constant name="EMBarrel_dz" value="(ECalBarrel_half_length)-CryoBarrelSide"/>
<!-- thickness of active volume between two absorber plates at EMBarrel_rmin, measuring perpendicular to the readout plate -->
<constant name="LAr_thickness" value="LArGapThickness"/>
<!-- passive layer consists of lead in the middle and steel on the outside, glued -->
<!-- When employing trapezoidal planes Pb_thickness corresponds to the minimum thickness, i.e at the front of the calo -->
<constant name="Pb_thickness" value="1.80*mm"/>
<constant name="planeLength" value="-EMBarrel_rmin*cos(InclinationAngle) + sqrt(EMBarrel_rmax*EMBarrel_rmax - EMBarrel_rmin*EMBarrel_rmin*sin(InclinationAngle)*sin(InclinationAngle))"/>
<constant name="phi" value="asin(planeLength / EMBarrel_rmax * sin(InclinationAngle))"/>
<!-- use a different value for Pb_thickness_max when employing trapezoidal planes -->
<!-- approximate constant sampling fraction: make the absorber grow linearly with the radius,
taking into account the angular projection effect -->
<!-- <constant name="Pb_thickness_max" value="1.3 * Pb_thickness * EMBarrel_rmax/EMBarrel_rmin *
cos(InclinationAngle - phi) / cos(InclinationAngle)" />-->
<constant name="Pb_thickness_max" value="Pb_thickness" />
<!-- total amount of steel in one passive plate: it is divided for the outside layer on top and bottom -->
<constant name="Steel_thickness" value="0.1*mm"/>
<!-- total amount of glue in one passive plate: it is divided for the outside layer on top and bottom -->
<constant name="Glue_thickness" value="0.1*mm"/>
<!-- readout in between two absorber plates -->
<constant name="readout_thickness" value="1.2*mm"/>
<constant name="ECal_cell_size" value="5.1*mm"/>
</define>

<display>
<vis name="ecal_envelope" r="0.1" g="0.2" b="0.6" alpha="1" showDaughers="false" visible="true" />
</display>

<readouts>
<!-- readout for the simulation -->
<!-- offset in eta is eta max value including cryostat -->
<readout name="ECalBarrelEta">
<!-- segmentation type="GridTheta" grid_size_theta="0.5625" offset_theta="-0.83"/ -->
<segmentation type="GridEta" grid_size_eta="0.01" offset_eta="-1.2"/>
<id>system:4,cryo:1,type:3,subtype:3,layer:8,module:11,eta:9</id>
</readout>
<!-- readout for the reconstruction -->
<!-- phi position is calculated based on the centre of volume (hence it cannot be done in the simulation from energy deposits position) -->
<readout name="ECalBarrelPhiEta">
<!-- segmentation type="FCCSWGridPhiTheta" grid_size_theta="0.5625" phi_bins="768" offset_theta="-0.83" offset_phi="-pi+(pi/768.)"/ -->
<segmentation type="FCCSWGridPhiEta" grid_size_eta="0.01" phi_bins="768" offset_eta="-1.2" offset_phi="-pi+(pi/768.)"/>
<id>system:4,cryo:1,type:3,subtype:3,layer:8,eta:9,phi:10</id>
</readout>
</readouts>

<detectors>
<detector id="DetID_ECal_Barrel" name="ECalBarrel" type="EmCaloBarrelInclined" readout="ECalBarrelEta">
<type_flags type=" DetType_CALORIMETER + DetType_ELECTROMAGNETIC + DetType_BARREL"/>
<sensitive type="SimpleCalorimeterSD"/>
<dimensions rmin="BarCryoECal_rmin" rmax="BarCryoECal_rmax" dz="BarCryoECal_dz" vis="ecal_envelope"/>
<cryostat name="ECAL_Cryo">
<material name="Aluminum"/>
<dimensions rmin1="BarCryoECal_rmin" rmin2="BarCryoECal_rmin+CryoBarrelFront" rmax1="BarCryoECal_rmax-CryoBarrelBack" rmax2="BarCryoECal_rmax" dz="BarCryoECal_dz"/>
<front sensitive="false"/> <!-- inner wall of the cryostat -->
<side sensitive="false"/> <!-- both sides of the cryostat -->
<back sensitive="false"/> <!-- outer wall of the cryostat -->
</cryostat>
<calorimeter name="EM_barrel">
<!-- offset defines the numbering of the modules: module==0 for phi=0 direction -->
<dimensions rmin="EMBarrel_rmin" rmax="EMBarrel_rmax" dz="EMBarrel_dz" offset="-InclinationAngle"/>
<active thickness="LAr_thickness">
<material name="LAr"/>
<!-- overlap offset is a specific feature of the construction; do not change! -->
<!-- one volume for a gap on both side of the readout) -->
<overlap offset="0.5"/>
</active>
<passive>
<rotation angle="InclinationAngle"/> <!-- inclination angle -->
<inner thickness="Pb_thickness" sensitive="false">
<material name="Lead"/>
</inner>
<innerMax thickness="Pb_thickness_max" sensitive="false">
<material name="Lead"/>
</innerMax>
<glue thickness="Glue_thickness" sensitive="false">
<material name="lArCaloGlue"/>
</glue>
<outer thickness="Steel_thickness" sensitive="false">
<material name="lArCaloSteel"/>
</outer>
</passive>
<readout thickness="readout_thickness" sensitive="false">
<material name="PCB"/>
</readout>
<layers> <!-- pcb electrode segmentation in the radial direction -->
<layer thickness="1.5*cm" repeat="1"/>
<layer thickness="3.5*cm" repeat="11"/>
</layers>
</calorimeter>
</detector>
</detectors>
</lccdd>
38 changes: 38 additions & 0 deletions FCCee/CLD/compact/CLD_o4_v05/materials.xml
Original file line number Diff line number Diff line change
@@ -1,6 +1,44 @@
<?xml version="1.0" encoding="UTF-8"?>
<materials>

<material name="FR4">
<D type="density" unit="g/cm3" value="1.7"/>
<fraction n="0.1808" ref="Si"/>
<fraction n="0.4056" ref="O"/>
<fraction n="0.2780" ref="C"/>
<fraction n="0.0684" ref="H"/>
<fraction n="0.0671" ref="Br"/>
</material>

<material name="LAr">
<D value="1.396" unit="g/cm3" />
<composite n="1" ref="Ar"/>
</material>

<material name="LKr">
<D value="2.416" unit="g/cm3" />
<composite n="1" ref="Kr"/>
</material>

<!-- prepreg glue (C5H8O4Si) between steel and lead in absorber of EM calo -->
<!-- from ATL-LARG-PUB-2009-001 -->
<material name="lArCaloGlue">
<D value="1.69" unit="g/cm3" />
<composite n="5" ref="C"/>
<composite n="8" ref="H" />
<composite n="4" ref="O" />
<composite n="1" ref="Si" />
</material>

<!-- stainless steel as in ATLAS EM calo -->
<!-- from ATL-LARG-PUB-2009-001 -->
<material name="lArCaloSteel">
<D value="7.84" unit="g/cm3" />
<fraction n="0.7175" ref="Fe"/>
<fraction n="0.19" ref="Cr" />
<fraction n="0.0925" ref="Ni" />
</material>

<!--
Air by weight from
Expand Down

0 comments on commit 1a708a9

Please sign in to comment.