genie::rew::GReWeightNuXSecCCQE Class Reference

Reweighting CCQE GENIE neutrino cross sections. More...

#include <GReWeightNuXSecCCQE.h>

Inheritance diagram for genie::rew::GReWeightNuXSecCCQE:

Public Member Functions
	GReWeightNuXSecCCQE ()
	~GReWeightNuXSecCCQE ()
bool	IsHandled (GSyst_t syst)
	does the current weight calculator handle the input nuisance param? More...
void	SetSystematic (GSyst_t syst, double val)
	update the value for the specified nuisance param More...
void	Reset (void)
	set all nuisance parameters to default values More...
void	Reconfigure (void)
	propagate updated nuisance parameter values to actual MC, etc More...
double	CalcWeight (const EventRecord &event)
	calculate a weight for the input event using the current nuisance param values More...
void	SetMode (int mode)
void	RewNue (bool tf)
void	RewNuebar (bool tf)
void	RewNumu (bool tf)
void	RewNumubar (bool tf)
void	SetMaPath (string p)
void	SetZExpPath (string p)
Public Member Functions inherited from genie::rew::GReWeightI
virtual	~GReWeightI ()

Static Public Attributes
static const int	kModeMa = 0
static const int	kModeNormAndMaShape = 1
static const int	kModeZExp = 2
static const int	fZExpMaxSyst = 4
	maximum number of systematics More...

Private Member Functions
void	Init (void)
double	CalcWeightNorm (const EventRecord &event)
double	CalcWeightMaShape (const EventRecord &event)
double	CalcWeightMa (const EventRecord &event)
double	CalcWeightZExp (const EventRecord &event)

Private Attributes
XSecAlgorithmI *	fXSecModelDef
	default model More...
XSecAlgorithmI *	fXSecModel
	tweaked model More...
Registry *	fXSecModelConfig
	config in tweaked model More...
string	fFFModel
int	fMode
	0: Ma, 1: Norm and MaShape, 2: Z-Expansion More...
bool	fRewNue
	reweight nu_e CC? More...
bool	fRewNuebar
	reweight nu_e_bar CC? More...
bool	fRewNumu
	reweight nu_mu CC? More...
bool	fRewNumubar
	reweight nu_mu_bar CC? More...
string	fMaPath
	M_{A} path in config Registry. More...
double	fNormTwkDial
double	fNormDef
double	fNormCurr
double	fMaTwkDial
double	fMaDef
double	fMaCurr
int	fZExpCurrIdx
	current coefficient index More...
int	fZExpMaxCoef
	max number of coefficients to use More...
string	fZExpPath
	algorithm path to get coefficients More...
double	fZExpTwkDial [fZExpMaxSyst]
double	fZExpDef [fZExpMaxSyst]
double	fZExpCurr [fZExpMaxSyst]
	array of current parameter values More...

Additional Inherited Members
Protected Member Functions inherited from genie::rew::GReWeightI
	GReWeightI ()

Detailed Description

Reweighting CCQE GENIE neutrino cross sections.

Author

Costas Andreopoulos <costas.andreopoulos stfc.ac.uk> University of Liverpool & STFC Rutherford Appleton Lab

Jim Dobson <J.Dobson07 imperial.ac.uk> Imperial College London

Aug 1, 2009

Copyright (c) 2003-2017, GENIE Neutrino MC Generator Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Definition at line 45 of file GReWeightNuXSecCCQE.h.

Constructor & Destructor Documentation

GReWeightNuXSecCCQE::GReWeightNuXSecCCQE

(

)

Definition at line 73 of file GReWeightNuXSecCCQE.cxx.

{
  this->Init();
}

GReWeightNuXSecCCQE::~GReWeightNuXSecCCQE

(

)

Definition at line 78 of file GReWeightNuXSecCCQE.cxx.

{
#ifdef _G_REWEIGHT_CCQE_DEBUG_
  fTestFile->cd();
  fTestNtp ->Write();
  fTestFile->Close();
  delete fTestFile;
#endif
}

Member Function Documentation

double GReWeightNuXSecCCQE::CalcWeight ( const EventRecord &  event )

virtual

calculate a weight for the input event using the current nuisance param values

Implements genie::rew::GReWeightI.

Definition at line 270 of file GReWeightNuXSecCCQE.cxx.

