C++ Interface for the Fortran Bhabha and exclusive two photon generator BABAYAGA.NLO. More...

#include <BabayagaNLO.h>

Public Member Functions
	BabayagaNLO ()
	Constructor.

	~BabayagaNLO ()
	Destrucotr.

void	setDefaultSettings ()
	Sets the default settings for the BhWide Fortran generator.

void	setScatAngle (std::pair< double, double > angleRange)
	Sets the theta scattering angle range for the scattered particles.

void	setNSearchMax (int nSearchMax)
	Sets the number of events used to search maximum.

void	setFMax (double fMax)
	Maximum differential cross section.

void	setMaxAcollinearity (double maxAcollinearity)
	Sets maximum acollinearity angle between finale state leptons/photons in degrees.

void	setCmsEnergyNominal (double cmsEnergyNominal)
	Sets nominal ECM.

void	setMinEnergy (double minEnergy)
	Sets minimum energy for leptons/photons in the final state, in GeV.

void	setEpsilon (double epsilon)
	Sets soft/hard photon energy separator.

void	setVacPol (const std::string &vacPol)
	Sets vacuum polarization.

void	setModel (const std::string &model)
	Sets model: matched or ps.

void	setMode (const std::string &mode)
	Sets mode: weighted or unweighted.

void	setOrder (const std::string &order)
	Sets Order: born, alpha or exp.

void	setFinalState (const std::string &finalState)
	Sets final state minimum energy.

void	setEnergySpread (double spread)
	TEMPORARY SOLUTION! Approximate energy spread per beam (CMS)

void	setVPUncertainty (bool vpuncertainty)
	Calculate VP uncertainty by internal reweighting (on/off)

void	setNPhotons (int nPhot)
	Sets the fixed number of nphot (hard) photons.

void	setUserMode (const std::string &usermode)
	Sets User mode similar to TEEGG: ETRON, EGAMMA, GAMMA or PRESCALE or NONE (default)

void	setNSKDataFile (const std::string &NSKDataFile)
	Sets NSK VP data file.

void	setEEMIN (double eemin)
	Sets the minimum CMS energy of the tagged e-/e+.

void	setTEMIN (double temin)
	Sets the minimum CMS angle between the tagged e-/e+ and -z axis.

void	setEGMIN (double egmin)
	Sets the minimum CMS energy of the gamma.

void	setTGMIN (double tgmin)
	Sets the minimum CMS angle between the gamma and -z axis.

void	setEEVETO (double eeveto)
	Sets the minimum CMS energy to veto e-/e+.

void	setTEVETO (double teveto)
	Sets the maximum CMS theta of e-/e+ in final state.

void	setEGVETO (double egveto)
	Sets the minimum CMS energy to veto gamma.

void	setTGVETO (double tgveto)
	Sets the maximum CMS angle between the gamma and -z axis.

void	setMaxPrescale (double maxprescale)
	Sets the maximum prescale value.

void	init ()
	Initializes the generator.

void	initExtraInfo ()
	Initializes the extra info.

void	generateEvent (MCParticleGraph &mcGraph, double ecm, ROOT::Math::XYZVector vertex, ROOT::Math::LorentzRotation boost)
	Generates one single event.

double	getSDif ()
	Returns weight.

void	term ()
	Terminates the generator.

Protected Member Functions
void	applySettings ()
	Apply the settings to the internal Fortran generator.

void	storeParticle (MCParticleGraph &mcGraph, const double *mom, int pdg, ROOT::Math::XYZVector vertex, ROOT::Math::LorentzRotation boost, bool isVirtual=false, bool isInitial=false, bool isISRFSR=false)
	Store a single generated particle into the MonteCarlo graph.

Protected Attributes
bool	m_applyBoost
	Apply a boost to the MCParticles.

double	m_pi
	pi=3.1415....

double	m_conversionFactor
	Conversion factor for hbarc to nb.

double	m_alphaQED0
	QED coupling constant at Q=0.

double	m_massElectron
	muon mass.

double	m_massMuon
	electron mass.

double	m_massW
	W mass [GeV] for on shell sin2theta and GF.

double	m_massZ
	Z mass [GeV].

double	m_widthZ
	Z width [GeV] (may be recalculated by EW library).

std::string	m_vacPol
	vacuum polarization: off, hadr5 (Jegerlehner) or hmnt (Teubner).

std::string	m_finalState
	final state: ee, gg or mm.

std::string	m_model
	model: matched or ps.

