Reweighting GENIE coherent neutrino-nucleus cross sections. More...

#include <GReWeightNuXSecCOH.h>

Inheritance diagram for genie::rew::GReWeightNuXSecCOH:

Public Member Functions
	GReWeightNuXSecCOH ()

	~GReWeightNuXSecCOH ()

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	RewNue (bool tf)

void	RewNuebar (bool tf)

void	RewNumu (bool tf)

void	RewNumubar (bool tf)

void	RewCC (bool tf)

void	RewNC (bool tf)

void	SetMaPath (string p)

void	SetR0Path (string p)

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

Private Member Functions
void	Init (void)

Private Attributes
XSecAlgorithmI *	fXSecModel

Registry *	fXSecModelConfig

bool	fRewNue
	reweight nu_e? More...

bool	fRewNuebar
	reweight nu_e_bar? More...

bool	fRewNumu
	reweight nu_mu? More...

bool	fRewNumubar
	reweight nu_mu_bar? More...

bool	fRewCC
	reweight CC? More...

bool	fRewNC
	reweight NC? More...

string	fMaPath
	M_{A} path in config Registry. More...

string	fR0Path
	R_{0} path in config Registry. More...

double	fMaTwkDial

double	fMaDef

double	fMaCurr

double	fR0TwkDial

double	fR0Def

double	fR0Curr

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

Detailed Description

Reweighting GENIE coherent neutrino-nucleus 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

Definition at line 40 of file GReWeightNuXSecCOH.h.

Constructor & Destructor Documentation

GReWeightNuXSecCOH::GReWeightNuXSecCOH ( )

Definition at line 47 of file GReWeightNuXSecCOH.cxx.

 {
   this->Init();
 }

GReWeightNuXSecCOH::~GReWeightNuXSecCOH ( )

Definition at line 52 of file GReWeightNuXSecCOH.cxx.

 {
 
 }

Member Function Documentation

double GReWeightNuXSecCOH::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 120 of file GReWeightNuXSecCOH.cxx.

 {
   Interaction * interaction = event.Summary();
 
   bool is_coh = interaction->ProcInfo().IsCoherent();
   if(!is_coh) return 1.;
 
   bool tweaked = 
       (TMath::Abs(fMaTwkDial) > controls::kASmallNum) ||
       (TMath::Abs(fR0TwkDial) > controls::kASmallNum);
   if(!tweaked) return 1.0;
 
   bool is_cc  = interaction->ProcInfo().IsWeakCC();
   bool is_nc  = interaction->ProcInfo().IsWeakNC();
   if(is_cc && !fRewCC) return 1.;
   if(is_nc && !fRewNC) return 1.;
 
   int nupdg = interaction->InitState().ProbePdg();
   if(nupdg==kPdgNuMu     && !fRewNumu   ) return 1.;
   if(nupdg==kPdgAntiNuMu && !fRewNumubar) return 1.;
   if(nupdg==kPdgNuE      && !fRewNue    ) return 1.;
   if(nupdg==kPdgAntiNuE  && !fRewNuebar ) return 1.;
 
   interaction->KinePtr()->UseSelectedKinematics();
 
   double old_xsec   = event.DiffXSec();
   double old_weight = event.Weight();
   double new_xsec   = fXSecModel->XSec(interaction, kPSxyfE);
   double new_weight = old_weight * (new_xsec/old_xsec);
 
   interaction->KinePtr()->ClearRunningValues();
 
   return new_weight;
 }

void GReWeightNuXSecCOH::Init ( void )

private

Definition at line 155 of file GReWeightNuXSecCOH.cxx.

 {
   AlgId id("genie::ReinSehgalCOHPiPXSec","Default");
 
   AlgFactory * algf = AlgFactory::Instance();
   Algorithm * alg = algf->AdoptAlgorithm(id);
 
   fXSecModel = dynamic_cast<XSecAlgorithmI*>(alg);
   fXSecModel->AdoptSubstructure();
 
   fXSecModelConfig = new Registry(fXSecModel->GetConfig());
 //LOG("ReW", pNOTICE) << *fXSecModelConfig;
 
   this->RewNue    (true);
   this->RewNuebar (true);
   this->RewNumu   (true);
   this->RewNumubar(true);
   this->RewCC     (true);
   this->RewNC     (true);
 
   this->SetMaPath("Ma");
   this->SetR0Path("Ro");
 
   fMaTwkDial   = 0.; 
   fMaDef       = fXSecModelConfig->GetDouble(fMaPath);
   fMaCurr      = fMaDef;
   fR0TwkDial   = 0.; 
   fR0Def       = fXSecModelConfig->GetDouble(fR0Path);
   fR0Curr      = fR0Def;
 }

