genie::rew::GReWeightNuXSecCCQEaxial Class Reference

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

#include <GReWeightNuXSecCCQEaxial.h>

Inheritance diagram for genie::rew::GReWeightNuXSecCCQEaxial:

Public Member Functions
	GReWeightNuXSecCCQEaxial ()
	~GReWeightNuXSecCCQEaxial ()
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...
double	CalcChisq (void)
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_dpl
	CCQE model with dipole f/f (default) More...
XSecAlgorithmI *	fXSecModel_zexp
	CCQE model with z-expansion 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...
TFile *	fTestFile
TNtupleD *	fTestNtp

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> STFC, Rutherford Appleton Laboratory

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 43 of file GReWeightNuXSecCCQEaxial.h.

Constructor & Destructor Documentation

GReWeightNuXSecCCQEaxial::GReWeightNuXSecCCQEaxial

(

)

Definition at line 51 of file GReWeightNuXSecCCQEaxial.cxx.

{
  this->Init();
}

GReWeightNuXSecCCQEaxial::~GReWeightNuXSecCCQEaxial

(

)

Definition at line 56 of file GReWeightNuXSecCCQEaxial.cxx.

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

Member Function Documentation

double GReWeightNuXSecCCQEaxial::CalcChisq

(

void

)

Definition at line 165 of file GReWeightNuXSecCCQEaxial.cxx.

{
  double chisq = TMath::Power(fFFTwkDial, 2.);
  return chisq;
}

double GReWeightNuXSecCCQEaxial::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 102 of file GReWeightNuXSecCCQEaxial.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/dipole, twk dial: 1 -> zexp).
  // 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 zexp_xsec            = fXSecModel_zexp->XSec(interaction, kPSQ2fE);
  //double def_integrated_xsec  = fXSecModel_dpl->Integral(interaction);
  //double zexp_integrated_xsec = fXSecModel_zexp->Integral(interaction);

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

  //double def_ratio  = def_xsec  / def_integrated_xsec;
  //double zexp_ratio = zexp_xsec / zexp_integrated_xsec;
  double def_ratio  = def_xsec ;
  double zexp_ratio = zexp_xsec;

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

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

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


  return weight;
}

void GReWeightNuXSecCCQEaxial::Init ( void  )

private

Definition at line 171 of file GReWeightNuXSecCCQEaxial.cxx.

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

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

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

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

  fFFTwkDial = 0.;

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

}

bool GReWeightNuXSecCCQEaxial::IsHandled ( GSyst_t  syst )

virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 66 of file GReWeightNuXSecCCQEaxial.cxx.

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

void GReWeightNuXSecCCQEaxial::Reconfigure ( void  )

virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 97 of file GReWeightNuXSecCCQEaxial.cxx.

{

}

void GReWeightNuXSecCCQEaxial::Reset ( void  )

virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 92 of file GReWeightNuXSecCCQEaxial.cxx.

{
  fFFTwkDial = 0.;
}

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

inline

Definition at line 58 of file GReWeightNuXSecCCQEaxial.h.

{ fRewNue     = tf;   }

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

inline

Definition at line 59 of file GReWeightNuXSecCCQEaxial.h.

{ fRewNuebar  = tf;   }

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

inline

Definition at line 60 of file GReWeightNuXSecCCQEaxial.h.

{ fRewNumu    = tf;   }

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

inline

Definition at line 61 of file GReWeightNuXSecCCQEaxial.h.

{ fRewNumubar = tf;   }

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

virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 79 of file GReWeightNuXSecCCQEaxial.cxx.

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

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

Member Data Documentation

double genie::rew::GReWeightNuXSecCCQEaxial::fFFTwkDial

private

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

Definition at line 70 of file GReWeightNuXSecCCQEaxial.h.

bool genie::rew::GReWeightNuXSecCCQEaxial::fRewNue

private

reweight nu_e CC?

Definition at line 72 of file GReWeightNuXSecCCQEaxial.h.

bool genie::rew::GReWeightNuXSecCCQEaxial::fRewNuebar

private

reweight nu_e_bar CC?

Definition at line 73 of file GReWeightNuXSecCCQEaxial.h.

bool genie::rew::GReWeightNuXSecCCQEaxial::fRewNumu

private

reweight nu_mu CC?

Definition at line 74 of file GReWeightNuXSecCCQEaxial.h.

bool genie::rew::GReWeightNuXSecCCQEaxial::fRewNumubar

private

reweight nu_mu_bar CC?

Definition at line 75 of file GReWeightNuXSecCCQEaxial.h.

TFile* genie::rew::GReWeightNuXSecCCQEaxial::fTestFile

private

Definition at line 77 of file GReWeightNuXSecCCQEaxial.h.

TNtupleD* genie::rew::GReWeightNuXSecCCQEaxial::fTestNtp

private

Definition at line 78 of file GReWeightNuXSecCCQEaxial.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQEaxial::fXSecModel_dpl

private

CCQE model with dipole f/f (default)

Definition at line 67 of file GReWeightNuXSecCCQEaxial.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecCCQEaxial::fXSecModel_zexp

private

CCQE model with z-expansion f/f ("maximally" tweaked)

Definition at line 68 of file GReWeightNuXSecCCQEaxial.h.

---

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

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