470 likes | 761 Views
Combined Atmospheric Analysis: Study of 6 Month’s Data. Andy Blake Cambridge University Wednesday June 13 th 2007. Overview. Developing a combined atmospheric neutrino analysis: – Contained vertex muon neutrinos (John Chapman). – Upward-going muons (Brian Rebel, John Chapman).
E N D
Combined Atmospheric Analysis: Study of 6 Month’s Data. Andy Blake Cambridge University Wednesday June 13th 2007
Overview • Developing a combined atmospheric neutrino analysis: • – Contained vertex muon neutrinos (John Chapman). • – Upward-going muons (Brian Rebel, John Chapman). • – Contained vertex electron neutrinos (Ben Speakman). • Currently developing event selection code: • – All track and shower selection variables now calculated • in AtNuOutput classes and stored in AtNuEvent classes. • – Event selection implemented in AtNuAna package. • – Analysis variables written to short AtNuAna Ntuples. • Have now completed study using six month data sample. Andy Blake, Cambridge University Combined Analysis, slide 2
Event Selection Outline of event selection: DATA / MC Data: HV, coil, GPS, LI, bad chips etc. MC: large scatters. DATA QUALITY CHECKS LI rejection, clean planes, “fiducial” event GOOD EVENTS John and Brian’s track selection Ben’s shower selection GOOD TRACKS GOOD SHOWERS UPWARD-GOING MUONS CONTAINED VERTEX ELECTRON NEUTRINOS Andy Blake, Cambridge University Combined Analysis, slide 3
Event Types typedef enum EEventType { kNothing = 0x00, kGoodSlice = 0x01, kGoodEvent = 0x02, kGoodTrack = 0x04, kGoodShower = 0x08, kCV = 0x10, kCE = 0x20, kUP = 0x40, kFC = 0x80, kPC = 0x100, kPCDN = 0x200, kPCUP = 0x400, kUPMU = 0x800, kCosmic = 0x1000, kStoppingMuon = 0x2000, kThroughGoingMuon = 0x4000, kMultipleMuon = 0x8000, kVeto = 0x10000, kSpill = 0x20000, kBeamSpill = 0x40000, kFakeSpill = 0x80000, kGoodDirection = 0x100000, kGoodCharge = 0x200000, kPositiveCharge = 0x400000, kNegativeCharge = 0x800000, kAtmosNumu = 0x1000000, kAtmosNumuCV = 0x2000000, kAtmosNumuUP = 0x4000000, kAtmosNue = 0x8000000, kLIevent = 0x10000000, kMCscatter = 0x20000000 } EventType_t; typedef enum EAtmosEventType { kNothing = 0x00, kGoodEvent = 0x01, kSpill = 0x02, kBeamSpill = 0x04, kVeto = 0x08, kCV = 0x10, kUP = 0x20, kGoodTrack = 0x40, kGoodShower = 0x80, kAtmosNumu = 0x100, kAtmosNumuCV = 0x200, kAtmosNumuUP = 0x400, kAtmosNue = 0x800, kAtmosNumuCVcut1 = 0x1000, kAtmosNumuCVcut2 = 0x2000, kAtmosNumuCVcut3 = 0x4000, kAtmosNumuCVcut4 = 0x8000, kAtmosNumuUPcut1 = 0x10000, kAtmosNumuUPcut2 = 0x20000, kAtmosNumuUPcut3 = 0x40000, kAtmosNumuUPcut4 = 0x80000, kAtmosNueCVcut1 = 0x100000, kAtmosNueCVcut2 = 0x200000, kAtmosNueCVcut3 = 0x400000, kAtmosNueCVcut4 = 0x800000 } AtmosEventType_t; containment Atmospheric Neutrino Event Types containment cuts cosmics Event Types veto shield, beam spills upward going cuts good direction, good charge electron neutrino cuts Above event types defined by methods in AtNuAna class. light injection MC scatter Andy Blake, Cambridge University Combined Analysis, slide 4
6 Month Study • Data. • – runs 31812/0-33480/6. • (Jun 1st 2005 - Dec 31st 2005). • – 206 days live time. • Atmospheric Neutrino MC. • – runs 233-242 (Cambridge). • – 10x 178.3 kT-Yrs. • Upward Muon MC. • – runs 900-999 (Brian). • – 100x 2 yrs . • Cosmic Muon MC. • – runs 651-750 (Cambridge). • – 100x 17 hrs (70 days) . Andy Blake, Cambridge University Combined Analysis, slide 5
Data Quality Checks Good Data Bad Data Andy Blake, Cambridge University Combined Analysis, slide 6
Data Quality Checks BAD HV Integrated Live Time BAD COIL 100% line Andy Blake, Cambridge University Combined Analysis, slide 7
Data/MC Checks Stopping Muons: Reconstructed Zenith Angle Andy Blake, Cambridge University Combined Analysis, slide 8
Data/MC Checks Stopping Muons: RMS of Down-Going Timing Fit Timing fitting code under further study Andy Blake, Cambridge University Combined Analysis, slide 9
Data/MC Checks Stopping Muons: Track Fit Relative Error Andy Blake, Cambridge University Combined Analysis, slide 10
Reconstruction Checks Selected Contained Atmospheric Muon Neutrinos: Reconstructed Muon Momentum From Range. Andy Blake, Cambridge University Combined Analysis, slide 11
Reconstruction Checks Selected Contained Atmospheric Muon Neutrinos: Reconstructed Muon Momentum From Curvature. Andy Blake, Cambridge University Combined Analysis, slide 12
Reconstruction Checks Selected Contained Atmospheric Muon Neutrinos: Reconstructed Hadronic Shower Energy. Andy Blake, Cambridge University Combined Analysis, slide 13
Contained Event Selection • Good Event Selection • – Data Quality Cuts. • Good Track Selection. • – Track Planes. • – Track Like Planes. • – Track Pulse Height Fraction. • – “Fiducial” Track. • Containment Cuts. • – Contained Vertex. • – Trace Cut. • – Topology Cuts (I). • – Topology Cuts (II). • – Veto Shield Cut. John C’s Selection Cuts Andy Blake, Cambridge University Combined Analysis, slide 14
Trace Cut Contained Event Selection: Trace cut Andy Blake, Cambridge University Combined Analysis, slide 15
Topology Cuts (I) Contained Event Selection: Mean/RMS position of strips around vertex Andy Blake, Cambridge University Combined Analysis, slide 16
Topology Cuts (I) Contained Event Selection: Maximum plane charge, Strip displacement at vertex Andy Blake, Cambridge University Combined Analysis, slide 17
Topology Cuts (II) Contained Event Selection: Vertex Charge vs Track Direction Region removed by selection cut Andy Blake, Cambridge University Combined Analysis, slide 18
Veto Shield Cut Time of Closest Veto Shield Hit (Cut at ±50 ns) Andy Blake, Cambridge University Combined Analysis, slide 19
Results Contained Event Selection expectation = 86 ± 12 events (no oscillations) Andy Blake, Cambridge University Combined Analysis, slide 20
Up-Going Event Selection • Good Event Selection. • – Data Quality Cuts. • Good Track Selection. • – Track Planes. • – Track Like Planes. • – Track Pulse Height Fraction. • – “Fiducial” Track. • Up Going Cuts. • – Topology Cuts. • – Timing Cuts (I). • – Timing Cuts (II). John C’s Selection Cuts Andy Blake, Cambridge University Combined Analysis, slide 21
Timing Cuts (I) Upward Event Selection: RMS timing deviations (up going, down going timing fits) Andy Blake, Cambridge University Combined Analysis, slide 22
Timing Cuts (I) Upward Event Selection: RMS(up) / Range cut Andy Blake, Cambridge University Combined Analysis, slide 23
Timing Cuts (I) Upward Event Selection: time slope cut Andy Blake, Cambridge University Combined Analysis, slide 24
Timing Cuts (II) Upward Event Selection: RMS(up)-RMS(down) Andy Blake, Cambridge University Combined Analysis, slide 25
Results Upward Event Selection expectation = 69 ± 10 events (no oscillations) Andy Blake, Cambridge University Combined Analysis, slide 26
Electron Neutrino Selection • Good Event Selection. • – Data Quality Cuts. • Good Shower Selection. • – Clean Planes. • – Shower Planes. • – Shower Pulse Height Fraction. • – “Fiducial” Shower. • Up Going Cuts. • – Contained Vertex. • – Shower Trace. • – Mean/RMS strips/charge per plane. • – Moment of Inertia Tensor. Ben’s Selection Cuts Andy Blake, Cambridge University Combined Analysis, slide 27
Shower Trace Cut Electron Neutrino Selection: Shower Trace. SHORT SHOWERS LONG SHOWERS Andy Blake, Cambridge University Combined Analysis, slide 28
Topology Cuts Electron Neutrino Selection: RMS Strips Per Plane. SHORT SHOWERS LONG SHOWERS Andy Blake, Cambridge University Combined Analysis, slide 29
MOI Cuts Electron Neutrino Selection: Moment of Inertia Tensor. SHORT SHOWERS LONG SHOWERS Andy Blake, Cambridge University Combined Analysis, slide 30
Veto Shield Cut Time of Closest Veto Shield Hit (Cut at ±50 ns) Andy Blake, Cambridge University Combined Analysis, slide 31
Results Electron Neutrino Selection expectation = 40 ± 5 events (no oscillations) Andy Blake, Cambridge University Combined Analysis, slide 32
Summary Muon Neutrino Selection Electron Neutrino Selection Andy Blake, Cambridge University Combined Analysis, slide 33
Selected FC/PC Events Selected FC/PC Events: Vertex Distance from Edge of Detector Andy Blake, Cambridge University Combined Analysis, slide 34
Selected FC/PC Events Selected FC/PC Events: Reconstructed Neutrino Energy Andy Blake, Cambridge University Combined Analysis, slide 35
Selected FC/PC Events Selected FC/PC Events: Reconstructed L/E Andy Blake, Cambridge University Combined Analysis, slide 36
Vetoed FC/PC Events Vetoed FC/PC Events: Vertex Distance from Edge of Detector Andy Blake, Cambridge University Combined Analysis, slide 37
Vetoed FC/PC Events Vetoed FC/PC Events: Reconstructed Zenith Angle (from top of track) Andy Blake, Cambridge University Combined Analysis, slide 38
Vetoed FC/PC Events Vetoed FC/PC Events: Track Planes. Andy Blake, Cambridge University Combined Analysis, slide 39
Selected UPMU Events Selected UPMU Events: Reconstructed Muon Momentum Andy Blake, Cambridge University Combined Analysis, slide 40
Selected UPMU Events Selected UPMU Events: Reconstructed Zenith Angle Need to add the horizontal events Andy Blake, Cambridge University Combined Analysis, slide 41
Selected NUE Events Selected NUE Events: Vertex Distance from Edge of Detctor Andy Blake, Cambridge University Combined Analysis, slide 42
Selected NUE Events Selected NUE Events: Shower Planes Andy Blake, Cambridge University Combined Analysis, slide 43
Selected NUE Events Selected NUE Events: Reconstructed Shower Energy. Reconstructed using hadronic energy scale, So needs re-tuning Andy Blake, Cambridge University Combined Analysis, slide 44
Vetoed NUE Events Vetoed NUE Events: Vertex Distance from Edge of Detector Andy Blake, Cambridge University Combined Analysis, slide 45
Vetoed NUE Events Vetoed NUE Events: Reconstructed Shower Energy Large disagreement between data and MC Andy Blake, Cambridge University Combined Analysis, slide 46
Summary • Much of Analysis Implemented. • – Data Quality Checks. • – Data/MC Comparisons. • – MC Checks. • – Event Selection. • Some disagreements. • – Timing Simulation. • – Shower Energy Reconstruction. • – Electron Neutrino Background. • Future Work. • – Timing Calibration for Latest Data. • – Study Veto Shield. • – Bayesian L/E Resolution. Andy Blake, Cambridge University Combined Analysis, slide 47