Changes between Version 2 and Version 3 of HowToUseNUISCOMP-GENIE

Timestamp:

Jun 12, 2017, 6:04:20 PM (7 years ago)

Author:

Patrick Stowell

Comment:

--

HowToUseNUISCOMP-GENIE

@@ -240 +240 @@
 }}}
 
-This will produce a lot of logging output , but eventually, if there are no errors, it should finish with:
-{{{
+In total this will produce a lot of logging output, so only some snippets are included below for comparison.
+{{{
+[LOG Fitter]: Starting nuiscomp.exe
+
+...
+
+[LOG Minmzr]:- Loading Sample : MiniBooNE_CC1pip_XSec_1DTpi_nu
+[LOG Sample]:-- Loading Sample : MiniBooNE_CC1pip_XSec_1DTpi_nu
+[LOG Sample]:-- Creating GENIEInputHandler : MiniBooNE_CC1pip_XSec_1DTpi_nu
+                |-> Total Entries    : 2500000
+                |-> Event Integral   : 1.47491e-28 events/nucleon
+                |-> Flux Integral    : 1.67753e+10 /cm2
+                |-> Event/Flux       : 8.79217e-39 cm2/nucleon
+
+
+...
+
+[LOG Fitter]: Generating Comparison.
+[LOG Reconf]:--- Starting Reconfigure iter. 0
+[LOG Reconf]:--- Event Manager Reconfigure
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 0 events. [M, W] = [33, 1]
+[LOG Reconf]:--- MiniBooNE_CC1pip_XSec_1DTpi_nu : Processed 500000 events. [M, W] = [13, 1]
+
+...
+
 [LOG Fitter]: Saving current full FCN predictions
 [LOG Minmzr]:- Writing each of the data classes...
@@ -249 +272 @@
 [LOG Fitter]: Comparison Complete.
 [LOG Fitter]: ------------------------------------ -
-}}}
+
+}}}
+
+
 
 To check that the comparison definitely finished successfully, lets open the root file and check the file is not empty.
@@ -262 +288 @@
 If you see something similar to this then the comparisons should have ran successfully.
 
+== Analysing the output ==
+
+== Providing reweight parameters ==
+The comparisons application also allows different GENIE reweight parameters to be provided in our card file.
+
+The format for this is:
+{{{
+genie_parameter  NAME   VALUE   STATE
+}}}
+
+- **NAME** : Specifies the name of the reweight dial. These can be found in '$GENIE/src/ReWeight/GSyst.cxx'
+- **VALUE**: Value of the reweight dial in units of 1-sigma variations (1-sigma defined by GENIE)
+- **STATE**: State of this parameter dial, for most cases this should be left as 'FIX'.
+
+In our example we shall generate another set of comparisons this time with two parameters shifted.
+- Charged Current Resonant Axial Mass : Axial mass parameter used in the resonant form factor
+- Charged Current Resonant Normalisation : Total normalisation of CCRES events.
+
+First we look for the possible name in GENIE reweight:
+
+**$GENIE/src/ReWeight/GSyst.h**
+{{{
+class GSyst {
+public:
+ //......................................................................................
+ static string AsString(GSyst_t syst)
+ {
+   switch(syst) {
+     case ( kXSecTwkDial_MaNCEL           ) : return "MaNCEL";               break;
+     
+     ...
+
+     case ( kXSecTwkDial_NormCCRES        ) : return "NormCCRES";            break;
+     case ( kXSecTwkDial_MaCCRESshape     ) : return "MaCCRESshape";         break;
+     case ( kXSecTwkDial_MvCCRESshape     ) : return "MvCCRESshape";         break;
+     case ( kXSecTwkDial_MaCCRES          ) : return "MaCCRES";              break;
+     case ( kXSecTwkDial_MvCCRES          ) : return "MvCCRES";              break;
+}}}
+
+We can see the dials we are interested in this list, specified by the strings: 'MaCCRES' and 'NormCCRES' respectively.
+
+Now we want to change these dials to +1 for MaCCES and -1 sigma for NormCCRES. Therefore we add the following lines to our card file:
+
+{{{
+genie_parameter MaCCRES       +1.0  FIX
+genie_parameter NormCCRES   -1.0   FIX
+}}}
+
+so that it now looks like the following
+
+**genie_tutorial.card**
+{{{
+genie_parameter MaCCRES       +1.0  FIX
+genie_parameter NormCCRES   -1.0   FIX
+
+sample MiniBooNE_CC1pip_XSec_1DTpi_nu  GENIE: gntp.2063030.ghep.root
+sample MiniBooNE_CC1pip_XSec_1DTu_nu GENIE: gntp.2063030.ghep.root  
+}}}
+
+With this new card file we can then run our comparisons again, 
+