genie::Kinematics Class Reference

Generated/set kinematical variables for an event. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/genie/v3_00_06_p01/Linux64bit+3.10-2.17-e17-debug/GENIE-Generator/src/Framework/Interaction/Kinematics.h"

Inheritance diagram for genie::Kinematics:

Public Member Functions
	Kinematics ()
	Kinematics (const Kinematics &kv)
	Kinematics (TRootIOCtor *)
	~Kinematics ()
double	x (bool selected=false) const
double	y (bool selected=false) const
double	Q2 (bool selected=false) const
double	q2 (bool selected=false) const
double	W (bool selected=false) const
double	t (bool selected=false) const
double	Logx (bool selected=false) const
double	Logy (bool selected=false) const
double	LogQ2 (bool selected=false) const
double	LogW (bool selected=false) const
double	Log10x (bool selected=false) const
double	Log10y (bool selected=false) const
double	Log10Q2 (bool selected=false) const
double	Log10W (bool selected=false) const
const TLorentzVector &	FSLeptonP4 (void) const
const TLorentzVector &	HadSystP4 (void) const
void	Setx (double x, bool selected=false)
void	Sety (double y, bool selected=false)
void	SetQ2 (double Q2, bool selected=false)
void	Setq2 (double q2, bool selected=false)
void	SetW (double W, bool selected=false)
void	Sett (double t, bool selected=false)
void	SetFSLeptonP4 (const TLorentzVector &p4)
void	SetFSLeptonP4 (double px, double py, double pz, double E)
void	SetHadSystP4 (const TLorentzVector &p4)
void	SetHadSystP4 (double px, double py, double pz, double E)
bool	KVSet (KineVar_t kv) const
double	GetKV (KineVar_t kv) const
void	SetKV (KineVar_t kv, double value)
void	ClearRunningValues (void)
void	UseSelectedKinematics (void)
void	Reset (void)
void	Copy (const Kinematics &kine)
void	Print (ostream &stream) const
Kinematics &	operator= (const Kinematics &kine)

Private Member Functions
void	Init (void)
	initialize More...
void	CleanUp (void)
	clean-up More...

Private Attributes
map< KineVar_t, double >	fKV
	selected kinematics More...
TLorentzVector *	fP4Fsl
	generated final state primary lepton 4-p (LAB) More...
TLorentzVector *	fP4HadSyst
	generated final state hadronic system 4-p (LAB) More...

Friends
ostream &	operator<< (ostream &stream, const Kinematics &kine)

Detailed Description

Generated/set kinematical variables for an event.

Author

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

May 08, 2004

Copyright (c) 2003-2019, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Definition at line 40 of file Kinematics.h.

Constructor & Destructor Documentation

Kinematics::Kinematics

(

)

Definition at line 43 of file Kinematics.cxx.

References Init().

                       :
TObject()
{
  this->Init();
}

Kinematics::Kinematics

(

const Kinematics &

kv

)

Definition at line 49 of file Kinematics.cxx.

References Copy(), and Init().

                                                    :
TObject()
{
  this->Init();
  this->Copy(kinematics);
}

Kinematics::Kinematics

(

TRootIOCtor *

)

Definition at line 56 of file Kinematics.cxx.

                                   :
TObject(),
fP4Fsl(0),
fP4HadSyst(0)
{

}

Kinematics::~Kinematics

(

)

Definition at line 64 of file Kinematics.cxx.

References CleanUp().

{
  this->CleanUp();
}

Member Function Documentation

void Kinematics::CleanUp ( void  )

private

clean-up

Definition at line 77 of file Kinematics.cxx.

References fKV, fP4Fsl, and fP4HadSyst.

Referenced by HadSystP4(), and ~Kinematics().

{
  fKV.clear();

  delete fP4Fsl;
  delete fP4HadSyst;
}

void Kinematics::ClearRunningValues

(

void

)

Definition at line 357 of file Kinematics.cxx.

