Compute the matrix element for a generic factorised process defined in a Fortran weighting function. More...

#include <include/CepGen/Process/FortranFactorisedProcess.h>

Inheritance diagram for FortranFactorisedProcess:

Collaboration diagram for FortranFactorisedProcess:

Public Types
enum class	Mapping { linear = 0 , exponential , square , power_law }
	Type of mapping to apply on the variable. More...

Public Member Functions
	FortranFactorisedProcess (const ParametersList &, const std::function< double()> &func)
	Construct a Fortran-CepGen interface object using a double precision argument-less F77 function.

ProcessPtr	clone () const override
	Copy all process attributes into a new object.

double	computeWeight () override
	Compute the phase space point weight.

void	fillKinematics () final
	Fill the Event object with the particles' kinematics.

bool	validatedBeamKinematics ()

utils::RandomGenerator &	randomGenerator () const
	Accessor for this process' random number generator.

void	clear ()
	Reset process prior to the phase space and variables definition.

void	clearEvent ()
	Restore the event object to its initial state.

void	initialise ()
	Initialise the process once the kinematics has been set.

const Kinematics &	kinematics () const
	Constant reference to the process kinematics.

Kinematics &	kinematics ()
	Reference to the process kinematics.

void	setKinematics ()

double	weight (const std::vector< double > &)
	Compute the weight for a phase-space point.

void	dumpPoint (std::ostream *=nullptr) const
	Dump the coordinate of the phase-space point being evaluated.

void	dumpVariables (std::ostream *=nullptr) const
	List all variables handled by this generic process.

size_t	ndim () const
	Number of dimensions to perform integration.

bool	hasEvent () const
	Does the process contain (and hold) an event?

const Event &	event () const
	Handled particles objects and their relationships.

Event &	event ()
	Event object read/write accessor.

Event *	eventPtr ()
	Event pointer read/write accessor.

const Momentum &	pA () const
	Positive-z incoming beam particle's 4-momentum.

Momentum &	pA ()
	Positive-z incoming beam particle's 4-momentum.

double	mA2 () const
	Positive-z incoming beam particle's squared mass.

double	mA () const
	Positive-z incoming beam particle's mass.

const Momentum &	pB () const
	Negative-z incoming beam particle's 4-momentum.

Momentum &	pB ()
	Negative-z incoming beam particle's 4-momentum.

double	mB2 () const
	Negative-z incoming beam particle's squared mass.

double	mB () const
	Negative-z incoming beam particle's mass.

const Momentum &	pX () const
	Positive-z outgoing beam particle's 4-momentum.

Momentum &	pX ()
	Positive-z outgoing beam particle's 4-momentum.

double	mX2 () const
	Positive-z outgoing beam particle's squared mass.

double &	mX2 ()
	Positive-z outgoing beam particle's squared mass.

double	mX () const
	Positive-z outgoing beam particle's mass.

const Momentum &	pY () const
	Negative-z outgoing beam particle's 4-momentum.

Momentum &	pY ()
	Negative-z outgoing beam particle's 4-momentum.

double	mY2 () const
	Negative-z outgoing beam particle's squared mass.

double &	mY2 ()
	Negative-z outgoing beam particle's squared mass.

double	mY () const
	Negative-z outgoing beam particle's mass.

double	x1 () const
	Positive-z incoming parton's fractional momentum.

double &	x1 ()
	Positive-z incoming parton's fractional momentum.

double	t1 () const
	Positive-z incoming parton's squared mass.

double &	t1 ()
	Positive-z incoming parton's squared mass.

const Momentum &	q1 () const
	Positive-z incoming parton's 4-momentum.

Momentum &	q1 ()
	Positive-z incoming parton's 4-momentum.

double	x2 () const
	Negative-z incoming parton's fractional momentum.

double &	x2 ()
	Negative-z incoming parton's fractional momentum.

double	t2 () const
	Negative-z incoming parton's squared mass.

double &	t2 ()
	Negative-z incoming parton's squared mass.

const Momentum &	q2 () const
	Negative-z incoming parton's 4-momentum.

Momentum &	q2 ()
	Negative-z incoming parton's 4-momentum.

double	wCM () const
	Two-parton centre of mass energy.

Momentum &	pc (size_t)
	Central particle's 4-momentum.

const Momentum &	pc (size_t) const
	Central particle's 4-momentum.

double	s () const
	Two-beam squared centre of mass energy.

double	sqrtS () const
	Two-beam centre of mass energy.

double	inverseS () const
	Inverse two-beam squared centre of mass energy.

double	inverseSqrtS () const
	Inverse two-beam centre of mass energy.

Process &	defineVariable (double &out, const Mapping &type, const Limits &lim, const std::string &name, const std::string &description="")