std::string	m_mode
	mode: weighted or unweighted.

std::string	m_order
	order: born, alpha or exp.

double	m_EnergySpread
	TEMPORARY SOLUTION! Approximate energy spread per beam (CMS).

double	m_cmsEnergyNominal
	Nominal CMS Energy = 2*Ebeam [GeV].

std::pair< double, double >	m_ScatteringAngleRange
	Min and Max val.

std::pair< double, double >	m_ScatteringAngleRangePhoton
	Min and Max val.

double	m_epsilon
	Soft/hard photon separator in units of CMS/2.

double	m_eMin
	[GeV].

double	m_maxAcollinearity
	Maximum acollinearity of the electron-positron pair.

int	m_nSearchMax
	Events used to search maximum of differential cross section.

int	m_nPhot
	fixed number of nphot (hard) photons are generated.

bool	m_VPUncertainty
	vary all VP related parameters and extracted total uncertainty.

std::string	m_NSKDataFile
	data file for NSK VP.

std::string	m_userMode
	User mode similar to TEEGG: ETRON, EGAMMA, GAMMA or PRESCALE.

double	m_eemin
	Minimum CMS energy of the tagged e-/e+ (GeV).

double	m_temin
	Minimum CMS angle between the tagged e-/e+ and -z axis (deg).

double	m_egmin
	Minimum CMS energy of the gamma (GeV).

double	m_tgmin
	Minimum CMS angle between the gamma and -z axis (deg).

double	m_eeveto
	Minimum CMS energy to veto e-/e+ (GeV).

double	m_teveto
	Maximum CMS theta of e-/e+ in final state (deg).

double	m_egveto
	Minimum CMS energy to veto gamma (GeV).

double	m_tgveto
	Maximum CMS angle between the gamma and -z axis (deg).

double	m_maxprescale
	Maximum prescale value.

double	m_fMax
	Maximum of differential cross section.

double	m_sDif
	Differential xsec/weight used for event.

ROOT::Math::LorentzRotation	m_boostVector
	The Lorentz boost vector for the transformation CMS to LAB frame.

Private Attributes
int	m_npar [100]
	Integer parameters for BabayagaNLO.

double	m_xpar [100]
	Double parameters for BabayagaNLO.

Detailed Description

C++ Interface for the Fortran Bhabha and exclusive two photon generator BABAYAGA.NLO.

Definition at line 27 of file BabayagaNLO.h.

Constructor & Destructor Documentation

◆ BabayagaNLO()

BabayagaNLO ( )

Constructor.

Sets the default settings.

Definition at line 257 of file BabayagaNLO.cc.

{
  for (int i = 0; i < 100; ++i) {
    m_npar[i] = 0;
    m_xpar[i] = 0.0;
  }
 
  m_sDif = 0.;
 
  setDefaultSettings();
}

◆ ~BabayagaNLO()

~BabayagaNLO ( )

Destrucotr.

Definition at line 269 of file BabayagaNLO.cc.

{
 
}

Member Function Documentation

◆ applySettings()

void applySettings ( )

protected

Apply the settings to the internal Fortran generator.

Definition at line 403 of file BabayagaNLO.cc.

