#include <ForwardProton.h>

Public Types
enum	{ dimension = 5 }
	parameter dimension More...

typedef math::ErrorF< dimension >::type	CovarianceMatrix
	dimension-parameter covariance matrix More...

enum	Index { Index::xi, Index::th_x, Index::vtx_x, Index::th_y, Index::vtx_y, Index::num_indices = dimension }
	indices to the covariance matrix More...

enum	LHCSector { LHCSector::invalid = -1, LHCSector::sector45, LHCSector::sector56 }
	LHC sector. More...

typedef math::XYZPointF	Point
	point in the space More...

enum	ReconstructionMethod { ReconstructionMethod::invalid = -1, ReconstructionMethod::singleRP, ReconstructionMethod::multiRP }
	type of reconstruction applied for this track More...

typedef math::XYZVectorF	Vector
	spatial vector More...

Public Member Functions
float	chi2 () const
	chi-squared of the fit More...

const CTPPSLocalTrackLiteRefVector &	contributingLocalTracks () const
	list of RP tracks that contributed to this global track More...

double	error (Index i) const
	return the uncertainty on a given component More...

	ForwardProton ()
	default constructor More...

	ForwardProton (double chi2, double ndof, const Point &vtx, const Vector &momentum, float xi, const CovarianceMatrix &cov, ReconstructionMethod method, const CTPPSLocalTrackLiteRefVector &local_tracks, bool valid, const float time=0., const float time_err=0.)
	constructor from refit parameters, fitted vertex and momentum, and longitudinal fractional momentum loss More...

LHCSector	lhcSector () const

ReconstructionMethod	method () const
	reconstruction method for this track More...

const Vector &	momentum () const
	fitted track direction More...

unsigned int	ndof () const
	number of degrees of freedom for the track fit More...

float	normalizedChi2 () const
	chi-squared divided by ndof (or chi-squared * 1e6 if ndof is zero) More...

float	p () const
	scalar norm of fitted track momentum More...

float	pt () const
	scalar fitted track transverse momentum More...

float	px () const
	fitted track momentum horizontal component More...

float	py () const
	fitted track momentum vertical component More...

float	pz () const
	fitted track momentum longitudinal component More...

void	setContributingLocalTracks (const CTPPSLocalTrackLiteRefVector &v)
	store the list of RP tracks that contributed to this global track More...

void	setMethod (const ReconstructionMethod &method)
	set the reconstruction method for this track More...

void	setTime (float time)

void	setTimeError (float time_err)

void	setValidFit (bool valid=true)
	set the flag for the fit validity More...

float	t () const
	four-momentum transfer squared, in GeV^2 More...

float	thetaX () const
	vertical scattering angle, in rad More...

float	thetaXError () const
	uncertainty on fitted momentum horizontal angle opening More...

float	thetaY () const
	horizontal scattering angle, in rad More...

float	thetaYError () const
	uncertainty on fitted momentum vertical angle opening More...

float	time () const
	time of proton arrival at forward stations More...

float	timeError () const
	uncertainty on time of proton arrival at forward stations More...

bool	validFit () const
	flag for the fit validity More...

const Point &	vertex () const
	fitted vertex position More...

float	vx () const
	fitted vertex horizontal position More...

float	vxError () const
	uncertainty on fitted vertex horizontal position More...

float	vy () const
	fitted vertex vertical position More...

float	vyError () const
	uncertainty on fitted vertex vertical position More...

float	vz () const
	vertex longitudinal position (conventionally set to 0) More...

float	xi () const
	longitudinal fractional momentum loss More...

float	xiError () const
	uncertainty on longitudinal fractional momentum loss More...

Static Public Member Functions
static float	calculateT (double beam_mom, double proton_mom, double theta)
	compute the squared four-momentum transfer from incident and scattered momenta, and angular information More...

static float	mass ()
	proton mass in GeV More...

Private Attributes
float	chi2_
	chi-squared More...

CTPPSLocalTrackLiteRefVector	contributing_local_tracks_
	collection of references to tracks contributing to this object definition More...

CovarianceMatrix	covariance_
	5x5 covariance matrix More...

ReconstructionMethod	method_
	type of reconstruction applied More...

Vector	momentum_
	reconstructed momentum vector More...

unsigned int	ndof_
	number of degrees of freedom More...

float	time_
	reconstructed time at forward detectors More...

