Event Energy Reconstruction Beatriz Tapia 03/26/2019 Big Thanks to - - PowerPoint PPT Presentation
Event Energy Reconstruction Beatriz Tapia 03/26/2019 Big Thanks to - - PowerPoint PPT Presentation
Event Energy Reconstruction Beatriz Tapia 03/26/2019 Big Thanks to Michel and Jose! 3 GeV nueCC events distributed through the entire volume Jose generated events with the latest photon library and energy deposits From Edep [MeV] in
❖ 3 GeV nueCC events distributed through the entire volume ❖ Jose generated events with the latest photon library and energy deposits From Edep [MeV] in one voxel to detected by one pmt: Edep(voxel) * 24000 photons/MeV * 0.12 * visibility(voxel, pmt) = PE(pmt) Ereco(voxel) = PE(pmt) / (24000 photons/MeV * 0.12 * visibility(voxel, pmt))
Ereco = Σ_voxel [Edep(voxel)] = Σ_pmt [Σ_voxel [PE(pmt) * weight(voxel) / (24000 photons/MeV * 0.12 * visibility(voxel, pmt))]] Sum over PMTs with PEs and voxels with Edep, with weight(voxel) = Edep(voxel) / Total Edep
2
- To 1st approximation -> assume all the energy deposit happened at int-vertex-voxel
- 770 events, only 561 events that pass fiducial cut (int-vertex 0.5 m away from boundary)
- 1st sum and then divide -> treat all PMTs as being only one
Ereco = (Σ_pmt [PE(pmt)] / Σ_pmt [visibility(int-vertex-voxel, pmt)])/(24000 photons/MeV * 0.12)
3
Ereco [GeV] Ereco [GeV] TrueX [cm]
Effect of Fiducial Cut on Edep
4
Edep [GeV] Edep [GeV]
Backups
5
6
Used reconstructed information: OpHits related to MC Truth via EventID
7
- To 2nd approximation -> assume all the energy deposit happened at int-vertex-voxel
- 770 events, only 561 events that pass fiducial cut (0.5 m away from boundary)
- 1st divide then sum
Ereco = (Σ_pmt [PE(pmt) / visibility(int-vertex-voxel, pmt)])/(24000 photons/MeV * 0.12)
Ereco [GeV] TrueX [cm] Ereco [GeV]