Changes between Version 2 and Version 3 of HowToUseNUISMIN-GENIE

Timestamp:

Jun 13, 2017, 3:09:24 PM (8 years ago)

Author:

Patrick Stowell

Comment:

--

HowToUseNUISMIN-GENIE

@@ -248 +248 @@
 If you want to check this is valid you could run this card file with 'nuiscomp' as shown in the tutorial here [wiki:HowToUseNUISCOMP-GENIE] and check the nominal prediction looks okay.
 
+
+=== Reweighting GENIE Comparisons ===
+The comparisons application also allows different GENIE reweight parameters to be provided in our card file.
+
+The format for fixed dial values is:
+{{{
+genie_parameter  NAME   VALUE   STATE
+}}}
+or for free dial variations
+{{{
+genie_parameter  NAME   VALUE  LOW HIGH STEP 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)
+- **LOW** : Lower limit that the dial is allowed to float to.
+- **HIGH** : Upper limit that the dial is allowed to float to.
+- **STEP** : Step size passed to Minuit. 
+- **STATE**: State of this parameter dial, for most cases this should be left as 'FIX'.
+
+For further examples on how to include these structures in card files please see [wiki:CardFileExamples Card File Examples].
+
+
+In our example we shall generate another set of comparisons this time with two parameters shifted.
+- Charged Current Quasielastic Axial Mass : Axial mass parameter used in the quasielastic form factor
+- Charged Current Quasielastic Normalisation : Total normalisation of CCQE 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: 'MaCCQE' and 'NormCCQE' 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  
+}}}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
 === Setting up our reweight dials ===
 
+Now that we have our samples, the next thing we want to do is choose some dials to vary. If we look inside 
+
+
 ==== Defining Free dials ====