float	time_err_
	uncertainty on reconstructed time at forward detectors More...

bool	valid_fit_
	fit validity flag More...

Point	vertex_
	reconstructed vertex position at z/s = 0 More...

float	xi_
	fractional momentum loss (positive for diffractive protons) More...

Static Private Attributes
static constexpr float	mass_ = 0.938272046
	proton mass, GeV More...

static constexpr float	massSquared_ = 0.88035443
	proton mass squared, GeV^2 More...

Detailed Description

Definition at line 21 of file ForwardProton.h.

Member Typedef Documentation

◆ CovarianceMatrix

typedef math::ErrorF<dimension>::type reco::ForwardProton::CovarianceMatrix

dimension-parameter covariance matrix

Definition at line 31 of file ForwardProton.h.

◆ Point

typedef math::XYZPointF reco::ForwardProton::Point

point in the space

Definition at line 35 of file ForwardProton.h.

◆ Vector

typedef math::XYZVectorF reco::ForwardProton::Vector

spatial vector

Definition at line 33 of file ForwardProton.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

parameter dimension

Enumerator
dimension

Definition at line 27 of file ForwardProton.h.

◆ Index

enum reco::ForwardProton::Index

strong

indices to the covariance matrix

Enumerator
xi
th_x
vtx_x
th_y
vtx_y
num_indices

Definition at line 29 of file ForwardProton.h.

◆ LHCSector

enum reco::ForwardProton::LHCSector

strong

LHC sector.

Enumerator
invalid
sector45
sector56

Definition at line 139 of file ForwardProton.h.

◆ ReconstructionMethod

enum reco::ForwardProton::ReconstructionMethod

strong

type of reconstruction applied for this track

Enumerator
invalid
singleRP
multiRP

Definition at line 37 of file ForwardProton.h.

51 { return vertex_; }

Constructor & Destructor Documentation

◆ ForwardProton() [1/2]

ForwardProton::ForwardProton ( )

default constructor

Definition at line 12 of file ForwardProton.cc.

     : time_(0.),
       time_err_(0.),
       xi_(0.),
       chi2_(0.),
       ndof_(0),
       valid_fit_(false),

◆ ForwardProton() [2/2]

ForwardProton::ForwardProton	(	double	chi2,
		double	ndof,
		const Point &	vtx,
		const Vector &	momentum,
		float	xi,
		const CovarianceMatrix &	cov,
		ReconstructionMethod	method,
		const CTPPSLocalTrackLiteRefVector &	local_tracks,
		bool	valid,
		const float	time = `0.`,
		const float	time_err = `0.`
	)

constructor from refit parameters, fitted vertex and momentum, and longitudinal fractional momentum loss

Definition at line 21 of file ForwardProton.cc.

     : vertex_(vtx),
       momentum_(momentum),
       time_(time),
       time_err_(time_err),
       xi_(xi),
       covariance_(cov),
       chi2_(chi2),
       ndof_(ndof),
       valid_fit_(valid),
       method_(method),

Member Function Documentation

◆ calculateT()

float ForwardProton::calculateT	(	double	beam_mom,
		double	proton_mom,
		double	theta
	)

static

compute the squared four-momentum transfer from incident and scattered momenta, and angular information

