Reweight GENIE NC resonance neutrino-production cross section. Basically a clone of the corresponding CC code. More...

#include <GReWeightNuXSecNCRES.h>

Inheritance diagram for genie::rew::GReWeightNuXSecNCRES:

Public Member Functions
	GReWeightNuXSecNCRES ()

	~GReWeightNuXSecNCRES ()

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	SetMaPath (string p)

void	SetMvPath (string p)

void	RewNue (bool tf)

void	RewNuebar (bool tf)

void	RewNumu (bool tf)

void	RewNumubar (bool tf)

Public Member Functions inherited from genie::rew::GReWeightI
virtual	~GReWeightI ()

Static Public Attributes
static const int	kModeMaMv = 0

static const int	kModeNormAndMaMvShape = 1

Private Member Functions
void	Init (void)

double	CalcWeightNorm (const EventRecord &event)

double	CalcWeightMaMvShape (const EventRecord &event)

double	CalcWeightMaMv (const EventRecord &event)

Private Attributes
XSecAlgorithmI *	fXSecModelDef
	default model More...

XSecAlgorithmI *	fXSecModel
	tweaked model More...

Registry *	fXSecModelConfig
	config in tweaked model More...

int	fMode
	0: Ma/Mv, 1: Norm and MaShape/MvShape More...

string	fMaPath
	M_{A} path in configuration. More...

string	fMvPath
	M_{V} path in configuration. More...

bool	fRewNue
	reweight nu_e NC? More...

bool	fRewNuebar
	reweight nu_e_bar NC? More...

bool	fRewNumu
	reweight nu_mu NC? More...

bool	fRewNumubar
	reweight nu_mu_bar NC? More...

double	fNormTwkDial

double	fNormDef

double	fNormCurr

double	fMaTwkDial

double	fMaDef

double	fMaCurr

double	fMvTwkDial

double	fMvDef

double	fMvCurr

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

Detailed Description

Reweight GENIE NC resonance neutrino-production cross section. Basically a clone of the corresponding CC code.

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

Definition at line 41 of file GReWeightNuXSecNCRES.h.

Constructor & Destructor Documentation

GReWeightNuXSecNCRES::GReWeightNuXSecNCRES ( )

Definition at line 52 of file GReWeightNuXSecNCRES.cxx.

 {
   this->Init();
 }

GReWeightNuXSecNCRES::~GReWeightNuXSecNCRES ( )

Definition at line 57 of file GReWeightNuXSecNCRES.cxx.

 {
 
 }

Member Function Documentation

double GReWeightNuXSecNCRES::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 161 of file GReWeightNuXSecNCRES.cxx.

 {
   bool is_res = event.Summary()->ProcInfo().IsResonant();
   bool is_nc  = event.Summary()->ProcInfo().IsWeakNC();
   if(!is_res || !is_nc) 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==kModeMaMv) {
      double wght = this->CalcWeightMaMv(event);
      return wght;
   } 
   else 
   if(fMode==kModeNormAndMaMvShape) {
      double wght = 
          this->CalcWeightNorm      (event) *
          this->CalcWeightMaMvShape (event);
      return wght;
   }
 
   return 1.;
 }

double GReWeightNuXSecNCRES::CalcWeightMaMv ( const EventRecord & event )

private

Definition at line 235 of file GReWeightNuXSecNCRES.cxx.

 {
   bool tweaked = 
      (TMath::Abs(fMaTwkDial) > controls::kASmallNum) ||
      (TMath::Abs(fMvTwkDial) > controls::kASmallNum);
   if(!tweaked) return 1.0;
 
   Interaction * interaction = event.Summary();
 
   interaction->KinePtr()->UseSelectedKinematics();
 
   double old_xsec   = event.DiffXSec();
   double old_weight = event.Weight();
   double new_xsec   = fXSecModel->XSec(interaction, kPSWQ2fE);
   double new_weight = old_weight * (new_xsec/old_xsec);
 
   interaction->KinePtr()->ClearRunningValues();
 
   return new_weight;
 }

