Shower Characterisation James Pillow 19/12/2019 EM Task Force - PowerPoint PPT Presentation

Shower Characterisation James Pillow 19/12/2019 EM Task Force

Events • MC events from ProtoDUNE SP Production 2 with space charge. • Using 3D position from shower calorimetry object. • Using energy estimation from method I’ve talked about before (in back up slides) James Pillow | 19/12/2019 | EM Task Force 2

Previously on EM Task Force… Peak too high James Pillow | 19/12/2019 | EM Task Force 3

dE/dx • Move to using the shower calorimetry 3D points (still using my energy estimator for each point though) • Peak still slightly too high… James Pillow | 19/12/2019 | EM Task Force 4

dE/dx • Aaron pointed out that the expected approximate recomb correction should be 0.7. • I had been using 0.63. James Pillow | 19/12/2019 | EM Task Force 5

dE/dx • Change recombination correction to 0.7: • Get peak at 2.1MeV/cm! James Pillow | 19/12/2019 | EM Task Force 6

Energy Estimation • Increasing recombination makes total energy estimation worse. James Pillow | 19/12/2019 | EM Task Force 7

Energy Estimation from ADC • Made a SCE On sample of ~0.5- >~7.5 GeV electrons (~100,000). • Relationship between ADC and energy seems very linear James Pillow | 19/12/2019 | EM Task Force 8

Energy Estimation from ADC • Use a robust linear estimator to do ADC->Energy • Also take gradient from fit, but force an (0,0) intercept • Pushes peak to 0 on the right. • Less symmetrical distribution though. James Pillow | 19/12/2019 | EM Task Force 9

Longitudinal Profile • Longitudinal profile of EM shower described by: !# = 𝐹 & 𝑐 (# )*+ , *-. • !" /(1) From last talk. Plot on top was simply plotting the wrong thing James Pillow | 19/12/2019 | EM Task Force 10

Transverse Profile • Transverse profile looks as I expect • First bin is low simply 90% Energy from the fact the 3D points won’t necessarily be right next to the primary axis of the shower. James Pillow | 19/12/2019 | EM Task Force 11

Back-up James Pillow 19/12/2019 EM Task Force

Energy Estimation – ‘Traditional’ • I estimate the energy of each hit using: ( caloAlg.ElectronsFromADCArea( hit->Integral(), plane) * caloAlg.LifetimeCorrection(hit->PeakTime()) ) / kGeVtoElectrons * recombination Recombination correction CalorimetryAlg.cxx function to estimate Lifetime correction Convert number of number of electrons based on hit charge electrons to GeV (from (average value from modbox from Nick) physical constant) • I then total the energy estimation of each hit. • I also then add together the total energy estimation for all particles in the PFParticle hierarchy. James Pillow | 19/12/2019 | EM Task Force 13 13

Energy Estimator ADC Fits James Pillow | 19/12/2019 | EM Task Force 14