1 / 20

A. Margiotta on behalf of Bologna group Dipartimento di Fisica – Università di Bologna e

MUPAGE and its parameterization MU on PA rametric GE nerator A generator of atmospheric muon events STATUS. A. Margiotta on behalf of Bologna group Dipartimento di Fisica – Università di Bologna e Sezione INFN Bologna. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009.

yolandamunoz
Télécharger la présentation

A. Margiotta on behalf of Bologna group Dipartimento di Fisica – Università di Bologna e

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. MUPAGEand its parameterizationMUon PArametric GEneratorA generator of atmospheric muon eventsSTATUS A. Margiotta on behalf of Bologna group Dipartimento di Fisica – Università di Bologna e Sezione INFN Bologna

  2. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 MUPAGE is back again! • Some bugs found and corrected • In v2r1 showers hitting lateral surface leave the can volume, • instead of entering  the number of multimuon events entering • the lateral surface is suppressed • BUT • Unfortunately a new bug was introduced in v3r1  a problem • with a variable name • the corrected MUPAGE v3r2 is now available • on CVS Lyon Repository • cvs co mupage 2

  3. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 A short story of MUPAGE • Standalone c++ code based on the parameterization of the flux of single and multiple atmospheric muons (Astropart.Phys.25:1-13,2006) • Authors: Giada Carminati – Annarita Margiotta – Maurizio Spurio • INPUT: geometrical and physical information: • detector geometry, min and max muon energy, max lateral distance, min and max multiplicity……… • OUTPUT: ASCII file in table format or ANTARES format (I/O ANTARES library required) C++ code C++ code MUPAGE v3r2 parametric formulas parametric formulas

  4. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 The main steps rnd • (uniform) random generation of • bundle direction () • impact point on a can surrounding the active volume (X, Y, Z) • bundle multiplicity M • Hit-or-Miss method to reproduce the dependence of the flux on depth, zenith and multiplicity • M=1 (single muons): muon direction and impact point coincident with bundle direction and impact point. Energy according to the parametric formula • M>1 : radial and energy distribution of muons in a bundle extracted according to the parametric formulas H-M

  5. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 CAN for a km3 detector: Radius ≈ 500 m Height ≈ 1000 m Max.depth≈3.35 kmw.e. Given Ngenevents the program evaluates the equivalent livetime T • 8x106 events • LT ~ 1 day • CPU time ~ 50 minutes (in table format) (biprocessor Intel Xeon 3.06 GHz)

  6. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Zenith distribution of reconstructed tracks (Aart Strategy - ANTARES) BEFORE AFTER −5% of reconstructed events

  7. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Consistency check – 1 – flux vs depth Fexp = 6.0 x 10-5 m-2 s-1 sr-1 Fsim = (6.0 ± 0.1) x 10-5m-2 s-1 sr-1 Parametric formula

  8. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Consistency check – 2 – energy distribution MUPAGEv3r2

  9. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 ANTARES results with 5 lines (QC) data MUPAGE CORSIKA (QGSJET01) + NSU model CORSIKA (QGSJET01) + Horandel model MC uncert.

  10. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Ingredients of the parameterization • Primary CR flux • unpublished MACRO model constrained by measurements of underground muons • Primary CR interaction + shower propagation in atmosphere • HEMAS + DMPJET • Muon propagation in water • MUSIC MUPAGE

  11. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Underwater muons: data

  12. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Underwater muon: parameterizations

  13. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Flux vs. zenith angle

  14. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Comparison of parameterization with MACRO data Not only you can use MUPAGE, also its parameterization can be quoted agreement within 5%

  15. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Conclusions • MUPAGE is now OK • (dis-)agreement with experimental data (ANTARES) is within 25-40%, depending on QC • parameterization • can reproduce very well old calculations (Okada and Bugaev, f.e.) • excellent agreement with underwater/ice and underground data

  16. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009

  17. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 MUPAGE

  18. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 45 deg 20 deg impact point of the bundle axis on the can uniform random generation 70 deg

  19. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 Hit or Miss to model the flux dependence on h,  and m according to the parametric formulas; . • Event (h*, *, m*)is stored if, after generating arandom number u: 0< u < max=(hmin,0o,m=1) h*, *, m*

  20. A.Margiotta, KM3NeT meeting – Paris, 23/ 02/ 2009 data MUPAGE CORSIKA (QGSJET01) + NSU model CORSIKA (QGSJET01) + Horandel model MC uncert. ANTARES results with 5 lines (no QC)

More Related