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 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) Ereco [GeV] Ereco [GeV] TrueX [cm] 3
Effect of Fiducial Cut on Edep Edep [GeV] Edep [GeV] 4
Backups 5
Used reconstructed information: OpHits related to MC Truth via EventID 6
● 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] Ereco [GeV] TrueX [cm] 7
Recommend
More recommend