genie::rew::GReWeightNuXSecCCQEvec Class Reference

Reweighting vector form factors in GENIE CCQE neutrino cross section calculations. More...

#include <GReWeightNuXSecCCQEvec.h>

Inheritance diagram for genie::rew::GReWeightNuXSecCCQEvec:

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

Private Member Functions
void	Init (void)

Private Attributes
XSecAlgorithmI *	fXSecModel_bba
	CCQE model with BBA05 f/f (default) More...
XSecAlgorithmI *	fXSecModel_dpl
	CCQE model with dipole f/f ("maximally" tweaked) More...
double	fFFTwkDial
	tweaking dial (0: bba/default, +1: dipole) 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...

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

Detailed Description

Reweighting vector form factors in GENIE CCQE neutrino cross section calculations.

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

May 24, 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 45 of file GReWeightNuXSecCCQEvec.h.

Constructor & Destructor Documentation

GReWeightNuXSecCCQEvec::GReWeightNuXSecCCQEvec

(

)

Definition at line 46 of file GReWeightNuXSecCCQEvec.cxx.

{
  this->Init();
}

GReWeightNuXSecCCQEvec::~GReWeightNuXSecCCQEvec

(

)

Definition at line 51 of file GReWeightNuXSecCCQEvec.cxx.

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

Member Function Documentation

double GReWeightNuXSecCCQEvec::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 97 of file GReWeightNuXSecCCQEvec.cxx.

{
  bool tweaked = (TMath::Abs(fFFTwkDial) > controls::kASmallNum);
  if(!tweaked) return 1.;

  Interaction * interaction = event.Summary();

  bool is_qe = interaction->ProcInfo().IsQuasiElastic();
  bool is_cc = interaction->ProcInfo().IsWeakCC();
  if(!is_qe || !is_cc) return 1.;

  bool charm = interaction->ExclTag().IsCharmEvent(); // skip CCQE charm
  if(charm) 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.;
 
  //
  // Calculate weight
  // Input tweaking dial changes elastic nucleon form factors
  // (twk dial: 0 -> default/BBA, twk dial: 1 -> dipole).
  // Calculated weight includes `shape only effect in dsigma/dQ2
  // (normalized to constant integrated cross section)
  //

  interaction->KinePtr()->UseSelectedKinematics();
  interaction->SetBit(kIAssumeFreeNucleon);

  double dial                = fFFTwkDial;
  double old_weight          = event.Weight();
  double def_xsec            = event.DiffXSec();
  double dpl_xsec            = fXSecModel_dpl->XSec(interaction, kPSQ2fE);
  double def_integrated_xsec = fXSecModel_bba->Integral(interaction);
  double dpl_integrated_xsec = fXSecModel_dpl->Integral(interaction);

  assert(def_integrated_xsec > 0.);
  assert(dpl_integrated_xsec > 0.);
//  if(def_integrated_xsec <= 0 || dpl_integrated_xsec <= 0) return 1.;

  double def_ratio = def_xsec / def_integrated_xsec;
  double dpl_ratio = dpl_xsec / dpl_integrated_xsec;

  assert(def_ratio > 0.);
//  if(def_ratio <= 0) return 1.;

  double weight = old_weight * (dial * dpl_ratio + (1-dial)*def_ratio) / def_ratio;

#ifdef _G_REWEIGHT_CCQE_VEC_DEBUG_
  double E  = interaction->InitState().ProbeE(kRfHitNucRest);
  double Q2 = interaction->Kine().Q2(true);
  fTestNtp->Fill(
    E,Q2,weight,def_integrated_xsec,dpl_integrated_xsec,def_xsec,dpl_xsec);
#endif


  return weight;
}

void GReWeightNuXSecCCQEvec::Init ( void  )

private

Definition at line 158 of file GReWeightNuXSecCCQEvec.cxx.

{
  AlgFactory * algf = AlgFactory::Instance();

  AlgId id0("genie::LwlynSmithQELCCPXSec","Default");
  Algorithm * alg0 = algf->AdoptAlgorithm(id0);
  fXSecModel_bba = dynamic_cast<XSecAlgorithmI*>(alg0);
  fXSecModel_bba->AdoptSubstructure();

  AlgId id1("genie::LwlynSmithQELCCPXSec","DipoleELFF");
  Algorithm * alg1 = algf->AdoptAlgorithm(id1);
  fXSecModel_dpl = dynamic_cast<XSecAlgorithmI*>(alg1);
  fXSecModel_dpl->AdoptSubstructure();

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

  fFFTwkDial = 0.;

#ifdef _G_REWEIGHT_CCQE_VEC_DEBUG_
  fTestFile = new TFile("./ccqevec_reweight_test.root","recreate");
  fTestNtp  = new TNtupleD("testntp","","E:Q2:wght:sig0:sig:dsig0:dsig");
#endif

}

bool GReWeightNuXSecCCQEvec::IsHandled ( GSyst_t  syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 61 of file GReWeightNuXSecCCQEvec.cxx.

{
   switch(syst) {
    case ( kXSecTwkDial_VecFFCCQEshape ) : 
       return true;
       break;
    default:
       return false;
       break;
   }
   return false;
}

void GReWeightNuXSecCCQEvec::Reconfigure ( void  )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 92 of file GReWeightNuXSecCCQEvec.cxx.

{

}

void GReWeightNuXSecCCQEvec::Reset ( void  )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 87 of file GReWeightNuXSecCCQEvec.cxx.

{
  fFFTwkDial = 0.;
}

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

inline

Definition at line 59 of file GReWeightNuXSecCCQEvec.h.

{ fRewNue     = tf;   }

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

inline

Definition at line 60 of file GReWeightNuXSecCCQEvec.h.

{ fRewNuebar  = tf;   }

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

inline

Definition at line 61 of file GReWeightNuXSecCCQEvec.h.

{ fRewNumu    = tf;   }

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

inline

Definition at line 62 of file GReWeightNuXSecCCQEvec.h.

{ fRewNumubar = tf;   }

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 74 of file GReWeightNuXSecCCQEvec.cxx.

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

  switch(syst) {
    case ( kXSecTwkDial_VecFFCCQEshape ) : 
       fFFTwkDial = twk_dial;
       break;
    default:
       return;
  }
}

Member Data Documentation

double genie::rew::GReWeightNuXSecCCQEvec::fFFTwkDial

private

tweaking dial (0: bba/default, +1: dipole)

Definition at line 71 of file GReWeightNuXSecCCQEvec.h.

bool genie::rew::GReWeightNuXSecCCQEvec::fRewNue

private

reweight nu_e CC?

Definition at line 73 of file GReWeightNuXSecCCQEvec.h.

bool genie::rew::GReWeightNuXSecCCQEvec::fRewNuebar

private

reweight nu_e_bar CC?

Definition at line 74 of file GReWeightNuXSecCCQEvec.h.

bool genie::rew::GReWeightNuXSecCCQEvec::fRewNumu

private

reweight nu_mu CC?

Definition at line 75 of file GReWeightNuXSecCCQEvec.h.

bool genie::rew::GReWeightNuXSecCCQEvec::fRewNumubar

private

reweight nu_mu_bar CC?

Definition at line 76 of file GReWeightNuXSecCCQEvec.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQEvec::fXSecModel_bba

private

CCQE model with BBA05 f/f (default)

Definition at line 68 of file GReWeightNuXSecCCQEvec.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQEvec::fXSecModel_dpl

private

CCQE model with dipole f/f ("maximally" tweaked)

Definition at line 69 of file GReWeightNuXSecCCQEvec.h.

---

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

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