Definition at line 44 of file ForwardProton.cc.

                                                                                 {
   const double t0 = 2. * (massSquared_ + beam_mom * proton_mom -
                           sqrt((massSquared_ + beam_mom * beam_mom) * (massSquared_ + proton_mom * proton_mom)));
   const double S = sin(theta / 2.);
   return t0 - 4. * beam_mom * proton_mom * S * S;

References massSquared_, S(), funct::sin(), mathSSE::sqrt(), FrontierCondition_GT_autoExpress_cfi::t0, and theta().

Referenced by CTPPSProtonReconstructionSimulationValidator::fillPlots(), and t().

◆ chi2()

float reco::ForwardProton::chi2 ( void ) const

inline

chi-squared of the fit

Definition at line 77 of file ForwardProton.h.

79 { return xi_; }

◆ contributingLocalTracks()

const CTPPSLocalTrackLiteRefVector& reco::ForwardProton::contributingLocalTracks ( ) const

inline

list of RP tracks that contributed to this global track

Definition at line 136 of file ForwardProton.h.

Referenced by CTPPSProtonReconstructionPlotter::analyze(), CTPPSProtonReconstructionSimulationValidator::analyze(), and CTPPSProtonReconstructionPlotter::CalculateTimingTrackingDistance().

◆ error()

double reco::ForwardProton::error ( Index i ) const

inline

return the uncertainty on a given component

Definition at line 93 of file ForwardProton.h.

93 { return error(Index::th_x); }

References error(), and th_x.

Referenced by Page1Parser.Page1Parser::check_for_whole_start_tag(), and error().

◆ lhcSector()

LHCSector reco::ForwardProton::lhcSector ( ) const

inline

Definition at line 140 of file ForwardProton.h.

          :
     static constexpr float mass_ = 0.938272046;        
     static constexpr float massSquared_ = 0.88035443;  
  

Referenced by CTPPSProtonReconstructionSimulationValidator::analyze().

◆ mass()

static float reco::ForwardProton::mass ( )

inlinestatic

proton mass in GeV

Definition at line 107 of file ForwardProton.h.

111 { return time_; }

Referenced by Particle.Particle::__str__(), and DiObject.DiMuon::__str__().

◆ method()

ReconstructionMethod reco::ForwardProton::method ( ) const

inline

reconstruction method for this track

Definition at line 131 of file ForwardProton.h.

131 { return contributing_local_tracks_; }

References contributing_local_tracks_.

Referenced by CTPPSProtonReconstructionSimulationValidator::analyze(), and setValidFit().

◆ momentum()

const Vector& reco::ForwardProton::momentum ( ) const

inline

fitted track direction

Definition at line 64 of file ForwardProton.h.

65 { return momentum_.x(); }

◆ ndof()

unsigned int reco::ForwardProton::ndof ( ) const

inline

number of degrees of freedom for the track fit

Definition at line 79 of file ForwardProton.h.

79 { return xi_; }

References xi_.

◆ normalizedChi2()

float reco::ForwardProton::normalizedChi2 ( ) const

inline

chi-squared divided by ndof (or chi-squared * 1e6 if ndof is zero)

Definition at line 81 of file ForwardProton.h.

81 { return px() / p(); }

References p(), and px().

◆ p()

float reco::ForwardProton::p ( ) const

inline

scalar norm of fitted track momentum

Definition at line 66 of file ForwardProton.h.

67 { return momentum_.y(); }

Referenced by normalizedChi2(), Electron.Electron::ptErr(), and t().

◆ pt()

float reco::ForwardProton::pt ( ) const

inline

scalar fitted track transverse momentum

Definition at line 68 of file ForwardProton.h.

69 { return momentum_.z(); }

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

◆ px()

float reco::ForwardProton::px ( ) const

inline

fitted track momentum horizontal component

Definition at line 70 of file ForwardProton.h.

72 { return chi2_; }

Referenced by normalizedChi2().

◆ py()

float reco::ForwardProton::py ( ) const

inline

fitted track momentum vertical component

Definition at line 72 of file ForwardProton.h.

72 { return chi2_; }

References chi2_.

◆ pz()

float reco::ForwardProton::pz ( ) const

inline

fitted track momentum longitudinal component

Definition at line 74 of file ForwardProton.h.

74 { return ndof_; }

References ndof_.

◆ setContributingLocalTracks()

void reco::ForwardProton::setContributingLocalTracks ( const CTPPSLocalTrackLiteRefVector & v )

inline

store the list of RP tracks that contributed to this global track

Definition at line 134 of file ForwardProton.h.

134 { invalid = -1, sector45, sector56 };

References invalid, sector45, and sector56.

Referenced by CTPPSProtonProducer::produce().

◆ setMethod()

void reco::ForwardProton::setMethod ( const ReconstructionMethod & method )

inline

set the reconstruction method for this track

Definition at line 129 of file ForwardProton.h.

129 { contributing_local_tracks_ = v; }

References contributing_local_tracks_, and findQualityFiles::v.

◆ setTime()

void reco::ForwardProton::setTime ( float time )

inline

Definition at line 117 of file ForwardProton.h.

119 { valid_fit_ = valid; }

References RunInfoPI::valid, and valid_fit_.

Referenced by CTPPSProtonProducer::produce().

◆ setTimeError()

void reco::ForwardProton::setTimeError ( float time_err )

inline

Definition at line 121 of file ForwardProton.h.

121 { return valid_fit_; }

References valid_fit_.

Referenced by CTPPSProtonProducer::produce().

◆ setValidFit()

void reco::ForwardProton::setValidFit ( bool valid = true )

inline

set the flag for the fit validity

Definition at line 124 of file ForwardProton.h.

124 { method_ = method; }

References method(), and method_.

◆ t()

float ForwardProton::t ( ) const

four-momentum transfer squared, in GeV^2

Definition at line 51 of file ForwardProton.cc.

                              {
   const double beam_mom = p() / (1. - xi());
   const double theta = std::hypot(thetaX(), thetaY());
   return calculateT(beam_mom, p(), theta);

References calculateT(), p(), theta(), thetaX(), thetaY(), and xi().

Referenced by CTPPSProtonReconstructionSimulationValidator::fillPlots().

◆ thetaX()

float reco::ForwardProton::thetaX ( ) const

inline

vertical scattering angle, in rad

Definition at line 86 of file ForwardProton.h.

88 { return sqrt(covariance_((unsigned int)i, (unsigned int)i)); }

Referenced by CTPPSProtonReconstructionSimulationValidator::fillPlots(), and t().

◆ thetaXError()

float reco::ForwardProton::thetaXError ( ) const

inline

uncertainty on fitted momentum horizontal angle opening

Definition at line 98 of file ForwardProton.h.

99 { return error(Index::vtx_y); }

◆ thetaY()

float reco::ForwardProton::thetaY ( ) const

inline

horizontal scattering angle, in rad

Definition at line 88 of file ForwardProton.h.

88 { return sqrt(covariance_((unsigned int)i, (unsigned int)i)); }

References covariance_, mps_fire::i, and mathSSE::sqrt().

Referenced by CTPPSProtonReconstructionPlotter::SingleMultiCorrelationPlots::fill(), CTPPSProtonReconstructionSimulationValidator::fillPlots(), and t().

◆ thetaYError()

float reco::ForwardProton::thetaYError ( ) const

inline

uncertainty on fitted momentum vertical angle opening

Definition at line 100 of file ForwardProton.h.

102 { return mass_; }

◆ time()

float reco::ForwardProton::time ( ) const

inline

time of proton arrival at forward stations

Definition at line 116 of file ForwardProton.h.

116 { time_err_ = time_err; }

References time_err_.

Referenced by CTPPSProtonReconstructionSimulationValidator::analyze().

◆ timeError()

float reco::ForwardProton::timeError ( ) const

inline

uncertainty on time of proton arrival at forward stations

Definition at line 120 of file ForwardProton.h.

121 { return valid_fit_; }

Referenced by CTPPSProtonReconstructionSimulationValidator::analyze().

◆ validFit()

bool reco::ForwardProton::validFit ( ) const

inline

flag for the fit validity

Definition at line 126 of file ForwardProton.h.

126 { return method_; }

References method_.

Referenced by CTPPSProtonReconstructionSimulationValidator::analyze().

◆ vertex()

const Point& reco::ForwardProton::vertex ( ) const

inline

fitted vertex position

Definition at line 56 of file ForwardProton.h.

57 { return vertex_.z(); }

Referenced by Tau.Tau::dxy(), and CTPPSProtonReconstructionSimulationValidator::fillPlots().

◆ vx()

float reco::ForwardProton::vx ( ) const

inline

fitted vertex horizontal position

Definition at line 58 of file ForwardProton.h.

59 { return momentum_; }

◆ vxError()

float reco::ForwardProton::vxError ( ) const

inline

uncertainty on fitted vertex horizontal position

Definition at line 102 of file ForwardProton.h.

102 { return mass_; }

References mass_.

◆ vy()

float reco::ForwardProton::vy ( ) const

inline

fitted vertex vertical position

Definition at line 60 of file ForwardProton.h.

61 { return momentum_.r(); }

◆ vyError()

float reco::ForwardProton::vyError ( ) const

inline

uncertainty on fitted vertex vertical position

Definition at line 104 of file ForwardProton.h.

111 { return time_; }

◆ vz()

float reco::ForwardProton::vz ( ) const

inline

vertex longitudinal position (conventionally set to 0)

Definition at line 62 of file ForwardProton.h.

63 { return momentum_.rho(); }

◆ xi()

float reco::ForwardProton::xi ( ) const

inline

longitudinal fractional momentum loss

Definition at line 84 of file ForwardProton.h.

88 { return sqrt(covariance_((unsigned int)i, (unsigned int)i)); }

Referenced by CTPPSProtonReconstructionPlotter::SingleMultiCorrelationPlots::fill(), CTPPSProtonReconstructionSimulationValidator::fillPlots(), and t().

◆ xiError()

float reco::ForwardProton::xiError ( ) const

inline

uncertainty on longitudinal fractional momentum loss

Definition at line 96 of file ForwardProton.h.

97 { return error(Index::vtx_x); }

Member Data Documentation

◆ chi2_

float reco::ForwardProton::chi2_

private

chi-squared

Definition at line 165 of file ForwardProton.h.

Referenced by py().

◆ contributing_local_tracks_

CTPPSLocalTrackLiteRefVector reco::ForwardProton::contributing_local_tracks_

private

collection of references to tracks contributing to this object definition

Definition at line 173 of file ForwardProton.h.

Referenced by method(), and setMethod().

◆ covariance_

CovarianceMatrix reco::ForwardProton::covariance_

private

5x5 covariance matrix

Definition at line 163 of file ForwardProton.h.

Referenced by thetaY().

◆ mass_

constexpr float reco::ForwardProton::mass_ = 0.938272046

staticconstexprprivate

proton mass, GeV

Definition at line 149 of file ForwardProton.h.

Referenced by vxError().

◆ massSquared_

constexpr float reco::ForwardProton::massSquared_ = 0.88035443

staticconstexprprivate

proton mass squared, GeV^2

Definition at line 150 of file ForwardProton.h.

Referenced by calculateT().

◆ method_

ReconstructionMethod reco::ForwardProton::method_

private

type of reconstruction applied

Definition at line 171 of file ForwardProton.h.

Referenced by setValidFit(), and validFit().

◆ momentum_

Vector reco::ForwardProton::momentum_

private

reconstructed momentum vector

Definition at line 155 of file ForwardProton.h.

◆ ndof_

unsigned int reco::ForwardProton::ndof_

private

number of degrees of freedom

Definition at line 167 of file ForwardProton.h.

Referenced by pz().

◆ time_

float reco::ForwardProton::time_

private

reconstructed time at forward detectors

Definition at line 157 of file ForwardProton.h.

◆ time_err_

float reco::ForwardProton::time_err_

private

uncertainty on reconstructed time at forward detectors

Definition at line 159 of file ForwardProton.h.

Referenced by time().

◆ valid_fit_

bool reco::ForwardProton::valid_fit_

private

fit validity flag

Definition at line 169 of file ForwardProton.h.

Referenced by setTime(), and setTimeError().

◆ vertex_

Point reco::ForwardProton::vertex_

private

reconstructed vertex position at z/s = 0

Definition at line 153 of file ForwardProton.h.

◆ xi_

float reco::ForwardProton::xi_

private

fractional momentum loss (positive for diffractive protons)

Definition at line 161 of file ForwardProton.h.

Referenced by ndof().

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ CovarianceMatrix

◆ Point

◆ Vector

Member Enumeration Documentation

◆ anonymous enum

◆ Index

◆ LHCSector

◆ ReconstructionMethod

Constructor & Destructor Documentation

◆ ForwardProton() [1/2]

◆ ForwardProton() [2/2]

Member Function Documentation

◆ calculateT()

◆ chi2()

◆ contributingLocalTracks()

◆ error()

◆ lhcSector()

◆ mass()

◆ method()

◆ momentum()

◆ ndof()

◆ normalizedChi2()

◆ p()

◆ pt()

◆ px()

◆ py()

◆ pz()

◆ setContributingLocalTracks()

◆ setMethod()

◆ setTime()

◆ setTimeError()

◆ setValidFit()

◆ t()

◆ thetaX()

◆ thetaXError()

◆ thetaY()

◆ thetaYError()

◆ time()

◆ timeError()

◆ validFit()

◆ vertex()

◆ vx()

◆ vxError()

◆ vy()

◆ vyError()

◆ vz()

◆ xi()

◆ xiError()

Member Data Documentation

◆ chi2_

◆ contributing_local_tracks_

◆ covariance_

◆ mass_

◆ massSquared_

◆ method_

◆ momentum_

◆ ndof_

◆ time_

◆ time_err_

◆ valid_fit_

◆ vertex_

◆ xi_