Skip to content

Add v15 ROC calculation script for issue #100 - #123

Open
1SelenaRao wants to merge 1 commit into
IncandelaLab:masterfrom
1SelenaRao:roc-studies-v15
Open

Add v15 ROC calculation script for issue #100#123
1SelenaRao wants to merge 1 commit into
IncandelaLab:masterfrom
1SelenaRao:roc-studies-v15

Conversation

@1SelenaRao

Copy link
Copy Markdown

Adds roc_calcs_v15.py for issue #100 — computes ECal radius-of-containment
for v15 geometry (400μm Si), at both 68% and 95% containment.

  • Updated ecal_layerZs to v15 values from ldmx-sw EcalGeometry.v15()
  • Points at the v4.5.4 signal samples under /home/vamitamas/Samples8GeV/v4.5.4_signal/
  • Fixed pass name (sim, not signal) to match actual branch names in the v15 files
  • roc_frac is now a CLI arg; ran both 0.68 and 0.95 via the included slurm job scripts
  • Full run: ~11M events, ~80% yield after standard cuts (no clean SP hit, angle overflow, missing energy/distance data)

Closes #100

- Copy roc_calcs_v14.py to roc_calcs_v15.py
- Update ecal_layerZs to v15 geometry (400um Si) per
  ldmx-sw EcalGeometry.v15()
- Point file_paths at v4.5.4_signal samples
  (/home/vamitamas/Samples8GeV/v4.5.4_signal/)
- Fix pass_name from 'signal' to 'sim' to match actual
  branch names in v15 production files
- Make roc_frac a CLI arg to support both 0.68 and 0.95
  containment runs
- Add slurm job scripts for both containment fractions
dist_energy[layer_index].append( (dist_ele_traj, edep) )
debugPrint(f' Hit appended to dist_energy at {layer_index = }; {dist_ele_traj = }, {edep = }')

roc_values = np.zeros(34)

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

we actually have 32 layers in the v14 geom

debugPrint(f' Event has {bin_index = }', debug)

dist_energy = defaultdict(list)
energy_total = np.zeros(34)

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

we actually have 32 layers in the v14 geom


ang_bins = [0, 10, 15, 25, 30, 40, 50, 60, 70]
num_ang_bins = len(ang_bins) - 1
num_ecal_layers = 34

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

same here

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Have a RoC with 95% containment

2 participants