nugen/NuReweight/GENIEReweight.cxx
///<Set the nominal values for the reweight parameters.
<Set the nominal values for the reweight parameters.
<Set the nominal values for the reweight parameters .
void GENIEReweight::SetNominalValues() {
    fNominalParameters[(int)rwgt::fReweightMaNCEL] = 0.99;
    fNominalParameters[(int)rwgt::fReweightEtaNCEL] = 0.12;
    fNominalParameters[(int)rwgt::fReweightNormCCQE] = 1.0;
    fNominalParameters[(int)rwgt::fReweightNormCCQEenu] = 1.0;
    fNominalParameters[(int)rwgt::fReweightMaCCQEshape] = 0.99;
    fNominalParameters[(int)rwgt::fReweightMaCCQE] = 0.99;
    fNominalParameters[(int)rwgt::fReweightVecCCQEshape] = 0.84;
    fNominalParameters[(int)rwgt::fReweightNormCCRES] = 1.0;
    fNominalParameters[(int)rwgt::fReweightMaCCRESshape] = 1.12;
    fNominalParameters[(int)rwgt::fReweightMvCCRESshape]  = 0.84;
    fNominalParameters[(int)rwgt::fReweightMaCCRES] = 1.12;
    fNominalParameters[(int)rwgt::fReweightMvCCRES] = 0.84;
    fNominalParameters[(int)rwgt::fReweightNormNCRES] = 1.0;
    fNominalParameters[(int)rwgt::fReweightMaNCRESshape] = 1.12;
    fNominalParameters[(int)rwgt::fReweightMvNCRESshape] = 0.84;
    fNominalParameters[(int)rwgt::fReweightMaNCRES] = 1.12;
    fNominalParameters[(int)rwgt::fReweightMvNCRES] = 0.84;
    fNominalParameters[(int)rwgt::fReweightMaCOHpi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightR0COHpi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvpCC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvpCC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvpNC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvpNC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvnCC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvnCC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvnNC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvnNC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarpCC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarpCC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarpNC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarpNC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarnCC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarnCC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarnNC1pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRvbarnNC2pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightAhtBY] = 0.538;
    fNominalParameters[(int)rwgt::fReweightBhtBY] = 0.305;
    fNominalParameters[(int)rwgt::fReweightCV1uBY] = 0.291;
    fNominalParameters[(int)rwgt::fReweightCV2uBY] = 0.189;
    fNominalParameters[(int)rwgt::fReweightAhtBYshape] = 0.538;
    fNominalParameters[(int)rwgt::fReweightBhtBYshape] = 0.305;
    fNominalParameters[(int)rwgt::fReweightCV1uBYshape] = 0.291;
    fNominalParameters[(int)rwgt::fReweightCV2uBYshape] = 0.189;
    fNominalParameters[(int)rwgt::fReweightNormDISCC] = 1.0;
    fNominalParameters[(int)rwgt::fReweightRnubarnuCC] = 0.0;
    fNominalParameters[(int)rwgt::fReweightDISNuclMod] = 0.0;
    fNominalParameters[(int)rwgt::fReweightNC] = 1.0;
     Hadronization [free nucleon target]
    fNominalParameters[(int)rwgt::fReweightAGKY_xF1pi] = 0.385;
    fNominalParameters[(int)rwgt::fReweightAGKY_pT1pi] = 1./6.625;
     Medium-effects to hadronization
    fNominalParameters[(int)rwgt::fReweightFormZone] = 1.0;
     Intranuclear rescattering systematics.
    fNominalParameters[(int)rwgt::fReweightMFP_pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightMFP_N] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrCEx_pi] = 1.0;
#ifdef GENIE_PRE_R3
    fNominalParameters[(int)rwgt::fReweightFrElas_pi] = 1.0;
#endif
    fNominalParameters[(int)rwgt::fReweightFrInel_pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrAbs_pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrPiProd_pi] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrCEx_N] = 1.0;
#ifdef GENIE_PRE_R3
    fNominalParameters[(int)rwgt::fReweightFrElas_N] = 1.0;
#endif
    fNominalParameters[(int)rwgt::fReweightFrInel_N] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrAbs_N] = 1.0;
    fNominalParameters[(int)rwgt::fReweightFrPiProd_N] = 1.0;
    RFG Nuclear model
    fNominalParameters[(int)rwgt::fReweightCCQEPauliSupViaKF] = 1.0;
    fNominalParameters[(int)rwgt::fReweightCCQEMomDistroFGtoSF] = 0.0;
    Continous "switch" at 0.0 full FG model.  At 1.0 full spectral function model.
    From genie code it looks like weird stuff may happen for <0 and >1.
    This parameter does not have an "uncertainty" value associated with it.
    The tweaked dial value gets passed all the way through unchanged to the weight calculator
     Resonance decays
    fNominalParameters[(int)rwgt::fReweightBR1gamma] = 1.0;
    fNominalParameters[(int)rwgt::fReweightBR1eta] = 1.0;
    fNominalParameters[(int)rwgt::fReweightTheta_Delta2Npi] = 0.0;
    Continous "switch" at 0.0 full isotropic pion angular distribution.  At 1.0 full R/S pion angular distribtuion.
    This parameter does not have an "uncertainty" value associated with it.
    The tweaked dial value gets passed all the way through unchanged to the weight calculator
  }
void GENIEReweight :: SetNominalValues ( ) { fNominalParameters [ ( int ) rwgt :: fReweightMaNCEL ] = 0.99 ; fNominalParameters [ ( int ) rwgt :: fReweightEtaNCEL ] = 0.12 ; fNominalParameters [ ( int ) rwgt :: fReweightNormCCQE ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightNormCCQEenu ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightMaCCQEshape ] = 0.99 ; fNominalParameters [ ( int ) rwgt :: fReweightMaCCQE ] = 0.99 ; fNominalParameters [ ( int ) rwgt :: fReweightVecCCQEshape ] = 0.84 ; fNominalParameters [ ( int ) rwgt :: fReweightNormCCRES ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightMaCCRESshape ] = 1.12 ; fNominalParameters [ ( int ) rwgt :: fReweightMvCCRESshape ] = 0.84 ; fNominalParameters [ ( int ) rwgt :: fReweightMaCCRES ] = 1.12 ; fNominalParameters [ ( int ) rwgt :: fReweightMvCCRES ] = 0.84 ; fNominalParameters [ ( int ) rwgt :: fReweightNormNCRES ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightMaNCRESshape ] = 1.12 ; fNominalParameters [ ( int ) rwgt :: fReweightMvNCRESshape ] = 0.84 ; fNominalParameters [ ( int ) rwgt :: fReweightMaNCRES ] = 1.12 ; fNominalParameters [ ( int ) rwgt :: fReweightMvNCRES ] = 0.84 ; fNominalParameters [ ( int ) rwgt :: fReweightMaCOHpi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightR0COHpi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvpCC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvpCC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvpNC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvpNC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvnCC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvnCC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvnNC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvnNC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarpCC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarpCC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarpNC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarpNC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarnCC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarnCC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarnNC1pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRvbarnNC2pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightAhtBY ] = 0.538 ; fNominalParameters [ ( int ) rwgt :: fReweightBhtBY ] = 0.305 ; fNominalParameters [ ( int ) rwgt :: fReweightCV1uBY ] = 0.291 ; fNominalParameters [ ( int ) rwgt :: fReweightCV2uBY ] = 0.189 ; fNominalParameters [ ( int ) rwgt :: fReweightAhtBYshape ] = 0.538 ; fNominalParameters [ ( int ) rwgt :: fReweightBhtBYshape ] = 0.305 ; fNominalParameters [ ( int ) rwgt :: fReweightCV1uBYshape ] = 0.291 ; fNominalParameters [ ( int ) rwgt :: fReweightCV2uBYshape ] = 0.189 ; fNominalParameters [ ( int ) rwgt :: fReweightNormDISCC ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightRnubarnuCC ] = 0.0 ; fNominalParameters [ ( int ) rwgt :: fReweightDISNuclMod ] = 0.0 ; fNominalParameters [ ( int ) rwgt :: fReweightNC ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightAGKY_xF1pi ] = 0.385 ; fNominalParameters [ ( int ) rwgt :: fReweightAGKY_pT1pi ] = 1. / 6.625 ; fNominalParameters [ ( int ) rwgt :: fReweightFormZone ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightMFP_pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightMFP_N ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrCEx_pi ] = 1.0 ; #ifdef GENIE_PRE_R3 fNominalParameters [ ( int ) rwgt :: fReweightFrElas_pi ] = 1.0 ; #endif fNominalParameters [ ( int ) rwgt :: fReweightFrInel_pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrAbs_pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrPiProd_pi ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrCEx_N ] = 1.0 ; #ifdef GENIE_PRE_R3 fNominalParameters [ ( int ) rwgt :: fReweightFrElas_N ] = 1.0 ; #endif fNominalParameters [ ( int ) rwgt :: fReweightFrInel_N ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrAbs_N ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightFrPiProd_N ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightCCQEPauliSupViaKF ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightCCQEMomDistroFGtoSF ] = 0.0 ; fNominalParameters [ ( int ) rwgt :: fReweightBR1gamma ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightBR1eta ] = 1.0 ; fNominalParameters [ ( int ) rwgt :: fReweightTheta_Delta2Npi ] = 0.0 ; }
<Set the nominal values for the reweight parameters.
<Set the nominal values for the reweight parameters .
//CCQE Nominal Values //CCQE Nominal Values //Resonance Nominal Values //Coherent pion nominal values // Non-resonance background tweaking parameters: // DIS tweaking parameters - applied for DIS events with [Q2>Q2o, W>Wo], typically Q2okReweight =1GeV^2, WokReweight =1.7-2.0GeV //v to vbar ratio reweighting is not working in GENIE at the moment //The DIS nuclear model switch is not working in GENIE at the moment // // // Hadronization [free nucleon target] // // // Medium-effects to hadronization // // // Intranuclear rescattering systematics. // //RFG Nuclear model // //Continous "switch" at 0.0 full FG model.  At 1.0 full spectral function model. //From genie code it looks like weird stuff may happen for <0 and >1. //This parameter does not have an "uncertainty" value associated with it. //The tweaked dial value gets passed all the way through unchanged to the weight calculator // // Resonance decays // //Continous "switch" at 0.0 full isotropic pion angular distribution.  At 1.0 full R/S pion angular distribtuion. //This parameter does not have an "uncertainty" value associated with it. //The tweaked dial value gets passed all the way through unchanged to the weight calculator