genie::rew::GReWeightNuXSecNCEL Class Reference

Reweighting NCEL GENIE neutrino cross sections. More...

#include <GReWeightNuXSecNCEL.h>

Inheritance diagram for genie::rew::GReWeightNuXSecNCEL:

Public Member Functions
	GReWeightNuXSecNCEL ()
	~GReWeightNuXSecNCEL ()
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	SetMaPath (string p)
void	SetEtaPath (string p)
Public Member Functions inherited from genie::rew::GReWeightI
virtual	~GReWeightI ()

Private Member Functions
void	Init (void)

Private Attributes
XSecAlgorithmI *	fXSecModelDef
	default model More...
XSecAlgorithmI *	fXSecModel
	tweaked model More...
Registry *	fXSecModelConfig
	config in tweaked model 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...
string	fEtaPath
	eta path in config Registry More...
double	fMaTwkDial
double	fMaDef
double	fMaCurr
double	fEtaTwkDial
double	fEtaDef
double	fEtaCurr

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

Detailed Description

Reweighting NCEL GENIE neutrino cross sections.

Author

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

Nov 25, 2010

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 42 of file GReWeightNuXSecNCEL.h.

Constructor & Destructor Documentation

GReWeightNuXSecNCEL::GReWeightNuXSecNCEL

(

)

Definition at line 43 of file GReWeightNuXSecNCEL.cxx.

{
  this->Init();
}

GReWeightNuXSecNCEL::~GReWeightNuXSecNCEL

(

)

Definition at line 48 of file GReWeightNuXSecNCEL.cxx.

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

Member Function Documentation

double GReWeightNuXSecNCEL::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 122 of file GReWeightNuXSecNCEL.cxx.

