EHE MC Check

1. EHE MC Check K. Mase* on behalf of EHE penguins, Chiba university*

2. The purpose of this study To check the MC for the Extremely High Energy (EHE) nutrinos (muons). • Test data generated at Madison (Aya, Paolo) and Chiba (Keiichi, Shigeru) • Simulated condition @ Chiba Muons are simulated with JULIeT. JULIeT → photonics → ROMEO → DOM simulator and calibrator → weighting-module 6 < log10 E < 9 [GeV] dF/dE = E-1 The direction is uniform (default) The position is uniform within 880 m radius. The interaction of pair creation, brems. and photo-nuclear, and the decay is switched on. (default) bulk ice photonics table Collaboration Meeting Oxford

3. 1015 eV 1016 eV 1017 eV 1018 eV 1019 eV 1020 eV

4. Primary Particle # # Log10(Primary energy) [GeV] cos(theta) [rad.] # This is fine! Collaboration Meeting Oxford Azimuth [rad.]

5. The highest energy cascade along a muon track # # X [m] Y [m] # # Log10(Primary energy) [GeV] Collaboration Meeting Oxford Z [m]

6. The highest energy cascade along a muon track Pair creation+brems.+photo-nuclear Only pair creation Log10(highest cascade energy) Log10(Primary energy) [GeV] Collaboration Meeting Oxford

7. The z position of the highest energy cascade downgoing upgoing uniform Strange even only pair creation case Need debug Collaboration Meeting Oxford

8. upgoing downgoing Log10(dF/dLog10E) [/cm2/s/sr] Log10(dF/dLog10E) [/cm2/s/sr] Log10(Primary energy) [GeV] Log10(Primary energy) [GeV] upgoing downgoing Log10(Simulated Flux) [/cm2/sr] Log10(Simulated Flux) [/cm2/sr] Log10(Primary energy) [GeV] Log10(Primary energy) [GeV] • Weight modulefor JULIeT(Aya) JULIeT numerical calc. → weight Simulated flux Collaboration Meeting Oxford

9. Detector response log10(OM Hits) log10(Total Np.e.) Log10(Primary energy) [GeV] Log10(Primary energy) [GeV] There are correlation. Collaboration Meeting Oxford

10. 5.36 s/event +284 289 s/event +94 383 s/event +118 501 s/event • Consuming time • Xeon 2.4GHz • Used Project • I3DefaultPhysicsSource • I3DefaultDetectorStatusSource • I3DefaultCalibrationSource • I3TextFileGeometrySource • I3MCCalibrationFiller • I3MCRawDOMStatusFiller • I3CalibrateStatusModule • I3JulietPrimaryParticleSource • I3Juliet • I3HitConstructor • I3NoiseGeneratorModule • I3RomeoInterfaceModule • I3DOMsimulator • I3DOMcalibrator • I3BoostPhysicsWriter • I3BoostGeometryWriter • I3BoostCalibrationWriter Only 1 big event ← 1 hour Acceptable level (1000 events/week) Averaged over 100 events It takes too much time for high energy (log10E > 9 GeV). Collaboration Meeting Oxford

11. Summary and future plan • The MC software to study EHE neutrinos are almost ready. • However, there are several problems. • strange distribution of the maximum cascade z position • It takes too much time to simulate high energy events (log10E > 9 GeV). • Other parameters and the waveform should be checked. • We will start mass production of MC events after fixing the above problems. • We will develop a reconstruction method for the EHE neutrino events. Collaboration Meeting Oxford

12. Backup Collaboration Meeting Oxford

13. Technical problem • When I run large number of events, I got a error. > private/dataclasses/I3DOMCalibration.cxx:84 Invalid ATWD channel in I3DOMCalibration > private/root-icetray/RootI3Tray.cxx:443 Invalid ATWD channel in I3DOMCalibration > private/root-icetray/RootI3Tray.cxx:457 Fatal exception caught: Invalid ATWD channel in I3DOMCalibration Collaboration Meeting Oxford

14. S Neutrino-Induced Cascades: Signal and Background • Signature ofne and nt are hadronic and electro-magnetic cascades. • Neutral Current interactions of all neutrino flavors produce hadronic cascades • Background consists of atmospheric muons, emitting energetic secondaries ~ 5 m Collaboration Meeting Oxford

15. 11000m 2800 m 1400 m Angular Spectrum Up Down Collaboration Meeting Oxford Bartol 2004

16. Events in Energy-direction plane By Monte Carlo GZK domain Atm m domain Collaboration Meeting Oxford

17. Energy deposit profile in km3 • By the JULIeT Monte Carlo function Collaboration Meeting Oxford

18. IceCube EHE n Sensitivity 90% C.L. for 10 year observation Collaboration Meeting Oxford

19. Particle identification IceCube • Solid: particle ID, direction, and energy • Shaded: energy only • (Supernova detection) Collaboration Meeting Oxford

20. JULIeT Photonics ROMEO DOM Sim • Software for EHE neutrinos EHE propagator Cherenkov Speed improved 　→ need more. 2D absolute sensitivity included saturation included Detector Electronics Now, the software is almost ready. Collaboration Meeting Oxford

