NuGen III. Friday, May 10, 13

Transcription

1 NuGen III 1

2 NuGen III 1

3 History of neutrino simulation NUSIM (Gary) ANIS (Marek) Gary s fast neutrino simulation for AMANDA Next generation of AMANDA neutrino simulation Cross section tables, final interaction tables are provided (used in NuGen too) NeutrinoGenerator (Aya) originally it was I3 version of ANIS, but renamed when Aya implemented neutrino re-generation in NuTau CC channel 2

4 Now NuGen is the third generation! Before release V ( NuGen 1) From V to now ( NuGen II) Used up to IC40 standard production Contains bugs in weight calculation Weight bug fixed, added weighted calculation for propagation inside the Earth Code maintainer is changed from Aya to Kotoyo Will be included in the first release of IceSim V4 Coming soon! (NuGen III) Older versions are still available with new IceSim, and will be maintained as long as I m the maintainer... 3

5 Why NuGen III? 4

6 Why NuGen III? I want to use Genie at detector! 4

7 Why NuGen III? I want to use Genie at detector! I want to use CORSIKA as neutrino source and want to simulate muon - neutrino coincidence events! 4

8 Why NuGen III? I want to use Genie at detector! I want to use CORSIKA as neutrino source and want to simulate muon - neutrino coincidence events! I need only neardetector simulation because I use my own (oscillation) calculation for in Earth propagation! 4

9 Why NuGen III? I want to use Genie at detector! I want to compare effect of In-Earth propagation with various EarthModel or neutrino cross section, while keeping neardetector simulation unchanged! I want to use CORSIKA as neutrino source and want to simulate muon - neutrino coincidence events! I need only neardetector simulation because I use my own (oscillation) calculation for in Earth propagation! 4

