genie::rew::GSyst Class Reference

#include <GSyst.h>

Static Public Member Functions
static string	AsString (GSyst_t syst)
static GSyst_t	FromString (string syst)
static bool	IsINukePionFateSystematic (GSyst_t syst)
static bool	IsINukeNuclFateSystematic (GSyst_t syst)
static bool	IsINukeFateSystematic (GSyst_t syst)
static bool	IsINukePionMeanFreePathSystematic (GSyst_t syst)
static bool	IsINukeNuclMeanFreePathSystematic (GSyst_t syst)
static bool	IsINukeMeanFreePathSystematic (GSyst_t syst)
static GSyst_t	NextPionFateSystematic (int i)
static GSyst_t	NextNuclFateSystematic (int i)
static GSyst_t	INukeFate2GSyst (INukeFateHA_t fate, int pdgc)
static GSyst_t	RBkg (InteractionType_t itype, int probe, int hitnuc, int npi)

Detailed Description

Definition at line 172 of file GSyst.h.

Member Function Documentation

static string genie::rew::GSyst::AsString ( GSyst_t  syst )

inlinestatic

Definition at line 175 of file GSyst.h.

 {
   switch(syst) {
     case ( kXSecTwkDial_MaNCEL           ) : return "MaNCEL";               break;
     case ( kXSecTwkDial_EtaNCEL          ) : return "EtaNCEL";              break;
     case ( kXSecTwkDial_NormCCQE         ) : return "NormCCQE";             break;
     case ( kXSecTwkDial_NormCCQEenu      ) : return "NormCCQEenu";          break;
     case ( kXSecTwkDial_MaCCQE           ) : return "MaCCQE";               break;
     case ( kXSecTwkDial_MaCCQEshape      ) : return "MaCCQEshape";          break;
     case ( kXSecTwkDial_ZNormCCQE        ) : return "ZNormCCQE";            break;
     case ( kXSecTwkDial_ZExpA1CCQE       ) : return "ZExpA1CCQE";           break;
     case ( kXSecTwkDial_ZExpA2CCQE       ) : return "ZExpA2CCQE";           break;
     case ( kXSecTwkDial_ZExpA3CCQE       ) : return "ZExpA3CCQE";           break;
     case ( kXSecTwkDial_ZExpA4CCQE       ) : return "ZExpA4CCQE";           break;
     case ( kXSecTwkDial_AxFFCCQEshape    ) : return "AxFFCCQEshape";        break;
     case ( kXSecTwkDial_VecFFCCQEshape   ) : return "VecFFCCQEshape";       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;
     case ( kXSecTwkDial_NormNCRES        ) : return "NormNCRES";            break;
     case ( kXSecTwkDial_MaNCRESshape     ) : return "MaNCRESshape";         break;
     case ( kXSecTwkDial_MvNCRESshape     ) : return "MvNCRESshape";         break;
     case ( kXSecTwkDial_MaNCRES          ) : return "MaNCRES";              break;
     case ( kXSecTwkDial_MvNCRES          ) : return "MvNCRES";              break;
     case ( kXSecTwkDial_MaCOHpi          ) : return "MaCOHpi";              break;
     case ( kXSecTwkDial_R0COHpi          ) : return "R0COHpi";              break;
     case ( kXSecTwkDial_RvpCC1pi         ) : return "NonRESBGvpCC1pi";      break;
     case ( kXSecTwkDial_RvpCC2pi         ) : return "NonRESBGvpCC2pi";      break;
     case ( kXSecTwkDial_RvpNC1pi         ) : return "NonRESBGvpNC1pi";      break;
     case ( kXSecTwkDial_RvpNC2pi         ) : return "NonRESBGvpNC2pi";      break;
     case ( kXSecTwkDial_RvnCC1pi         ) : return "NonRESBGvnCC1pi";      break;
     case ( kXSecTwkDial_RvnCC2pi         ) : return "NonRESBGvnCC2pi";      break;
     case ( kXSecTwkDial_RvnNC1pi         ) : return "NonRESBGvnNC1pi";      break;
     case ( kXSecTwkDial_RvnNC2pi         ) : return "NonRESBGvnNC2pi";      break;
     case ( kXSecTwkDial_RvbarpCC1pi      ) : return "NonRESBGvbarpCC1pi";   break;
     case ( kXSecTwkDial_RvbarpCC2pi      ) : return "NonRESBGvbarpCC2pi";   break;
     case ( kXSecTwkDial_RvbarpNC1pi      ) : return "NonRESBGvbarpNC1pi";   break;
     case ( kXSecTwkDial_RvbarpNC2pi      ) : return "NonRESBGvbarpNC2pi";   break;
     case ( kXSecTwkDial_RvbarnCC1pi      ) : return "NonRESBGvbarnCC1pi";   break;
     case ( kXSecTwkDial_RvbarnCC2pi      ) : return "NonRESBGvbarnCC2pi";   break;
     case ( kXSecTwkDial_RvbarnNC1pi      ) : return "NonRESBGvbarnNC1pi";   break;
     case ( kXSecTwkDial_RvbarnNC2pi      ) : return "NonRESBGvbarnNC2pi";   break;
     case ( kXSecTwkDial_AhtBY            ) : return "AhtBY";                break;
     case ( kXSecTwkDial_BhtBY            ) : return "BhtBY";                break;
     case ( kXSecTwkDial_CV1uBY           ) : return "CV1uBY";               break;
     case ( kXSecTwkDial_CV2uBY           ) : return "CV2uBY";               break;
     case ( kXSecTwkDial_AhtBYshape       ) : return "AhtBYshape";           break;
     case ( kXSecTwkDial_BhtBYshape       ) : return "BhtBYshape";           break;
     case ( kXSecTwkDial_CV1uBYshape      ) : return "CV1uBYshape";          break;
     case ( kXSecTwkDial_CV2uBYshape      ) : return "CV2uBYshape";          break;
     case ( kXSecTwkDial_NormDISCC        ) : return "NormDISCC";            break;
     case ( kXSecTwkDial_RnubarnuCC       ) : return "RnubarnuCC";           break;
     case ( kXSecTwkDial_DISNuclMod       ) : return "DISNuclMod";           break;
     case ( kHadrAGKYTwkDial_xF1pi        ) : return "AGKYxF1pi";            break;
     case ( kHadrAGKYTwkDial_pT1pi        ) : return "AGKYpT1pi";            break;
     case ( kHadrNuclTwkDial_FormZone     ) : return "FormZone";             break;
     case ( kINukeTwkDial_MFP_pi          ) : return "MFP_pi";               break;
     case ( kINukeTwkDial_MFP_N           ) : return "MFP_N";                break;
     case ( kINukeTwkDial_FrCEx_pi        ) : return "FrCEx_pi";             break;
     case ( kINukeTwkDial_FrElas_pi       ) : return "FrElas_pi";            break;
     case ( kINukeTwkDial_FrInel_pi       ) : return "FrInel_pi";            break;
     case ( kINukeTwkDial_FrAbs_pi        ) : return "FrAbs_pi";             break;
     case ( kINukeTwkDial_FrPiProd_pi     ) : return "FrPiProd_pi";          break;
     case ( kINukeTwkDial_FrCEx_N         ) : return "FrCEx_N";              break;
     case ( kINukeTwkDial_FrElas_N        ) : return "FrElas_N";             break;
     case ( kINukeTwkDial_FrInel_N        ) : return "FrInel_N";             break;
     case ( kINukeTwkDial_FrAbs_N         ) : return "FrAbs_N";              break;
     case ( kINukeTwkDial_FrPiProd_N      ) : return "FrPiProd_N";           break;
     case ( kSystNucl_CCQEPauliSupViaKF   ) : return "CCQEPauliSupViaKF";    break;
     case ( kSystNucl_CCQEMomDistroFGtoSF ) : return "CCQEMomDistroFGtoSF";  break;
     case ( kRDcyTwkDial_BR1gamma         ) : return "RDecBR1gamma";         break;
     case ( kRDcyTwkDial_BR1eta           ) : return "RDecBR1eta";           break;
     case ( kRDcyTwkDial_Theta_Delta2Npi  ) : return "Theta_Delta2Npi";      break;

     default: 
       return "-";
   }
   return "";
 }

