Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ECalInclinedFCCee updates #60

Open
wants to merge 13 commits into
base: main
Choose a base branch
from
Open

ECalInclinedFCCee updates #60

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
385ee39
HCAL endcaps in IDEA-LAr
faltovaj Jan 25, 2023
13b114e
CMakeFile fixed
faltovaj Feb 2, 2023
925b154
Merge remote-tracking branch 'origin/main' into fullLarDetector
faltovaj Mar 3, 2023
443d7dd
Merge remote-tracking branch 'origin/main' into fullLarDetector
faltovaj Apr 4, 2023
21b5774
ECal inclined src for FCCee
faltovaj Apr 13, 2023
6e8d700
ECAL barrel double-vessel front cryo
faltovaj May 23, 2023
eb1c999
ECAL barrel double-vessel back cryo
faltovaj May 23, 2023
4cc4627
double vessel
faltovaj Aug 28, 2023
ddec6ab
merge
faltovaj Aug 28, 2023
4e6d030
xml file fixed
faltovaj Aug 29, 2023
8e30cdd
double vessel cryo
faltovaj Aug 30, 2023
973e99f
info added
faltovaj Aug 31, 2023
9c47afe
Merge remote-tracking branch 'origin' into cryostatUpdate
faltovaj Aug 31, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 6 additions & 0 deletions Detector/DetFCCeeECalInclined/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,3 +3,9 @@
################################################################################

install(DIRECTORY ${CMAKE_CURRENT_LIST_DIR}/compact DESTINATION ${CMAKE_INSTALL_DATADIR}/${CMAKE_PROJECT_NAME}/Detector/DetFCCeeECalInclined)

file(GLOB sources src/*.cpp)
add_dd4hep_plugin(DetFCCeeECalInclined SHARED ${sources})
target_link_libraries(DetFCCeeECalInclined DD4hep::DDCore)

set(LIBRARY_OUTPUT_PATH ${CMAKE_LIBRARY_OUTPUT_DIRECTORY})
6 changes: 4 additions & 2 deletions Detector/DetFCCeeECalInclined/compact/FCCee_ECalBarrel_upstream.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,9 @@
</comment>
</info>

<define><constant name="BarECal_rmax" value="3840*mm"/>
<define>
<!-- Redefinition of a constant from Detector/DetFCCeeIDEA-LAr/compact/FCCee_DectDimensions.xml -->
<constant name="BarECal_rmax" value="3840*mm"/>
<!-- 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 -->
Expand Down Expand Up @@ -132,4 +134,4 @@
</calorimeter>
</detector>
</detectors>
</lccdd>
</lccdd>
163 changes: 163 additions & 0 deletions Detector/DetFCCeeECalInclined/compact/FCCee_ECalBarrel_v2.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,163 @@
<?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 using FCCee factory (EmCaloBarrelInclinedFCCee).
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="4*mm"/>
<constant name="CryoAirGapThickness" value="50*mm"/> <!-- 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+CryoAirGapThickness"/>

<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+CryoAirGapThickness"/>

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

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

<!-- air margin around calorimeter -->
<constant name="BarCryoECal_rmin" value="BarECal_rmin+AirMarginThickness"/>
<constant name="BarCryoECal_rmax" value="BarECal_rmax-AirMarginThickness"/>
<constant name="BarCryoECalWarm_dz" value="BarECal_dz"/>
<constant name="BarCryoECalCold_dz" value="BarECal_dz-CryoBarrelSideWarm-CryoAirGapThickness"/>
<!-- 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="BarECal_dz-CryoBarrelSide"/>
<constant name="EMBarrel_warm_dz" value="BarECal_dz-CryoBarrelSideWarm"/>
<!-- 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"/>
</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:4,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:4,subtype:3,layer:8,eta:9,phi:10</id>
</readout>
</readouts>

<detectors>
<detector id="BarECal_id" name="ECalBarrel" type="EmCaloBarrelInclinedFCCee" readout="ECalBarrelEta">
<type_flags type=" DetType_CALORIMETER + DetType_ELECTROMAGNETIC + DetType_BARREL"/>
<sensitive type="SimpleCalorimeterSD"/>
<dimensions rmin="BarCryoECal_rmin" rmax="BarCryoECal_rmax" dz="BarCryoECalWarm_dz" vis="ecal_envelope"/>
<cryostat name="ECAL_Cryo">
<material name="Aluminum"/>
<cryo_front_warm sentive="false">
<dimensions rmin="BarCryoECal_rmin" rmax = "BarCryoECal_rmin+CryoBarrelFrontWarm" dz="BarCryoECalWarm_dz"/>
</cryo_front_warm>
<cryo_front_cold sentive="false">
<dimensions rmin="BarCryoECal_rmin+CryoBarrelFrontWarm+CryoAirGapThickness" rmax="BarCryoECal_rmin+CryoBarrelFrontWarm+CryoAirGapThickness+CryoBarrelFrontCold" dz="BarCryoECalCold_dz"/>
</cryo_front_cold>
<cryo_back_cold sentive="false">
<dimensions rmin="BarCryoECal_rmax-CryoBarrelBackWarm-CryoAirGapThickness-SolenoidBarrel-CryoBarrelBackCold" rmax="BarCryoECal_rmax-CryoBarrelBackWarm-CryoAirGapThickness-SolenoidBarrel" dz="BarCryoECalCold_dz"/>
</cryo_back_cold>
<cryo_solenoid sentive="false">
<dimensions rmin="BarCryoECal_rmax-CryoBarrelBackWarm-0.5*CryoAirGapThickness-SolenoidBarrel" rmax="BarCryoECal_rmax-CryoBarrelBackWarm-0.5*CryoAirGapThickness" dz="BarCryoECalCold_dz"/>
</cryo_solenoid>
<cryo_back_warm sentive="false">
<dimensions rmin="BarCryoECal_rmax-CryoBarrelBackWarm" rmax = "BarCryoECal_rmax" dz="BarCryoECalWarm_dz"/>
</cryo_back_warm>
<cryo_side_cold sensitive="false">
<dimensions rmin="BarCryoECal_rmin+CryoBarrelFrontWarm+CryoAirGapThickness+CryoBarrelFrontCold" rmax="BarCryoECal_rmax-CryoBarrelBackWarm-CryoBarrelBackCold-SolenoidBarrel-CryoAirGapThickness" dz="BarCryoECalCold_dz"/>
</cryo_side_cold>
<cryo_side_warm sensitive="false">
<dimensions rmin="BarCryoECal_rmin+CryoBarrelFrontWarm" rmax="BarCryoECal_rmax-CryoBarrelBackWarm" dz="BarCryoECalWarm_dz"/>
</cryo_side_warm>
</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" dim_z="EMBarrel_warm_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>
Loading