{
  //--------------------
  // Integer parameters
  //--------------------
  m_npar[0] = m_nPhot;
  m_npar[1] = m_nSearchMax;
  if (m_VPUncertainty) m_npar[5] = 1;
  else m_npar[5] = 0;
 
  //--------------------
  // Double parameters
  //--------------------
  m_xpar[0]  = m_cmsEnergyNominal;
  m_xpar[1]  = m_pi;
  m_xpar[2]  = m_conversionFactor;
  m_xpar[3]  = m_alphaQED0;
  m_xpar[4]  = m_massElectron;
  m_xpar[5]  = m_massMuon;
  m_xpar[6]  = m_massW;
  m_xpar[7]  = m_massZ;
  m_xpar[8]  = m_widthZ;
  m_xpar[9]  = m_eMin;
  m_xpar[10] = m_maxAcollinearity;
  m_xpar[11] = m_epsilon;
  m_xpar[20] = m_ScatteringAngleRange.first;
  m_xpar[21] = m_ScatteringAngleRange.second;
  m_xpar[30] = m_fMax;
  m_xpar[40] = m_EnergySpread; // only for some internal checks, not the actual smaering!
 
  //user cuts
  m_xpar[50] = m_eemin;
  m_xpar[51] = m_temin;
  m_xpar[52] = m_egmin;
  m_xpar[53] = m_tgmin;
  m_xpar[54] = m_eeveto;
  m_xpar[55] = m_teveto;
  m_xpar[56] = m_egveto;
  m_xpar[57] = m_tgveto;
  m_xpar[58] = m_maxprescale;
 
  //--------------------
  // string parameters, wrapped to integers
  //--------------------
  if (m_finalState == "ee") m_npar[20] = 1;
  else if (m_finalState == "gg") m_npar[20] = 2;
  else if (m_finalState == "mm") m_npar[20] = 3;
  else B2FATAL("Invalid final state: " << m_finalState);
 
  if (m_model == "matched") m_npar[21] = 1;
  else if (m_model == "ps") m_npar[21] = 2;
  else B2FATAL("Invalid matching model: " << m_model);
 
  if (m_order == "born") m_npar[22] = 1;
  else if (m_order == "alpha") m_npar[22] = 2;
  else if (m_order == "exp") m_npar[22] = 3;
  else B2FATAL("Invalid QED order: " << m_order);
 
  if (m_vacPol == "off") m_npar[23] = 1;
  else if (m_vacPol == "hadr5") m_npar[23] = 2;
  else if (m_vacPol == "hlmnt") m_npar[23] = 3;
  else B2FATAL("Invalid vacuum polarization code: " << m_vacPol);
 
  if (m_mode == "unweighted" || m_mode == "uw") m_npar[24] = 1;
  else if (m_mode == "weighted" || m_mode == "w") m_npar[24] = 2;
  else B2FATAL("Invalid mode: " << m_mode);
 
  if (m_userMode == "NONE") m_npar[25] = 1;
  else if (m_userMode == "GAMMA") m_npar[25] = 2;
  else if (m_userMode == "EGAMMA") m_npar[25] = 3;
  else if (m_userMode == "ETRON") m_npar[25] = 4;
  else if (m_userMode == "PRESCALE") m_npar[25] = 5;
  else B2FATAL("Invalid user mode: " << m_userMode);
 
  // set the data file for the novosibirsk routine (not supported yet, but needs to be set)
  size_t fileLength = m_NSKDataFile.size();
  babayaga_setvpolnsk_(m_NSKDataFile.c_str(), &fileLength);
 
  //use mode to control init/generation/finalize in FORTRAN code
  int mode = -1;
  double ecm = -1.;
  main_belle2_(&mode, &ecm, m_xpar, m_npar);
}

◆ generateEvent()

void generateEvent	(	MCParticleGraph &	mcGraph,
		double	ecm,
		ROOT::Math::XYZVector	vertex,
		ROOT::Math::LorentzRotation	boost
	)

Generates one single event.

Parameters

mcGraph	Reference to the MonteCarlo graph into which the generated particles will be stored.
ecm	Centre of mass energy.
vertex	Production vertex.
boost	Lorentz boost vector.

Definition at line 331 of file BabayagaNLO.cc.

{
  //Generate event
  int mode = 1;
  main_belle2_(&mode, &ecm, m_xpar, m_npar);
 
  //Store the initial particles as initial (non-virtual) particles into the MCParticleGraph
  storeParticle(mcGraph, momset_.bq1, 11, vertex, boost, false, true, false);
  storeParticle(mcGraph, momset_.bp1, -11, vertex, boost, false, true, false);
 
  //Store the final state fermions or photons (for 'gg') as real particle into the MCParticleGraph
  int pdg = 11;
  int antipdg = -11;
 
  if (m_finalState == "gg") {
    pdg = 22;
    antipdg = 22;
  } else  if (m_finalState == "mm") {
    pdg = 13;
    antipdg = -13;
  }
 
  storeParticle(mcGraph, momset_.bq2, pdg, vertex, boost, false, false, false);
  storeParticle(mcGraph, momset_.bp2, antipdg, vertex, boost, false, false, false);
 
  //Store the real ISR and FSR photons into the MCParticleGraph
  for (int iPhot = 0; iPhot <  momset_.bnphot; ++iPhot) {
    double photMom[4] = {momset_.bphot[0][iPhot], momset_.bphot[1][iPhot], momset_.bphot[2][iPhot], momset_.bphot[3][iPhot]};
    storeParticle(mcGraph, photMom, 22, vertex, boost, false, false, true);
  }
 
  //set event weight
  StoreObjPtr<EventMetaData> eventMetaDataPtr("EventMetaData", DataStore::c_Event);
  eventMetaDataPtr->setGeneratedWeight(momset_.bweight);
 
  //set event prescale
  StoreObjPtr<EventExtraInfo> eventExtraInfo;
  if (not eventExtraInfo.isValid())
    eventExtraInfo.create();
  if (eventExtraInfo->hasExtraInfo("GeneratedPrescale")) {
    B2WARNING("EventExtraInfo with given name is already set! I won't set it again!");
  } else {
    float targetValue = prescale_.bprescale;
    eventExtraInfo->addExtraInfo("GeneratedPrescale", targetValue);
  }
 
}

◆ getSDif()

double getSDif ( )

inline

Returns weight.

Definition at line 190 of file BabayagaNLO.h.

190{return m_sDif;};

◆ init()

void init ( )

Initializes the generator.

Definition at line 325 of file BabayagaNLO.cc.

{
  applySettings();
}

◆ initExtraInfo()

void initExtraInfo ( )

Initializes the extra info.

Definition at line 319 of file BabayagaNLO.cc.

{
  StoreObjPtr<EventExtraInfo> extrainfo;
  extrainfo.registerInDataStore();
}

◆ setCmsEnergyNominal()

void setCmsEnergyNominal ( double cmsEnergyNominal )

inline

Sets nominal ECM.

Parameters

cmsEnergyNominal nominal ecm in GeV

Definition at line 64 of file BabayagaNLO.h.

64{ m_cmsEnergyNominal = cmsEnergyNominal; }

◆ setDefaultSettings()

void setDefaultSettings ( )

Sets the default settings for the BhWide Fortran generator.

Definition at line 274 of file BabayagaNLO.cc.

{
//   m_cmsEnergy = -1.;
  m_cmsEnergyNominal = -1.;
 
  m_applyBoost = true;
 
  m_finalState = "ee";
  m_vacPol     = "hlmnt";
  m_order      = "exp";
  m_model      = "matched";
  m_mode       = "unweighted";
  m_userMode   = "NONE";
 
  m_pi               = 3.1415926535897932384626433832795029;
  m_conversionFactor = 0.389379660e6;
  m_alphaQED0        = 1.0 / 137.0359895;
  m_massElectron     = 0.51099906 * Unit::MeV;
  m_massMuon         = 105.65836900 * Unit::MeV;
  m_massW            = 80.385 * Unit::GeV;
  m_massZ            = 91.1882 * Unit::GeV;
  m_widthZ           = 2.4952 * Unit::GeV;
  m_eMin             = 0.1 * Unit::GeV;
  m_epsilon          = 5.e-7;
  m_maxAcollinearity = 180.0;// * Unit::Deg;
 
  m_ScatteringAngleRange = make_pair(15.0, 165.0); //in [deg]
  m_ScatteringAngleRangePhoton = make_pair(15.0, 165.0); //in [deg] (unused in original babayaga)
 
  m_nPhot      = -1;
  m_nSearchMax = 50000;
  m_fMax       = -1.;
  m_EnergySpread = 5e-3;
  m_VPUncertainty = false;
  m_eemin = 0.0;
  m_temin = 0.0;
  m_egmin = 0.0;
  m_tgmin = 0.0;
  m_eeveto = 0.0;
  m_teveto = 0.0;
  m_egveto = 0.0;
  m_tgveto = 0.0;
  m_maxprescale = 1.0;
}

◆ setEEMIN()

void setEEMIN ( double eemin )

inline

Sets the minimum CMS energy of the tagged e-/e+.

Parameters

eemin minimum CMS energy of the tagged e-/e+

Definition at line 129 of file BabayagaNLO.h.

129{ m_eemin = eemin; }

◆ setEEVETO()

void setEEVETO ( double eeveto )

inline

Sets the minimum CMS energy to veto e-/e+.

Parameters

eeveto minimum CMS energy to veto e-/e+

Definition at line 149 of file BabayagaNLO.h.

149{ m_eeveto = eeveto; }

◆ setEGMIN()

void setEGMIN ( double egmin )

inline

Sets the minimum CMS energy of the gamma.

Parameters

egmin minimum CMS energy of the gamma

Definition at line 139 of file BabayagaNLO.h.

139{ m_egmin = egmin; }

◆ setEGVETO()

void setEGVETO ( double egveto )

inline

Sets the minimum CMS energy to veto gamma.

Parameters

egveto minimum CMS energy to veto gamma

Definition at line 159 of file BabayagaNLO.h.

159{ m_egveto = egveto; }

◆ setEnergySpread()

void setEnergySpread ( double spread )

inline

TEMPORARY SOLUTION! Approximate energy spread per beam (CMS)

Parameters

spread energy spread per beam (CMS)

Definition at line 104 of file BabayagaNLO.h.

104{ m_EnergySpread = spread; }

◆ setEpsilon()

void setEpsilon ( double epsilon )

inline

Sets soft/hard photon energy separator.

Parameters

epsilon soft/hard photon energy separator, in units of ECMS/2

Definition at line 74 of file BabayagaNLO.h.

74{ m_epsilon = epsilon; }

◆ setFinalState()

void setFinalState ( const std::string & finalState )

inline

Sets final state minimum energy.

Parameters

finalState Minimum energy for leptons/photons in the final state, in GeV

Definition at line 99 of file BabayagaNLO.h.

99{ m_finalState = finalState; }

◆ setFMax()

void setFMax ( double fMax )

inline

Maximum differential cross section.

Parameters

fMax	Maximum differential cross section

Definition at line 54 of file BabayagaNLO.h.

54{ m_fMax = fMax; }

◆ setMaxAcollinearity()

void setMaxAcollinearity ( double maxAcollinearity )

inline

Sets maximum acollinearity angle between finale state leptons/photons in degrees.

Parameters

maxAcollinearity Maximum acollinearity angle between finale state leptons/photons in degrees

Definition at line 59 of file BabayagaNLO.h.

59{ m_maxAcollinearity = maxAcollinearity; }

◆ setMaxPrescale()

void setMaxPrescale ( double maxprescale )

inline

Sets the maximum prescale value.

Parameters

maxprescale maximum prescale value (at maxmin theta)

Definition at line 169 of file BabayagaNLO.h.

169{ m_maxprescale = maxprescale; }

◆ setMinEnergy()

void setMinEnergy ( double minEnergy )

inline

Sets minimum energy for leptons/photons in the final state, in GeV.

Parameters

minEnergy Minimum energy for leptons/photons in the final state, in GeV

Definition at line 69 of file BabayagaNLO.h.

69{ m_eMin = minEnergy; }

◆ setMode()

void setMode ( const std::string & mode )

inline

Sets mode: weighted or unweighted.

Parameters

mode	weighted or unweighted

Definition at line 89 of file BabayagaNLO.h.

89{ m_mode = mode; }

◆ setModel()

void setModel ( const std::string & model )

inline

Sets model: matched or ps.

Parameters

model matched or ps

Definition at line 84 of file BabayagaNLO.h.

84{ m_model = model; }

◆ setNPhotons()

void setNPhotons ( int nPhot )

inline

Sets the fixed number of nphot (hard) photons.

Parameters

nPhot number of nphot (hard) photons

Definition at line 114 of file BabayagaNLO.h.

114{ m_nPhot = nPhot; }

◆ setNSearchMax()

void setNSearchMax ( int nSearchMax )

inline

Sets the number of events used to search maximum.

Parameters

nSearchMax Number of events for maximum search

Definition at line 49 of file BabayagaNLO.h.

49{ m_nSearchMax = nSearchMax; }

◆ setNSKDataFile()

void setNSKDataFile ( const std::string & NSKDataFile )

inline

Sets NSK VP data file.

Parameters

NSKDataFile (default)

Definition at line 124 of file BabayagaNLO.h.

124{ m_NSKDataFile = NSKDataFile; }

◆ setOrder()

void setOrder ( const std::string & order )

inline

Sets Order: born, alpha or exp.

Parameters

order born, alpha or exp

Definition at line 94 of file BabayagaNLO.h.

94{ m_order = order; }

◆ setScatAngle()

void setScatAngle ( std::pair< double, double > angleRange )

inline

Sets the theta scattering angle range for the scattered particles.

Parameters

angleRange A pair of values, representing the min and max theta angle of the scattered particles in [deg].

Definition at line 44 of file BabayagaNLO.h.

44{ m_ScatteringAngleRange = angleRange; }

◆ setTEMIN()

void setTEMIN ( double temin )

inline

Sets the minimum CMS angle between the tagged e-/e+ and -z axis.

Parameters

temin minimum CMS energy of the tagged e-/e+

Definition at line 134 of file BabayagaNLO.h.

134{ m_temin = temin; }

◆ setTEVETO()

void setTEVETO ( double teveto )

inline

Sets the maximum CMS theta of e-/e+ in final state.

Parameters

teveto maximum CMS theta of e-/e+ in final state

Definition at line 154 of file BabayagaNLO.h.