static GSyst_t genie::rew::GSyst::FromString ( string  syst )

inlinestatic

Definition at line 257 of file GSyst.h.

 {
   const GSyst_t systematics[] =
   {
       kXSecTwkDial_MaNCEL,
       kXSecTwkDial_EtaNCEL,
       kXSecTwkDial_NormCCQE,   
       kXSecTwkDial_NormCCQEenu,   
       kXSecTwkDial_MaCCQE,        
       kXSecTwkDial_MaCCQEshape,   
       kXSecTwkDial_ZNormCCQE,   
       kXSecTwkDial_ZExpA1CCQE,        
       kXSecTwkDial_ZExpA2CCQE,        
       kXSecTwkDial_ZExpA3CCQE,        
       kXSecTwkDial_ZExpA4CCQE,        
       kXSecTwkDial_AxFFCCQEshape,
       kXSecTwkDial_VecFFCCQEshape,
       kXSecTwkDial_NormCCRES,    
       kXSecTwkDial_MaCCRESshape, 
       kXSecTwkDial_MvCCRESshape, 
       kXSecTwkDial_MaCCRES,      
       kXSecTwkDial_MvCCRES,      
       kXSecTwkDial_NormNCRES,    
       kXSecTwkDial_MaNCRESshape, 
       kXSecTwkDial_MvNCRESshape, 
       kXSecTwkDial_MaNCRES,      
       kXSecTwkDial_MvNCRES,      
       kXSecTwkDial_MaCOHpi,      
       kXSecTwkDial_R0COHpi,    
       kXSecTwkDial_RvpCC1pi,   
       kXSecTwkDial_RvpCC2pi,   
       kXSecTwkDial_RvpNC1pi,    
       kXSecTwkDial_RvpNC2pi,     
       kXSecTwkDial_RvnCC1pi,     
       kXSecTwkDial_RvnCC2pi,     
       kXSecTwkDial_RvnNC1pi,     
       kXSecTwkDial_RvnNC2pi,     
       kXSecTwkDial_RvbarpCC1pi,  
       kXSecTwkDial_RvbarpCC2pi,  
       kXSecTwkDial_RvbarpNC1pi,  
       kXSecTwkDial_RvbarpNC2pi,  
       kXSecTwkDial_RvbarnCC1pi,  
       kXSecTwkDial_RvbarnCC2pi,  
       kXSecTwkDial_RvbarnNC1pi,  
       kXSecTwkDial_RvbarnNC2pi,  
       kXSecTwkDial_AhtBY,        
       kXSecTwkDial_BhtBY,        
       kXSecTwkDial_CV1uBY,       
       kXSecTwkDial_CV2uBY,       
       kXSecTwkDial_AhtBYshape,   
       kXSecTwkDial_BhtBYshape,   
       kXSecTwkDial_CV1uBYshape,  
       kXSecTwkDial_CV2uBYshape,  
       kXSecTwkDial_NormDISCC,    
       kXSecTwkDial_RnubarnuCC,   
       kXSecTwkDial_DISNuclMod,   
       kHadrAGKYTwkDial_xF1pi,    
       kHadrAGKYTwkDial_pT1pi,    
       kHadrNuclTwkDial_FormZone, 
       kINukeTwkDial_MFP_pi,      
       kINukeTwkDial_MFP_N,       
       kINukeTwkDial_FrCEx_pi,    
       kINukeTwkDial_FrElas_pi,   
       kINukeTwkDial_FrInel_pi,   
       kINukeTwkDial_FrAbs_pi,    
       kINukeTwkDial_FrPiProd_pi, 
       kINukeTwkDial_FrCEx_N,    
       kINukeTwkDial_FrElas_N,   
       kINukeTwkDial_FrInel_N,   
       kINukeTwkDial_FrAbs_N,    
       kINukeTwkDial_FrPiProd_N, 
       kSystNucl_CCQEPauliSupViaKF,   
       kSystNucl_CCQEMomDistroFGtoSF,
       kRDcyTwkDial_BR1gamma,       
       kRDcyTwkDial_BR1eta,         
       kRDcyTwkDial_Theta_Delta2Npi,
       kNullSystematic
   };

   int isyst=0;
   while(systematics[isyst]!=kNullSystematic) {
     if( AsString(systematics[isyst]) == syst ) {
        return systematics[isyst];
     }
     isyst++;
   }

   return kNullSystematic;
 }

static GSyst_t genie::rew::GSyst::INukeFate2GSyst ( INukeFateHA_t  fate, int  pdgc  )

inlinestatic

Definition at line 469 of file GSyst.h.

 {
  // get the corresponding GSyst_t systematic parameter enumeration from the
  // input intranuke fate enumeration and PDG code
  //
  if(pdg::IsPion(pdgc)) {
     switch (fate) {
      case kIHAFtUndefined : return kNullSystematic;            break;
      case kIHAFtCEx       : return kINukeTwkDial_FrCEx_pi;     break;
      case kIHAFtElas      : return kINukeTwkDial_FrElas_pi;    break;
      case kIHAFtInelas    : return kINukeTwkDial_FrInel_pi;    break;
      case kIHAFtAbs       : return kINukeTwkDial_FrAbs_pi;     break;
      case kIHAFtPiProd    : return kINukeTwkDial_FrPiProd_pi;  break;
      default              : return kNullSystematic;            break;
     }
  } else
  if(pdg::IsNucleon(pdgc)) {
     switch (fate) {
      case kIHAFtUndefined : return kNullSystematic;           break;
      case kIHAFtCEx       : return kINukeTwkDial_FrCEx_N;     break;
      case kIHAFtElas      : return kINukeTwkDial_FrElas_N;    break;
      case kIHAFtInelas    : return kINukeTwkDial_FrInel_N;    break;
      case kIHAFtAbs       : return kINukeTwkDial_FrAbs_N;     break;
      case kIHAFtPiProd    : return kINukeTwkDial_FrPiProd_N;  break;
      default              : return kNullSystematic;           break;
     }
  }
  return kNullSystematic;
 }

