Skip to content

Include APA mesh atteenuation factor in the PMT efficiencies - #435

Merged
fjnicolas merged 2 commits into
developfrom
feature/pmteff_withapamesh
Mar 22, 2024
Merged

Include APA mesh atteenuation factor in the PMT efficiencies#435
fjnicolas merged 2 commits into
developfrom
feature/pmteff_withapamesh

Conversation

@fjnicolas

@fjnicolas fjnicolas commented Mar 20, 2024 •
edited
Loading

Copy link
Copy Markdown
Contributor

Updating the Photon Detection Efficiencies after the new measurements for the X-ARAPUCAs (#412) we realized that the light attenuation due to the APA grounding mesh was not accounted for the PMTs. This component is not present in the GDML file for the SBND geometry (i.e. not accounted by the semi-analytic model corrections), so it needs to be included in the overal detection efficiency. This PR is to include its effect in the simulation of the PMT signals.

The mesh transparency is 85% [ref], but the light attenuation depends on the photons incident angle, as can be seen in the figure below

The distribution of incident angles for the two components VUV and visible are shown below

We see that the maximum of the distribution is at 42º for the VUV and at 37º for the visible light. A zoom of the attenuation factor around those values can be seen below

According to these results we will use in our simulations an attenuation factor of 0.8

@fjnicolas
fjnicolas marked this pull request as draft March 20, 2024 15:01
@fjnicolas fjnicolas changed the title Feature/pmteff withapamesh Include APA mesh atteenuation factor in the PMT efficiencies Mar 20, 2024
@fjnicolas

fjnicolas commented Mar 20, 2024 •
edited
Loading

Copy link
Copy Markdown
Contributor Author

As for validation, the current PDS MC paper draft reflects the current performance up to Reco1 using the efficiencies introduced by this PR. As for Reco2, since the CI will not really pick any possible drop in performaces due to this change, I made some validation plots for the relevant higher level variables that rely on PDS information (OpT0Finder, SimpleFlashMatch and CRUMBS). The plots compare current (black histograms) with the updated version including the mesh attenuation (red crosses)

  • CRUMBS score (neutrinos left, cosmics right):
  • OpT0Finder chi2 (neutrinos left, cosmics right):
  • OpT0Finder PEFractional (neutrinos left, cosmics right):
  • OpT0Finder PEMeasured (neutrinos left, cosmics right):

This one is obviously the one that is affected the most, since we are indeed reducing the amount of deteceted of photons by an overall scaling factor of 0.8. Since the plot above may be a bit hard to read, some numbers:
- Mean MeasPE (standard) : 58781.8
- Mean MeasPE (with mesh): 48628.4

  • SimpleFlashMatch score (neutrinos left, cosmics right):

The ratio plots look sensible to me, so we are probably good to go without having to re-train the input weights for the tools. I'm asking the main developers of these tools to review this! @henrylay97 @lynnt20 @Robert-Darby

@fjnicolas
fjnicolas marked this pull request as ready for review March 20, 2024 16:22
lynnt20
lynnt20 approved these changes Mar 20, 2024
View reviewed changes
Comment thread sbndcode/OpT0Finder/job/opt0finder_sbnd.fcl
# below is the UPDATED PMT EFFICIENCIES
sbnd_opt0_finder_one_to_many.FlashMatchConfig.PhotonLibHypothesis.VUVEfficiency: [0., 0., 0., 0., 0., 0., 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.12 , 0.12 , 0.12 ,0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 , 0.12 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021]
sbnd_opt0_finder_one_to_many.FlashMatchConfig.PhotonLibHypothesis.VISEfficiency: [0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0., 0., 0., 0., 0., 0., 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.19 , 0.19 , 0.19 , 0.19 , 0.19 , 0.19 , 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0., 0., 0., 0., 0., 0., 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0., 0., 0.,0., 0., 0., 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.19 , 0.19 , 0.19 ,0.19 , 0.19 , 0.19 , 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0., 0., 0.,0., 0., 0., 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0., 0., 0., 0., 0., 0.,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.19 , 0.19 , 0.19 , 0.19 , 0.19 , 0.19 ,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0., 0., 0., 0., 0., 0.,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0., 0., 0., 0., 0., 0., 0.014, 0.014, 0.014,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.014,0.19 , 0.19 , 0.19 , 0.19 , 0.19 , 0.19 , 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0.014, 0.014, 0.014, 0.014, 0.014, 0.014, 0.13 , 0.13 , 0.13 ,0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 , 0.13 ,0., 0., 0., 0., 0., 0.]
sbnd_opt0_finder_one_to_many.FlashMatchConfig.PhotonLibHypothesis.VUVEfficiency: [0., 0., 0., 0., 0., 0., 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0., 0., 0., 0., 0., 0.,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0., 0., 0.,0., 0., 0., 0., 0., 0., 0., 0., 0.,0., 0., 0., 0., 0., 0., 0.021, 0.021, 0.021,0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021, 0.021,0., 0., 0., 0., 0., 0., 0.096 , 0.096 , 0.096 ,0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 , 0.096 ,0.021, 0.021, 0.021, 0.021, 0.021, 0.021]

@lynnt20 lynnt20 Mar 20, 2024

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

thank you for going through all of these and changing them 🥰

henrylay97
henrylay97 approved these changes Mar 22, 2024
View reviewed changes

@henrylay97 henrylay97 left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Nice one Fran, thanks for all the validation. As discussed offline I'm happy with this PR going in. The only direct impact on CRUMBS will be through the flash measured PE variable from OpT0Finder, but this is not a crazy change. Doesn't require an immediate update, can just be rolled into the next training automatically :D

@fjnicolas

fjnicolas commented Mar 22, 2024

Copy link
Copy Markdown
Contributor Author

CI build and sequential tests can be found in this link. Expected changes in the product sizes for all the light related data products (generally a drop given the lower ScintPreScale and detection efficiency). The lower efficiency can be also seen in the OpT0 PE related variables:

473:   rec.slc[0].opt0.hypoPE differs: 109979.6328125 vs 87983.703125
474:   rec.slc[0].opt0.measPE differs: 123453.671875 vs 99095.1484375

@fjnicolas
fjnicolas merged commit dad921a into develop Mar 22, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

@ggamezdiego ggamezdiego ggamezdiego approved these changes
@henrylay97 henrylay97 henrylay97 approved these changes
@lynnt20 lynnt20 lynnt20 approved these changes
@Robert-Darby Robert-Darby Awaiting requested review from Robert-Darby

Assignees

No one assigned

Labels

None yet

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

4 participants

@fjnicolas @ggamezdiego @henrylay97 @lynnt20