Skip to content

Multi sensitive volume mpgd #1014

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

Open
ybedfer wants to merge 8 commits into
base: main
Choose a base branch
from
+845 −326
Open

Multi sensitive volume mpgd #1014

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
58ac412
MPGD Barrel: Add "5SensVolumes" versions of XML. Modified Detector Co…
Nov 13, 2025
0b0b444
Merge branch 'main' into Multi-SensitiveVolume-MPGD
Nov 14, 2025
84febdd
MPGD Barrel: Simplification. "5SensVolumes" renamed "2DStrip". Add de…
Nov 16, 2025
22eec1e
Merge branch 'main' into Multi-SensitiveVolume-MPGD
Nov 26, 2025
5c3e53f
MPGD Barrel: Backward compatibility via cmake option "-DMPGD_INSTALL_…
Nov 26, 2025
f30c59f
MPGD Barrel: Install "_pixel" flavors of XMLs by default. Erase pre-e…
Dec 8, 2025
53a2b0b
[pre-commit.ci] auto fixes from pre-commit.com hooks
pre-commit-ci[bot] Dec 8, 2025
9e66406
Merge branch 'main' into Multi-SensitiveVolume-MPGD
ybedfer Dec 8, 2025
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
52 changes: 50 additions & 2 deletions CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -92,11 +92,59 @@ install(FILES ${CMAKE_CURRENT_BINARY_DIR}/${PROJECT_NAME}_full.xml
)

#-----------------------------------------------------------------------------------
# Install the detector description files.
# Install the detector description files...
# ...excluding MPGD pixel/2DStrip, which are handled later on
install(DIRECTORY compact/
DESTINATION ${CMAKE_INSTALL_DATADIR}/${PROJECT_NAME}/compact
FILES_MATCHING PATTERN "*.xml" PATTERN "*.gdml"
FILES_MATCHING
PATTERN "*.xml"
PATTERN "*.gdml"
PATTERN "*_pixel.xml" EXCLUDE
PATTERN "*_2DStrip.xml" EXCLUDE
)
# Install MPGD files:
# - pixel or 2DStrip version depending upon option "MPGD_INSTALL_PIXEL"
# - The "install" command turns out to not be able to overwrite pre-existing.
# => We use the "file" command to work around that. Which later command takes
# effect immediately, at cmake execution time. => Emit a warning message.
set(MPGD_INSTALL_PIXEL TRUE CACHE BOOL "Use \"pixel\" flavour of MPGD XMLs instead of \"2DStrip\". It's temporarily the default.")
IF(${MPGD_INSTALL_PIXEL})
file(GLOB_RECURSE MPGD_FILES
RELATIVE ${CMAKE_CURRENT_SOURCE_DIR}/compact/tracking
"compact/tracking/*_pixel.xml"
)
foreach(FILE IN ITEMS ${MPGD_FILES})
execute_process(COMMAND bash -c "echo ${FILE} | sed -e 's/_pixel//'" OUTPUT_VARIABLE INSTALL_FILE OUTPUT_STRIP_TRAILING_WHITESPACE)
set(PRE_EXISTING_FILE ${CMAKE_INSTALL_PREFIX}/${CMAKE_INSTALL_DATADIR}/${PROJECT_NAME}/compact/tracking/${INSTALL_FILE})
if(EXISTS ${PRE_EXISTING_FILE})
message(STATUS "MPGD: Remove ${PRE_EXISTING_FILE}")
file(REMOVE_RECURSE ${PRE_EXISTING_FILE})
endif()
install(FILES compact/tracking/${FILE}
DESTINATION ${CMAKE_INSTALL_DATADIR}/${PROJECT_NAME}/compact/tracking
RENAME ${INSTALL_FILE}
)
endforeach()
message(STATUS "MPGD: \"pixel\" flavour of XMLs \"${MPGD_FILES}\" to be installed.")
ELSE()
file(GLOB_RECURSE MPGD_FILES
RELATIVE ${CMAKE_CURRENT_SOURCE_DIR}/compact/tracking
"compact/tracking/*_2DStrip.xml"
)
foreach(FILE IN ITEMS ${MPGD_FILES})
execute_process(COMMAND bash -c "echo ${FILE} | sed -e 's/_2DStrip//'" OUTPUT_VARIABLE INSTALL_FILE OUTPUT_STRIP_TRAILING_WHITESPACE)
set(PRE_EXISTING_FILE ${CMAKE_INSTALL_PREFIX}/${CMAKE_INSTALL_DATADIR}/${PROJECT_NAME}/compact/tracking/${INSTALL_FILE})
if(EXISTS ${PRE_EXISTING_FILE})
message(STATUS "MPGD: Remove ${PRE_EXISTING_FILE}")
file(REMOVE_RECURSE ${PRE_EXISTING_FILE})
endif()
install(FILES compact/tracking/${FILE}
DESTINATION ${CMAKE_INSTALL_DATADIR}/${PROJECT_NAME}/compact/tracking
RENAME ${INSTALL_FILE}
)
endforeach()
MESSAGE(STATUS "MPGD: \"2DStrip\" flavour of XMLs \"${MPGD_FILES}\" to be installed.")
ENDIF()

#-----------------------------------------------------------------------------------
# Install the detector calibration files.
Expand Down
151 changes: 151 additions & 0 deletions compact/tracking/mpgd_barrel_2DStrip.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,151 @@
<!-- SPDX-License-Identifier: LGPL-3.0-or-later -->
<!-- Copyright (C) 2022 Sylvester Joosten, Wouter Deconinck, Matt Posik, Yann Bedfer -->

<lccdd>
<info name="mpgd_barrel.xml"
title="Micro Pattern Gas Detectors"
author="mposik1983"
url="https://github.com/mposik1983"
status="development"
version="2.0"
><comment/>
</info>

<define>
<comment>
Inner MPGD tracking layer = CyMBaL
Note: MPGDs come in two distinct flavours: pixel and 2DStrip.
</comment>
<constant name="InnerMPGDBarrel_2DStrip" value="1"/>

<constant name="InnerMPGDBarrelLayer_length" value="InnerMPGDBarrel_length" />

<comment> Intrinsic parameters... </comment>
<comment> ...Thicknesses </comment>
<constant name="MMKaptonOverlay_thickness" value="50*um"/>
<constant name="MMCuGround_thickness" value="1.58*um"/>
<constant name="MMPCB_thickness" value="100*um"/>
<constant name="MMCuStrip_thickness" value="12*um"/>
<constant name="MMKaptonStrip_thickness" value="75*um"/>
<constant name="MMResistiveStrip_thickness" value="20*um"/>
<constant name="MMGas_thickness" value="128*um"/>
<constant name="MMMesh_thickness" value="18*um"/>
<constant name="MMGasGap_thickness" value="3000*um"/>
<constant name="MMDriftCuElectrode_thickness" value="5*um"/>
<constant name="MMDriftKapton_thickness" value="250*um"/>
<constant name="MMDriftCuGround_thickness" value="0.41*um"/>
<comment> FIXME: No support material is here, so fudge factor used to bring material budget to ~0.5% for barrel. </comment>
<constant name="MMFudgeInnerMPGDBarrel_thickness" value="570*um"/>
<comment> FIXME: No definite plan for connections/services to the inner sectors yet, so guess value. </comment>
<constant name="MMInnerService_thickness" value="0.05*cm"/>
<comment> 2DStrip: Multiple Sensitive Volumes. => Subdivide GasGap: 3 thin slices (pixel and (p|n)strips) at the centre, enclosed by 2 thick radiators. </comment>
<constant name="MMThinGap_thickness" value="10*um"/>
<constant name="MMThickGap_thickness" value="(MMGasGap_thickness - MMThinGap_thickness * 3) / 2"/>
<comment> ...Dimensions </comment>
<constant name="MMModuleWidth" value="46.0*cm"/>
<constant name="MMModuleLength" value="61.0*cm"/>
<constant name="MMFrameWidth" value="1.0*cm"/>
<constant name="MMOutwardFrameWidth" value="5.0*cm"/>
<constant name="MMOuterSector_R" value="57.7*cm"/>
<constant name="MMInnerSector_R" value="55.5*cm"/>
<comment> ...Positioning </comment>
<constant name="MMRadial_offset" value="1.0*cm"/>
<constant name="MMCentral_gap" value="0.2*cm"/>

<comment> In principle, stave count should not be independenty specified: it is fixed by above intrinsic parameters </comment>
<constant name="MPGDBarrelStave_count" value="8"/>
<constant name="InnerMPGDBarrelLayer_rmin" value="InnerMPGDBarrel_rmin"/>
<constant name="InnerMPGDBarrelLayer_rmax" value="MMOuterSector_R + MMRadial_offset / 2 + 1*cm"/>
<constant name="rmin2Sensor" value="MMDriftCuGround_thickness + MMDriftKapton_thickness + MMDriftCuElectrode_thickness + MMGasGap_thickness / 2" />
<constant name="MMInnerSensor_R" value="MMInnerSector_R + rmin2Sensor" />
<constant name="MMOuterSensor_R" value="MMOuterSector_R + rmin2Sensor" />
<constant name="MMInnerAperture" value="MMModuleWidth/MMInnerSensor_R" />
<constant name="MMOuterAperture" value="MMModuleWidth/MMOuterSensor_R" />
</define>

<detectors>
<detector
id="TrackerBarrel_2_ID"
name="InnerMPGDBarrel"
type="epic_CylinderMPGDBarrel"
readout="MPGDBarrelHits"
insideTrackingVolume="true"
vis="TrackerVis">
<type_flags type="DetType_TRACKER + DetType_BARREL"/>
<dimensions
width="MMModuleWidth"
length="MMModuleLength"/>

<module name="InnerMPGDBarrel_Mod1" vis="TrackerLayerVis" >
<comment> Models describe the overlapping of staves along phi, either (I) via offsets or (II) via distinct radii ("rmin1/2"). Solution (II) is only partially implemented: would require FOUR distinct "radius" in the readout segmentation. Note that in any case "rsensor" must be set equal the "radius" of the corresponding segmentation and is used internally to double-check the consistency between the stack of "module_component" and the segmentation "radius". </comment>
<model name="MMInnerSector" rmin1="MMInnerSector_R" rmin2="MMInnerSector_R" offset="MMRadial_offset" rsensor="MMInnerSensor_R"/>
<model name="MMOuterSector" rmin1="MMOuterSector_R" rmin2="MMOuterSector_R" offset="MMRadial_offset" rsensor="MMOuterSensor_R"/>
<frame name="MMFrame" material="CarbonFiber" width="MMFrameWidth" vis="TrackerSupportVis" />
<frame name="MMOutwardFrame" material="CarbonFiber" width="MMOutwardFrameWidth" vis="TrackerServiceVis" />
<service name="MMInnerService" material="Kapton" thickness="MMInnerService_thickness" vis="TrackerServiceVis" /> <comment> This includes cables => Kapton as a material is overoptimistic. </comment>
<comment> Going from the inside (sensitive) side to the readout side </comment>
<module_component name="DriftCuGround" thickness="MMDriftCuGround_thickness" material="Copper" vis="TrackerModuleVis"/>
<module_component name="DriftKapton" thickness="MMDriftKapton_thickness" material="Kapton" vis="TrackerVis" />
<module_component name="DriftCuElectrode" thickness="MMDriftCuElectrode_thickness" material="Copper" vis="TrackerVis"/>
<module_component name="InnerRadiator" thickness="MMThickGap_thickness" material="Ar90IsoButane" sensitive="True" vis="TrackerMPGDGasVis"/>
<module_component name="pStripThinGap" thickness="MMThinGap_thickness" material="Ar90IsoButane" sensitive="True" vis="TrackerMPGDGasVis"/>
<module_component name="pixelThinGap" thickness="MMThinGap_thickness" material="Ar90IsoButane" sensitive="True" vis="TrackerMPGDGasVis"/>
<module_component name="nStripThinGap" thickness="MMThinGap_thickness" material="Ar90IsoButane" sensitive="True" vis="TrackerMPGDGasVis"/>
<module_component name="OuterRadiator" thickness="MMThickGap_thickness" material="Ar90IsoButane" sensitive="True" vis="TrackerMPGDGasVis"/>
<module_component name="Mesh" thickness="MMMesh_thickness" material="MMGAS_InoxForMesh" vis="TrackerVis"/>
<module_component name="Fudge" thickness="MMFudgeInnerMPGDBarrel_thickness" material="Kapton" vis="TrackerVis"/>
<module_component name="Gas" thickness="MMGas_thickness" material="Ar90IsoButane" vis="TrackerMPGDGasVis"/>
<module_component name="ResistiveStrips" thickness="MMResistiveStrip_thickness" material="MMGAS_ResistivePaste" vis="TrackerVis"/>
<module_component name="KaptonStrips" thickness="MMKaptonStrip_thickness" material="Kapton" vis="TrackerVis" />
<module_component name="CuStrips" thickness="MMCuStrip_thickness" material="Copper" vis="TrackerVis" />
<module_component name="PCB" thickness="MMPCB_thickness" material="Fr4" vis="TrackerVis"/>
<module_component name="CuGround" thickness="MMCuGround_thickness" material="Copper" vis="TrackerVis"/>
<module_component name="KaptonOverlay" thickness="MMKaptonOverlay_thickness" material="Kapton" vis="TrackerMPGDVis"/> <!-- Instead of "TrackerSupportVis" in order to highligh frames -->
</module>