static bool genie::rew::GSyst::IsINukeFateSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 381 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_FrCEx_pi    ) : 
     case ( kINukeTwkDial_FrElas_pi   ) :
     case ( kINukeTwkDial_FrInel_pi   ) :
     case ( kINukeTwkDial_FrAbs_pi    ) :
     case ( kINukeTwkDial_FrPiProd_pi ) :
     case ( kINukeTwkDial_FrCEx_N     ) : 
     case ( kINukeTwkDial_FrElas_N    ) :
     case ( kINukeTwkDial_FrInel_N    ) :
     case ( kINukeTwkDial_FrAbs_N     ) :
     case ( kINukeTwkDial_FrPiProd_N  ) :
       return true;
       break;
     
     default:
       return false;
       break;
   }
   return false;
 }

static bool genie::rew::GSyst::IsINukeMeanFreePathSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 432 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_MFP_pi ) : 
     case ( kINukeTwkDial_MFP_N  ) : 
       return true;
       break;
     
     default:
       return false;
       break;
   }
   return false;
 }

static bool genie::rew::GSyst::IsINukeNuclFateSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 364 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_FrCEx_N   ) : 
     case ( kINukeTwkDial_FrElas_N  ) : 
     case ( kINukeTwkDial_FrInel_N  ) : 
     case ( kINukeTwkDial_FrAbs_N   ) : 
     case ( kINukeTwkDial_FrPiProd_N) : 
        return true;
        break;
     default: 
        return false;
        break;
   }
   return false;
 }

static bool genie::rew::GSyst::IsINukeNuclMeanFreePathSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 418 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_MFP_N  ) : 
       return true;
       break;
     
     default:
       return false;
       break;
   }
   return false;
 }

