Reweight GENIE CC resonance neutrino-production. More...

#include <GReWeightNuXSecCCRES.h>

Inheritance diagram for genie::rew::GReWeightNuXSecCCRES:

Public Member Functions
	GReWeightNuXSecCCRES ()

	~GReWeightNuXSecCCRES ()

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 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...

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 CC resonance neutrino-production.

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 45 of file GReWeightNuXSecCCRES.h.

Constructor & Destructor Documentation

GReWeightNuXSecCCRES::GReWeightNuXSecCCRES ( )

Definition at line 56 of file GReWeightNuXSecCCRES.cxx.

 {
   this->Init();
 }

GReWeightNuXSecCCRES::~GReWeightNuXSecCCRES ( )

Definition at line 61 of file GReWeightNuXSecCCRES.cxx.

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

Member Function Documentation

double GReWeightNuXSecCCRES::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 170 of file GReWeightNuXSecCCRES.cxx.

 {
   bool is_res = event.Summary()->ProcInfo().IsResonant();
   bool is_cc  = event.Summary()->ProcInfo().IsWeakCC();
   if(!is_res || !is_cc) 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 GReWeightNuXSecCCRES::CalcWeightMaMv ( const EventRecord & event )

private

Definition at line 249 of file GReWeightNuXSecCCRES.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 GReWeightNuXSecCCRES::CalcWeightMaMvShape ( const EventRecord & event )

private

Definition at line 270 of file GReWeightNuXSecCCRES.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();
 
 #ifdef _G_REWEIGHT_CCRES_DEBUG_
   double E  = interaction->InitState().ProbeE(kRfHitNucRest);
   double Q2 = interaction->Kine().Q2(true);
   double W  = interaction->Kine().W(true);
   fTestNtp->Fill(E,Q2,W,new_weight);
 #endif
 
   return new_weight;
 }

double GReWeightNuXSecCCRES::CalcWeightNorm ( const EventRecord & event )

private

Definition at line 240 of file GReWeightNuXSecCCRES.cxx.

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

void GReWeightNuXSecCCRES::Init ( void )

private

Definition at line 197 of file GReWeightNuXSecCCRES.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;
 
 #ifdef _G_REWEIGHT_CCRES_DEBUG_
   fTestFile = new TFile("./ccres_reweight_test.root","recreate");
   fTestNtp  = new TNtupleD("testntp","","E:Q2:W:wght");
 #endif
 }