double GReWeightNuXSecNCRES::CalcWeightMaMvShape ( const EventRecord & event )

private

Definition at line 256 of file GReWeightNuXSecNCRES.cxx.

 {
   bool tweaked = 
      (TMath::Abs(fMaTwkDial) > controls::kASmallNum) ||
      (TMath::Abs(fMvTwkDial) > controls::kASmallNum);
   if(!tweaked) return 1.0;
 
   Interaction * interaction = event.Summary();
 
   interaction->KinePtr()->UseSelectedKinematics();
 
   double old_xsec   = event.DiffXSec();
   double old_weight = event.Weight();
   double new_xsec   = fXSecModel->XSec(interaction, kPSWQ2fE);
   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 twk_integrated_xsec = fXSecModel    -> Integral(interaction);   
   assert(twk_integrated_xsec > 0);
   new_weight *= (old_integrated_xsec/twk_integrated_xsec);
 
 //LOG("ReW", pDEBUG) << "integrated cross section (old) = " << old_integrated_xsec;
 //LOG("ReW", pDEBUG) << "integrated cross section (twk) = " << twk_integrated_xsec;
 //LOG("ReW", pDEBUG) << "new weight (normalized to const integral) = " << new_weight;
 
   interaction->KinePtr()->ClearRunningValues();
 
   return new_weight;
 }

double GReWeightNuXSecNCRES::CalcWeightNorm ( const EventRecord & event )

private

Definition at line 226 of file GReWeightNuXSecNCRES.cxx.

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

void GReWeightNuXSecNCRES::Init ( void )

private

Definition at line 188 of file GReWeightNuXSecNCRES.cxx.

 {
   AlgId id("genie::ReinSehgalRESPXSec","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());
 //LOG("ReW", pNOTICE) << *fXSecModelConfig;
 
   this->SetMode(kModeNormAndMaMvShape);
 
   this->RewNue    (true);
   this->RewNuebar (true);
   this->RewNumu   (true);
   this->RewNumubar(true);
 
   this->SetMaPath("Ma");
   this->SetMvPath("Mv");
 
   fNormTwkDial = 0.;
   fNormDef     = 1.;
   fNormCurr    = fNormDef;
   fMaTwkDial   = 0.; 
   fMaDef       = fXSecModelConfig->GetDouble(fMaPath);
   fMaCurr      = fMaDef;
   fMvTwkDial   = 0.; 
   fMvDef       = fXSecModelConfig->GetDouble(fMvPath);
   fMvCurr      = fMvDef;
 }

bool GReWeightNuXSecNCRES::IsHandled ( GSyst_t syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 62 of file GReWeightNuXSecNCRES.cxx.

 {
    bool handle;
 
    switch(syst) {
 
      case ( kXSecTwkDial_NormNCRES    ) :
      case ( kXSecTwkDial_MaNCRESshape ) :
      case ( kXSecTwkDial_MvNCRESshape ) :
        if(fMode==kModeNormAndMaMvShape) { 
           handle = true;  
        } else { 
           handle = false; 
        }
        break;
 
      case ( kXSecTwkDial_MaNCRES ) :
      case ( kXSecTwkDial_MvNCRES ) :
        if(fMode==kModeMaMv) { 
           handle = true;  
        } else { 
           handle = false; 
        }
        break;
 
      default:
           handle = false;
           break;
    }
 
    return handle;
 }

void GReWeightNuXSecNCRES::Reconfigure ( void )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 128 of file GReWeightNuXSecNCRES.cxx.

 {
   GSystUncertainty * fracerr = GSystUncertainty::Instance();
 
   if(fMode==kModeMaMv) {   
      double fracerr_ma = fracerr->OneSigmaErr(kXSecTwkDial_MaNCRES);
      double fracerr_mv = fracerr->OneSigmaErr(kXSecTwkDial_MvNCRES);
      fMaCurr = fMaDef * (1. + fMaTwkDial * fracerr_ma);
      fMvCurr = fMvDef * (1. + fMvTwkDial * fracerr_mv);
   }
   else
   if(fMode==kModeNormAndMaMvShape) { 
      double fracerr_norm = fracerr->OneSigmaErr(kXSecTwkDial_NormNCRES);
      double fracerr_mash = fracerr->OneSigmaErr(kXSecTwkDial_MaNCRESshape);
      double fracerr_mvsh = fracerr->OneSigmaErr(kXSecTwkDial_MvNCRESshape);
      fNormCurr = fNormDef * (1. + fNormTwkDial * fracerr_norm);
      fMaCurr   = fMaDef   * (1. + fMaTwkDial   * fracerr_mash);
      fMvCurr   = fMvDef   * (1. + fMvTwkDial   * fracerr_mvsh);
   }
 
   fNormCurr = TMath::Max(0., fNormCurr);
   fMaCurr   = TMath::Max(0., fMaCurr  );
   fMvCurr   = TMath::Max(0., fMvCurr  );
 
   Registry & r = *fXSecModelConfig;
 
   r.Set(fMaPath, fMaCurr); 
   r.Set(fMvPath, fMvCurr); 
   fXSecModel->Configure(r);
 
 //LOG("ReW, pDEBUG) << *fXSecModel;
 }

void GReWeightNuXSecNCRES::Reset ( void )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 116 of file GReWeightNuXSecNCRES.cxx.

 {
   fNormTwkDial = 0.;
   fNormCurr    = fNormDef;
   fMaTwkDial   = 0.; 
   fMaCurr      = fMaDef;
   fMvTwkDial   = 0.; 
   fMvCurr      = fMvDef;
 
   this->Reconfigure();
 }

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

inline

Definition at line 61 of file GReWeightNuXSecNCRES.h.

61 { fRewNue = tf; }

genie::rew::GReWeightNuXSecNCRES::fRewNue

bool fRewNue

reweight nu_e NC?

Definition: GReWeightNuXSecNCRES.h:80

tf

Definition: tf_graph.h:23

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

inline

Definition at line 62 of file GReWeightNuXSecNCRES.h.

62 { fRewNuebar = tf; }

genie::rew::GReWeightNuXSecNCRES::fRewNuebar

bool fRewNuebar

reweight nu_e_bar NC?

Definition: GReWeightNuXSecNCRES.h:81

tf

Definition: tf_graph.h:23

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

inline

Definition at line 63 of file GReWeightNuXSecNCRES.h.

63 { fRewNumu = tf; }

genie::rew::GReWeightNuXSecNCRES::fRewNumu

bool fRewNumu

reweight nu_mu NC?

Definition: GReWeightNuXSecNCRES.h:82

tf

Definition: tf_graph.h:23

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

inline

Definition at line 64 of file GReWeightNuXSecNCRES.h.

64 { fRewNumubar = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecNCRES::fRewNumubar

bool fRewNumubar

reweight nu_mu_bar NC?

Definition: GReWeightNuXSecNCRES.h:83

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

inline

Definition at line 59 of file GReWeightNuXSecNCRES.h.

59 { fMaPath = p; }

test.p

p

Definition: test.py:228

genie::rew::GReWeightNuXSecNCRES::fMaPath

string fMaPath

M_{A} path in configuration.

Definition: GReWeightNuXSecNCRES.h:78

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

inline

Definition at line 58 of file GReWeightNuXSecNCRES.h.

58 { fMode = mode; }

genie::rew::GReWeightNuXSecNCRES::fMode

int fMode

0: Ma/Mv, 1: Norm and MaShape/MvShape

Definition: GReWeightNuXSecNCRES.h:77

makeVarFluxes.mode

mode

Definition: makeVarFluxes.py:9

void genie::rew::GReWeightNuXSecNCRES::SetMvPath ( string p )

inline

Definition at line 60 of file GReWeightNuXSecNCRES.h.

60 { fMvPath = p; }

genie::rew::GReWeightNuXSecNCRES::fMvPath

string fMvPath

M_{V} path in configuration.

Definition: GReWeightNuXSecNCRES.h:79

test.p

p

Definition: test.py:228

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 95 of file GReWeightNuXSecNCRES.cxx.

 {
   if(!this->IsHandled(syst)) return;
 
   switch(syst) {
     case ( kXSecTwkDial_NormNCRES ) :
       fNormTwkDial = twk_dial;
       break;
     case ( kXSecTwkDial_MaNCRESshape ) :
     case ( kXSecTwkDial_MaNCRES ) :
       fMaTwkDial = twk_dial;
       break;
     case ( kXSecTwkDial_MvNCRESshape ) :
     case ( kXSecTwkDial_MvNCRES ) :
       fMvTwkDial = twk_dial;
       break;
     default:
       break;
   }
 }

Member Data Documentation

double genie::rew::GReWeightNuXSecNCRES::fMaCurr

private

Definition at line 89 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fMaDef

private

Definition at line 88 of file GReWeightNuXSecNCRES.h.

string genie::rew::GReWeightNuXSecNCRES::fMaPath

private

M_{A} path in configuration.

Definition at line 78 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fMaTwkDial

private

Definition at line 87 of file GReWeightNuXSecNCRES.h.

int genie::rew::GReWeightNuXSecNCRES::fMode

private

0: Ma/Mv, 1: Norm and MaShape/MvShape

Definition at line 77 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fMvCurr

private

Definition at line 92 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fMvDef

private

Definition at line 91 of file GReWeightNuXSecNCRES.h.

string genie::rew::GReWeightNuXSecNCRES::fMvPath

private

M_{V} path in configuration.

Definition at line 79 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fMvTwkDial

private

Definition at line 90 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fNormCurr

private

Definition at line 86 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fNormDef

private

Definition at line 85 of file GReWeightNuXSecNCRES.h.

double genie::rew::GReWeightNuXSecNCRES::fNormTwkDial

private

Definition at line 84 of file GReWeightNuXSecNCRES.h.

bool genie::rew::GReWeightNuXSecNCRES::fRewNue

private

reweight nu_e NC?

Definition at line 80 of file GReWeightNuXSecNCRES.h.

bool genie::rew::GReWeightNuXSecNCRES::fRewNuebar

private

reweight nu_e_bar NC?

Definition at line 81 of file GReWeightNuXSecNCRES.h.

bool genie::rew::GReWeightNuXSecNCRES::fRewNumu

private

reweight nu_mu NC?

Definition at line 82 of file GReWeightNuXSecNCRES.h.

bool genie::rew::GReWeightNuXSecNCRES::fRewNumubar

private

reweight nu_mu_bar NC?

Definition at line 83 of file GReWeightNuXSecNCRES.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecNCRES::fXSecModel

private

tweaked model

Definition at line 74 of file GReWeightNuXSecNCRES.h.

Registry* genie::rew::GReWeightNuXSecNCRES::fXSecModelConfig

private

config in tweaked model

Definition at line 75 of file GReWeightNuXSecNCRES.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecNCRES::fXSecModelDef

private

default model

Definition at line 73 of file GReWeightNuXSecNCRES.h.

const int GReWeightNuXSecNCRES::kModeMaMv = 0

static

Definition at line 44 of file GReWeightNuXSecNCRES.h.

const int GReWeightNuXSecNCRES::kModeNormAndMaMvShape = 1

static

Definition at line 45 of file GReWeightNuXSecNCRES.h.

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

GENIE/src/ReWeight/GReWeightNuXSecNCRES.h
GENIE/src/ReWeight/GReWeightNuXSecNCRES.cxx

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation