120 likes | 209 Views
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.
E N D
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…)