21. Consuming time • Xeon 2.4GHz • Used Project • I3DefaultPhysicsSource • I3DefaultDetectorStatusSource • I3DefaultCalibrationSource • I3TextFileGeometrySource • I3MCCalibrationFiller • I3MCRawDOMStatusFiller • I3CalibrateStatusModule • I3JulietPrimaryParticleSource • I3Juliet • I3HitConstructor • I3NoiseGeneratorModule • I3RomeoInterfaceModule • I3DOMsimulator • I3DOMcalibrator • I3BoostPhysicsWriter • I3BoostGeometryWriter • I3BoostCalibrationWriter Only 1 big event 100 events 21.45s/4events 1154s/4events 1533s/4events 2003s/4events 3015s/4events It takes too much time for high energy (log10E > 9 GeV). Collaboration Meeting Oxford

22. amanda-core http://code.icecube.wisc.edu/svn/projects/amanda-core/releases/V01-00-02 dataclasses http://code.icecube.wisc.edu/svn/projects/dataclasses/releases/V01-08-02 dataio http://code.icecube.wisc.edu/svn/projects/dataio/releases/V02-00-01 examples-offline http://code.icecube.wisc.edu/svn/projects/examples-offline/releases/V01-04-01 icetray http://code.icecube.wisc.edu/svn/projects/icetray/releases/V01-05-01 offline-mk http://code.icecube.wisc.edu/svn/projects/offline-mk/releases/V01-06-01 phys-services http://code.icecube.wisc.edu/svn/projects/phys-services/releases/V01-07-01 root-icetray http://code.icecube.wisc.edu/svn/projects/root-icetray/releases/V01-05-00 ithon http://code.icecube.wisc.edu/svn/projects/ithon/releases/V01-00-01 eventviewer http://code.icecube.wisc.edu/svn/projects/eventviewer/releases/V01-05-00 eventviewer-modules http://code.icecube.wisc.edu/svn/projects/eventviewer-modules/releases/V01-05-00 I3Db http://code.icecube.wisc.edu/svn/projects/I3Db/releases/V01-04-01 analysis-tree http://code.icecube.wisc.edu/svn/projects/analysis-tree/releases/V02-03-00 monolith-reader http://code.icecube.wisc.edu/svn/projects/monolith-reader/releases/V02-00-02 triggerUtil-C++ http://code.icecube.wisc.edu/svn/projects/triggerUtil-C++/releases/V00-01-03 DOMcalibrator http://code.icecube.wisc.edu/svn/projects/DOMcalibrator/releases/V01-00-01/ ICETRAY-APP http://code.icecube.wisc.edu/svn/projects/ICETRAY-APP/releases/V01-05-00 OFFLINE-SOFTWARE http://code.icecube.wisc.edu/svn/projects/OFFLINE-SOFTWARE/releases/V01-06-03 DOMsimulator http://code.icecube.wisc.edu/svn/projects/DOMsimulator/releases/V01-02-01 PSInterface http://code.icecube.wisc.edu/svn/sandbox/ehe_sandbox/PSInterface c2j-icetray http://code.icecube.wisc.edu/svn/projects/c2j-icetray/releases/V02-01-01 corsika++ http://code.icecube.wisc.edu/svn/projects/corsika++/releases/V01-05-02 hit-constructor http://code.icecube.wisc.edu/svn/sandbox/ehe_sandbox/hit-constructor mmc http://code.icecube.wisc.edu/svn/projects/mmc/releases/V00-02-04 mmc-icetray http://code.icecube.wisc.edu/svn/projects/mmc-icetray/releases/V01-01-01 noise-generator http://code.icecube.wisc.edu/svn/projects/noise-generator/releases/V00-01-00 romeo http://code.icecube.wisc.edu/svn/projects/romeo/trunk romeo-amanda http://code.icecube.wisc.edu/svn/projects/romeo-amanda/releases/V00-00-04 romeo-interface http://code.icecube.wisc.edu/svn/projects/romeo-interface/trunk i3anis http://code.icecube.wisc.edu/svn/projects/i3anis/releases/V00-00-02 simple-muon http://code.icecube.wisc.edu/svn/projects/simple-muon/releases/V00-00-02/ juliet http://code.icecube.wisc.edu/svn/projects/juliet/trunk juliet-interface http://code.icecube.wisc.edu/svn/projects/juliet-interface/trunk ucr-icetray http://code.icecube.wisc.edu/svn/projects/ucr-icetray/releases/V00-01-00/ TWRSimulator http://code.icecube.wisc.edu/svn/projects/TWRSimulator/releases/V00-00-02/ cascade-generator http://code.icecube.wisc.edu/svn/projects/cascade-generator/releases/V00-00-01/ examples-simulation http://code.icecube.wisc.edu/svn/projects/examples-simulation/releases/V00-10-00/ topsim http://code.icecube.wisc.edu/svn/projects/topsim/releases/V00-00-01/ weighting-module http://code.icecube.wisc.edu/svn/sandbox/weighting-module FeatureExtractor http://code.icecube.wisc.edu/svn/projects/FeatureExtractor/trunk atmo-flux http://code.icecube.wisc.edu/svn/sandbox/atmo-flux fat-reader http://code.icecube.wisc.edu/svn/projects/fat-reader/trunk Collaboration Meeting Oxford