double	variableValue (size_t i, double x) const
	Retrieve the physical value for one variable.

const std::string &	name () const
	Module unique indexing name.

bool	operator== (const SteeredObject &oth) const
	Equality operator.

bool	operator!= (const SteeredObject &oth) const
	Inequality operator.

const ParametersList &	parameters () const override
	Module user-defined parameters.

void	setParameters (const ParametersList &params) override
	Set module parameters.

void	setDescribedParameters (const ParametersList &params_orig)
	Set (documented) module parameters.

Static Public Member Functions
static ParametersDescription	description ()

Static Public Attributes
static ParametersList	kProcParameters

Protected Member Functions
void	addEventContent () override
	Set the incoming and outgoing state to be expected in the process.

void	prepareKinematics () final
	Compute the incoming state kinematics.

virtual void	computeBeamKinematics ()
	Compute the outgoing protons (or remnants) kinematics.

double	that () const
	\(\hat t=\frac{1}{2}\left[(p_1-p_3)^2+(p_2-p_4)^2\right]\)

double	uhat () const
	\(\hat u=\frac{1}{2}\left[(p_1-p_4)^2+(p_2-p_3)^2\right]\)

double	shat () const
	\(\hat s=(p_1+p_2)^2=(p_3+...)^2\)

double	generateVariables () const
	Generate and initialise all variables handled by this process.

void	setEventContent (const std::unordered_map< Particle::Role, spdgids_t > &)
	Set the incoming and outgoing states to be defined in this process (and prepare the Event object accordingly)

double	alphaEM (double q) const
	Compute the electromagnetic running coupling algorithm at a given scale.

double	alphaS (double q) const
	Compute the strong coupling algorithm at a given scale.

const std::vector< double > &	lastCoordinates () const
	Last coordinates fed.

template<typename T >
T	steer (const std::string &key) const
	Retrieve a parameters as previously steered.

template<typename T , typename U >
U	steerAs (const std::string &key) const
	Retrieve a recast parameters as previously steered.

std::string	steerName () const
	Retrieve module name from parameters.

std::string	steerPath (const std::string &key) const
	Retrieve a path from common search paths.

Protected Attributes
const std::unique_ptr< PhaseSpaceGenerator >	phase_space_generator_
	Kinematic variables generator for the phase space coverage.

const bool	symmetrise_

const bool	store_alphas_

double	mp_
	Proton mass, in GeV/c \(^2\).

double	mp2_
	Squared proton mass, in GeV \(^2\)/c \(^4\).

const std::string	name_
	Module unique indexing name.

ParametersList	params_
	Module parameters.

Detailed Description

Definition at line 26 of file FortranFactorisedProcess.h.

Member Enumeration Documentation

◆ Mapping

enum class Mapping

stronginherited

Enumerator
linear	a linear \({\rm d}x\) mapping
exponential	an exponential \(\frac{\dot{x}}{x} = \dot{\log x}\) mapping
square	a square \({\rm d}x^2=2x\cdot\dot{x}\) mapping
power_law	a power-law mapping inherited from LPAIR

Definition at line 125 of file Process.h.

Constructor & Destructor Documentation

◆ FortranFactorisedProcess()

FortranFactorisedProcess	(	const ParametersList &	,
		const std::function< double()> &	func )

explicit

Parameters

[in] func a double precision argument-less Fortran function returning the event weight

Member Function Documentation

◆ addEventContent()

void addEventContent ( )

overrideprotectedvirtualinherited

Implements Process.

◆ alphaEM()

double alphaEM ( double q ) const

protectedinherited

◆ alphaS()

double alphaS ( double q ) const

protectedinherited

◆ clear()

void clear ( )

inherited

◆ clearEvent()

void clearEvent ( )

inherited

◆ clone()

ProcessPtr clone ( ) const

overridevirtual

Reimplemented from Process.

Definition at line 31 of file FortranFactorisedProcess.h.

◆ computeBeamKinematics()

virtual void computeBeamKinematics ( )

protectedvirtualinherited

◆ computeWeight()

double computeWeight ( )

overridevirtualinherited

Implements Process.

◆ defineVariable()

Process & defineVariable	(	double &	out,
		const Mapping &	type,
		const Limits &	lim,
		const std::string &	name,
		const std::string &	description = "" )

inherited

◆ description()

static ParametersDescription description ( )

staticinherited

◆ dumpPoint()

void dumpPoint ( std::ostream * = nullptr ) const

inherited

◆ dumpVariables()

void dumpVariables ( std::ostream * = nullptr ) const

inherited

◆ event() [1/2]

Event & event ( )

inherited

◆ event() [2/2]

const Event & event ( ) const

inherited

◆ eventPtr()

Event * eventPtr ( )

inherited

◆ fillKinematics()

void fillKinematics ( )

finalvirtualinherited

Implements Process.

◆ generateVariables()

double generateVariables ( ) const

protectedinherited

Returns: Phase space point-dependent component of the Jacobian weight of the point in the phase space for integration

Note: To be run at each point computation (therefore, to be optimised!)

◆ hasEvent()

bool hasEvent ( ) const

inherited

Definition at line 62 of file Process.h.

◆ initialise()

void initialise ( )

inherited

◆ inverseS()

double inverseS ( ) const

inherited

Definition at line 109 of file Process.h.

◆ inverseSqrtS()

double inverseSqrtS ( ) const

inherited

Definition at line 110 of file Process.h.

◆ kinematics() [1/2]

Kinematics & kinematics ( )

inherited

Definition at line 51 of file Process.h.

◆ kinematics() [2/2]

const Kinematics & kinematics ( ) const

inherited

Definition at line 50 of file Process.h.

◆ lastCoordinates()

const std::vector< double > & lastCoordinates ( ) const

protectedinherited

Definition at line 159 of file Process.h.

◆ mA()

double mA ( ) const

inherited

Definition at line 69 of file Process.h.

◆ mA2()

double mA2 ( ) const

inherited

Definition at line 68 of file Process.h.

◆ mB()

double mB ( ) const

inherited

Definition at line 72 of file Process.h.

◆ mB2()

double mB2 ( ) const

inherited

Definition at line 71 of file Process.h.

◆ mX()

double mX ( ) const

inherited

Definition at line 75 of file Process.h.

◆ mX2() [1/2]

double & mX2 ( )

inherited

Definition at line 97 of file Process.h.

◆ mX2() [2/2]

double mX2 ( ) const

inherited

Definition at line 74 of file Process.h.

◆ mY()

double mY ( ) const

inherited

Definition at line 78 of file Process.h.

◆ mY2() [1/2]

double & mY2 ( )

inherited

Definition at line 98 of file Process.h.

◆ mY2() [2/2]

double mY2 ( ) const

inherited

Definition at line 77 of file Process.h.

◆ name()

const std::string & name ( ) const

inherited

Definition at line 42 of file NamedModule.h.

◆ ndim()

size_t ndim ( ) const

inherited

Definition at line 59 of file Process.h.

◆ operator!=()

template<typename T >

bool operator!= ( const SteeredObject< T > & oth ) const

inherited

Definition at line 50 of file SteeredObject.h.

◆ operator==()

template<typename T >

bool operator== ( const SteeredObject< T > & oth ) const

inherited

Definition at line 48 of file SteeredObject.h.

◆ pA() [1/2]

Momentum & pA ( )

inherited

◆ pA() [2/2]

const Momentum & pA ( ) const

inherited

◆ parameters()

template<typename T >

const ParametersList & parameters ( ) const

overridevirtualinherited

Reimplemented from Steerable.

Definition at line 53 of file SteeredObject.h.

◆ pB() [1/2]

Momentum & pB ( )

inherited

◆ pB() [2/2]

const Momentum & pB ( ) const

inherited

◆ pc() [1/2]

Momentum & pc ( size_t )

inherited

◆ pc() [2/2]

const Momentum & pc ( size_t ) const

inherited

◆ prepareKinematics()

void prepareKinematics ( )

finalprotectedvirtualinherited

Reimplemented from Process.

◆ pX() [1/2]

Momentum & pX ( )

inherited

◆ pX() [2/2]

const Momentum & pX ( ) const

inherited

◆ pY() [1/2]

Momentum & pY ( )

inherited

◆ pY() [2/2]

const Momentum & pY ( ) const

inherited

◆ q1() [1/2]

Momentum & q1 ( )

inherited

◆ q1() [2/2]

const Momentum & q1 ( ) const

inherited

◆ q2() [1/2]

Momentum & q2 ( )

inherited

◆ q2() [2/2]

const Momentum & q2 ( ) const

inherited

◆ randomGenerator()

utils::RandomGenerator & randomGenerator ( ) const

inherited

◆ s()

double s ( ) const

inherited

Definition at line 107 of file Process.h.

◆ setDescribedParameters()

template<typename T >

void setDescribedParameters ( const ParametersList & params_orig )

inherited

Definition at line 72 of file SteeredObject.h.

◆ setEventContent()

void setEventContent ( const std::unordered_map< Particle::Role, spdgids_t > & )

protectedinherited

◆ setKinematics()

void setKinematics ( )

inherited

◆ setParameters()

template<typename T >

void setParameters ( const ParametersList & )

overridevirtualinherited

Reimplemented from Steerable.

Definition at line 61 of file SteeredObject.h.

◆ shat()

double shat ( ) const