bool GReWeightNuXSecCOH::IsHandled ( GSyst_t syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 57 of file GReWeightNuXSecCOH.cxx.

 {
    bool handle;
 
    switch(syst) {
      case ( kXSecTwkDial_MaCOHpi ) : 
      case ( kXSecTwkDial_R0COHpi ) : 
        handle = true;  
        break;
      default:
        handle = false;
        break;
    }
    return handle;
 }

void GReWeightNuXSecCOH::Reconfigure ( void )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 97 of file GReWeightNuXSecCOH.cxx.

 {
   GSystUncertainty * fracerr = GSystUncertainty::Instance();
 
   double fracerr_ma = fracerr->OneSigmaErr(kXSecTwkDial_MaCOHpi);
   double fracerr_r0 = fracerr->OneSigmaErr(kXSecTwkDial_R0COHpi);
 
   fMaCurr = fMaDef * (1. + fMaTwkDial * fracerr_ma);
   fR0Curr = fR0Def * (1. + fR0TwkDial * fracerr_r0);
 
   fMaCurr = TMath::Max(0., fMaCurr  );
   fR0Curr = TMath::Max(0., fR0Curr  );
 
   Registry & r = *fXSecModelConfig;
 
   r.Set(fMaPath, fMaCurr); 
   r.Set(fR0Path, fR0Curr); 
 
   fXSecModel->Configure(r);
 
 //LOG("ReW", pDEBUG) << *fXSecModel;
 }

void GReWeightNuXSecCOH::Reset ( void )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 87 of file GReWeightNuXSecCOH.cxx.

 {
   fMaTwkDial   = 0.; 
   fMaCurr      = fMaDef;
   fR0TwkDial   = 0.; 
   fR0Curr      = fR0Def;
 
   this->Reconfigure();
 }

void genie::rew::GReWeightNuXSecCOH::RewCC ( bool tf )

inline

Definition at line 58 of file GReWeightNuXSecCOH.h.

58 { fRewCC = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCOH::fRewCC

bool fRewCC

reweight CC?

Definition: GReWeightNuXSecCOH.h:74

void genie::rew::GReWeightNuXSecCOH::RewNC ( bool tf )

inline

Definition at line 59 of file GReWeightNuXSecCOH.h.

59 { fRewNC = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCOH::fRewNC

bool fRewNC

reweight NC?

Definition: GReWeightNuXSecCOH.h:75

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

inline

Definition at line 54 of file GReWeightNuXSecCOH.h.

54 { fRewNue = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCOH::fRewNue

bool fRewNue

reweight nu_e?

Definition: GReWeightNuXSecCOH.h:70

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

inline

Definition at line 55 of file GReWeightNuXSecCOH.h.

55 { fRewNuebar = tf; }

genie::rew::GReWeightNuXSecCOH::fRewNuebar

bool fRewNuebar

reweight nu_e_bar?

Definition: GReWeightNuXSecCOH.h:71

tf

Definition: tf_graph.h:23

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

inline

Definition at line 56 of file GReWeightNuXSecCOH.h.

56 { fRewNumu = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCOH::fRewNumu

bool fRewNumu

reweight nu_mu?

Definition: GReWeightNuXSecCOH.h:72

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

inline

Definition at line 57 of file GReWeightNuXSecCOH.h.

57 { fRewNumubar = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCOH::fRewNumubar

bool fRewNumubar

reweight nu_mu_bar?

Definition: GReWeightNuXSecCOH.h:73

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

inline

Definition at line 60 of file GReWeightNuXSecCOH.h.

60 { fMaPath = p; }

test.p

p

Definition: test.py:228

genie::rew::GReWeightNuXSecCOH::fMaPath

string fMaPath

M_{A} path in config Registry.

Definition: GReWeightNuXSecCOH.h:76

void genie::rew::GReWeightNuXSecCOH::SetR0Path ( string p )

inline

Definition at line 61 of file GReWeightNuXSecCOH.h.

61 { fR0Path = p; }

test.p

p

Definition: test.py:228

genie::rew::GReWeightNuXSecCOH::fR0Path

string fR0Path

R_{0} path in config Registry.

Definition: GReWeightNuXSecCOH.h:77

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 73 of file GReWeightNuXSecCOH.cxx.

 {
    switch(syst) {
      case ( kXSecTwkDial_MaCOHpi ) : 
        fMaTwkDial = twk_dial;
        break;
      case ( kXSecTwkDial_R0COHpi ) : 
        fR0TwkDial = twk_dial;
        break;
      default:
        break;
    }
 }

Member Data Documentation

double genie::rew::GReWeightNuXSecCOH::fMaCurr

private

Definition at line 80 of file GReWeightNuXSecCOH.h.

double genie::rew::GReWeightNuXSecCOH::fMaDef

private

Definition at line 79 of file GReWeightNuXSecCOH.h.

string genie::rew::GReWeightNuXSecCOH::fMaPath

private

M_{A} path in config Registry.

Definition at line 76 of file GReWeightNuXSecCOH.h.

double genie::rew::GReWeightNuXSecCOH::fMaTwkDial

private

Definition at line 78 of file GReWeightNuXSecCOH.h.

double genie::rew::GReWeightNuXSecCOH::fR0Curr

private

Definition at line 83 of file GReWeightNuXSecCOH.h.

double genie::rew::GReWeightNuXSecCOH::fR0Def

private

Definition at line 82 of file GReWeightNuXSecCOH.h.

string genie::rew::GReWeightNuXSecCOH::fR0Path

private

R_{0} path in config Registry.

Definition at line 77 of file GReWeightNuXSecCOH.h.

double genie::rew::GReWeightNuXSecCOH::fR0TwkDial

private

Definition at line 81 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewCC

private

reweight CC?

Definition at line 74 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewNC

private

reweight NC?

Definition at line 75 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewNue

private

reweight nu_e?

Definition at line 70 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewNuebar

private

reweight nu_e_bar?

Definition at line 71 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewNumu

private

reweight nu_mu?

Definition at line 72 of file GReWeightNuXSecCOH.h.

bool genie::rew::GReWeightNuXSecCOH::fRewNumubar

private

reweight nu_mu_bar?

Definition at line 73 of file GReWeightNuXSecCOH.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCOH::fXSecModel

private

Definition at line 67 of file GReWeightNuXSecCOH.h.

Registry* genie::rew::GReWeightNuXSecCOH::fXSecModelConfig

private

Definition at line 68 of file GReWeightNuXSecCOH.h.

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

GENIE/src/ReWeight/GReWeightNuXSecCOH.h
GENIE/src/ReWeight/GReWeightNuXSecCOH.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation