First look at new GENIE syst CAFs LBL meeting 30 Jul 2018 Chris - PowerPoint PPT Presentation

First look at new GENIE syst CAFs LBL meeting – 30 Jul 2018 Chris Backhouse – University College London

Introduction ◮ Using example file /dune/data/users/marshalc/CAF FHC.root ◮ Weights now readable in CAFAna as dune.genie wgt[][] ◮ Systematic names/indices from GetGenieWeightNames() ◮ May not be final interface ◮ Current files contain 50 random universes per syst ◮ Making very crude plots – Ev reco and no selection C. Backhouse (UCL) Calibrator 1 / 10

Systs 1/8 RvpCC1pi MaCCQE × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvpCC1pi MaCCQE 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 2 / 10

Systs 2/8 RvpCC2pi RvpNC1pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvpCC2pi RvpNC1pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 3 / 10

Systs 3/8 RvnNC2pi RvnNC1pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvnNC2pi RvnNC1pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 4 / 10

Systs 4/8 RvbarpCC1pi RvbarpCC2pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvbarpCC1pi RvbarpCC2pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 5 / 10

Systs 5/8 RvbarnCC1pi RvbarnCC2pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvbarnCC1pi RvbarnCC2pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 6 / 10

Systs 6/8 RvbarpNC1pi RvbarpNC2pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvbarpNC1pi RvbarpNC2pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 7 / 10

Systs 7/8 RvbarnNC1pi RvbarnNC2pi × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) RvbarnNC1pi RvbarnNC2pi 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 8 / 10

Systs 8/8 BR1gamma BR1eta × 3 × 3 10 10 Events Events 300 300 200 200 100 100 0 0 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) BR1gamma BR1eta 1.1 1.1 Ratio to nominal Ratio to nominal 1.05 1.05 1 1 0.95 0.95 0.9 0.9 0 2 4 6 8 10 0 2 4 6 8 10 Reconstructed energy (GeV) Reconstructed energy (GeV) C. Backhouse (UCL) Calibrator 9 / 10

Conclusion ◮ Can succesfully read weights ◮ Why are they > 1 for all systs except BR1eta ? ◮ To fit pull terms for individual systs need systematic stride through the parameter, not random sampling C. Backhouse (UCL) Calibrator 10 / 10