190 likes | 300 Views
Status of the IceTop Shower and Primary Energy Reconstruction. Stefan Klepser for the IceTop people DESY, Zeuthen, Humboldt-Universität zu Berlin April 2007. Outline. What ' s new What ' s up next. Standard Reconstruction Chain. Merge & Calibrate Waveforms DOMcalibrator.
E N D
Status of the IceTop Showerand Primary EnergyReconstruction Stefan Klepser for the IceTop people DESY, Zeuthen, Humboldt-Universität zu Berlin April 2007
Outline • What's new • What's up next Stefan Klepser: Status of the IceTop Reconstruction
Standard Reconstruction Chain Merge & CalibrateWaveforms DOMcalibrator Pulsify Waveforms→ Charge, Time topwaveprocessor 1011001... PEDOM, tDOM, wDOM Tank-Calibration,Grouping Pulses,HG/LG Treatment topeventbuilder PhysicsAnalysis Reconstruct Shower toprec E0, X, Y, θ, ... (dE/dx)TANK, tTANK, wTANK Stefan Klepser: Status of the IceTop Reconstruction
Set ofTank Pulses( ) lg(q) shower size Rref lg(R) IceTop Reconstruction Core First Guess:Center of Gravity (COG) Direction First Guess:Plane Wave Fit Core & Shower Size Log-Llh:Lateral Fit Curved Shower Front& Core & Showersize Fit:Combined Fit Stefan Klepser: Status of the IceTop Reconstruction
I3 simulation data VEM Pulse Spectra VEM = Vertical Equivalent Muon ~ Charge = O(150 PE) ~ dE/dx = 207 MeV ~ 1 PeV events Stefan Klepser: Status of the IceTop Reconstruction
Event Rate Event Rate ≥ 5 Stations = 0.74 Hz (un-pre-scaled) Deadtime Fraction < 0.1% data Stefan Klepser: Status of the IceTop Reconstruction
Event Survival Fractions Dataset: October 2006 Events ≥ 5 Stations: 513480 100% First Guesses OK: 513126 99.993% LateralFit OK: 461261 89.8% Core contained: 370407 72.1% Core 50m contained: 232426 45.2% Zenith < 40° 193251 37.6% Stefan Klepser: Status of the IceTop Reconstruction
Plot classification datameans real data CORSIKA simulationmeans shower (no array) simulation(particles weighted with average detector responses taken from GEANT4 simulation) I3 simulationmeans full icetray simulations(tanks parametrised (no GEANT4), and using a special branch of romeo to have pulse height conservation) fakedatameans „this is just a symbolic sketch“ Stefan Klepser: Status of the IceTop Reconstruction
CORSIKA simulation plot: F. Kislat Stefan Klepser: Status of the IceTop Reconstruction
159 PeV,age: 1.47 data S100 NEW: Lateral Distribution Function → Internal Report icecube/200702001 plot: F. Kislat CORSIKA simulation DLP function: Stefan Klepser: Status of the IceTop Reconstruction
NEW: Pulse Height Fluctuation Parametrisation • Needed in the Lateral Pulse Height fit • parametrised with power laws • compared in data and simulation • will be further investigated → Internal Report icecube/200702001 plot: F. Kislat,A. Van Overloop Stefan Klepser: Status of the IceTop Reconstruction
12 PeV,age: 1.75 data NEW: Pnohit in LLH fakedata Stefan Klepser: Status of the IceTop Reconstruction
Core Resolution For 1 PeV vertical showers: σ = 13.56 m (with Pnohit, 14.86 m before) a.u. data I3 simulation 50m-inside-cut Stefan Klepser: Status of the IceTop Reconstruction
Direction Reconstruction 1-2 PeV showers I3 simulation σplane = 3.8° σcombined = 2.0° plot: A. Van Overloop Stefan Klepser: Status of the IceTop Reconstruction
Smean_lgR S100(E0, θ) → E0(S100, θ) data NEW: Energy Extraction CORSIKA simulation NEW!(to be checked) Stefan Klepser: Status of the IceTop Reconstruction
Energy Resolution For 1 PeV vertical showers: • σ(log10E) = 0.0905 (with Pnohit, 0.1013 without) • Δ(log10E) = -0.082 (with Pnohit, -0.05 without) • Reconstruction Efficiency @ 1 PeV: 87 % I3 simulation Stefan Klepser: Status of the IceTop Reconstruction
Preliminary Energy Spectrum fitted Knee position: 2.6 PeV γ = -1.74 Zenith Angles Other Experiments γ = -2.05 preliminary! (not acceptance-corrected,not simulated, nothing!) Stefan Klepser: Status of the IceTop Reconstruction
Status Summary • Direction LLH fits work well • Core/Shower Size LLH fits work well • Energy extraction qualitatively good, but needs more quantitative investigation • Everything’s ready for physics analyses! Stefan Klepser: Status of the IceTop Reconstruction
What's up next • Extended Reconstruction Efficiency Studies • Shower Curvature Studies • Combined Fit Development • Deconvolution of the Energy Spectrum • Composition Analyses • Verify, verify, verify (simulation & data) • Things are improving every week → stay tuned! Stefan Klepser: Status of the IceTop Reconstruction