GENIEGen.fcl
Go to the documentation of this file.
1 BEGIN_PROLOG
2 
3 standard_genie:
4 {
5  module_type: "GENIEGen"
6  PassEmptySpills: false
7  BeamName: "numi" # numi or booster at this point - really for bookkeeping
8 #
9 # Monochromatic flux
10 #FluxType: "mono" # mono, histogram, ntuple, or simple_flux
11 #MonoEnergy: 2. # in GEV
12 #
13 # Functional flux options
14 #FunctionalBinning: 1000 # use with "function" option: number of flux bins
15 #FunctionalFlux: 1. # use with "function" option: functional form of flux
16 #FluxEmax: 2.01 # use with "function" option: what you think it is
17 #FluxEmin: 2. # use with "function" option: what you think it is
18 #
19 # Optimized flux, Nov 2017:
20  FluxType: "dk2nu"
21  FluxFiles: ["g4lbne_v3r5p4_QGSP_BERT_OptimizedEngineeredNov2017_neutrino_?????.dk2nu.root"]
22  FluxCopyMethod: "IFDH"
23  FluxSearchPaths: "/pnfs/dune/persistent/users/ljf26/fluxfiles/g4lbne/v3r5p4/QGSP_BERT/OptimizedEngineeredNov2017/neutrino/flux/dk2nu"
24 #
25 # Select entire MPD or TPC only volume:
26  TopVolume: "volWorld" # volume in which to produce interactions: entire MPD
27 #TopVolume: "TPCGas_vol" # volume in which to produce interactions: gas TPC only
28 #
29  EventsPerSpill: 1. # set != 0 to get n events per spill
30 #POTPerSpill: 7.5e13 # 1.2MWatt@120GeV, 56% uptime,1sec cycle is 6.24e13. Peer pressure is 7.5e13
31  BeamDirection: [0., 0., 1.] # all in the z direction
32  BeamCenter: [0.1, 0, -2.0] # center of the beam relative to detector coordinate origin, in meters for GENIE
33  BeamRadius: 3. # in meters for GENIE. n.b. ignored for "mono" FluxType
34  SurroundingMass: 0.0 # mass surrounding the detector to use
35 #GlobalTimeOffset: 10000. # in ns - 10000 means the spill appears 10 us into the readout window
36 #RandomTimeOffset: 10000. # length of spill in ns
37  GlobalTimeOffset: 0. # in ns - 10000 means the spill appears 10 us into the readout window
38  RandomTimeOffset: 0. # length of spill in ns
39  FiducialCut: "none" # fiducial cut, see https://cdcvs.fnal.gov/redmine/projects/nusoft/wiki/GENIEHelper
40 #GenFlavors: [14] # Make nu_mu only
41  GenFlavors: [12,14,-12,-14] # pdg codes of flux generator neutrino flavors
42 # Environment: ["GPRODMODE", "YES"] # deprecated
43  Environment: [ ] # obsolete
44  ProductionMode: "yes" # turn off the GENIE verbosity
45  EventGeneratorList: "Default"
46  DetectorLocation: "LBNE-NearDet_02" # location name for flux window.
47  MixerConfig: "none" # no flux mixing by default
48  #MixerConfig: "swap 12:16 14:16 -12:-16 -14:-16" # example flavor swapping
49  MixerBaseline: 0. # distance from tgt to flux window needs to be set if using histogram flx
50  DebugFlags: 0 # no debug flags on by default
51 #RandomSeed: 12345 # for parallelization, this can be set manually for each job, option 0 seems to produce strange distributions of generated verticies when used with "function"
52  XSecTable: "gxspl-FNALsmall.xml" # default cross section
53 }
54 
55 standard_genie_atmo_flux: @local::standard_genie
56 standard_genie_atmo_flux.FluxType: "atmo_BARTOL" # use the BARTOL fluxes in GENIE
57 standard_genie_atmo_flux.AtmoEmin: 0.01 # minimum neutrino energy, in GeV
58 standard_genie_atmo_flux.AtmoEmax: 100. # maximum neutrino energy, in GeV
59 standard_genie_atmo_flux.Rl: 20. # radius of sphere on which to generate neutrinos
60 standard_genie_atmo_flux.Rt: 20. # radius of transverse area on sphere
61 standard_genie_atmo_flux.TopVolume: "TPCGas_vol"
62 standard_genie_atmo_flux.GenFlavors: [12, 14, -12, -14]
63 standard_genie_atmo_flux.GeomScan: "flux: 10000 1.1"
64 
65 standard_genie_uniform_flux: @local::standard_genie
66 standard_genie_uniform_flux.TopVolume: "TPCGas_vol"
67 standard_genie_uniform_flux.FluxFiles: ["dune_scan_flux.root"]
68 standard_genie_uniform_flux.GenFlavors: [12,14]
69 standard_genie_uniform_flux.GeomScan: "flux: 10000 1.1"
70 
71 standard_genie_uniform_flux_nc: @local::standard_genie_uniform_flux
72 standard_genie_uniform_flux_nc.EventGeneratorList: "NC"
73 #standard_genie_uniform_flux_nc.Environment: ["GPRODMODE",
74 # "YES",
75 # "GEVGL",
76 # "NC"]
77 
78 standard_genie_uniform_flux_cc: @local::standard_genie_uniform_flux
79 standard_genie_uniform_flux_cc.EventGeneratorList: "CC"
80 #standard_genie_uniform_flux_cc.Environment: ["GPRODMODE",
81 # "YES",
82 # "GEVGL",
83 # "CC"]
84 
85 #nucleon decay
86 standard_ndk:
87 {
88  module_type: "NucleonDecay"
89  DecayMode: 10
90 }
91 # Nucleon decay mode ID:
92 # ---------------------------------------------------------
93 # ID | Decay Mode | Current Limit
94 # | | (1E+34 yrs)
95 # ---------------------------------------------------------
96 # 0 | p --> e^{+} + \pi^{0} | 1.3
97 # 1 | p --> \mu^{+} + \pi^{0} | 1.1
98 # 2 | p --> e^{+} + \eta^{0} | 0.42
99 # 3 | p --> \mu^{+} + \eta^{0} | 0.13
100 # 4 | p --> e^{+} + \rho^{0} | 0.07
101 # 5 | p --> \mu^{+} + \rho^{0} | 0.02
102 # 6 | p --> e^{+} + \omega^{0} | 0.03
103 # 7 | p --> \mu^{+} + \omega^{0} | 0.08
104 # 8 | n --> e^{+} + \pi^{-} | 0.2
105 # 9 | n --> \mu^{+} + \pi^{-} | 0.1
106 # 10 | p --> \bar{\nu}} + K^{+} | 0.4
107 # ---------------------------------------------------------
108 
109 
110 END_PROLOG
Generated by   doxygen 1.8.11