Search/Navigation:
Related:
SLAC/EPP/HPS Public
Jefferson Lab/Hall B/HPS Run Wiki
S30XL-LESA/LDMX
Instructions on how to make cluster amplitude plots and hit time plots.
For a data file, run the following command on an EVIO calibration file. You can configure the following cleanup cuts in the steering file (default is on) - <trackQualityCut> requires tracks to have a Chi2 better than <maxChi2> (default is 10) and 6-hit tracks; <feeCut> requires electrons of momentum between ebeam*<minFee> (default is 0.85) and ebeam*<maxFee> (default is 1.15).
Code Block |
---|
java -jar <jar file> hps-java/steering-files/src/main/resources/org/hps/steering/analysis/SvtClusterPlots.lcsim <input recon file> -DoutputFile=<output file> |
These commands will output a root file with both cluster charge (amplitude) and signal to noise for single hit clusters, multiple hit clusters, and hits on tracks that pass the desired quality cuts. There are also cluster hit times, single hit clusters on track hit times (used in the analysis), and several other plots.
In order to analyze, you need to clone the HPS-CODE repository.
Code Block |
---|
git clone https://github.com/JeffersonLab/HPS-CODE |
Run the following command on the output of the previous java command to produce plots.
Use the -m option to fit the multiple hit cluster plots (default is hits on track plots). Use the -c option to fit cluster amplitude and signal to noise plots with a landau convoluted with a Gaussian (default is just a landau).
Code Block |
---|
python HPS-CODE/ANALYSIS/SVT/MakeClusterPlots.py <output file base name> <input root file> |
Warning: If you try to fit using the landau-Gaussian convolution (which is the correct way), you may have to manually adjust fitting ranges in the code or the fitting may get stuck.
Output:
<output file base name>_sigtonoise.pdf - fitted signal to noise at each sensor and the most probable value at each sensor
<output file base name>_clusteramplitude.pdf - fitted cluster amplitude at each sensor and the most probable value at each sensor
<output file base name>_hittime.pdf - fitted hit times (for single hit clusters on track) at each sensor and the fitted sigma of the hit time at each sensor
<output file base name>.root - same plots as the three pdf files above but as root histograms
Order of "Sensor Number" is as follows (starting at 0):
0. module_L1b_halfmodule_axial_sensor0
1. module_L1b_halfmodule_stereo_sensor0
2. module_L1t_halfmodule_axial_sensor0
3. module_L1t_halfmodule_stereo_sensor0
4. module_L2b_halfmodule_axial_sensor0
5. module_L2b_halfmodule_stereo_sensor0
6. module_L2t_halfmodule_axial_sensor0
7. module_L2t_halfmodule_stereo_sensor0
8. module_L3b_halfmodule_axial_sensor0
9. module_L3b_halfmodule_stereo_sensor0
10. module_L3t_halfmodule_axial_sensor0
11. module_L3t_halfmodule_stereo_sensor0
12. module_L4b_halfmodule_axial_hole_sensor0
13. module_L4b_halfmodule_axial_slot_sensor0
14. module_L4b_halfmodule_stereo_hole_sensor0
15. module_L4b_halfmodule_stereo_slot_sensor0
16. module_L4t_halfmodule_axial_hole_sensor0
17. module_L4t_halfmodule_axial_slot_sensor0
18. module_L4t_halfmodule_stereo_hole_sensor0
19. module_L4t_halfmodule_stereo_slot_sensor0
20. module_L5b_halfmodule_axial_hole_sensor0
21. module_L5b_halfmodule_axial_slot_sensor0
22. module_L5b_halfmodule_stereo_hole_sensor0
23. module_L5b_halfmodule_stereo_slot_sensor0
24. module_L5t_halfmodule_axial_hole_sensor0
25. module_L5t_halfmodule_axial_slot_sensor0
26. module_L5t_halfmodule_stereo_hole_sensor0
27. module_L5t_halfmodule_stereo_slot_sensor0
28. module_L6b_halfmodule_axial_hole_sensor0
29. module_L6b_halfmodule_axial_slot_sensor0
30. module_L6b_halfmodule_stereo_hole_sensor0
31. module_L6b_halfmodule_stereo_slot_sensor0
32. module_L6t_halfmodule_axial_hole_sensor0
33. module_L6t_halfmodule_axial_slot_sensor0
34. module_L6t_halfmodule_stereo_hole_sensor0
35. module_L6t_halfmodule_stereo_slot_sensor0
For SVT NIM Paper
Files:
Cuts: For hits on track - 6-hit tracks, trackChi2 < 10, 0.85*ebeam < trackP < 1.15*ebeam, and trackOmega > 0