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

▶
shower characterisation
SMART_READER_LITE
LIVE PREVIEW

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

Nov 07, 2022 413 likes •567 views

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 Ive talked

slide-1
SLIDE 1

James Pillow

Shower Characterisation

19/12/2019 EM Task Force

slide-2
SLIDE 2

James Pillow

  • 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)

| 19/12/2019 | 2 EM Task Force

Events

slide-3
SLIDE 3

James Pillow | 19/12/2019 | 3 EM Task Force

Previously on EM Task Force…

Peak too high

slide-4
SLIDE 4

James Pillow | 19/12/2019 | 4 EM Task Force

  • Move to using the shower

calorimetry 3D points (still using my energy estimator for each point though)

  • Peak still slightly too high…

dE/dx

slide-5
SLIDE 5

James Pillow | 19/12/2019 | 5 EM Task Force

  • Aaron pointed out that the

expected approximate recomb correction should be 0.7.

  • I had been using 0.63.

dE/dx

slide-6
SLIDE 6

James Pillow | 19/12/2019 | 6 EM Task Force

  • Change recombination

correction to 0.7:

  • Get peak at 2.1MeV/cm!

dE/dx

slide-7
SLIDE 7

James Pillow | 19/12/2019 | 7 EM Task Force

  • Increasing recombination makes

total energy estimation worse.

Energy Estimation

slide-8
SLIDE 8

James Pillow

  • Made a SCE On sample of ~0.5-

>~7.5 GeV electrons (~100,000).

  • Relationship between ADC and

energy seems very linear

| 19/12/2019 | 8 EM Task Force

Energy Estimation from ADC

slide-9
SLIDE 9

James Pillow

  • 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.

| 19/12/2019 | 9 EM Task Force

Energy Estimation from ADC

slide-10
SLIDE 10

James Pillow

  • Longitudinal profile of EM

shower described by:

  • !"

!# = 𝐹&𝑐 (#)*+,*-. /(1)

| 19/12/2019 | 10 EM Task Force

Longitudinal Profile

From last talk. Plot on top was simply plotting the wrong thing

slide-11
SLIDE 11

James Pillow

  • Transverse profile

looks as I expect

  • First bin is low simply

from the fact the 3D points won’t necessarily be right next to the primary axis of the shower.

| 19/12/2019 | 11 EM Task Force

Transverse Profile

90% Energy

slide-12
SLIDE 12

James Pillow

Back-up

19/12/2019 EM Task Force

slide-13
SLIDE 13

James Pillow | 19/12/2019 | 13 EM Task Force

Energy Estimation – ‘Traditional’

  • I estimate the energy of each hit using:
  • 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.

13

( caloAlg.ElectronsFromADCArea( hit->Integral(), plane) * caloAlg.LifetimeCorrection(hit->PeakTime()) ) / kGeVtoElectrons * recombination

CalorimetryAlg.cxx function to estimate number of electrons based on hit charge Lifetime correction Convert number of electrons to GeV (from physical constant) Recombination correction (average value from modbox from Nick)

slide-14
SLIDE 14

James Pillow | 19/12/2019 | 14 EM Task Force

Energy Estimator ADC Fits