{
  Interaction * interaction = event.Summary();

  bool is_qe = interaction->ProcInfo().IsQuasiElastic();
  bool is_nc = interaction->ProcInfo().IsWeakNC();
  if(!is_qe || !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.;

  bool tweaked_ma  = (TMath::Abs(fMaTwkDial ) > controls::kASmallNum);
  bool tweaked_eta = (TMath::Abs(fEtaTwkDial) > controls::kASmallNum);
  bool tweaked     = tweaked_ma || tweaked_eta;
  if(!tweaked) return 1.0;

  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;

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

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

  return new_weight;
}

void GReWeightNuXSecNCEL::Init ( void  )

private

Definition at line 166 of file GReWeightNuXSecNCEL.cxx.

{
  AlgId id("genie::AhrensNCELPXSec","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", pDEBUG) << *fXSecModelConfig;

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

  this->SetMaPath ("Ma");
  this->SetEtaPath("Eta");

  fMaTwkDial   = 0.; 
  fMaDef       = fXSecModelConfig->GetDouble(fMaPath);
  fMaCurr      = fMaDef;
  fEtaTwkDial  = 0.; 
  fEtaDef      = fXSecModelConfig->GetDouble(fEtaPath);
  fEtaCurr     = fEtaDef;

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

bool GReWeightNuXSecNCEL::IsHandled ( GSyst_t  syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 58 of file GReWeightNuXSecNCEL.cxx.

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

void GReWeightNuXSecNCEL::Reconfigure ( void  )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 97 of file GReWeightNuXSecNCEL.cxx.

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

  int sign_matwk  = utils::rew::Sign(fMaTwkDial );
  int sign_etatwk = utils::rew::Sign(fEtaTwkDial);
  
  double fracerr_ma  = fracerr->OneSigmaErr(kXSecTwkDial_MaNCEL,  sign_matwk );
  double fracerr_eta = fracerr->OneSigmaErr(kXSecTwkDial_EtaNCEL, sign_etatwk);
  
  fMaCurr  = fMaDef  * (1. + fMaTwkDial  * fracerr_ma);
  fEtaCurr = fEtaDef * (1. + fEtaTwkDial * fracerr_eta);

  fMaCurr  = TMath::Max(0., fMaCurr  );
  fEtaCurr = TMath::Max(0., fEtaCurr );

  Registry & r = *fXSecModelConfig;

  r.Set(fMaPath,  fMaCurr ); 
  r.Set(fEtaPath, fEtaCurr); 
  fXSecModel->Configure(r);

//LOG("ReW, pDEBUG) << *fXSecModel;
}

void GReWeightNuXSecNCEL::Reset ( void  )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 87 of file GReWeightNuXSecNCEL.cxx.

{
  fMaTwkDial   = 0.; 
  fMaCurr      = fMaDef;
  fEtaTwkDial  = 0.; 
  fEtaCurr     = fEtaDef;

  this->Reconfigure();
}

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

inline

Definition at line 56 of file GReWeightNuXSecNCEL.h.

{ fRewNue     = tf;   }

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

inline

Definition at line 57 of file GReWeightNuXSecNCEL.h.

{ fRewNuebar  = tf;   }

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

inline

Definition at line 58 of file GReWeightNuXSecNCEL.h.

{ fRewNumu    = tf;   }

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

inline

Definition at line 59 of file GReWeightNuXSecNCEL.h.

{ fRewNumubar = tf;   }

void genie::rew::GReWeightNuXSecNCEL::SetEtaPath ( string  p )

inline

Definition at line 61 of file GReWeightNuXSecNCEL.h.

{ fEtaPath    = p;    }

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

inline

Definition at line 60 of file GReWeightNuXSecNCEL.h.

{ fMaPath     = p;    }

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 71 of file GReWeightNuXSecNCEL.cxx.

{
  if(!this->IsHandled(syst)) return;

  switch(syst) {
    case ( kXSecTwkDial_MaNCEL ) :
      fMaTwkDial = twk_dial;
      break;
    case ( kXSecTwkDial_EtaNCEL ) :
      fEtaTwkDial = twk_dial;
      break;
    default:
      break;
  }
}

Member Data Documentation

double genie::rew::GReWeightNuXSecNCEL::fEtaCurr

private

Definition at line 82 of file GReWeightNuXSecNCEL.h.

double genie::rew::GReWeightNuXSecNCEL::fEtaDef

private

Definition at line 81 of file GReWeightNuXSecNCEL.h.

string genie::rew::GReWeightNuXSecNCEL::fEtaPath

private

eta path in config Registry

Definition at line 76 of file GReWeightNuXSecNCEL.h.

double genie::rew::GReWeightNuXSecNCEL::fEtaTwkDial

private

Definition at line 80 of file GReWeightNuXSecNCEL.h.

double genie::rew::GReWeightNuXSecNCEL::fMaCurr

private

Definition at line 79 of file GReWeightNuXSecNCEL.h.

double genie::rew::GReWeightNuXSecNCEL::fMaDef

private

Definition at line 78 of file GReWeightNuXSecNCEL.h.

string genie::rew::GReWeightNuXSecNCEL::fMaPath

private

M_{A} path in config Registry.

Definition at line 75 of file GReWeightNuXSecNCEL.h.

double genie::rew::GReWeightNuXSecNCEL::fMaTwkDial

private

Definition at line 77 of file GReWeightNuXSecNCEL.h.

bool genie::rew::GReWeightNuXSecNCEL::fRewNue

private

reweight nu_e CC?

Definition at line 71 of file GReWeightNuXSecNCEL.h.

bool genie::rew::GReWeightNuXSecNCEL::fRewNuebar

private

reweight nu_e_bar CC?

Definition at line 72 of file GReWeightNuXSecNCEL.h.

bool genie::rew::GReWeightNuXSecNCEL::fRewNumu

private

reweight nu_mu CC?

Definition at line 73 of file GReWeightNuXSecNCEL.h.

bool genie::rew::GReWeightNuXSecNCEL::fRewNumubar

private

reweight nu_mu_bar CC?

Definition at line 74 of file GReWeightNuXSecNCEL.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecNCEL::fXSecModel

private

tweaked model

Definition at line 68 of file GReWeightNuXSecNCEL.h.

Registry* genie::rew::GReWeightNuXSecNCEL::fXSecModelConfig

private

config in tweaked model

Definition at line 69 of file GReWeightNuXSecNCEL.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecNCEL::fXSecModelDef

private

default model

Definition at line 67 of file GReWeightNuXSecNCEL.h.

---

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

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