1 / 22

Status of Sirene

Status of Sirene. Maarten de Jong. What?. Sirene is yet another program that simulates the detector response to muons and showers It uses a general purpose collections framework for PDF tables allows for optimisation of accuracy and speed of interpolations facilitates I/O. JTools.

piera
Download Presentation

Status of Sirene

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. Status of Sirene Maarten de Jong

  2. What? • Sirene is yet another program that simulates the detector response to muons and showers • It uses a general purpose collections framework for PDF tables • allows for optimisation of accuracy and speed of interpolations • facilitates I/O

  3. JTools • A Java inspired collections framework in C++ with multi-dimensional interpolation functionality • different collection types • each collection has map functionality • variable bin size O(log(N)) access time • equidistant bins O(1) access time • allows for any number of dimensions and bins • binning in dimension n may depend on bin in dimension n-1 • allows for different interpolation algorithms • polynomial • spline • allows for a mixture of different collection types (new) • use of type lists

  4. JPhysics • Formalism for light production, propagation and detection • semi-analytical PDFs, based on single scattering approximation • Custom classes for PDF/CDF tables • treatment of weight functions to improve accuracy of interpolations (or equivalently reduce number of bins) • may serve as input to Sirene • Custom muon energy loss cross sections and shower energy generation (Paul Kooijman) • Bremmstrahlung (threshold 10 MeV) • pair production (threshold 100 MeV) • implementation photo-nuclear interactions pending

  5. Procedure • Make PDF tables • at present, JPhysics table generators are used • other applications can be interfaced readily • Convert PDFs to CDFs • allows to make various checks • Run event processor • read detector geometry • setup radiation tables (to speed up determination of cross section and shower energy for pair production) • setup tables with maximal light yield (eliminate dependence on orientation of PMT) • event loop

  6. Standard Antares Monte Carlo production “km3 + geasim”

  7. Processing steps • Read Event from MonteCarloEventWriter output file • Remove existing hits, if any • Propagate muon(s) • simulate energy loss and EM-showers • generate hits • Process shower particles from primary vertex • simulate hadronic showers as an EM shower, but assign different weights to the energies of different particle species • generate hits as one shower but use correct longitudinal profile of individual showers • Merge hits (to speed-up TriggerEfficiency) • DTmax typically 15 ps • Write Event to MonteCarloEventWriter compatible output file

  8. Trigger effective volume Volume [m3] km3 (V3.7) Sirene ± 10% Ratio En [GeV]

  9. Muon light yield scattered direct km3 (V3.7) Sirene light yield [n.p.e] En [GeV] En [GeV]

  10. Muon PDF scattered direct km3 (V3.7) Sirene light yield [n.p.e] Dt [ns] Dt [ns]

  11. EM-shower light yield scattered direct km3 (V3.7) Sirene light yield [n.p.e] En [GeV] En [GeV]

  12. EM-shower PDF scattered direct km3 (V3.7) Sirene light yield [n.p.e] Dt [ns] Dt [ns]

  13. Performance Average time per event ¶ time [ms] En [GeV] ¶ Events with ≥ 1 hits

  14. First tests KM3NeTMonte Carlo production

  15. KM3NeT • Detector • 154 towers, each with 40 optical modules, each with 31 PMTs • number of PMTs 154 x 40 x 31 = 190,960 • geometry ref154_3inch31pm180_204006.det • PMT • type 3 inch, standard Bialkali • light collector ring acceptance += 20% • Generation • particles muon neutrinos • Energy spectrum E-1.4 • Energy range [102,107] GeV • Running time • km3 (V4.3) ~135 minutes / 35,000 events • Sirene~12 minutes / 35,000 events

  16. muon light yield scattered direct km3 (V4.3) Sirene light yield [n.p.e] En [GeV] En [GeV]

  17. Muon PDF scattered direct km3 (V4.3) Sirene light yield [n.p.e] Dt [ns] Dt [ns]

  18. EM-shower light yield direct + scattered km3 (V4.3) In km3 V4.3 direct light from EM-showers is tabulated together with scattered light from EM-showers Sirene light yield [n.p.e] En [GeV]

  19. EM-shower PDF direct + scattered km3 (V4.3) In km3 V4.3 direct light from EM-showers is tabulated together with scattered light from EM-showers Sirene light yield [n.p.e] Dt [ns]

  20. Performance Average time per event ¶ time [ms] En [GeV] ¶ Events with ≥ 1 hits

  21. Summary & Outlook • Sirene is a new program that simulates detector response to muons and showers • Antares • trigger effective volume is similar • distributions of arrival times are slightly different¶ • Sirene is ready for use... • KM3NeT • distributions of arrival times in km3 V4.3 seem to be more reliable than V3.7 (thanks to Clancy James) • Sirene is ready for tests... • ¶ problem in Sirene with light from EM-showers fixed (shower energy)

  22. Summary & Outlook (II) • Possible improvements • include photo-nuclear cross sections • include [multiple] scattering of muon • Possible enhancements • generation of neutrino vertex and simulation of detector response, possibly including real-time trigger, could operate in a single step for the run-by-run Monte Carlo simulation

More Related