Skip to content

Low-Q2 reconstruction benchmark improvements #239

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

Draft
wants to merge 11 commits into
base: master
Choose a base branch
from
Draft

Low-Q2 reconstruction benchmark improvements #239

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
a603c6f
Add tagger acceptance plots
simonge Nov 21, 2025
3a4b554
Add and fix some kinematics and add basic script to plot x-q2 accepta…
simonge Nov 21, 2025
03321e7
Swap global angle with scattering angle
simonge Nov 21, 2025
1222a29
Add another formatting script and tweek histograms
simonge Nov 21, 2025
63699e7
Tag formatting scripts as such.
simonge Nov 21, 2025
b753b58
Add Energy-Q2 formatting
simonge Nov 21, 2025
6c380b1
Lots of formatting changes, swapping x/y axis adding new plots, movin…
simonge Nov 23, 2025
4e5c252
Formatting and workflow improvement
simonge Nov 24, 2025
4526da5
Add line showing smallest resolvable Q2 region
simonge Nov 25, 2025
5a985f6
remove grid
simonge Nov 25, 2025
e7be324
Merge remote-tracking branch 'origin/master' into low-q2-reconstructi…
simonge Nov 27, 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
153 changes: 103 additions & 50 deletions benchmarks/lowq2_reconstruction/Snakefile
View file
Open in desktop

Large diffs are not rendered by default.

138 changes: 138 additions & 0 deletions benchmarks/lowq2_reconstruction/epic_ip6_extended_10x100.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,138 @@
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">

<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>

<documentation level="-1">
# EPIC Detector
- https://github.com/eic/epic
- https://github.com/eic/ip6
</documentation>

<!-- Some information about detector -->
<info name="EPIC Detector" title="EPIC Detector"
author="EPIC Collaboration"
url="https://github.com/eic/epic"
status="development"
version="v1 2021-03-16">
<comment> EPIC </comment>
</info>

<!-- Standard environmental conditions: 20 °C, 1 atm -->
<std_conditions type="NTP"/>

<define>
<documentation level="2">
## Main Constant Definitions

The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="${DETECTOR_PATH}/compact/fields/beamline_10x100.xml" />
<include ref="${DETECTOR_PATH}/compact/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/far_forward/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/far_backward/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/version.xml"/>
<include ref="${DD4hepINSTALL}/DDDetectors/compact/detector_types.xml"/>
</define>


<includes>
<gdmlFile ref="${DD4hepINSTALL}/DDDetectors/compact/elements.xml"/>
<gdmlFile ref="${DETECTOR_PATH}/compact/materials.xml"/>
</includes>
<include ref="${DETECTOR_PATH}/compact/optical_materials.xml"/>

<include ref="${DETECTOR_PATH}/compact/calibrations.xml"/>

<limits>
<limitset name="world_limits">
</limitset>
<limitset name="beamline_limits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="kill_limits">
<limit name="time_max" particles="*" value="0.0" unit="ns"/>
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>

<regions>
<region name="world_region" eunit="MeV" lunit="mm" cut="0.001" threshold="0.001">
<limitsetref name="world_limits"/>
</region>
</regions>

<display>
<include ref="${DETECTOR_PATH}/compact/colors.xml"/>
<include ref="${DETECTOR_PATH}/compact/display.xml"/>
</display>

<documentation level="0">
## World Volume

The world is a simple box, but could be a union of multiple regions.
</documentation>
<world material="Air">
<shape type="Box" dx="world_dx" dy="world_dy" dz="world_dz"/>
<!--regionref name="world_region"/-->
<!--limitsetref name="world_limits"/-->
</world>

<documentation level="0">
## Detector Subsystems

### IP Subsystems

The interaction point subsystems are included before the central detector subsystems.
This is becuase the IP subsystems, for example the beampipe, will define paramters
which are subsquently used in the central detector construction -- e.g. the vertex tracker
uses the beampipe OD to help define its placement.

The IP subsystems include the Far forward and backward regions. The list of subsystem includes:
- Interaction region beampipe
- B0 tracker
- Off-momentum tracker
- Far forward roman pots
- Zero Degree Calorimeter
- Beam line magnets.
- and more...
</documentation>


<documentation level="10">
## Central tracking detectors
</documentation>

<include ref="${DETECTOR_PATH}/compact/tracking/definitions_craterlake.xml"/>

<documentation level="11">
## Central beam pipe
</documentation>
<include ref="${DETECTOR_PATH}/compact/central_beampipe.xml"/>

<documentation level="11">
## Far foward detectors
</documentation>
<include ref="${DETECTOR_PATH}/compact/far_forward/default.xml"/>