10 Why NuGen III? I want to use Genie at detector! I want to compare effect of In-Earth propagation with various EarthModel or neutrino cross section, while keeping neardetector simulation unchanged! Well... yes, I can implement them, but future NuGen maintainers will cry... :( (NuGen is complicated enough already!) I want to use CORSIKA as neutrino source and want to simulate muon - neutrino coincidence events! I need only neardetector simulation because I use my own (oscillation) calculation for in Earth propagation! 4

11 How structures of code is changed? v3 internal classes are in namespace nugen (some of them are renamed from v2) v1, v2 classes are in legacy directory Source modules put one primary neutrino into frame with the name of NuGPrimary 5

12 NuGen III consist of three modules 6

13 NuGen III consist of three modules Source module put a primary neutrino 6

14 NuGen III consist of three modules Source module put a primary neutrino Propagator module propagates the neutrino up to detection volume 6

15 NuGen III consist of three modules Source module put a primary neutrino Propagator module propagates the neutrino up to detection volume VTXGenerator module makes secondaries within detection volume 6

16 NuGen III consist of three modules Source module put a primary neutrino Propagator module propagates the neutrino up to detection volume VTXGenerator module makes secondaries within detection volume You may skip Propagator if you want 6

17 NuGen Sources Three modules exist now and easy to add more I3NuGDiffuseSource I3NuGPointSource I3NuGSourceSelector selects a neutrino from I3MCTree (use Nathan s Corsika-reader) Currently select a neutrino randomly from all candidate Stores NuGPrimary and generation weight as I3MCWeightDict in frame μ ν 7

18 Added new option for Injection area Current setting of Injection Radius : 1200m or tau extention for mu for straight downgoing (upgoing) tracks, more than 80% of injection area is outside of the detector Some tracks are even out of range of photonics table (radius 680m) 8

19 NUSIM style injection area Define cylindrical detector area outside of IceCube Inject neutrinos which pass through the detector area the injection area varies with injection zenith angle Injection area (r14 in Nusim) is recorded in MCWeightDict OneWeight will compensate the effect extention for mu or tau Square injection area θ 2Z sin(θ) 2Z sin(θ) + 2R cos(θ) 2R cos(θ) 9

20 Propagation inside the Earth I3NuGPropagator Take a primary neutrino NuGPrimary from frame Propagate it up to near by the detector Store one final-candidate neutrino in frame as NuGInIceNeutrino and add propagation weight to I3MCWeightDict Will support table output in near future (Similar as JULIeT table but contains only neutrinos) 10

21 Vertex Generator Module Take a particle NuGInIceNeutrino (if the neutrino does not exist take NuGPrimary instead) and calculate detection volume Then make interaction vertex inside the detection volume Three interaction modes will be available using NuGen s final interaction tables using Genie s cross sections and vertex generator Hybrid (use Genie for low energy, NuGen for high energy) Coding is not finished yet, but first integration of Genie to NuGen is done! 11

22 Genie Integration to NuGen GMCJDriver is a large program Not easy to understand, IceCube-type experiment is not supposed by developers, thus may contain problematic change in any update All we need is Good bye genie::gmcjdriver! 1. Cross Section and 2. Functions to make Final State Particles Such core codes will be relatively stable, then we don t have to worry about version dependence of Genie too much (hopefully). 12

23 NuGen rock(sio2) ice(h2o) Psurv = exp(- X * N * σ) L g/cm^2 1/L Simulation Wpos Pint = 1-Psurv(L) X (VertexPosition, g/cm^2) Genie Si O H Ptot = Pint_Si+ Pint_Si = L_Si * N_Si * σ_si Pint_O = L_O * N_O * σ_o Pint_H = L_H * N_H * σ_h Pint_O + Pint_H Pint = Ptot / GlobScale GlobScale Ptot 1/L VertexPosition(g/cm^2) 13

24 Step 1 : Genie Stand-alone test 14

25 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 14

26 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 14

27 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 2. Calculate total cross section for the injected neutrino and volume material 14

28 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 2. Calculate total cross section for the injected neutrino and volume material 3. Calculate survival probability Psurv then get interaction probability Pint = Psurv 14

29 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 2. Calculate total cross section for the injected neutrino and volume material 3. Calculate survival probability Psurv then get interaction probability Pint = Psurv 4. Select a vertex randomly, and calculate position weight 14

30 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 2. Calculate total cross section for the injected neutrino and volume material 3. Calculate survival probability Psurv then get interaction probability Pint = Psurv 4. Select a vertex randomly, and calculate position weight 5. Make a vertex and generate secondaries 14

31 Step 1 : Genie Stand-alone test Make an stand-alone genie program that follows exactly same procedure as NuGen 1. Define detection cylinder volume, and inject a neutrino from the forward cap of the cylinder 2. Calculate total cross section for the injected neutrino and volume material 3. Calculate survival probability Psurv then get interaction probability Pint = Psurv 4. Select a vertex randomly, and calculate position weight 5. Make a vertex and generate secondaries Then compare it to standard Genie! 14

32 events simulation results (no weight applied) Blue : GMCJDriver with weight opt (ice) Green : GMCJDriver default (ice) Red : Genie-NuG driver (ice) Violet : Genie-NuG driver (ice + rock) Low statistics at Low E 15

33 events simulation results (no weight applied) Secondary = secondary lepton (Muon for CC, Neutrino for NC) 16

34 events simulation results (weighted) 17

35 events simulation results 18

36 Step 2 : I3NuGGenie (interaction module) 1. Port stand-alone Genie-NuG program to nugen as an interaction module (I3NuGGenie) with correct geometry (Ice + Rock) Energy range : 1GeV ~ 200GeV Zenith range : 0deg ~ 180deg PrimarySource : I3NuGDiffuseSource 2. Compare the result with current NuGen II simulation (with InjectOnCylinder option) 19

37 Comparison : NuGen and NuGGenie 20

38 Step 3 : Make a unified module Thinking about Hybrid option, a unified VertexGenerator will be easier to manage Core functions like genie::getxsec (get cross section) and genie::generateevent (making secondary particles) could be implemented as a service in genie-icetray project But for now, we may use I3NuGGenie (with Ice + Bedrock) as a substitute of genie-icetray, especially for HE muon background study 21

39 Plan towards next release Send request to genie-icetray maintainer (Ken? Claudio?) to move following headers to public directory ConvertToGST.h ConvertPDGEncodingToI3ParticleType.h If LowEnergy group want to use or check I3NuGGenie very soon, make new release as NuGen 2.5 after cleaning source codes Add table writer for propagator Add unified vertex generator Then release NuGen 3! 22

40 Backup 23

41 Benchmark mode New event class dedicated for generating benchmark dataset may be useful to get muon / cascade effective area skips InEarth propagation, interaction weight is always 1.0 currently three options are supported 1. generate charged lepton at front edge of the detector area, with primary s energy and direction 2. generate charged lepton inside the detector area, with primary s energy and direction 3. generate charged leptons (and hadron cascades) inside the detector area with standard process (energy of secondary particles are smaller than primary) 24

42 Benchmark Mode 25

43 gevgen_simpledriver 26

44 gevgen_i3 Basically same as gevgen (with GMCJDriver) : Primary flux is always E-1 Interaction vertex is uniformly distributed in detection volume (no position weight) Should perform exactly same as genie-icetray Weighted option is available (Word Weight used in Genie is a little different from what we usually use : it s the interaction weight but doesn t mean 1 - Psurv.) 27

45 Weight difference IntWeight : 1 - exp (- n*x*sigma) IntWeight_lowE : n*x*sigma IceOnly IceRock 28

46 Vertex positions (weighted) 29