<layer module="InnerMPGDBarrel_Mod1" id="1" vis="InvisibleWithDaughters">
<barrel_envelope
inner_r="InnerMPGDBarrelLayer_rmin"
outer_r="InnerMPGDBarrelLayer_rmax"
z_length="InnerMPGDBarrelLayer_length"
z0="InnerMPGDBarrel_zoffset"/>
<layer_material surface="inner" binning="binPhi,binZ" bins0="MPGDBarrelStave_count*24" bins1="100" />
<layer_material surface="outer" binning="binPhi,binZ" bins0="MPGDBarrelStave_count*24" bins1="100" />
<rphi_layout phi0="0.0"/>
<z_layout gap="MMCentral_gap"/>
</layer>
</detector>
</detectors>

<plugins>
<plugin name="DD4hep_ParametersPlugin">
<argument value="InnerMPGDBarrel"/>
<argument value="layer_pattern: str=InnerMPGDBarrel_layer\d"/>
</plugin>
</plugins>

<readouts>
<readout name="MPGDBarrelHits">
<comment> 768 strips in phi yield a segmentation of 2pi/8/768 ~= 1mrad. Which in turn, yields ~160 um in resolution.
</comment>
<segmentation type="MultiSegmentation" key="sensor">
<segmentation name="Inner" type="MultiSegmentation" key_value="0" key="strip">
<segmentation name="InnerPix" type="CylindricalGridPhiZ" key_value="0" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R" />
<segmentation name="InnerPhi" type="CylindricalGridPhiZ" key_value="1" grid_size_phi="1*mrad" grid_size_z="MMModuleLength" radius="MMInnerSensor_R - MMThinGap_thickness" />
<segmentation name="InnerZ" type="CylindricalGridPhiZ" key_value="2" grid_size_phi="MMInnerAperture" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R + MMThinGap_thickness" />
<segmentation name="InnerBot" type="CylindricalGridPhiZ" key_value="3" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R - 3 * MMThinGap_thickness / 2 - MMThickGap_thickness / 2" />
<segmentation name="InnerTop" type="CylindricalGridPhiZ" key_value="4" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R + 3 * MMThinGap_thickness / 2 + MMThickGap_thickness / 2" />
</segmentation>
<segmentation name="Outer" type="MultiSegmentation" key_value="1" key="strip">
<segmentation name="OuterPix" type="CylindricalGridPhiZ" key_value="0" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R" />
<segmentation name="OuterPhi" type="CylindricalGridPhiZ" key_value="1" grid_size_phi="1*mrad" grid_size_z="MMModuleLength" radius="MMOuterSensor_R - MMThinGap_thickness" />
<segmentation name="OuterZ" type="CylindricalGridPhiZ" key_value="2" grid_size_phi="MMOuterAperture" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R + MMThinGap_thickness" />
<segmentation name="OuterTop" type="CylindricalGridPhiZ" key_value="3" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R - 3 * MMThinGap_thickness / 2 - MMThickGap_thickness / 2" />
<segmentation name="OuterBot" type="CylindricalGridPhiZ" key_value="4" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R + 3 * MMThinGap_thickness / 2 + MMThickGap_thickness / 2" />
</segmentation>
</segmentation>
<id>system:8,layer:4,module:12,sensor:2,strip:28:4,phi:-16,z:-16</id>
</readout>
</readouts>

</lccdd>
34 changes: 13 additions & 21 deletions compact/tracking/mpgd_barrel.xml → compact/tracking/mpgd_barrel_pixel.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,12 +11,10 @@

<define>
<comment>
Inner MPGD tracking layer(s)

Note: the inner and outer layers are implemented as separate detectors, as they
belong to different ACTS tracking volumes. If this restriction goes away
in the future they could be put together in a single tag.
Inner MPGD tracking layer = CyMBaL
Note: MPGDs come in two distinct flavours: pixel and 2DStrip.
</comment>
<constant name="InnerMPGDBarrel_2DStrip" value="0"/>

<constant name="InnerMPGDBarrelLayer_length" value="InnerMPGDBarrel_length" />

Expand All @@ -34,9 +32,9 @@
<constant name="MMDriftCuElectrode_thickness" value="5*um"/>
<constant name="MMDriftKapton_thickness" value="250*um"/>
<constant name="MMDriftCuGround_thickness" value="0.41*um"/>
<comment> FIXME: No support material is here, so fudge factor used to bring material budget to ~0.5% for barrel </comment>
<comment> FIXME: No support material is here, so fudge factor used to bring material budget to ~0.5% for barrel. </comment>
<constant name="MMFudgeInnerMPGDBarrel_thickness" value="570*um"/>
<comment> FIXME: No definite plan for connections/services to the inner sectors yet, so guess value </comment>
<comment> FIXME: No definite plan for connections/services to the inner sectors yet, so guess value. </comment>
<constant name="MMInnerService_thickness" value="0.05*cm"/>
<comment> ...Dimensions </comment>
<constant name="MMModuleWidth" value="46.0*cm"/>
Expand All @@ -51,6 +49,8 @@

<comment> In principle, stave count should not be independenty specified: it is fixed by above intrinsic parameters </comment>
<constant name="MPGDBarrelStave_count" value="8"/>
<constant name="InnerMPGDBarrelLayer_rmin" value="InnerMPGDBarrel_rmin"/>
<constant name="InnerMPGDBarrelLayer_rmax" value="MMOuterSector_R + MMRadial_offset / 2 + 1*cm"/>
<constant name="rmin2Sensor" value="MMDriftCuGround_thickness + MMDriftKapton_thickness + MMDriftCuElectrode_thickness + MMGasGap_thickness / 2" />
<constant name="MMInnerSensor_R" value="MMInnerSector_R + rmin2Sensor" />
<constant name="MMOuterSensor_R" value="MMOuterSector_R + rmin2Sensor" />
Expand Down Expand Up @@ -96,8 +96,8 @@

<layer module="InnerMPGDBarrel_Mod1" id="1" vis="InvisibleWithDaughters">
<barrel_envelope
inner_r="InnerMPGDBarrel_rmin"
outer_r="InnerMPGDBarrel_rmax"
inner_r="InnerMPGDBarrelLayer_rmin"
outer_r="InnerMPGDBarrelLayer_rmax"
z_length="InnerMPGDBarrelLayer_length"
z0="InnerMPGDBarrel_zoffset"/>
<layer_material surface="inner" binning="binPhi,binZ" bins0="MPGDBarrelStave_count*24" bins1="100" />
Expand All @@ -120,19 +120,11 @@
<comment> 768 strips in phi yield a segmentation of 2pi/8/768 ~= 1mrad. Which in turn, yields ~160 um in resolution.
</comment>
<segmentation type="MultiSegmentation" key="sensor">
<comment> Strip segmentations ("strip" field !=0) are not used at simulation time but only in digitization at reconstruction time, see digitization class, "MPGDTrackerDigi". </comment>
<segmentation name="Inner" type="MultiSegmentation" key_value="0" key="strip">
<segmentation name="InnerPix" type="CylindricalGridPhiZ" key_value="0" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R" />
<segmentation name="InnerPhi" type="CylindricalGridPhiZ" key_value="1" grid_size_phi="1*mrad" grid_size_z="MMModuleLength" radius="MMInnerSensor_R" />
<segmentation name="InnerZ" type="CylindricalGridPhiZ" key_value="2" grid_size_phi="MMInnerAperture" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R" />
</segmentation>
<segmentation name="Outer" type="MultiSegmentation" key_value="1" key="strip">
<segmentation name="OuterPix" type="CylindricalGridPhiZ" key_value="0" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R" />
<segmentation name="OuterPhi" type="CylindricalGridPhiZ" key_value="1" grid_size_phi="1*mrad" grid_size_z="MMModuleLength" radius="MMOuterSensor_R" />
<segmentation name="OuterZ" type="CylindricalGridPhiZ" key_value="2" grid_size_phi="MMOuterAperture" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R" />
</segmentation>
<segmentation name="InnerPix" type="CylindricalGridPhiZ" key_value="0" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMInnerSensor_R" />
<segmentation name="OuterPix" type="CylindricalGridPhiZ" key_value="1" grid_size_phi="1*mrad" grid_size_z="0.150*mm*sqrt(12)" radius="MMOuterSensor_R" />
</segmentation>
<id>system:8,layer:4,module:12,sensor:2,strip:30:2,phi:-16,z:-16</id>
<comment> "strip" field is ineffective... but required. </comment>
<id>system:8,layer:4,module:12,sensor:2,strip:28:4,phi:-16,z:-16</id>
</readout>
</readouts>

Expand Down
Loading
Loading