154{ m_teveto = teveto; }

◆ setTGMIN()

void setTGMIN ( double tgmin )

inline

Sets the minimum CMS angle between the gamma and -z axis.

Parameters

tgmin minimum CMS angle between the gamma and -z axis

Definition at line 144 of file BabayagaNLO.h.

144{ m_tgmin = tgmin; }

◆ setTGVETO()

void setTGVETO ( double tgveto )

inline

Sets the maximum CMS angle between the gamma and -z axis.

Parameters

tgveto maximum CMS angle between the gamma and -z axis

Definition at line 164 of file BabayagaNLO.h.

164{ m_tgveto = tgveto; }

◆ setUserMode()

void setUserMode ( const std::string & usermode )

inline

Sets User mode similar to TEEGG: ETRON, EGAMMA, GAMMA or PRESCALE or NONE (default)

Parameters

usermode ETRON, EGAMMA, GAMMA or NONE (default)

Definition at line 119 of file BabayagaNLO.h.

119{ m_userMode = usermode; }

◆ setVacPol()

void setVacPol ( const std::string & vacPol )

inline

Sets vacuum polarization.

Parameters

vacPol Vacuum polarization option

Definition at line 79 of file BabayagaNLO.h.

79{ m_vacPol = vacPol; }

◆ setVPUncertainty()

void setVPUncertainty ( bool vpuncertainty )

inline

Calculate VP uncertainty by internal reweighting (on/off)

Parameters

vpuncertainty boolean on/off

Definition at line 109 of file BabayagaNLO.h.

109{ m_VPUncertainty = vpuncertainty; }

◆ storeParticle()

void storeParticle	(	MCParticleGraph &	mcGraph,
		const double *	mom,
		int	pdg,
		ROOT::Math::XYZVector	vertex,
		ROOT::Math::LorentzRotation	boost,
		bool	isVirtual = `false`,
		bool	isInitial = `false`,
		bool	isISRFSR = `false`
	)

protected

Store a single generated particle into the MonteCarlo graph.

Parameters

mcGraph	Reference to the MonteCarlo graph into which the particle should be stored.
mom	The 3-momentum of the particle in [GeV].
pdg	The PDG code of the particle.
vertex	Production vertex.
boost	Lorentz boost vector.
isVirtual	If the particle is a virtual particle, such as the incoming particles, set this to true.
isInitial	If the particle is a initial particle for ISR, set this to true.
isISRFSR	If the particle is from ISR or FSR, set this to true.

Definition at line 488 of file BabayagaNLO.cc.

{
 
//   Create particle
  MCParticleGraph::GraphParticle& part = mcGraph.addParticle();
  if (isVirtual) {
    part.setStatus(MCParticle::c_IsVirtual);
  } else if (isInitial) {
    part.setStatus(MCParticle::c_Initial);
  }
 
  // all particle of a generator are primary
  part.addStatus(MCParticle::c_PrimaryParticle);
 
  // all particles produced by BABAYAGA are stable
  part.addStatus(MCParticle::c_StableInGenerator);
 
  // add ISR/FSR flags to all photons that are not the primary gg pair
  if (isISRFSR) {
    part.addStatus(MCParticle::c_IsISRPhoton);
    part.addStatus(MCParticle::c_IsFSRPhoton);
  }
 
  part.setPDG(pdg);
  part.setFirstDaughter(0);
  part.setLastDaughter(0);
  part.setMomentum(ROOT::Math::XYZVector(mom[0], mom[1], mom[2]));
  part.setMass(TDatabasePDG::Instance()->GetParticle(pdg)->Mass());
  part.setEnergy(mom[3]);
 
  //boost
  ROOT::Math::PxPyPzEVector p4 = part.get4Vector();
  p4 = boost * p4;
  part.set4Vector(p4);
 
  //set vertex
  if (!isInitial) {
    ROOT::Math::XYZVector v3 = part.getProductionVertex();
    v3 = v3 + vertex;
    part.setProductionVertex(v3);
    part.setValidVertex(true);
  }
}

◆ term()

void term ( )

Terminates the generator.

Closes the internal Fortran generator.

Definition at line 381 of file BabayagaNLO.cc.

