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 - - 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[][]
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
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events MaCCQE
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvpCC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal MaCCQE
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvpCC1pi
- C. Backhouse (UCL)
Calibrator 2 / 10
Systs 2/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvpCC2pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvpNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvpCC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvpNC1pi
- C. Backhouse (UCL)
Calibrator 3 / 10
Systs 3/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvnNC2pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvnNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvnNC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvnNC1pi
- C. Backhouse (UCL)
Calibrator 4 / 10
Systs 4/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarpCC1pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarpCC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarpCC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarpCC2pi
- C. Backhouse (UCL)
Calibrator 5 / 10
Systs 5/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarnCC1pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarnCC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarnCC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarnCC2pi
- C. Backhouse (UCL)
Calibrator 6 / 10
Systs 6/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarpNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarpNC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarpNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarpNC2pi
- C. Backhouse (UCL)
Calibrator 7 / 10
Systs 7/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarnNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events RvbarnNC2pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarnNC1pi
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal RvbarnNC2pi
- C. Backhouse (UCL)
Calibrator 8 / 10
Systs 8/8
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events BR1gamma
2 4 6 8 10
Reconstructed energy (GeV)
100 200 300
3
10 ×
Events BR1eta
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal BR1gamma
2 4 6 8 10
Reconstructed energy (GeV)
0.9 0.95 1 1.05 1.1
Ratio to nominal BR1eta
- 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