<documentation level="11">
## Far backward detectors
</documentation>
<include ref="${DETECTOR_PATH}/compact/far_backward/extended.xml"/>
</lccdd>
138 changes: 138 additions & 0 deletions benchmarks/lowq2_reconstruction/epic_ip6_extended_18x275.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,138 @@
<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
xmlns:xs="http://www.w3.org/2001/XMLSchema"
xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">

<debug>
<type name="surface" value="0"/>
<type name="material" value="0"/>
<type name="readout" value="0"/>
<type name="segmentation" value="0"/>
<type name="limits" value="0"/>
<type name="region" value="0"/>
<type name="includes" value="0"/>
</debug>

<documentation level="-1">
# EPIC Detector
- https://github.com/eic/epic
- https://github.com/eic/ip6
</documentation>

<!-- Some information about detector -->
<info name="EPIC Detector" title="EPIC Detector"
author="EPIC Collaboration"
url="https://github.com/eic/epic"
status="development"
version="v1 2021-03-16">
<comment> EPIC </comment>
</info>

<!-- Standard environmental conditions: 20 °C, 1 atm -->
<std_conditions type="NTP"/>

<define>
<documentation level="2">
## Main Constant Definitions

The ip6 (or other ip) defines should be included first.
These files have only a define tags.
</documentation>
<include ref="${DETECTOR_PATH}/compact/fields/beamline_18x275.xml" />
<include ref="${DETECTOR_PATH}/compact/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/far_forward/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/far_backward/definitions.xml" />
<include ref="${DETECTOR_PATH}/compact/version.xml"/>
<include ref="${DD4hepINSTALL}/DDDetectors/compact/detector_types.xml"/>
</define>


<includes>
<gdmlFile ref="${DD4hepINSTALL}/DDDetectors/compact/elements.xml"/>
<gdmlFile ref="${DETECTOR_PATH}/compact/materials.xml"/>
</includes>
<include ref="${DETECTOR_PATH}/compact/optical_materials.xml"/>

<include ref="${DETECTOR_PATH}/compact/calibrations.xml"/>

<limits>
<limitset name="world_limits">
</limitset>
<limitset name="beamline_limits">
<limit name="step_length_max" particles="*" value="1.0" unit="mm" />
<limit name="track_length_max" particles="*" value="1.0" unit="mm" />
<limit name="time_max" particles="*" value="0.1" unit="ns" />
<limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
<limit name="range_min" particles="*" value="0.1" unit="mm" />
</limitset>
<limitset name="kill_limits">
<limit name="time_max" particles="*" value="0.0" unit="ns"/>
</limitset>
<limitset name="cal_limits">
<limit name="step_length_max" particles="*" value="5.0" unit="mm"/>
</limitset>
</limits>

<regions>
<region name="world_region" eunit="MeV" lunit="mm" cut="0.001" threshold="0.001">
<limitsetref name="world_limits"/>
</region>
</regions>

<display>
<include ref="${DETECTOR_PATH}/compact/colors.xml"/>
<include ref="${DETECTOR_PATH}/compact/display.xml"/>
</display>

<documentation level="0">
## World Volume

The world is a simple box, but could be a union of multiple regions.
</documentation>
<world material="Air">
<shape type="Box" dx="world_dx" dy="world_dy" dz="world_dz"/>
<!--regionref name="world_region"/-->
<!--limitsetref name="world_limits"/-->
</world>

<documentation level="0">
## Detector Subsystems

### IP Subsystems

The interaction point subsystems are included before the central detector subsystems.
This is becuase the IP subsystems, for example the beampipe, will define paramters
which are subsquently used in the central detector construction -- e.g. the vertex tracker
uses the beampipe OD to help define its placement.

The IP subsystems include the Far forward and backward regions. The list of subsystem includes:
- Interaction region beampipe
- B0 tracker
- Off-momentum tracker
- Far forward roman pots
- Zero Degree Calorimeter
- Beam line magnets.
- and more...
</documentation>


<documentation level="10">
## Central tracking detectors
</documentation>

<include ref="${DETECTOR_PATH}/compact/tracking/definitions_craterlake.xml"/>

<documentation level="11">
## Central beam pipe
</documentation>
<include ref="${DETECTOR_PATH}/compact/central_beampipe.xml"/>

<documentation level="11">
## Far foward detectors
</documentation>
<include ref="${DETECTOR_PATH}/compact/far_forward/default.xml"/>

<documentation level="11">
## Far backward detectors
</documentation>
<include ref="${DETECTOR_PATH}/compact/far_backward/extended.xml"/>
</lccdd>
Loading