{
 
  B2RESULT("Babayaga.NLO: Final state:    " << m_finalState);
  B2RESULT("Babayaga.NLO: Mode:           " << m_mode);
  B2RESULT("Babayaga.NLO: Order:          " << m_order);
  B2RESULT("Babayaga.NLO: Model:          " << m_model);
  B2RESULT("Babayaga.NLO: Vac. pol. (VP): " << m_vacPol);
  B2RESULT("Babayaga.NLO: Usercuts:       " << m_userMode);
 
  // all other results are displayed from fortran using a callback to C extern logging
 
  int mode = 2;
  double ecm = -1.;
  main_belle2_(&mode, &ecm, m_xpar, m_npar);
 
}

Member Data Documentation

◆ m_alphaQED0

double m_alphaQED0

protected

QED coupling constant at Q=0.

Definition at line 202 of file BabayagaNLO.h.

◆ m_applyBoost

bool m_applyBoost

protected

Apply a boost to the MCParticles.

Definition at line 199 of file BabayagaNLO.h.

◆ m_boostVector

ROOT::Math::LorentzRotation m_boostVector

protected

The Lorentz boost vector for the transformation CMS to LAB frame.

Definition at line 240 of file BabayagaNLO.h.

◆ m_cmsEnergyNominal

double m_cmsEnergyNominal

protected

Nominal CMS Energy = 2*Ebeam [GeV].

Definition at line 215 of file BabayagaNLO.h.

◆ m_conversionFactor

double m_conversionFactor

protected

Conversion factor for hbarc to nb.

Definition at line 201 of file BabayagaNLO.h.

◆ m_eemin

double m_eemin

protected

Minimum CMS energy of the tagged e-/e+ (GeV).

Definition at line 227 of file BabayagaNLO.h.

◆ m_eeveto

double m_eeveto

protected

Minimum CMS energy to veto e-/e+ (GeV).

Definition at line 231 of file BabayagaNLO.h.

◆ m_egmin

double m_egmin

protected

Minimum CMS energy of the gamma (GeV).

Definition at line 229 of file BabayagaNLO.h.

◆ m_egveto

double m_egveto

protected

Minimum CMS energy to veto gamma (GeV).

Definition at line 233 of file BabayagaNLO.h.

◆ m_eMin

double m_eMin

protected

[GeV].

Definition at line 219 of file BabayagaNLO.h.

◆ m_EnergySpread

double m_EnergySpread

protected

TEMPORARY SOLUTION! Approximate energy spread per beam (CMS).

Definition at line 214 of file BabayagaNLO.h.

◆ m_epsilon

double m_epsilon

protected

Soft/hard photon separator in units of CMS/2.

Definition at line 218 of file BabayagaNLO.h.

◆ m_finalState

std::string m_finalState

protected

final state: ee, gg or mm.

Definition at line 210 of file BabayagaNLO.h.

◆ m_fMax

double m_fMax

protected

Maximum of differential cross section.

Definition at line 237 of file BabayagaNLO.h.

◆ m_massElectron

double m_massElectron

protected

muon mass.

Definition at line 203 of file BabayagaNLO.h.

◆ m_massMuon

double m_massMuon

protected

electron mass.

Definition at line 204 of file BabayagaNLO.h.

◆ m_massW

double m_massW

protected

W mass [GeV] for on shell sin2theta and GF.

Definition at line 205 of file BabayagaNLO.h.

◆ m_massZ

double m_massZ

protected

Z mass [GeV].

Definition at line 206 of file BabayagaNLO.h.

◆ m_maxAcollinearity

double m_maxAcollinearity

protected

Maximum acollinearity of the electron-positron pair.

Definition at line 220 of file BabayagaNLO.h.

◆ m_maxprescale

double m_maxprescale

protected

Maximum prescale value.

Definition at line 235 of file BabayagaNLO.h.

◆ m_mode

std::string m_mode

protected

mode: weighted or unweighted.

Definition at line 212 of file BabayagaNLO.h.

◆ m_model

std::string m_model

protected

model: matched or ps.

Definition at line 211 of file BabayagaNLO.h.

◆ m_npar

int m_npar[100]

private

Integer parameters for BabayagaNLO.

Definition at line 261 of file BabayagaNLO.h.

◆ m_nPhot

int m_nPhot

protected

fixed number of nphot (hard) photons are generated.

A negative value means all possible photons.

Definition at line 222 of file BabayagaNLO.h.

◆ m_nSearchMax

int m_nSearchMax

protected

Events used to search maximum of differential cross section.

Definition at line 221 of file BabayagaNLO.h.

◆ m_NSKDataFile

std::string m_NSKDataFile

protected

data file for NSK VP.

Definition at line 224 of file BabayagaNLO.h.

◆ m_order

std::string m_order

protected

