1 / 22

EHE MC Check

EHE MC Check. K. Mase* on behalf of EHE penguins, Chiba university*. 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.

leane
Download Presentation

EHE MC Check

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  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

More Related