static bool genie::rew::GSyst::IsINukePionFateSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 347 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_FrCEx_pi   ) : 
     case ( kINukeTwkDial_FrElas_pi  ) : 
     case ( kINukeTwkDial_FrInel_pi  ) : 
     case ( kINukeTwkDial_FrAbs_pi   ) : 
     case ( kINukeTwkDial_FrPiProd_pi) : 
        return true;
        break;
     default: 
        return false;
        break;
   }
   return false;
 }

static bool genie::rew::GSyst::IsINukePionMeanFreePathSystematic ( GSyst_t  syst )

inlinestatic

Definition at line 404 of file GSyst.h.

 {
   switch(syst) {
     case ( kINukeTwkDial_MFP_pi ) : 
       return true;
       break;
     
     default:
       return false;
       break;
   }
   return false;
 }

static GSyst_t genie::rew::GSyst::NextNuclFateSystematic ( int  i )

inlinestatic

Definition at line 458 of file GSyst.h.

 {
    if(i==0) return kINukeTwkDial_FrCEx_N;
    if(i==1) return kINukeTwkDial_FrElas_N;
    if(i==2) return kINukeTwkDial_FrInel_N;
    if(i==3) return kINukeTwkDial_FrAbs_N;
    if(i==4) return kINukeTwkDial_FrPiProd_N;

    return kNullSystematic;
 }

static GSyst_t genie::rew::GSyst::NextPionFateSystematic ( int  i )

inlinestatic

Definition at line 447 of file GSyst.h.

 {
    if(i==0) return kINukeTwkDial_FrCEx_pi;
    if(i==1) return kINukeTwkDial_FrElas_pi;
    if(i==2) return kINukeTwkDial_FrInel_pi;
    if(i==3) return kINukeTwkDial_FrAbs_pi;
    if(i==4) return kINukeTwkDial_FrPiProd_pi;

    return kNullSystematic;
 }

static GSyst_t genie::rew::GSyst::RBkg ( InteractionType_t  itype, int  probe, int  hitnuc, int  npi  )

inlinestatic

Definition at line 499 of file GSyst.h.

 {
   bool is_v    = pdg::IsNeutrino     (probe);
   bool is_vbar = pdg::IsAntiNeutrino (probe);
   bool is_p    = pdg::IsProton       (hitnuc);
   bool is_n    = pdg::IsNeutron      (hitnuc);
  
   // CC
   bool is_cc = (itype == kIntWeakCC);
   if(is_cc) {
     if(is_v && is_p) {
       if(npi==1) return kXSecTwkDial_RvpCC1pi;
       if(npi==2) return kXSecTwkDial_RvpCC2pi;   
     }
     if(is_v && is_n) {
       if(npi==1) return kXSecTwkDial_RvnCC1pi;
       if(npi==2) return kXSecTwkDial_RvnCC2pi;   
     }
     if(is_vbar && is_p) {
       if(npi==1) return kXSecTwkDial_RvbarpCC1pi;
       if(npi==2) return kXSecTwkDial_RvbarpCC2pi;   
     }
     if(is_vbar && is_n) {
       if(npi==1) return kXSecTwkDial_RvbarnCC1pi;
       if(npi==2) return kXSecTwkDial_RvbarnCC2pi;   
     }
   }//cc

   // NC
   bool is_nc = (itype == kIntWeakNC);
   if(is_nc) {
     if(is_v && is_p) {
       if(npi==1) return kXSecTwkDial_RvpNC1pi;
       if(npi==2) return kXSecTwkDial_RvpNC2pi;   
     }
     if(is_v && is_n) {
       if(npi==1) return kXSecTwkDial_RvnNC1pi;
       if(npi==2) return kXSecTwkDial_RvnNC2pi;   
     }
     if(is_vbar && is_p) {
       if(npi==1) return kXSecTwkDial_RvbarpNC1pi;
       if(npi==2) return kXSecTwkDial_RvbarpNC2pi;   
     }
     if(is_vbar && is_n) {
       if(npi==1) return kXSecTwkDial_RvbarnNC1pi;
       if(npi==2) return kXSecTwkDial_RvbarnNC2pi;   
     }
   }//nc

   return kNullSystematic;
 }

---

The documentation for this class was generated from the following file:

• GENIE/src/ReWeight/GSyst.h