genie::rew::GReWeightNonResonanceBkg Class Reference

Reweighting non-resonance background level. More...

#include <GReWeightNonResonanceBkg.h>

Inheritance diagram for genie::rew::GReWeightNonResonanceBkg:

Public Member Functions
	GReWeightNonResonanceBkg ()
	~GReWeightNonResonanceBkg ()
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	SetWminCut (double W)
Public Member Functions inherited from genie::rew::GReWeightI
virtual	~GReWeightI ()

Private Member Functions
void	Init (void)

Private Attributes
double	fWmin
	W_{min} cut. Reweight only events with W < W_{min}. More...
map< GSyst_t, double >	fRTwkDial
map< GSyst_t, double >	fRDef
map< GSyst_t, double >	fRCurr

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

Detailed Description

Reweighting non-resonance background level.

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 34 of file GReWeightNonResonanceBkg.h.

Constructor & Destructor Documentation

GReWeightNonResonanceBkg::GReWeightNonResonanceBkg

(

)

Definition at line 44 of file GReWeightNonResonanceBkg.cxx.

{
  this->Init();
}

GReWeightNonResonanceBkg::~GReWeightNonResonanceBkg

(

)

Definition at line 49 of file GReWeightNonResonanceBkg.cxx.

{

}

Member Function Documentation

double GReWeightNonResonanceBkg::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 138 of file GReWeightNonResonanceBkg.cxx.

{
  Interaction * interaction = event.Summary();

  bool is_dis = interaction->ProcInfo().IsDeepInelastic();
  if(!is_dis) return 1.;

  bool selected = true;
  double W = interaction->Kine().W(selected);
  bool in_transition = (W<fWmin);
  if(!in_transition) return 1.;

  int probe  = interaction->InitState().ProbePdg();
  int hitnuc = interaction->InitState().Tgt().HitNucPdg();
  InteractionType_t itype = interaction->ProcInfo().InteractionTypeId();

  int nhadmult = 0;
  int nnuc     = 0;
  int npi      = 0;

  GHepParticle * p = 0;
  TIter event_iter(&event);
  while ( (p = dynamic_cast<GHepParticle *>(event_iter.Next())) ) {

     int pdgc = p->Pdg();        
     int imom = p->FirstMother();
     if(imom == -1) continue;

     GHepParticle * mom = event.Particle(imom);
     if(!mom) continue;

     if(mom->Pdg() == kPdgHadronicSyst) 
     {
        nhadmult++;
        if ( pdg::IsNucleon(pdgc) ) { nnuc++; }
        if ( pdg::IsPion   (pdgc) ) { npi++;  }
     }
  }//p

  if(nhadmult < 2 || nhadmult > 3) return 1.;
  if(nnuc != 1) return 1.;

  GSyst_t syst = GSyst::RBkg(itype, probe, hitnuc, npi);
  if(syst == kNullSystematic) return 1.;

  double curr = fRCurr[syst];
  double def  = fRDef[syst];

  if(def>0. && curr>=0.) {
    double wght = curr / def;
    return wght;
  }

  return 1.;
}

void GReWeightNonResonanceBkg::Init ( void  )

private

Definition at line 194 of file GReWeightNonResonanceBkg.cxx.