References fKV, genie::kKVQ2, genie::kKVq2, genie::kKVt, genie::kKVW, genie::kKVx, and genie::kKVy.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), HadSystP4(), main(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::MECGenerator::SelectEmpiricalKinematics(), genie::RSPPResonanceSelector::SelectResonance(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode().

{
// clear the running values (leave the selected ones)
//
  fKV.erase( kKVx  );
  fKV.erase( kKVy  );
  fKV.erase( kKVQ2 );
  fKV.erase( kKVq2 );
  fKV.erase( kKVW  );
  fKV.erase( kKVt  );
}

void Kinematics::Copy

(

const Kinematics &

kine

)

Definition at line 93 of file Kinematics.cxx.

References fKV, fP4Fsl, fP4HadSyst, datagram_client::kv, Reset(), SetFSLeptonP4(), SetHadSystP4(), SetKV(), and febshutoff_auto::val.

Referenced by HadSystP4(), Kinematics(), operator=(), and genie::Interaction::SetKine().

{
  this->Reset();

  map<KineVar_t, double>::const_iterator iter;

  for(iter = kinematics.fKV.begin(); iter != kinematics.fKV.end(); ++iter) {
    KineVar_t kv  = iter->first;
    double    val = iter->second;
    this->SetKV(kv,val);
  }

  this->SetFSLeptonP4 (*kinematics.fP4Fsl);
  this->SetHadSystP4  (*kinematics.fP4HadSyst);
}

const TLorentzVector& genie::Kinematics::FSLeptonP4 ( void  ) const

inline

Definition at line 66 of file Kinematics.h.

References fP4Fsl.

Referenced by genie::BaryonResonanceDecayer::AcceptPionDecay(), genie::COHPrimaryLeptonGenerator::CalculatePrimaryLepton_AlvarezRuso(), genie::utils::EnergyDeltaFunctionSolutionQEL(), genie::LwlynSmithQELCCPXSec::FullDifferentialXSec(), genie::QELEventGenerator::ProcessEventRecord(), genie::AlvarezRusoCOHPiPXSec::XSec(), genie::BSKLNBaseRESPXSec2014::XSec(), and genie::NievesQELCCPXSec::XSec().

{ return *fP4Fsl;     }

double Kinematics::GetKV

(

KineVar_t

kv

)

const

Definition at line 333 of file Kinematics.cxx.

References genie::KineVar::AsString(), fKV, KVSet(), LOG, and pWARN.

Referenced by genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::SmithMonizQELCCPXSec::d2sQES_dQ2dv_SM(), genie::SmithMonizQELCCPXSec::d3sQES_dQ2dvdkF_SM(), HadSystP4(), Q2(), q2(), t(), genie::utils::mec::TensorContraction(), W(), x(), genie::AlamSimoAtharVacasSKPXSec2014::XSec(), genie::NievesSimoVacasMECPXSec2016::XSec(), and y().

{
  if(this->KVSet(kv)) {
     map<KineVar_t, double>::const_iterator iter = fKV.find(kv);
     return iter->second;
  } else {
    LOG("Interaction", pWARN)
        << "Kinematic variable: " << KineVar::AsString(kv) << " was not set";
  }
  return -99999;
}

const TLorentzVector& genie::Kinematics::HadSystP4 ( void  ) const

inline

Definition at line 67 of file Kinematics.h.

References CleanUp(), ClearRunningValues(), Copy(), E, fP4HadSyst, GetKV(), Init(), KVSet(), operator<<, operator=(), make_associated_cosmic_defs::p4, Print(), Reset(), SetFSLeptonP4(), SetHadSystP4(), SetKV(), SetQ2(), Setq2(), Sett(), SetW(), Setx(), Sety(), runNovaSAM::stream, and UseSelectedKinematics().

Referenced by genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::utils::EnergyDeltaFunctionSolutionQEL(), genie::LwlynSmithQELCCPXSec::FullDifferentialXSec(), genie::CharmHadronization::Hadronize(), genie::QELEventGenerator::ProcessEventRecord(), genie::AlvarezRusoCOHPiPXSec::XSec(), and genie::NievesQELCCPXSec::XSec().

{ return *fP4HadSyst; }

void Kinematics::Init ( void  )

private

initialize

Definition at line 69 of file Kinematics.cxx.

References fKV, fP4Fsl, and fP4HadSyst.

Referenced by HadSystP4(), and Kinematics().

{
  fKV.clear();

  fP4Fsl     = new TLorentzVector;
  fP4HadSyst = new TLorentzVector;
}

bool Kinematics::KVSet

(

KineVar_t

kv

)

const

Definition at line 327 of file Kinematics.cxx.

References fKV.

Referenced by GetKV(), HadSystP4(), genie::utils::ghep::NeutReactionCode(), Q2(), q2(), genie::QPMDISStrucFuncBase::Q2(), genie::utils::kinematics::Q2(), SetKV(), t(), genie::utils::kinematics::UpdateWQ2FromXY(), genie::utils::kinematics::UpdateXFromQ2Y(), genie::utils::kinematics::UpdateXYFromWQ2(), W(), genie::utils::kinematics::W(), x(), and y().

{
  if(fKV.count(kv) == 1) return true;
  else return false;
}

double Kinematics::Log10Q2

(

bool

selected = false

)

const

Definition at line 229 of file Kinematics.cxx.

References Q2().

{
  double Q2s = this->Q2(selected);
  return (Q2s>0) ? TMath::Log10(Q2s) : -99999;
}

double Kinematics::Log10W

(

bool

selected = false

)

const

Definition at line 235 of file Kinematics.cxx.

References W().

{
  double Ws = this->W(selected);
  return (Ws>0) ? TMath::Log10(Ws) : -99999;
}

double Kinematics::Log10x

(

bool

selected = false

)

const

Definition at line 217 of file Kinematics.cxx.

References x().

{
  double xs = this->x(selected);
  return (xs>0) ? TMath::Log10(xs) : -99999;
}

double Kinematics::Log10y

(

bool

selected = false

)

const

Definition at line 223 of file Kinematics.cxx.

References y().

{
  double ys = this->y(selected);
  return (ys>0) ? TMath::Log10(ys) : -99999;
}

double Kinematics::LogQ2

(

bool

selected = false

)

const

Definition at line 205 of file Kinematics.cxx.

References Q2().

{
  double Q2s = this->Q2(selected);
  return (Q2s>0) ? TMath::Log(Q2s) : -99999;
}

double Kinematics::LogW

(

bool

selected = false

)

const

Definition at line 211 of file Kinematics.cxx.

References W().

{
  double Ws = this->W(selected);
  return (Ws>0) ? TMath::Log(Ws) : -99999;
}

double Kinematics::Logx

(

bool

selected = false

)

const

Definition at line 193 of file Kinematics.cxx.

References x().

{
  double xs = this->x(selected);
  return (xs>0) ? TMath::Log(xs) : -99999;
}

double Kinematics::Logy

(

bool

selected = false

)

const

Definition at line 199 of file Kinematics.cxx.

References y().

{
  double ys = this->y(selected);
  return (ys>0) ? TMath::Log(ys) : -99999;
}

Kinematics & Kinematics::operator=

(

const Kinematics &

kine

)

Definition at line 401 of file Kinematics.cxx.

References Copy().

Referenced by HadSystP4().

{
  this->Copy(kinematics);
  return (*this);
}

void Kinematics::Print

(

ostream &

stream

)

const

Definition at line 388 of file Kinematics.cxx.

References genie::KineVar::AsString(), allTimeWatchdog::endl, fKV, datagram_client::kv, and febshutoff_auto::val.

Referenced by ClassImp(), and HadSystP4().

{
  stream << "[-] [Kinematics]" << endl;

  map<KineVar_t, double>::const_iterator iter;

  for(iter = fKV.begin(); iter != fKV.end(); ++iter) {
    KineVar_t kv  = iter->first;
    double    val = iter->second;
    stream << " |--> " << KineVar::AsString(kv) << " = " << val << endl;
  }
}

double Kinematics::Q2

(

bool

selected = false

)

const

Definition at line 135 of file Kinematics.cxx.

References GetKV(), genie::kKVQ2, genie::kKVq2, genie::kKVSelQ2, genie::kKVSelq2, KVSet(), LOG, and pWARN.

Referenced by genie::MECGenerator::AddFinalStateLepton(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), ConvertToGST(), genie::KovalenkoQELCharmPXSec::DR(), genie::BergerSehgalFMCOHPiPXSec2015::ExactKinematicTerm(), genie::BergerSehgalCOHPiPXSec2015::ExactKinematicTerm(), GenerateEvent(), genie::TransverseEnhancementFFModel::GetTransEnhMagFF(), genie::KPhaseSpace::IsAllowed(), genie::utils::kinematics::Jacobian(), Log10Q2(), LogQ2(), genie::utils::nuclear::NuclQELXSecSuppression(), genie::BergerSehgalFMCOHPiPXSec2015::PionCOMAbsMomentum(), genie::BergerSehgalCOHPiPXSec2015::PionCOMAbsMomentum(), genie::PrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QPMDISStrucFuncBase::Q2(), genie::utils::kinematics::Q2(), genie::KPhaseSpace::TLim(), genie::utils::kinematics::UpdateXFromQ2Y(), genie::utils::kinematics::UpdateXYFromWQ2(), genie::EmpiricalMECPXSec2015::XSec(), genie::COHElasticPXSec::XSec(), genie::AhrensNCELPXSec::XSec(), genie::RosenbluthPXSec::XSec(), genie::P33PaschosLalakulichPXSec::XSec(), genie::AhrensDMELPXSec::XSec(), genie::BergerSehgalFMCOHPiPXSec2015::XSec(), genie::BergerSehgalCOHPiPXSec2015::XSec(), genie::KovalenkoQELCharmPXSec::XSec(), and genie::KPhaseSpace::YLim().

{
// returns the running or selected value of momentum transfer Q2 (>0)

  if(selected) {
    if      (this->KVSet(kKVSelQ2) ) { return     this->GetKV(kKVSelQ2); }
    else if (this->KVSet(kKVSelq2) ) { return -1* this->GetKV(kKVSelq2); }
  } else {
    if      (this->KVSet(kKVQ2) )    { return     this->GetKV(kKVQ2); }
    else if (this->KVSet(kKVq2) )    { return -1* this->GetKV(kKVq2); }
  }

  LOG("Interaction", pWARN) << "Kinematic variable Q2 was not set";
  return -99999;
}

double Kinematics::q2

(

bool

selected = false

)

const

Definition at line 151 of file Kinematics.cxx.

References GetKV(), genie::kKVQ2, genie::kKVq2, genie::kKVSelQ2, genie::kKVSelq2, KVSet(), LOG, and pWARN.

Referenced by CalculateFormFactor(), genie::SmithMonizQELCCPXSec::dsQES_dQ2_SM(), genie::DipoleAxialFormFactorModel::FA(), genie::KuzminNaumov2016AxialFormFactorModel::FA(), genie::ZExpAxialFormFactorModel::FA(), genie::LwlynSmithFFNC::Fp(), genie::LwlynSmithFF::Fp(), genie::BBA03ELFormFactorsModel::Gen(), genie::BBA07ELFormFactorsModel::Gen(), genie::DipoleELFormFactorsModel::Gep(), genie::BBA03ELFormFactorsModel::Gep(), genie::BBA07ELFormFactorsModel::Gep(), genie::DipoleELFormFactorsModel::Gmn(), genie::BBA03ELFormFactorsModel::Gmn(), genie::BBA07ELFormFactorsModel::Gmn(), genie::DipoleELFormFactorsModel::Gmp(), genie::BBA03ELFormFactorsModel::Gmp(), genie::BBA07ELFormFactorsModel::Gmp(), genie::utils::nuclear::NuclQELXSecSuppression(), genie::BBA05ELFormFactorsModel::tau(), genie::LwlynSmithFF::tau(), genie::StrumiaVissaniIBDPXSec::XSec(), genie::LwlynSmithQELCCPXSec::XSec(), genie::PaisQELLambdaPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

{
// returns the running or selected value of momentum transfer q2 (<0)

  if(selected) {
    if      (this->KVSet(kKVSelQ2) ) { return -1* this->GetKV(kKVSelQ2); }
    else if (this->KVSet(kKVSelq2) ) { return     this->GetKV(kKVSelq2); }
  } else {
    if      (this->KVSet(kKVQ2) )    { return -1* this->GetKV(kKVQ2); }
    else if (this->KVSet(kKVq2) )    { return     this->GetKV(kKVq2); }
  }

  LOG("Interaction", pWARN) << "Kinematic variable q2 was not set";
  return -99999;
}

void Kinematics::Reset

(

void

)

Definition at line 85 of file Kinematics.cxx.

References fKV, SetFSLeptonP4(), and SetHadSystP4().

Referenced by Copy(), and HadSystP4().

{
  fKV.clear();

  this->SetFSLeptonP4 (0,0,0,0);
  this->SetHadSystP4  (0,0,0,0);
}

void Kinematics::SetFSLeptonP4

(

const TLorentzVector &

p4

)

Definition at line 307 of file Kinematics.cxx.

References fP4Fsl.

Referenced by genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::OutgoingDarkGenerator::AddToEventRecord(), genie::utils::ComputeFullQELPXSec(), Copy(), genie::utils::gsl::d5XSecAR::DoEval(), genie::utils::gsl::d5Xsec_dEldOmegaldOmegapi::DoEval(), genie::utils::gsl::d4Xsec_dEldThetaldOmegapi::DoEval(), genie::utils::gsl::d3Xsec_dOmegaldThetapi::DoEval(), HadSystP4(), Reset(), cafrwgt::CAFReweight::RetrieveGHEP(), genie::MECGenerator::SelectNSVLeptonKinematics(), and genie::COHKinematicsGenerator::SetKinematics().

{
  fP4Fsl->SetPxPyPzE(p4.Px(), p4.Py(), p4.Pz(), p4.E());
}

void Kinematics::SetFSLeptonP4	(	double	px,
		double	py,
		double	pz,
		double	E
	)

Definition at line 312 of file Kinematics.cxx.

References fP4Fsl.

{
  fP4Fsl->SetPxPyPzE(px,py,pz,E);
}

void Kinematics::SetHadSystP4

(

const TLorentzVector &

p4

)

Definition at line 317 of file Kinematics.cxx.

References fP4HadSyst.

Referenced by genie::HadronicSystemGenerator::AddFinalHadronicSyst(), genie::utils::ComputeFullQELPXSec(), Copy(), genie::utils::gsl::d5XSecAR::DoEval(), genie::utils::gsl::d5Xsec_dEldOmegaldOmegapi::DoEval(), genie::utils::gsl::d4Xsec_dEldThetaldOmegapi::DoEval(), genie::utils::gsl::d3Xsec_dOmegaldThetapi::DoEval(), genie::MECGenerator::GenerateNSVInitialHadrons(), HadSystP4(), Reset(), cafrwgt::CAFReweight::RetrieveGHEP(), and genie::COHKinematicsGenerator::SetKinematics().

{
  fP4HadSyst->SetPxPyPzE(p4.Px(), p4.Py(), p4.Pz(), p4.E());
}

void Kinematics::SetHadSystP4	(	double	px,
		double	py,
		double	pz,
		double	E
	)

Definition at line 322 of file Kinematics.cxx.

References fP4HadSyst.

{
  fP4HadSyst->SetPxPyPzE(px,py,pz,E);
}

void Kinematics::SetKV	(	KineVar_t	kv,
		double	value
	)

Definition at line 345 of file Kinematics.cxx.

References genie::KineVar::AsString(), fKV, datagram_client::kv, KVSet(), LOG, and pDEBUG.

Referenced by genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::SKKinematicsGenerator::ComputeMaxXSec(), genie::QELEventGeneratorSM::ComputeMaxXSec(), genie::QELEventGeneratorSM::ComputeMaxXSec2(), Copy(), genie::SmithMonizQELCCPXSec::d2sQES_dQ2dv_SM(), genie::utils::gsl::d2Xsec_dTCosth::DoEval(), genie::utils::gsl::d2Xsec_dQ2dv::DoEval(), HadSystP4(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::MECGenerator::SelectNSVLeptonKinematics(), SetQ2(), Setq2(), Sett(), SetW(), Setx(), and Sety().

{
  LOG("Interaction", pDEBUG)
            << "Setting " << KineVar::AsString(kv) << " to " << value;

  if(this->KVSet(kv)) {
     fKV[kv] = value;
  } else {
     fKV.insert( map<KineVar_t, double>::value_type(kv,value) );
  }
}

void Kinematics::SetQ2	(	double	Q2,
		bool	selected = false
	)

Definition at line 265 of file Kinematics.cxx.

References genie::kKVQ2, genie::kKVSelQ2, LOG, pWARN, and SetKV().

Referenced by BuildStdNtuple(), CalculateFormFactor(), genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::NievesQELCCPXSec::CompareNievesTensors(), genie::utils::ComputeFullQELPXSec(), genie::RESKinematicsGenerator::ComputeMaxXSec(), genie::IBDKinematicsGenerator::ComputeMaxXSec(), genie::QELKinematicsGenerator::ComputeMaxXSec(), genie::DMELKinematicsGenerator::ComputeMaxXSec(), genie::utils::gsl::dXSec_dQ2_E::DoEval(), genie::utils::gsl::d2XSecRESFast_dWQ2_E::DoEval(), genie::utils::gsl::d2XSec_dQ2dy_E::DoEval(), genie::utils::gsl::d2XSec_dQ2dydt_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_EW::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_EQ2::DoEval(), genie::LwlynSmithQELCCPXSec::FullDifferentialXSec(), GetCrossSection(), HadSystP4(), main(), MakePlots(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgal(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgalFM(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), cafrwgt::CAFReweight::RetrieveGHEP(), genie::MECGenerator::SelectEmpiricalKinematics(), genie::MECGenerator::SelectNSVLeptonKinematics(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::utils::kinematics::UpdateWQ2FromXY(), UseSelectedKinematics(), VerticalSlice(), and genie::NievesQELCCPXSec::XSec().

{
// sets the running or selected value of momentum transfer Q2 (>0)

  if(Qsqrd<0) {
     LOG("Interaction", pWARN)
                 << "Setting unphysical value for Q2 (Q2 = " << Qsqrd << ")";
  }
  KineVar_t kvar = (selected) ? kKVSelQ2 : kKVQ2;
  this->SetKV(kvar, Qsqrd);
}

void Kinematics::Setq2	(	double	q2,
		bool	selected = false
	)

Definition at line 277 of file Kinematics.cxx.

References genie::kKVq2, genie::kKVSelq2, LOG, pWARN, and SetKV().

Referenced by HadSystP4(), cafrwgt::CAFReweight::RetrieveGHEP(), and UseSelectedKinematics().

{
// sets the running or selected value of momentum transfer q2 (<0)

  if(qsqrd>0) {
     LOG("Interaction", pWARN)
                 << "Setting unphysical value for q2 (q2 = " << qsqrd << ")";
  }
  KineVar_t kvar = (selected) ? kKVSelq2 : kKVq2;
  this->SetKV(kvar, qsqrd);
}

void Kinematics::Sett	(	double	t,
		bool	selected = false
	)

Definition at line 301 of file Kinematics.cxx.

References genie::kKVSelt, genie::kKVt, and SetKV().

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), genie::utils::gsl::d2XSec_dQ2dydt_E::DoEval(), genie::utils::gsl::d3XSec_dxdydt_E::DoEval(), HadSystP4(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgalFM(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), cafrwgt::CAFReweight::RetrieveGHEP(), and UseSelectedKinematics().

{
  KineVar_t kvar = (selected) ? kKVSelt : kKVt;
  this->SetKV(kvar, tval);
}

void Kinematics::SetW	(	double	W,
		bool	selected = false
	)

Definition at line 289 of file Kinematics.cxx.

References genie::kKVSelW, genie::kKVW, LOG, pWARN, and SetKV().

Referenced by genie::DISHadronicSystemGenerator::AddFragmentationProducts(), genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::RESKinematicsGenerator::ComputeMaxXSec(), genie::QPMDISPXSec::DISRESJoinSuppressionFactor(), genie::QPMDMDISPXSec::DMDISRESJoinSuppressionFactor(), genie::utils::gsl::d2XSecRESFast_dWQ2_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_EW::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_EQ2::DoEval(), GetCrossSection(), HadSystP4(), main(), genie::utils::kinematics::PhaseSpaceVolume(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), cafrwgt::CAFReweight::RetrieveGHEP(), genie::MECGenerator::SelectEmpiricalKinematics(), genie::MECGenerator::SelectNSVLeptonKinematics(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::utils::kinematics::UpdateWQ2FromXY(), and UseSelectedKinematics().

{
// sets the running or selected value of invariant hadronic mass W

  if(hadr_mass_W<0) {
     LOG("Interaction", pWARN)
                   << "Setting unphysical value for W (W = " << hadr_mass_W << ")";
  }
  KineVar_t kvar = (selected) ? kKVSelW : kKVW;
  this->SetKV(kvar, hadr_mass_W);
}

void Kinematics::Setx	(	double	x,
		bool	selected = false
	)

Definition at line 241 of file Kinematics.cxx.

References genie::kKVSelx, genie::kKVx, LOG, pWARN, and SetKV().

Referenced by BuildStdNtuple(), genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), genie::DISKinematicsGenerator::ComputeMaxXSec(), genie::DMDISKinematicsGenerator::ComputeMaxXSec(), genie::utils::gsl::d2XSec_dxdy_E::DoEval(), genie::utils::gsl::d3XSec_dxdydt_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_E::DoEval(), genie::utils::gsl::d2XSec_dxdy_Ex::DoEval(), genie::utils::gsl::d2XSec_dxdy_Ey::DoEval(), genie::utils::gsl::d5XSecAR::DoEval(), genie::utils::gsl::d5Xsec_dEldOmegaldOmegapi::DoEval(), genie::utils::gsl::d4Xsec_dEldThetaldOmegapi::DoEval(), genie::utils::gsl::d3Xsec_dOmegaldThetapi::DoEval(), GetCrossSection(), HadSystP4(), main(), MakePlots(), genie::COHKinematicsGenerator::MaxXSec_ReinSehgal(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), cafrwgt::CAFReweight::RetrieveGHEP(), genie::MECGenerator::SelectEmpiricalKinematics(), genie::MECGenerator::SelectNSVLeptonKinematics(), genie::COHKinematicsGenerator::SetKinematics(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::utils::kinematics::UpdateXFromQ2Y(), genie::utils::kinematics::UpdateXYFromWQ2(), UseSelectedKinematics(), and VerticalSlice().

{
// sets the running or selected value of Bjorken scaling variable x

  if(xbj<0 || xbj>1) {
     LOG("Interaction", pWARN)
                      << "Setting unphysical value for x (x = " << xbj << ")";
  }
  KineVar_t kvar = (selected) ? kKVSelx : kKVx;
  this->SetKV(kvar, xbj);
}

void Kinematics::Sety	(	double	y,
		bool	selected = false
	)

Definition at line 253 of file Kinematics.cxx.

References genie::kKVSely, genie::kKVy, LOG, pWARN, and SetKV().

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), genie::NuEKinematicsGenerator::ComputeMaxXSec(), genie::COHElKinematicsGenerator::ComputeMaxXSec(), genie::DISKinematicsGenerator::ComputeMaxXSec(), genie::DMDISKinematicsGenerator::ComputeMaxXSec(), genie::utils::gsl::dXSec_dy_E::DoEval(), genie::utils::gsl::d2XSec_dxdy_E::DoEval(), genie::utils::gsl::d2XSec_dQ2dy_E::DoEval(), genie::utils::gsl::d2XSec_dQ2dydt_E::DoEval(), genie::utils::gsl::d3XSec_dxdydt_E::DoEval(), genie::utils::gsl::d2XSec_dWdQ2_E::DoEval(), genie::utils::gsl::d2XSec_dxdy_Ex::DoEval(), genie::utils::gsl::d2XSec_dxdy_Ey::DoEval(), genie::utils::gsl::d5XSecAR::DoEval(), genie::utils::gsl::d5Xsec_dEldOmegaldOmegapi::DoEval(), genie::utils::gsl::d4Xsec_dEldThetaldOmegapi::DoEval(), genie::utils::gsl::d3Xsec_dOmegaldThetapi::DoEval(), GetCrossSection(), HadSystP4(), main(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgal(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgalFM(), genie::COHKinematicsGenerator::MaxXSec_ReinSehgal(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::COHElKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), cafrwgt::CAFReweight::RetrieveGHEP(), genie::MECGenerator::SelectEmpiricalKinematics(), genie::MECGenerator::SelectNSVLeptonKinematics(), genie::COHKinematicsGenerator::SetKinematics(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::utils::kinematics::UpdateXYFromWQ2(), and UseSelectedKinematics().

{
// sets the running or selected value of inelasticity y

  if(inel_y<0 || inel_y>1) {
     LOG("Interaction", pWARN)
                     << "Setting unphysical value for y (y = " << inel_y << ")";
  }
  KineVar_t kvar = (selected) ? kKVSely : kKVy;
  this->SetKV(kvar, inel_y);
}

double Kinematics::t

(

bool

selected = false

)

const

Definition at line 180 of file Kinematics.cxx.

References GetKV(), genie::kKVSelt, genie::kKVt, KVSet(), LOG, and pWARN.

Referenced by genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), ConvertToGST(), genie::KPhaseSpace::IsAllowed(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::ReinDFRPXSec::XSec(), and genie::BergerSehgalFMCOHPiPXSec2015::XSec().

{
// returns the running or selected value of invariant hadronic mass W

  KineVar_t kvar = (selected) ? kKVSelt : kKVt;

  if(this->KVSet(kvar)) { return this->GetKV(kvar); }
  else {
    LOG("Interaction", pWARN) << "Kinematic variable t was not set";
  }
  return -99999;
}

void Kinematics::UseSelectedKinematics

(

void

)

Definition at line 369 of file Kinematics.cxx.

References fKV, genie::kKVSelQ2, genie::kKVSelq2, genie::kKVSelt, genie::kKVSelW, genie::kKVSelx, genie::kKVSely, SetQ2(), Setq2(), Sett(), SetW(), Setx(), and Sety().

Referenced by HadSystP4(), main(), and genie::RSPPResonanceSelector::SelectResonance().

{
// copy the selected kinematics into the running ones
//
  map<KineVar_t, double>::const_iterator iter;
  iter = fKV.find(kKVSelx);
  if(iter != fKV.end()) this->Setx(iter->second);
  iter = fKV.find(kKVSely);
  if(iter != fKV.end()) this->Sety(iter->second);
  iter = fKV.find(kKVSelQ2);
  if(iter != fKV.end()) this->SetQ2(iter->second);
  iter = fKV.find(kKVSelq2);
  if(iter != fKV.end()) this->Setq2(iter->second);
  iter = fKV.find(kKVSelW);
  if(iter != fKV.end()) this->SetW(iter->second);
  iter = fKV.find(kKVSelt);
  if(iter != fKV.end()) this->Sett(iter->second);
}

double Kinematics::W

(

bool

selected = false

)

const

Definition at line 167 of file Kinematics.cxx.

References GetKV(), genie::kKVSelW, genie::kKVW, KVSet(), LOG, and pWARN.

Referenced by genie::DFRKinematicsGenerator::ComputeMaxXSec(), ConvertToGHad(), ConvertToGST(), GenerateEvent(), genie::KNOPythiaHadronization::Hadronize(), genie::PythiaHadronization::Hadronize(), genie::CharmHadronization::Hadronize(), genie::KPhaseSpace::IsAllowed(), genie::utils::kinematics::Jacobian(), Log10W(), LogW(), genie::HadronizationModelBase::MaxMult(), genie::PythiaHadronization::MultiplicityProb(), genie::utils::ghep::NeutReactionCode(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::KNOPythiaHadronization::SelectHadronizer(), genie::utils::kinematics::UpdateXYFromWQ2(), genie::utils::kinematics::W(), genie::EmpiricalMECPXSec2015::XSec(), genie::P33PaschosLalakulichPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

{
// returns the running or selected value of invariant hadronic mass W

  KineVar_t kvar = (selected) ? kKVSelW : kKVW;

  if(this->KVSet(kvar)) { return this->GetKV(kvar); }
  else {
    LOG("Interaction", pWARN) << "Kinematic variable W was not set";
  }
  return -99999;
}

double Kinematics::x

(

bool

selected = false

)

const

Definition at line 109 of file Kinematics.cxx.

References GetKV(), genie::kKVSelx, genie::kKVx, KVSet(), LOG, and pWARN.

Referenced by genie::QPMDISStrucFuncBase::Calculate(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), ConvertToGHad(), ConvertToGST(), genie::QPMDISPXSec::DISRESJoinSuppressionFactor(), genie::QPMDMDISPXSec::DMDISRESJoinSuppressionFactor(), GenerateEvent(), genie::KPhaseSpace::IsAllowed(), genie::utils::kinematics::Jacobian(), Log10x(), Logx(), genie::QPMDISStrucFuncBase::NuclMod(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::QPMDISStrucFuncBase::Q2(), genie::utils::kinematics::Q2(), genie::QPMDISStrucFuncBase::R(), genie::BYStrucFunc::ScalingVar(), genie::QPMDISStrucFuncBase::ScalingVar(), genie::utils::kinematics::UpdateWQ2FromXY(), genie::utils::kinematics::W(), genie::SlowRsclCharmDISPXSecLO::XSec(), genie::ReinDFRPXSec::XSec(), genie::AivazisCharmPXSecLO::XSec(), genie::QPMDISPXSec::XSec(), genie::QPMDMDISPXSec::XSec(), genie::BergerSehgalCOHPiPXSec2015::XSec(), genie::ReinSehgalCOHPiPXSec::XSec(), and genie::KPhaseSpace::YLim_X().

{
// returns the running or selected value of Bjorken scaling variable x

  KineVar_t kvar = (selected) ? kKVSelx : kKVx;

  if(this->KVSet(kvar)) { return this->GetKV(kvar); }
  else {
    LOG("Interaction", pWARN) << "Kinematic variable x was not set";
  }
  return -99999;
}

double Kinematics::y

(

bool

selected = false

)

const

Definition at line 122 of file Kinematics.cxx.

References GetKV(), genie::kKVSely, genie::kKVy, KVSet(), LOG, and pWARN.

Referenced by genie::MECGenerator::AddFinalStateLepton(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), ConvertToGHad(), ConvertToGST(), genie::QPMDISPXSec::DISRESJoinSuppressionFactor(), genie::QPMDMDISPXSec::DMDISRESJoinSuppressionFactor(), genie::BergerSehgalFMCOHPiPXSec2015::ExactKinematicTerm(), genie::BergerSehgalCOHPiPXSec2015::ExactKinematicTerm(), GenerateEvent(), genie::KPhaseSpace::IsAllowed(), genie::utils::kinematics::Jacobian(), Log10y(), Logy(), genie::BergerSehgalFMCOHPiPXSec2015::PionCOMAbsMomentum(), genie::BergerSehgalCOHPiPXSec2015::PionCOMAbsMomentum(), genie::NuEPrimaryLeptonGenerator::ProcessEventRecord(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::PrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::ProcessEventRecord(), genie::QPMDISStrucFuncBase::Q2(), genie::utils::kinematics::Q2(), genie::KPhaseSpace::TLim(), genie::utils::kinematics::UpdateWQ2FromXY(), genie::utils::kinematics::UpdateXFromQ2Y(), genie::utils::kinematics::W(), genie::SlowRsclCharmDISPXSecLO::XSec(), genie::ReinDFRPXSec::XSec(), genie::AivazisCharmPXSecLO::XSec(), genie::IMDAnnihilationPXSec::XSec(), genie::QPMDISPXSec::XSec(), genie::BergerSehgalFMCOHPiPXSec2015::XSec(), genie::QPMDMDISPXSec::XSec(), genie::BergerSehgalCOHPiPXSec2015::XSec(), genie::NuElectronPXSec::XSec(), genie::ReinSehgalCOHPiPXSec::XSec(), and genie::BardinIMDRadCorPXSec::XSec().

{
// returns the running or selected value of inelasticity y

  KineVar_t kvar = (selected) ? kKVSely : kKVy;

  if(this->KVSet(kvar)) { return this->GetKV(kvar); }
  else {
    LOG("Interaction", pWARN) << "Kinematic variable y was not set";
  }
  return -99999;
}

Friends And Related Function Documentation

ostream& operator<< ( ostream &  stream, const Kinematics &  kine  )

friend

Referenced by HadSystP4().

Member Data Documentation

map<KineVar_t, double> genie::Kinematics::fKV

private

selected kinematics

Definition at line 103 of file Kinematics.h.

Referenced by CleanUp(), ClearRunningValues(), Copy(), GetKV(), Init(), KVSet(), Print(), Reset(), SetKV(), and UseSelectedKinematics().

TLorentzVector* genie::Kinematics::fP4Fsl

private

generated final state primary lepton 4-p (LAB)

Definition at line 104 of file Kinematics.h.

Referenced by CleanUp(), Copy(), FSLeptonP4(), Init(), and SetFSLeptonP4().

TLorentzVector* genie::Kinematics::fP4HadSyst

private

generated final state hadronic system 4-p (LAB)

Definition at line 105 of file Kinematics.h.

Referenced by CleanUp(), Copy(), HadSystP4(), Init(), and SetHadSystP4().

---

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

• Framework/Interaction/Kinematics.h
• Framework/Interaction/Kinematics.cxx