Systematic Errors in Shower Energy Scale Reconstruction for Neutrino Detectors
120 likes | 217 Vues
Explore systematic errors in shower energy scale reconstruction for neutrino detectors. Calculate total systematic error, evaluate PMT gain miscalibration, light level miscalibration, and more. Investigate remaining effects.
Systematic Errors in Shower Energy Scale Reconstruction for Neutrino Detectors
E N D
Presentation Transcript
Shower Energy Scale Systematics Philip Rodrigues Oxford Group Meeting 26 April 2007
Aims • Calculate total syst error on shower energy scale • Contributions: • PMT gain miscalibration • Light level miscalibration • Random miscalibrations • Dead channels, noise • Need ND and FD. Only ND so far
Method • Reconstruct reroot files with nominal and altered settings • Follow showers from each true nu in fid vol • Sum shower energy in mips
PMT gain • Only relevant to MC, which does this: • Divide by gain to do deposited energy PE • Multiplies by gain to do PE ADC • So it should cancel, but PE counting and randomness can change it
Sanity Check • Change gains by factor 5, to make sure this is doing what we think it is Seems reasonable: Gain calib should cancel in limit of high #photons/PEs (ie low gain regime), but quantization effects occur if low #photons/PE (ie high gain)
Effect of Random Seed • Alter random seed in PhotonTransport/DetSim • Reco same events
Rolling • Effect of random seed dominates effect of gain • To overcome this: • Take moderate # of fiducial events • Reroll each event 500 times • Average the shower energy over all rerollings • For error, use error on the mean (not RMS of distribution)
Result • Still no systematic shift • Low energy events have big spread – why?
Reconstruction Effects • Weird event:
Hybrid method • Reroll events as before • Now add binning – take weighted mean in 0.5 GeV bins in true shw E • Could almost believe there’s a systematic effect in gain up • Will average over all events to get one number
Light level • Should be simpler – 5% change in GeVmip calibration should give 5% change in shw energy • Lo and behold!
Still to do • Finish off gain systematic • Decide what to do with badly reco’d evts • Get one number to characterize the error • Repeat for FD • Check light level for FD • Investigate remaining effects (dead channels, noise, random miscalibrations…)