close
Warning:
Can't synchronize with repository "(default)" ("(default)" is not readable or not a Git repository.). Look in the Trac log for more information.
- Timestamp:
-
Jun 12, 2017, 7:13:53 PM (7 years ago)
- Author:
-
Patrick Stowell
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v12
|
v13
|
|
32 | 32 | The standard gevgen application options can be ran using |
33 | 33 | {{{ |
34 | | $ gevgen -f -n -e |
| 34 | $ gevgen -h |
| 35 | Syntax: |
| 36 | |
| 37 | gevgen [-h] |
| 38 | [-r run#] |
| 39 | -n nev |
| 40 | -e energy (or energy range) |
| 41 | -p neutrino_pdg |
| 42 | -t target_pdg |
| 43 | [-f flux_description] |
| 44 | [-o outfile_name] |
| 45 | [-w] |
| 46 | [--seed random_number_seed] |
| 47 | [--cross-sections xml_file] |
| 48 | [--event-generator-list list_name] |
| 49 | [--message-thresholds xml_file] |
| 50 | [--unphysical-event-mask mask] |
| 51 | [--event-record-print-level level] |
| 52 | [--mc-job-status-refresh-rate rate] |
| 53 | [--cache-file root_file] |
35 | 54 | }}} |
36 | 55 | |
… |
… |
|
41 | 60 | We want to run comparisons to MiniBooNE muon-neutrino scattering data on a mineral oil target, therefore to generate events we can run |
42 | 61 | {{{ |
43 | | $ gevgen -f MiniBooNE -t CH |
| 62 | $ source $GENIE_DIR/environment_setup.sh |
| 63 | $ export GXMLPATH=${GENIE_DIR}/genie_xsec/v2_12_0/NULL/DefaultPlusMECWithNC/data/ |
| 64 | $ gevgen -n 2500000 -t 1000060120[0.85714],1000010010[0.14286] -p 14 --cross-sections $GXMLPATH/gxspl-FNALsmall.xml --event-generator-list Default -f MiniBooNE_numu_flux.root,numu_mb -e 0,10 -o gntp.2063030.ghep.root |
44 | 65 | }}} |
45 | 66 | |
… |
… |
|
136 | 157 | }}} |
137 | 158 | |
138 | | === Running gevgen_nuisance (optional alternative) === |
139 | | |
140 | | Other than the standard gevgen application, if built against GENIE events later than GENIE 2.12.0, NUISANCE can also provide a gevgen_nuisance app, with additional features to support the creation of ‘NUISANCE-ready’ events. |
141 | | |
142 | | To build gevgen_nuisance, you will need to build with the option ‘-DBUILD_GEVGEN=ON’ during the cmake configure statement as shown in the build notes here: |
143 | | |
144 | | {{{ |
145 | | $ cmake ../ -DUSE_GENIE=1 -DBUILD_GEVGEN=1 |
146 | | }}} |
147 | | |
148 | | The gevgen_nuisance can be ran in a similar way to the normal gevgen application but flux inputs and target definitions can be provided in simpler formats |
149 | | |
150 | | In our example we want to generate events for a MiniBooNE muon neutrino flux on a CH2 target. Therefore we run |
151 | | |
152 | | {{{ |
153 | | $ gevgen_nuisance -n 250000 -f MiniBooNE_fhc_numu -t CH2 |
154 | | }}} |
155 | | |
156 | | We should then open the output event file to check it constrains the GENIE event tree ‘gtree’ object and the flux and event rate predictions. |
157 | | |
158 | | Notice that we no longer need to run PrepareGENIE on gevgen_nuisance events before using them nuisance. |
159 | | |
160 | | gevgen_nuisance also provides extra support for mixed target and flux specifications. For further notes on these see: |
161 | 159 | |
162 | 160 | == Making a GENIE Comparison == |