protectedinherited

◆ sqrtS()

double sqrtS ( ) const

inherited

Definition at line 108 of file Process.h.

◆ steer()

template<typename T >

T steer ( const std::string & key ) const

protectedinherited

Definition at line 39 of file Steerable.h.

◆ steerAs()

template<typename T , typename U >

U steerAs ( const std::string & key ) const

protectedinherited

Definition at line 44 of file Steerable.h.

◆ steerName()

std::string steerName ( ) const

protectedinherited

Definition at line 48 of file Steerable.h.

◆ steerPath()

std::string steerPath ( const std::string & key ) const

protectedinherited

◆ t1() [1/2]

double & t1 ( )

inherited

Definition at line 99 of file Process.h.

◆ t1() [2/2]

double t1 ( ) const

inherited

Definition at line 81 of file Process.h.

◆ t2() [1/2]

double & t2 ( )

inherited

Definition at line 101 of file Process.h.

◆ t2() [2/2]

double t2 ( ) const

inherited

Definition at line 84 of file Process.h.

◆ that()

double that ( ) const

protectedinherited

◆ uhat()

double uhat ( ) const

protectedinherited

◆ validatedBeamKinematics()

bool validatedBeamKinematics ( )

inherited

◆ variableValue()

double variableValue	(	size_t	i,
		double	x ) const

inherited

◆ wCM()

double wCM ( ) const

inherited

Definition at line 90 of file Process.h.

◆ weight()

double weight ( const std::vector< double > & )

inherited

◆ x1() [1/2]

double & x1 ( )

inherited

Definition at line 100 of file Process.h.

◆ x1() [2/2]

double x1 ( ) const

inherited

Definition at line 80 of file Process.h.

◆ x2() [1/2]

double & x2 ( )

inherited

Definition at line 102 of file Process.h.

◆ x2() [2/2]

double x2 ( ) const

inherited

Definition at line 83 of file Process.h.

Member Data Documentation

◆ kProcParameters

ParametersList kProcParameters

static

Definition at line 33 of file FortranFactorisedProcess.h.

◆ mp2_

double mp2_

protectedinherited

Definition at line 121 of file Process.h.

◆ mp_

double mp_

protectedinherited

Definition at line 120 of file Process.h.

◆ name_

const std::string name_

protectedinherited

Definition at line 45 of file NamedModule.h.

◆ params_

ParametersList params_

mutableprotectedinherited

Definition at line 50 of file Steerable.h.

◆ phase_space_generator_

const std::unique_ptr<PhaseSpaceGenerator> phase_space_generator_

protectedinherited

Definition at line 63 of file FactorisedProcess.h.

◆ store_alphas_

const bool store_alphas_

protectedinherited

Definition at line 65 of file FactorisedProcess.h.

◆ symmetrise_

const bool symmetrise_

protectedinherited

Definition at line 64 of file FactorisedProcess.h.

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ Mapping

Constructor & Destructor Documentation

◆ FortranFactorisedProcess()

Member Function Documentation

◆ addEventContent()

◆ alphaEM()

◆ alphaS()

◆ clear()

◆ clearEvent()

◆ clone()

◆ computeBeamKinematics()

◆ computeWeight()

◆ defineVariable()

◆ description()

◆ dumpPoint()

◆ dumpVariables()

◆ event() [1/2]

◆ event() [2/2]

◆ eventPtr()

◆ fillKinematics()

◆ generateVariables()

◆ hasEvent()

◆ initialise()

◆ inverseS()

◆ inverseSqrtS()

◆ kinematics() [1/2]

◆ kinematics() [2/2]

◆ lastCoordinates()

◆ mA()

◆ mA2()

◆ mB()

◆ mB2()

◆ mX()

◆ mX2() [1/2]

◆ mX2() [2/2]

◆ mY()

◆ mY2() [1/2]

◆ mY2() [2/2]

◆ name()

◆ ndim()

◆ operator!=()

◆ operator==()

◆ pA() [1/2]

◆ pA() [2/2]

◆ parameters()

◆ pB() [1/2]

◆ pB() [2/2]

◆ pc() [1/2]

◆ pc() [2/2]

◆ prepareKinematics()

◆ pX() [1/2]

◆ pX() [2/2]

◆ pY() [1/2]

◆ pY() [2/2]

◆ q1() [1/2]

◆ q1() [2/2]

◆ q2() [1/2]

◆ q2() [2/2]

◆ randomGenerator()

◆ s()

◆ setDescribedParameters()

◆ setEventContent()

◆ setKinematics()

◆ setParameters()

◆ shat()

◆ sqrtS()

◆ steer()

◆ steerAs()

◆ steerName()

◆ steerPath()

◆ t1() [1/2]

◆ t1() [2/2]