{
  this->SetWminCut(2.0*units::GeV);

  // Get the default  parameters 
  AlgConfigPool * conf_pool = AlgConfigPool::Instance();
  Registry * user_config = conf_pool->GlobalParameterList();

  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvpCC1pi,    user_config->GetDouble("DIS-HMultWgt-vp-CC-m2" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvpCC2pi,    user_config->GetDouble("DIS-HMultWgt-vp-CC-m3" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvpNC1pi,    user_config->GetDouble("DIS-HMultWgt-vp-NC-m2" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvpNC2pi,    user_config->GetDouble("DIS-HMultWgt-vp-NC-m3" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvnCC1pi,    user_config->GetDouble("DIS-HMultWgt-vn-CC-m2" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvnCC2pi,    user_config->GetDouble("DIS-HMultWgt-vn-CC-m3" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvnNC1pi,    user_config->GetDouble("DIS-HMultWgt-vn-NC-m2" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvnNC2pi,    user_config->GetDouble("DIS-HMultWgt-vn-NC-m3" )));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarpCC1pi, user_config->GetDouble("DIS-HMultWgt-vbp-CC-m2")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarpCC2pi, user_config->GetDouble("DIS-HMultWgt-vbp-CC-m3")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarpNC1pi, user_config->GetDouble("DIS-HMultWgt-vbp-NC-m2")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarpNC2pi, user_config->GetDouble("DIS-HMultWgt-vbp-NC-m3")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarnCC1pi, user_config->GetDouble("DIS-HMultWgt-vbn-CC-m2")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarnCC2pi, user_config->GetDouble("DIS-HMultWgt-vbn-CC-m3")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarnNC1pi, user_config->GetDouble("DIS-HMultWgt-vbn-NC-m2")));
  fRDef.insert(map<GSyst_t,double>::value_type(
     kXSecTwkDial_RvbarnNC2pi, user_config->GetDouble("DIS-HMultWgt-vbn-NC-m3")));

  map<GSyst_t,double>::const_iterator it = fRDef.begin();
  for( ; it != fRDef.end(); ++it) {
    fRCurr.insert(map<GSyst_t,double>::value_type(it->first, it->second));
  }
}

bool GReWeightNonResonanceBkg::IsHandled ( GSyst_t  syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 54 of file GReWeightNonResonanceBkg.cxx.

{
   bool handle;

   switch(syst) {
     case ( kXSecTwkDial_RvpCC1pi      ) : 
     case ( kXSecTwkDial_RvpCC2pi      ) : 
     case ( kXSecTwkDial_RvpNC1pi      ) : 
     case ( kXSecTwkDial_RvpNC2pi      ) : 
     case ( kXSecTwkDial_RvnCC1pi      ) : 
     case ( kXSecTwkDial_RvnCC2pi      ) : 
     case ( kXSecTwkDial_RvnNC1pi      ) : 
     case ( kXSecTwkDial_RvnNC2pi      ) : 
     case ( kXSecTwkDial_RvbarpCC1pi   ) : 
     case ( kXSecTwkDial_RvbarpCC2pi   ) :
     case ( kXSecTwkDial_RvbarpNC1pi   ) :
     case ( kXSecTwkDial_RvbarpNC2pi   ) : 
     case ( kXSecTwkDial_RvbarnCC1pi   ) : 
     case ( kXSecTwkDial_RvbarnCC2pi   ) : 
     case ( kXSecTwkDial_RvbarnNC1pi   ) : 
     case ( kXSecTwkDial_RvbarnNC2pi   ) : 

          handle = true;
          break;

     default:
          handle = false;
          break;
   }

   return handle;
}

void GReWeightNonResonanceBkg::Reconfigure ( void  )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 125 of file GReWeightNonResonanceBkg.cxx.

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

  map<GSyst_t,double>::iterator it = fRTwkDial.begin();
  for( ; it != fRTwkDial.end(); ++it) {
    GSyst_t syst = it->first;
    double twk_dial = it->second;
    double curr = fRDef[syst] * (1 + twk_dial * fracerr->OneSigmaErr(syst));
    fRCurr[syst] = TMath::Max(0.,curr);
  }
}

void GReWeightNonResonanceBkg::Reset ( void  )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 115 of file GReWeightNonResonanceBkg.cxx.

{
  map<GSyst_t,double>::iterator it = fRTwkDial.begin();
  for( ; it != fRTwkDial.end(); ++it) {
    it->second = 0.;
  }

  this->Reconfigure();
}

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 87 of file GReWeightNonResonanceBkg.cxx.

{
   switch(syst) {
     case ( kXSecTwkDial_RvpCC1pi      ) : 
     case ( kXSecTwkDial_RvpCC2pi      ) : 
     case ( kXSecTwkDial_RvpNC1pi      ) : 
     case ( kXSecTwkDial_RvpNC2pi      ) : 
     case ( kXSecTwkDial_RvnCC1pi      ) : 
     case ( kXSecTwkDial_RvnCC2pi      ) : 
     case ( kXSecTwkDial_RvnNC1pi      ) : 
     case ( kXSecTwkDial_RvnNC2pi      ) : 
     case ( kXSecTwkDial_RvbarpCC1pi   ) : 
     case ( kXSecTwkDial_RvbarpCC2pi   ) :
     case ( kXSecTwkDial_RvbarpNC1pi   ) :
     case ( kXSecTwkDial_RvbarpNC2pi   ) : 
     case ( kXSecTwkDial_RvbarnCC1pi   ) : 
     case ( kXSecTwkDial_RvbarnCC2pi   ) : 
     case ( kXSecTwkDial_RvbarnNC1pi   ) : 
     case ( kXSecTwkDial_RvbarnNC2pi   ) : 

          fRTwkDial[syst] = twk_dial;
          break;

     default:
          break;
   }
}

void genie::rew::GReWeightNonResonanceBkg::SetWminCut ( double  W )

inline

Definition at line 48 of file GReWeightNonResonanceBkg.h.

{ fWmin = W; }

Member Data Documentation

map<GSyst_t, double> genie::rew::GReWeightNonResonanceBkg::fRCurr

private

Definition at line 58 of file GReWeightNonResonanceBkg.h.

map<GSyst_t, double> genie::rew::GReWeightNonResonanceBkg::fRDef

private

Definition at line 57 of file GReWeightNonResonanceBkg.h.

map<GSyst_t, double> genie::rew::GReWeightNonResonanceBkg::fRTwkDial

private

Definition at line 56 of file GReWeightNonResonanceBkg.h.

double genie::rew::GReWeightNonResonanceBkg::fWmin

private

W_{min} cut. Reweight only events with W < W_{min}.

Definition at line 54 of file GReWeightNonResonanceBkg.h.

---

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

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