{
  bool is_qe = event.Summary()->ProcInfo().IsQuasiElastic();
  bool is_cc = event.Summary()->ProcInfo().IsWeakCC();
  if(!is_qe || !is_cc) return 1.;

  bool charm = event.Summary()->ExclTag().IsCharmEvent(); // skip CCQE charm
  if(charm) return 1.;

  int nupdg = event.Probe()->Pdg();
  if(nupdg==kPdgNuMu     && !fRewNumu   ) return 1.;
  if(nupdg==kPdgAntiNuMu && !fRewNumubar) return 1.;
  if(nupdg==kPdgNuE      && !fRewNue    ) return 1.;
  if(nupdg==kPdgAntiNuE  && !fRewNuebar ) return 1.;

  if(fMode==kModeMa && strcmp(fFFModel.c_str(),kModelDipole) == 0) {
     double wght = this->CalcWeightMa(event);
     return wght;
  } 
  else 
  if(fMode==kModeNormAndMaShape && strcmp(fFFModel.c_str(),kModelDipole) == 0) {
     double wght = 
         this->CalcWeightNorm    (event) *
         this->CalcWeightMaShape (event);
     return wght;
  }
  else
  if(fMode==kModeZExp && strcmp(fFFModel.c_str(),kModelZExp) == 0) {
     double wght = 
         this->CalcWeightNorm(event) *
         this->CalcWeightZExp(event);
     return wght;
  }

  return 1.;
}

double GReWeightNuXSecCCQE::CalcWeightMa ( const EventRecord &  event )

private

Definition at line 382 of file GReWeightNuXSecCCQE.cxx.

{
  bool tweaked = (TMath::Abs(fMaTwkDial) > controls::kASmallNum);
  if(!tweaked) return 1.0;

  Interaction * interaction = event.Summary();

  interaction->KinePtr()->UseSelectedKinematics();
  interaction->SetBit(kIAssumeFreeNucleon);

  double old_xsec   = event.DiffXSec();
  double old_weight = event.Weight();
  double new_xsec   = fXSecModel->XSec(interaction, kPSQ2fE);
  double new_weight = old_weight * (new_xsec/old_xsec);

//LOG("ReW", pDEBUG) << "differential cross section (old) = " << old_xsec;
//LOG("ReW", pDEBUG) << "differential cross section (new) = " << new_xsec;
//LOG("ReW", pDEBUG) << "event generation weight = " << old_weight;
//LOG("ReW", pDEBUG) << "new weight = " << new_weight;

  interaction->KinePtr()->ClearRunningValues();
  interaction->ResetBit(kIAssumeFreeNucleon);

  return new_weight;
}

double GReWeightNuXSecCCQE::CalcWeightMaShape ( const EventRecord &  event )

private

Definition at line 408 of file GReWeightNuXSecCCQE.cxx.

{
  bool tweaked = (TMath::Abs(fMaTwkDial) > controls::kASmallNum);
  if(!tweaked) return 1.0;

  Interaction * interaction = event.Summary();

  interaction->KinePtr()->UseSelectedKinematics();
  interaction->SetBit(kIAssumeFreeNucleon);

  double old_xsec   = event.DiffXSec();
  double old_weight = event.Weight();
  double new_xsec   = fXSecModel->XSec(interaction, kPSQ2fE);
  double new_weight = old_weight * (new_xsec/old_xsec);

//LOG("ReW", pDEBUG) << "differential cross section (old) = " << old_xsec;
//LOG("ReW", pDEBUG) << "differential cross section (new) = " << new_xsec;
//LOG("ReW", pDEBUG) << "event generation weight = " << old_weight;
//LOG("ReW", pDEBUG) << "new weight = " << new_weight;

//double old_integrated_xsec = event.XSec();
  double old_integrated_xsec = fXSecModelDef -> Integral(interaction);
  double new_integrated_xsec = fXSecModel    -> Integral(interaction);   
  assert(new_integrated_xsec > 0);
  new_weight *= (old_integrated_xsec/new_integrated_xsec);

//LOG("ReW", pDEBUG) << "integrated cross section (old) = " << old_integrated_xsec;
//LOG("ReW", pDEBUG) << "integrated cross section (new) = " << new_integrated_xsec;
//LOG("ReW", pDEBUG) << "new weight (normalized to const integral) = " << new_weight;

  interaction->KinePtr()->ClearRunningValues();
  interaction->ResetBit(kIAssumeFreeNucleon);

#ifdef _G_REWEIGHT_CCQE_DEBUG_
  double E  = interaction->InitState().ProbeE(kRfHitNucRest);
  double Q2 = interaction->Kine().Q2(true);
  fTestNtp->Fill(E,Q2,new_weight);
#endif

  return new_weight;
}