bool GReWeightNuXSecCCRES::IsHandled ( GSyst_t syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 71 of file GReWeightNuXSecCCRES.cxx.

 {
    bool handle;
 
    switch(syst) {
 
      case ( kXSecTwkDial_NormCCRES    ) :
      case ( kXSecTwkDial_MaCCRESshape ) :
      case ( kXSecTwkDial_MvCCRESshape ) :
        if(fMode==kModeNormAndMaMvShape) { 
           handle = true;  
        } else { 
           handle = false; 
        }
        break;
 
      case ( kXSecTwkDial_MaCCRES ) :
      case ( kXSecTwkDial_MvCCRES ) :
        if(fMode==kModeMaMv) { 
           handle = true;  
        } else { 
           handle = false; 
        }
        break;
 
      default:
           handle = false;
           break;
    }
 
    return handle;
 }

void GReWeightNuXSecCCRES::Reconfigure ( void )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 137 of file GReWeightNuXSecCCRES.cxx.

 {
   GSystUncertainty * fracerr = GSystUncertainty::Instance();
 
   if(fMode==kModeMaMv) {   
      double fracerr_ma = fracerr->OneSigmaErr(kXSecTwkDial_MaCCRES);
      double fracerr_mv = fracerr->OneSigmaErr(kXSecTwkDial_MvCCRES);
      fMaCurr = fMaDef * (1. + fMaTwkDial * fracerr_ma);
      fMvCurr = fMvDef * (1. + fMvTwkDial * fracerr_mv);
   }
   else
   if(fMode==kModeNormAndMaMvShape) { 
      double fracerr_norm = fracerr->OneSigmaErr(kXSecTwkDial_NormCCRES);
      double fracerr_mash = fracerr->OneSigmaErr(kXSecTwkDial_MaCCRESshape);
      double fracerr_mvsh = fracerr->OneSigmaErr(kXSecTwkDial_MvCCRESshape);
      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 GReWeightNuXSecCCRES::Reset ( void )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 125 of file GReWeightNuXSecCCRES.cxx.

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

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

inline

Definition at line 65 of file GReWeightNuXSecCCRES.h.

65 { fRewNue = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCCRES::fRewNue

bool fRewNue

reweight nu_e CC?

Definition: GReWeightNuXSecCCRES.h:84

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

inline

Definition at line 66 of file GReWeightNuXSecCCRES.h.

66 { fRewNuebar = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCCRES::fRewNuebar

bool fRewNuebar

reweight nu_e_bar CC?

Definition: GReWeightNuXSecCCRES.h:85

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

inline

Definition at line 67 of file GReWeightNuXSecCCRES.h.

67 { fRewNumu = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCCRES::fRewNumu

bool fRewNumu

reweight nu_mu CC?

Definition: GReWeightNuXSecCCRES.h:86

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

inline

Definition at line 68 of file GReWeightNuXSecCCRES.h.

68 { fRewNumubar = tf; }

tf

Definition: tf_graph.h:23

genie::rew::GReWeightNuXSecCCRES::fRewNumubar

bool fRewNumubar

reweight nu_mu_bar CC?

Definition: GReWeightNuXSecCCRES.h:87

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

inline

Definition at line 63 of file GReWeightNuXSecCCRES.h.

63 { fMaPath = p; }

test.p

p

Definition: test.py:228

genie::rew::GReWeightNuXSecCCRES::fMaPath

string fMaPath

M_{A} path in configuration.

Definition: GReWeightNuXSecCCRES.h:82

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

inline

Definition at line 62 of file GReWeightNuXSecCCRES.h.

62 { fMode = mode; }

makeVarFluxes.mode

mode

Definition: makeVarFluxes.py:9

genie::rew::GReWeightNuXSecCCRES::fMode

int fMode

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

Definition: GReWeightNuXSecCCRES.h:81

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

inline

Definition at line 64 of file GReWeightNuXSecCCRES.h.

64 { fMvPath = p; }

test.p

p

Definition: test.py:228

genie::rew::GReWeightNuXSecCCRES::fMvPath

string fMvPath

M_{V} path in configuration.

Definition: GReWeightNuXSecCCRES.h:83

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 104 of file GReWeightNuXSecCCRES.cxx.

 {
   if(!this->IsHandled(syst)) return;
 
   switch(syst) {
     case ( kXSecTwkDial_NormCCRES ) :
       fNormTwkDial = twk_dial;
       break;
     case ( kXSecTwkDial_MaCCRESshape ) :
     case ( kXSecTwkDial_MaCCRES ) :
       fMaTwkDial = twk_dial;
       break;
     case ( kXSecTwkDial_MvCCRESshape ) :
     case ( kXSecTwkDial_MvCCRES ) :
       fMvTwkDial = twk_dial;
       break;
     default:
       break;
   }
 }

Member Data Documentation

double genie::rew::GReWeightNuXSecCCRES::fMaCurr

private

Definition at line 93 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fMaDef

private

Definition at line 92 of file GReWeightNuXSecCCRES.h.

string genie::rew::GReWeightNuXSecCCRES::fMaPath

private

M_{A} path in configuration.

Definition at line 82 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fMaTwkDial

private

Definition at line 91 of file GReWeightNuXSecCCRES.h.

int genie::rew::GReWeightNuXSecCCRES::fMode

private

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

Definition at line 81 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fMvCurr

private

Definition at line 96 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fMvDef

private

Definition at line 95 of file GReWeightNuXSecCCRES.h.

string genie::rew::GReWeightNuXSecCCRES::fMvPath

private

M_{V} path in configuration.

Definition at line 83 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fMvTwkDial

private

Definition at line 94 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fNormCurr

private

Definition at line 90 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fNormDef

private

Definition at line 89 of file GReWeightNuXSecCCRES.h.

double genie::rew::GReWeightNuXSecCCRES::fNormTwkDial

private

Definition at line 88 of file GReWeightNuXSecCCRES.h.

bool genie::rew::GReWeightNuXSecCCRES::fRewNue

private

reweight nu_e CC?

Definition at line 84 of file GReWeightNuXSecCCRES.h.

bool genie::rew::GReWeightNuXSecCCRES::fRewNuebar

private

reweight nu_e_bar CC?

Definition at line 85 of file GReWeightNuXSecCCRES.h.

bool genie::rew::GReWeightNuXSecCCRES::fRewNumu

private

reweight nu_mu CC?

Definition at line 86 of file GReWeightNuXSecCCRES.h.

bool genie::rew::GReWeightNuXSecCCRES::fRewNumubar

private

reweight nu_mu_bar CC?

Definition at line 87 of file GReWeightNuXSecCCRES.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCRES::fXSecModel

private

tweaked model

Definition at line 78 of file GReWeightNuXSecCCRES.h.

Registry* genie::rew::GReWeightNuXSecCCRES::fXSecModelConfig

private

config in tweaked model

Definition at line 79 of file GReWeightNuXSecCCRES.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCRES::fXSecModelDef

private

default model

Definition at line 77 of file GReWeightNuXSecCCRES.h.

const int GReWeightNuXSecCCRES::kModeMaMv = 0

static

Definition at line 48 of file GReWeightNuXSecCCRES.h.

const int GReWeightNuXSecCCRES::kModeNormAndMaMvShape = 1

static

Definition at line 49 of file GReWeightNuXSecCCRES.h.

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

GENIE/src/ReWeight/GReWeightNuXSecCCRES.h
GENIE/src/ReWeight/GReWeightNuXSecCCRES.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