order: born, alpha or exp.

Definition at line 213 of file BabayagaNLO.h.

◆ m_pi

double m_pi

protected

pi=3.1415....

Definition at line 200 of file BabayagaNLO.h.

◆ m_ScatteringAngleRange

std::pair<double, double> m_ScatteringAngleRange

protected

Min and Max val.

for the scattering angle [deg] of the electron/positron.

Definition at line 216 of file BabayagaNLO.h.

◆ m_ScatteringAngleRangePhoton

std::pair<double, double> m_ScatteringAngleRangePhoton

protected

Min and Max val.

for the scattering angle [deg] of the g in gg mode.

Definition at line 217 of file BabayagaNLO.h.

◆ m_sDif

double m_sDif

protected

Differential xsec/weight used for event.

Definition at line 238 of file BabayagaNLO.h.

◆ m_temin

double m_temin

protected

Minimum CMS angle between the tagged e-/e+ and -z axis (deg).

Definition at line 228 of file BabayagaNLO.h.

◆ m_teveto

double m_teveto

protected

Maximum CMS theta of e-/e+ in final state (deg).

Definition at line 232 of file BabayagaNLO.h.

◆ m_tgmin

double m_tgmin

protected

Minimum CMS angle between the gamma and -z axis (deg).

Definition at line 230 of file BabayagaNLO.h.

◆ m_tgveto

double m_tgveto

protected

Maximum CMS angle between the gamma and -z axis (deg).

Definition at line 234 of file BabayagaNLO.h.

◆ m_userMode

std::string m_userMode

protected

User mode similar to TEEGG: ETRON, EGAMMA, GAMMA or PRESCALE.

Definition at line 226 of file BabayagaNLO.h.

◆ m_vacPol

std::string m_vacPol

protected

vacuum polarization: off, hadr5 (Jegerlehner) or hmnt (Teubner).

Definition at line 209 of file BabayagaNLO.h.

◆ m_VPUncertainty

bool m_VPUncertainty

protected

vary all VP related parameters and extracted total uncertainty.

Definition at line 223 of file BabayagaNLO.h.

◆ m_widthZ

double m_widthZ

protected

Z width [GeV] (may be recalculated by EW library).

Definition at line 207 of file BabayagaNLO.h.

◆ m_xpar

double m_xpar[100]

private

Double parameters for BabayagaNLO.

Definition at line 262 of file BabayagaNLO.h.

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

generators/babayaganlo/include/BabayagaNLO.h
generators/babayaganlo/src/BabayagaNLO.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BabayagaNLO()

◆ ~BabayagaNLO()

Member Function Documentation

◆ applySettings()

◆ generateEvent()

◆ getSDif()

◆ init()

◆ initExtraInfo()

◆ setCmsEnergyNominal()

◆ setDefaultSettings()

◆ setEEMIN()

◆ setEEVETO()

◆ setEGMIN()

◆ setEGVETO()

◆ setEnergySpread()

◆ setEpsilon()

◆ setFinalState()

◆ setFMax()

◆ setMaxAcollinearity()

◆ setMaxPrescale()

◆ setMinEnergy()

◆ setMode()

◆ setModel()

◆ setNPhotons()

◆ setNSearchMax()

◆ setNSKDataFile()

◆ setOrder()

◆ setScatAngle()

◆ setTEMIN()

◆ setTEVETO()

◆ setTGMIN()

◆ setTGVETO()

◆ setUserMode()

◆ setVacPol()

◆ setVPUncertainty()

◆ storeParticle()

◆ term()

Member Data Documentation

◆ m_alphaQED0

◆ m_applyBoost

◆ m_boostVector

◆ m_cmsEnergyNominal

◆ m_conversionFactor

◆ m_eemin

◆ m_eeveto

◆ m_egmin

◆ m_egveto

◆ m_eMin

◆ m_EnergySpread

◆ m_epsilon

◆ m_finalState

◆ m_fMax

◆ m_massElectron

◆ m_massMuon

◆ m_massW

◆ m_massZ

◆ m_maxAcollinearity

◆ m_maxprescale

◆ m_mode

◆ m_model

◆ m_npar

◆ m_nPhot

◆ m_nSearchMax

◆ m_NSKDataFile

◆ m_order

◆ m_pi

◆ m_ScatteringAngleRange

◆ m_ScatteringAngleRangePhoton

◆ m_sDif

◆ m_temin

◆ m_teveto

◆ m_tgmin

◆ m_tgveto

◆ m_userMode