double GReWeightNuXSecCCQE::CalcWeightNorm ( const EventRecord &  event )

private

Definition at line 373 of file GReWeightNuXSecCCQE.cxx.

{
  bool tweaked = (TMath::Abs(fNormTwkDial) > controls::kASmallNum);
  if(!tweaked) return 1.0;

  double wght = fNormCurr;
  return wght;
}

double GReWeightNuXSecCCQE::CalcWeightZExp ( const EventRecord &  event )

private

Definition at line 450 of file GReWeightNuXSecCCQE.cxx.

{
  // very similar to CalcWeightMa
  bool tweaked = false;
  for (int i=0;i<fZExpMaxCoef;i++)
  {
    tweaked = tweaked || (TMath::Abs(fZExpTwkDial[i]) > controls::kASmallNum);
  }
  if(!tweaked) { return 1.0; }

  Interaction * interaction = event.Summary();

  interaction->KinePtr()->UseSelectedKinematics();
  interaction->SetBit(kIAssumeFreeNucleon);

  double old_xsec   = event.DiffXSec();
  double old_weight = event.Weight();
  double new_xsec   = fXSecModel->XSec(interaction, kPSQ2fE);
  double new_weight = old_weight * (new_xsec/old_xsec);

  interaction->KinePtr()->ClearRunningValues();
  interaction->ResetBit(kIAssumeFreeNucleon);

  return new_weight;
}

void GReWeightNuXSecCCQE::Init ( void  )

private

Definition at line 307 of file GReWeightNuXSecCCQE.cxx.

{
  AlgId id("genie::LwlynSmithQELCCPXSec","Default");

  AlgFactory * algf = AlgFactory::Instance();

  Algorithm * algdef = algf->AdoptAlgorithm(id);
  fXSecModelDef = dynamic_cast<XSecAlgorithmI*>(algdef);
  fXSecModelDef->AdoptSubstructure();

  Algorithm * alg = algf->AdoptAlgorithm(id);
  fXSecModel = dynamic_cast<XSecAlgorithmI*>(alg);
  fXSecModel->AdoptSubstructure();

  fXSecModelConfig = new Registry(fXSecModel->GetConfig());
  fFFModel = fXSecModelConfig->GetAlg("FormFactorsAlg/AxialFormFactorModel").name;

  this->RewNue    (true);
  this->RewNuebar (true);
  this->RewNumu   (true);
  this->RewNumubar(true);

  //this->SetMaPath("FormFactorsAlg/Ma");
  this->SetMaPath("FormFactorsAlg/AxialFormFactorModel/QEL-Ma");
  this->SetZExpPath("FormFactorsAlg/AxialFormFactorModel/");

  if (strcmp(fFFModel.c_str(),kModelDipole) == 0)
  {
    this->SetMode(kModeNormAndMaShape);
    fMaDef       = fXSecModelConfig->GetDouble(fMaPath);
    fZExpMaxCoef = 0;
  } else
  if (strcmp(fFFModel.c_str(),kModelZExp) == 0)
  {
    this->SetMode(kModeZExp);
    fMaDef = 0.;
    fZExpMaxCoef =
      TMath::Min(fXSecModelConfig->GetInt(fZExpPath + "QEL-Kmax"),
       this->fZExpMaxSyst);
  }

  fNormTwkDial = 0.;
  fNormDef     = 1.;
  fNormCurr    = fNormDef;
  fMaTwkDial   = 0.; 
  fMaCurr      = fMaDef;

  ostringstream alg_key;
  for (int i=0;i<fZExpMaxSyst;i++)
  {
    alg_key.str("");
    alg_key << fZExpPath << "QEL-Z_A" << i+1;
    if (strcmp(fFFModel.c_str(),kModelZExp) == 0 && i < fZExpMaxCoef)
    { fZExpDef[i] = fXSecModelConfig->GetDouble(alg_key.str()); }
    else
    { fZExpDef[i] = 0.; }
    fZExpTwkDial[i] = 0.;
    fZExpCurr   [i] = fZExpDef[i];
  }

#ifdef _G_REWEIGHT_CCQE_DEBUG_
  fTestFile = new TFile("./ccqe_reweight_test.root","recreate");
  fTestNtp  = new TNtupleD("testntp","","E:Q2:wght");
#endif
}

