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:
-
Nov 15, 2016, 10:47:42 AM (8 years ago)
- Author:
-
Patrick Stowell
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v5
|
v6
|
|
3 | 3 | This page aims to outline how to use GENIE to create a sample of interactions on your favourite target. We'll also outline how to feed this generated sample through the NUISANCE framework and get good results. |
4 | 4 | |
5 | | = Using GENIE Inputs = |
| 5 | = Generating GENIE Events = |
| 6 | NUISANCE requires the GHepRecord format when reading in GENIE inputs. This format is required to apply the reweighting procedures to the GENIE events. The standard output of gevgen can be used, all that is required from the user is to ensure the correct target and flux is specified for a given sample. Some examples are below. |
6 | 7 | |
| 8 | ''' Generate 1E6 MiniBooNE FHC numu events ''' |
| 9 | {{{ |
| 10 | gevgen -f miniboone_fhc_numu.root,numu_mb -r 1 -n 1000000 -t 1000060120[0.857],1000010010[0.143] |
| 11 | }}} |
7 | 12 | |
8 | | = Generating GENIE Events = |
9 | | NUISANCE requires the GHepRecord format when reading in GENIE inputs. This format is required to apply the reweighting procedure to the events. |
| 13 | = Preparing GENIE Events = |
| 14 | Unfortunately because of the way the GHepRecord is setup an extra step is required to properly normalise the cross-section predictions. The PrepareGENIE app has been added to allow the splines used in the event generation to be fully reconstructed by NUISANCE. This app will loop over all events and build up the cross-section spline for each interaction channel. These are then added together to obtain the total cross-section per target in the sample. Finally the target specification is used to add these total splines together to obtain a single spline for the MC sample. |
10 | 15 | |
11 | | Unfortunately because of the way the GHepRecord is setup an extra step is required to properly normalise the cross-section predictions. The PrepareGENIE app has been added to allow the splines used in the event generation to be fully reconstructed by NUISANCE. This app will loop over all events and build up the cross-section spline for each interaction channel. These are then added together to obtain the total cross-section per target in the sample. Finally the target specification is used to add these total splines together to obtain a single spline for the MC sample. |
| 16 | This preparation step needs the target specification as well as the flux histogram passed to gevgen used to generate the events. |
| 17 | |
12 | 18 | {{{ |
13 | 19 | $ PrepareGENIE -h |
… |
… |
|
19 | 25 | [-i inputfile1.root,inputfile2.root,inputfile3.root,...] |
20 | 26 | [-f flux_root_file.root,flux_hist_name] |
21 | | [-t target1[frac1],target2[frac2],...] |
| 27 | [-t target1,target2,...] |
22 | 28 | |
23 | 29 | Prepare Mode [Default] : Takes a single GHep file, reconstructs the original GENIE splines, |