genie.fcl
Go to the documentation of this file.
1 BEGIN_PROLOG
2 
3 standard_genie:
4 {
5  module_type: "GENIEGen"
6 
7  DefinedVtxHistRange: false
8  VtxPosHistRange: [0. , 0., 0., 0., 0., 0.] #if DefinedVtxHistRange is set to true VtxPosHistRange sets the hist range of the vertex positions: [xmin, xmax, ymin, ymax, zmin, zmax].
9  #It is helpful for dual phase detector for which the range is asymmetric.
10  PassEmptySpills: false
11  FluxType: "mono" #mono, histogram, ntuple, or simple_flux
12  FluxFiles: ["flugg_L010z185i_neutrino_mode.root"] #name of file with flux histos
13  BeamName: "numi" #numi or booster at this point - really for bookkeeping
14  TopVolume: "volDetEnclosure" #volume in which to produce interactions
15  EventsPerSpill: 1. #set != 0 to get n events per spill
16  POTPerSpill: 5.e13 #should be obvious
17  MonoEnergy: 2. #in GEV
18  BeamCenter: [-1400., -350., 0.] #center of the beam in cm relative to detector coordinate origin, in meters for GENIE
19  BeamDirection: [0., 0., 1.] #all in the z direction
20  BeamRadius: 3. #in meters for GENIE
21  SurroundingMass: 0.0 #mass surrounding the detector to use
22  GlobalTimeOffset: 10000. #in ns - 10000 means the spill appears 10 us into the readout window
23  RandomTimeOffset: 10000. #length of spill in ns
24  FiducialCut: "none" #fiducial cut, see https://cdcvs.fnal.gov/redmine/projects/nusoft/wiki/GENIEHelper
25  GenFlavors: [12,14,-12,-14] #pdg codes of flux generator neutrino flavors
26  Environment: [ ] # obsolete
27  ProductionMode: "yes" #turn off the GENIE verbosity
28  EventGeneratorList: "Default"
29  DetectorLocation: "MINOS-NearDet" #location name for flux window
30  MixerConfig: "none" #no flux mixing by default
31  #MixerConfig: "swap 12:16 14:16 -12:-16 -14:-16" # example flavor swapping
32  MixerBaseline: 0. #distance from tgt to flux window needs to be set if using histogram flx
33  DebugFlags: 0 #no debug flags on by default
34  XSecTable: "gxspl-FNALsmall.xml" #default cross section
35 }
36 
37 standard_genie_atmo_flux: @local::standard_genie
38 standard_genie_atmo_flux.FluxType: "atmo_BARTOL" # use the BARTOL fluxes in GENIE
39 standard_genie_atmo_flux.AtmoEmin: 0.01 # minimum neutrino energy, in GeV
40 standard_genie_atmo_flux.AtmoEmax: 100. # maximum neutrino energy, in GeV
41 standard_genie_atmo_flux.Rl: 20. # radius of sphere on which to generate neutrinos
42 standard_genie_atmo_flux.Rt: 20. # radius of transverse area on sphere
43 standard_genie_atmo_flux.TopVolume: "volTPCActive"
44 standard_genie_atmo_flux.GenFlavors: [12, 14, -12, -14]
45 standard_genie_atmo_flux.GeomScan: "flux: 10000 1.1"
46 
47 standard_genie_uniform_flux: @local::standard_genie
48 standard_genie_uniform_flux.FluxType: "function"
49 standard_genie_uniform_flux.FunctionalFlux:"x"
50 standard_genie_uniform_flux.FunctionalBinning: 10000
51 standard_genie_uniform_flux.FluxEmin: 500
52 standard_genie_uniform_flux.FluxEmax: 600
53 standard_genie_uniform_flux.TopVolume: "volTPCActive"
54 standard_genie_uniform_flux.FluxFiles: ["dune_scan_flux.root"]
55 standard_genie_uniform_flux.GenFlavors: [12,14]
56 standard_genie_uniform_flux.GeomScan: "flux: 10000 1.1"
57 
58 standard_genie_uniform_flux_nc: @local::standard_genie_uniform_flux
59 standard_genie_uniform_flux_nc.EventGeneratorList: "NC"
60 
61 standard_genie_uniform_flux_cc: @local::standard_genie_uniform_flux
62 standard_genie_uniform_flux_cc.EventGeneratorList: "CC"
63 
64 argoneut_genie_hist_neutrino: @local::standard_genie
65 argoneut_genie_hist_neutrino.FluxType: "histogram"
66 argoneut_genie_hist_neutrino.FluxFiles: ["flugg_L010z185i_antineutrino_mode.root"]
67 argoneut_genie_hist_neutrino.BeamDirection: [0.0, -0.0583497, 0.998296]
68 
69 argoneut_genie_hist_antineutrino: @local::standard_genie
70 argoneut_genie_hist_antineutrino.FluxType: "histogram"
71 argoneut_genie_hist_antineutrino.FluxFiles: ["flugg_L010z185i_antineutrino_mode.root"]
72 argoneut_genie_hist_antineutrino.BeamDirection: [0.0, -0.0583497, 0.998296]
73 
74 argoneut_genie_simple_neutrino: @local::standard_genie
75 argoneut_genie_simple_neutrino.FluxType: "simple_flux"
76 argoneut_genie_simple_neutrino.FluxFiles: ["argoneut/gsimple_ArgoNeuT_le010z185i_run3_38l0-9r_00001.root"] #each file is 1e8 POT
77 argoneut_genie_simple_neutrino.TopVolume: "volTPCActive"
78 argoneut_genie_simple_neutrino.EventsPerSpill: 0
79 
80 
81 argoneut_genie_simple_antineutrino: @local::argoneut_genie_simple_neutrino
82 argoneut_genie_simple_antineutrino.FluxFiles: ["argoneut/gsimple_ArgoNeuT_le010z-185i_run4_38l0-9r_00001.root"] #each file is 1e8 POT
83 
84 #nucleon decay, default mode is p->nubarK+
85 #2016 PDG numbering scheme in pp.8-10 of http://www-pdg.lbl.gov/2016/listings/rpp2016-list-p.pdf (tau1 through tau60)
86 standard_ndk:
87 {
88  module_type: "NucleonDecay"
89  DecayMode: 19
90  DecayedNucleon: 2212
91 }
92 
93 #neutron-antineutron oscillation
94 standard_nosc:
95 {
96  module_type: "NeutronOsc"
97  DecayMode: 0
98 }
99 
100 END_PROLOG