bool GReWeightNuXSecCCQE::IsHandled ( GSyst_t  syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 88 of file GReWeightNuXSecCCQE.cxx.

{
// read form factor model and compare to mode
   bool handle;

   switch(syst) {

     case ( kXSecTwkDial_NormCCQE    ) :
       if(fMode==kModeNormAndMaShape && strcmp(fFFModel.c_str(),kModelDipole) == 0)
       {  
          handle = true;
       } else {  
          handle = false;
       }
       break;
     
     case ( kXSecTwkDial_MaCCQEshape ) :
       if(fMode==kModeNormAndMaShape && strcmp(fFFModel.c_str(),kModelDipole) == 0)
       {  
          handle = true;
       } else {  
          handle = false;
       }
       break;
     
     case ( kXSecTwkDial_MaCCQE ) :
       if(fMode==kModeMa && strcmp(fFFModel.c_str(),kModelDipole) == 0)
       {  
          handle = true;
       } else {  
          handle = false;
       }
       break;
     
     case ( kXSecTwkDial_ZNormCCQE    ) :
       if(fMode==kModeZExp && strcmp(fFFModel.c_str(),kModelZExp) == 0)
       {  
          handle = true;
       } else {  
          handle = false;
       }
       break;
     
     case ( kXSecTwkDial_ZExpA1CCQE ):
     case ( kXSecTwkDial_ZExpA2CCQE ):
     case ( kXSecTwkDial_ZExpA3CCQE ):
     case ( kXSecTwkDial_ZExpA4CCQE ):
       if(fMode==kModeZExp && strcmp(fFFModel.c_str(),kModelZExp) == 0)
       {  
          handle = true;
       } else {
          handle = false;
       }
       break;
     
     default:
          handle = false;
          break;
   }

   return handle;
}

void GReWeightNuXSecCCQE::Reconfigure ( void  )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 202 of file GReWeightNuXSecCCQE.cxx.

{
  GSystUncertainty * fracerr = GSystUncertainty::Instance();

  if(fMode==kModeMa && strcmp(fFFModel.c_str(),kModelDipole) == 0) {   
     int    sign_matwk = utils::rew::Sign(fMaTwkDial);
     double fracerr_ma = fracerr->OneSigmaErr(kXSecTwkDial_MaCCQE, sign_matwk);
     fMaCurr = fMaDef * (1. + fMaTwkDial * fracerr_ma);
     fMaCurr = TMath::Max(0., fMaCurr  );
  }
  else
  if(fMode==kModeNormAndMaShape && strcmp(fFFModel.c_str(),kModelDipole) == 0) {
     int    sign_normtwk = utils::rew::Sign(fNormTwkDial);
     int    sign_mashtwk = utils::rew::Sign(fMaTwkDial  );
     double fracerr_norm = fracerr->OneSigmaErr(kXSecTwkDial_NormCCQE,  sign_normtwk);
     double fracerr_mash = fracerr->OneSigmaErr(kXSecTwkDial_MaCCQEshape, sign_mashtwk);
     fNormCurr = fNormDef * (1. + fNormTwkDial * fracerr_norm);
     fMaCurr   = fMaDef   * (1. + fMaTwkDial   * fracerr_mash);
     fNormCurr = TMath::Max(0., fNormCurr);
     fMaCurr   = TMath::Max(0., fMaCurr  );
  }
  else
  if(fMode==kModeZExp && strcmp(fFFModel.c_str(),kModelZExp) == 0) {
     int     sign_twk = 0;
     int     sign_normtwk = utils::rew::Sign(fNormTwkDial);
     double  fracerr_zexp = 0.;
     double  fracerr_norm = fracerr->OneSigmaErr(kXSecTwkDial_ZNormCCQE, sign_normtwk);
     GSyst_t syst;
     // loop over all indices and update each
     for (int i=0;i<fZExpMaxCoef;i++)
     {
       switch(i){
         case 0: syst = kXSecTwkDial_ZExpA1CCQE; break;
         case 1: syst = kXSecTwkDial_ZExpA2CCQE; break;
         case 2: syst = kXSecTwkDial_ZExpA3CCQE; break;
         case 3: syst = kXSecTwkDial_ZExpA4CCQE; break;
         default: return; break;
       }
       sign_twk = utils::rew::Sign(fZExpTwkDial[i]);
       fracerr_zexp = fracerr->OneSigmaErr(syst, sign_twk);
       fZExpCurr[i] = fZExpDef[i] * (1. + fZExpTwkDial[i] * fracerr_zexp);
     }
     fNormCurr = fNormDef * (1. + fNormTwkDial * fracerr_norm);
     fNormCurr = TMath::Max(0., fNormCurr);
  }
  else {
    return;
  }

  Registry & r = *fXSecModelConfig;
  if (fMode==kModeMa || fMode==kModeNormAndMaShape)
  {
    r.Set(fMaPath, fMaCurr); 
  }
  else
  if (fMode==kModeZExp)
  {
    ostringstream alg_key;
    for (int i=0;i<fZExpMaxCoef;i++)
    {
      alg_key.str(""); // algorithm key for each coefficient
      alg_key << fZExpPath << "QEL-Z_A" << i+1;
      r.Set(alg_key.str(), fZExpCurr[i]);
    }
  }
  fXSecModel->Configure(r);
}

void GReWeightNuXSecCCQE::Reset ( void  )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 186 of file GReWeightNuXSecCCQE.cxx.

{
  fNormTwkDial = 0.;
  fNormCurr    = fNormDef;
  fMaTwkDial   = 0.; 
  fMaCurr      = fMaDef;

  for (int i=0;i<fZExpMaxSyst;i++)
  {
    fZExpTwkDial[i] = 0.;
    fZExpCurr   [i] = fZExpDef[i];
  } 

  this->Reconfigure();
}

void genie::rew::GReWeightNuXSecCCQE::RewNue ( bool  tf )

inline

Definition at line 66 of file GReWeightNuXSecCCQE.h.

{ fRewNue     = tf;   }

void genie::rew::GReWeightNuXSecCCQE::RewNuebar ( bool  tf )

inline

Definition at line 67 of file GReWeightNuXSecCCQE.h.

{ fRewNuebar  = tf;   }

void genie::rew::GReWeightNuXSecCCQE::RewNumu ( bool  tf )

inline

Definition at line 68 of file GReWeightNuXSecCCQE.h.

{ fRewNumu    = tf;   }

void genie::rew::GReWeightNuXSecCCQE::RewNumubar ( bool  tf )

inline

Definition at line 69 of file GReWeightNuXSecCCQE.h.

{ fRewNumubar = tf;   }

void genie::rew::GReWeightNuXSecCCQE::SetMaPath ( string  p )

inline

Definition at line 70 of file GReWeightNuXSecCCQE.h.

{ fMaPath     = p;    }

void genie::rew::GReWeightNuXSecCCQE::SetMode ( int  mode )

inline

Definition at line 65 of file GReWeightNuXSecCCQE.h.

{ fMode       = mode; }

void GReWeightNuXSecCCQE::SetSystematic ( GSyst_t  syst, double  val  )

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 151 of file GReWeightNuXSecCCQE.cxx.

{
  if(!this->IsHandled(syst))
  {
    LOG("ReW",pWARN) << "Systematic " << GSyst::AsString(syst) << " is not handled for algorithm "
      << fFFModel << " and mode " << fMode;
    return;
  }

    switch(syst) {
      case ( kXSecTwkDial_NormCCQE ) :
      case ( kXSecTwkDial_ZNormCCQE ) :
        fNormTwkDial = twk_dial;
        break;
      case ( kXSecTwkDial_MaCCQEshape ) :
      case ( kXSecTwkDial_MaCCQE ) :
        fMaTwkDial = twk_dial;
        break;
      case ( kXSecTwkDial_ZExpA1CCQE ) :
        if(fZExpMaxCoef>0){ fZExpTwkDial[0] = twk_dial; }
        break;
      case ( kXSecTwkDial_ZExpA2CCQE ) :
        if(fZExpMaxCoef>1){ fZExpTwkDial[1] = twk_dial; }
        break;
      case ( kXSecTwkDial_ZExpA3CCQE ) :
        if(fZExpMaxCoef>2){ fZExpTwkDial[2] = twk_dial; }
        break;
      case ( kXSecTwkDial_ZExpA4CCQE ) :
        if(fZExpMaxCoef>3){ fZExpTwkDial[3] = twk_dial; }
        break;
      default:
        break;
    }
}

void genie::rew::GReWeightNuXSecCCQE::SetZExpPath ( string  p )

inline

Definition at line 72 of file GReWeightNuXSecCCQE.h.

{ fZExpPath    = p;   }

Member Data Documentation

string genie::rew::GReWeightNuXSecCCQE::fFFModel

private

Definition at line 85 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fMaCurr

private

Definition at line 98 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fMaDef

private

Definition at line 97 of file GReWeightNuXSecCCQE.h.

string genie::rew::GReWeightNuXSecCCQE::fMaPath

private

M_{A} path in config Registry.

Definition at line 92 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fMaTwkDial

private

Definition at line 96 of file GReWeightNuXSecCCQE.h.

int genie::rew::GReWeightNuXSecCCQE::fMode

private

0: Ma, 1: Norm and MaShape, 2: Z-Expansion

Definition at line 87 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fNormCurr

private

Definition at line 95 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fNormDef

private

Definition at line 94 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fNormTwkDial

private

Definition at line 93 of file GReWeightNuXSecCCQE.h.

bool genie::rew::GReWeightNuXSecCCQE::fRewNue

private

reweight nu_e CC?

Definition at line 88 of file GReWeightNuXSecCCQE.h.

bool genie::rew::GReWeightNuXSecCCQE::fRewNuebar

private

reweight nu_e_bar CC?

Definition at line 89 of file GReWeightNuXSecCCQE.h.

bool genie::rew::GReWeightNuXSecCCQE::fRewNumu

private

reweight nu_mu CC?

Definition at line 90 of file GReWeightNuXSecCCQE.h.

bool genie::rew::GReWeightNuXSecCCQE::fRewNumubar

private

reweight nu_mu_bar CC?

Definition at line 91 of file GReWeightNuXSecCCQE.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQE::fXSecModel

private

tweaked model

Definition at line 83 of file GReWeightNuXSecCCQE.h.

Registry* genie::rew::GReWeightNuXSecCCQE::fXSecModelConfig

private

config in tweaked model

Definition at line 84 of file GReWeightNuXSecCCQE.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQE::fXSecModelDef

private

default model

Definition at line 82 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fZExpCurr[fZExpMaxSyst]

private

array of current parameter values

Definition at line 105 of file GReWeightNuXSecCCQE.h.

int genie::rew::GReWeightNuXSecCCQE::fZExpCurrIdx

private

current coefficient index

Definition at line 100 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fZExpDef[fZExpMaxSyst]

private

Definition at line 104 of file GReWeightNuXSecCCQE.h.

int genie::rew::GReWeightNuXSecCCQE::fZExpMaxCoef

private

max number of coefficients to use

Definition at line 101 of file GReWeightNuXSecCCQE.h.

const int GReWeightNuXSecCCQE::fZExpMaxSyst = 4

static

maximum number of systematics

Definition at line 52 of file GReWeightNuXSecCCQE.h.

string genie::rew::GReWeightNuXSecCCQE::fZExpPath

private

algorithm path to get coefficients

Definition at line 102 of file GReWeightNuXSecCCQE.h.

double genie::rew::GReWeightNuXSecCCQE::fZExpTwkDial[fZExpMaxSyst]

private

Definition at line 103 of file GReWeightNuXSecCCQE.h.

const int GReWeightNuXSecCCQE::kModeMa = 0

static

Definition at line 48 of file GReWeightNuXSecCCQE.h.

const int GReWeightNuXSecCCQE::kModeNormAndMaShape = 1

static

Definition at line 49 of file GReWeightNuXSecCCQE.h.

const int GReWeightNuXSecCCQE::kModeZExp = 2

static

Definition at line 50 of file GReWeightNuXSecCCQE.h.

---

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

• GENIE/src/ReWeight/GReWeightNuXSecCCQE.h
• GENIE/src/ReWeight/GReWeightNuXSecCCQE.cxx