nuisance is hosted by Hepforge, IPPP Durham
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.

Changes between Version 12 and Version 13 of HowToUseNUISCOMP-GENIE


Ignore:
Timestamp:
Jun 12, 2017, 7:13:53 PM (7 years ago)
Author:
Patrick Stowell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • HowToUseNUISCOMP-GENIE

    v12 v13  
    3232The standard gevgen application options can be ran using
    3333{{{
    34  $ gevgen -f  -n -e
     34 $ gevgen -h
     35Syntax:
     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]
    3554}}}
    3655
     
    4160We want to run comparisons to MiniBooNE muon-neutrino scattering data on a mineral oil target, therefore to generate events we can run
    4261{{{
    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
    4465}}}
    4566
     
    136157}}}
    137158
    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:
    161159
    162160== Making a GENIE Comparison ==