Add proton light quenching#879
Conversation
|
trigger build LArSoft/lar*@LARSOFT_SUITE_v10_12_02 SBNSoftware/sbndaq-artdaq-core@v1_10_06 SBNSoftware/sbn*@SBN_SUITE_v10_12_02 |
|
✔️ CI build for LArSoft Succeeded on slf7 for c14:prof -- details available through the CI dashboard |
|
✔️ CI build for LArSoft Succeeded on slf7 for e26:prof -- details available through the CI dashboard |
|
❌ CI build for SBND Failed at phase build SBND on slf7 for c14:prof -- details available through the CI dashboard 🚨 For more details about the failed phase, check the build SBND phase logs parent CI build details are available through the CI dashboard |
|
🚨 For more details about the warning phase, check the ci_tests SBND phase logs parent CI build details are available through the CI dashboard |
|
approved |
Description
Please provide a detailed description of the changes this pull request introduces.
This PR implements a light quenching factor for protons. The currently used modified box model is in agreement with experimental data (Ref Doke et al., Jpn. J. Appl. Phys. Vol. 41 (2002) pp. 1538–1545.) for MIPs but not for protons. The dotted line from the figure below shows the modified model prediction in the absence of electric field, which agrees with Figure 2 from the previous reference for MIPs but not for protons. Assuming the same light attenuation factor due to the electric field for protons and electrons (~0.6), the expected relative scintillation yield (L/L0) for protons can be calculated as the value quoted from the reference ~0.81 corrected by the electric field attenuation factor, that is L/L0 ~ 0.49.

From the figure below, we can see that the relative scintillation yield currently predicted by the modified box model is ~0.7, so we need an extra quenching factor of 0.7 to match the expected L/L0 ~ 0.49.
Checklist
Reviewers,AssigneesDevelopementRelevant PR links (optional)
Does this PR require merging another PR in a different repository (such as sbnanobj/sbnobj etc.)?
This PR needs LArSoft/larsim#163 to be merged.
Link(s) to docdb describing changes (optional)
Is there a docdb describing the issue this solves or the feature added?
More information can be found here https://sbn-docdb.fnal.gov/cgi-bin/sso/ShowDocument?docid=42966