reco::ForwardProton Class Reference

#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 19 of file ForwardProton.h.

Member Typedef Documentation

CovarianceMatrix

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

dimension-parameter covariance matrix

Definition at line 26 of file ForwardProton.h.

Point

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

point in the space

Definition at line 30 of file ForwardProton.h.

Vector

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

spatial vector

Definition at line 28 of file ForwardProton.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

parameter dimension

Enumerator
dimension

Definition at line 22 of file ForwardProton.h.

{ dimension = 5 };

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 24 of file ForwardProton.h.

{ xi, th_x, vtx_x, th_y, vtx_y, num_indices = dimension };

LHCSector

enum reco::ForwardProton::LHCSector

strong

LHC sector.

Enumerator
invalid
sector45
sector56

Definition at line 134 of file ForwardProton.h.

{ invalid = -1, sector45, sector56 };

ReconstructionMethod

enum reco::ForwardProton::ReconstructionMethod

strong

type of reconstruction applied for this track

Enumerator
invalid
singleRP
multiRP

Definition at line 32 of file ForwardProton.h.

{ invalid = -1, singleRP, multiRP };

Constructor & Destructor Documentation

ForwardProton() [1/2]

ForwardProton::ForwardProton

(

)

default constructor

Definition at line 13 of file ForwardProton.cc.

    : time_(0.),
      time_err_(0.),
      xi_(0.),
      chi2_(0.),
      ndof_(0),
      valid_fit_(false),
      method_(ReconstructionMethod::invalid) {}

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 22 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),
      contributing_local_tracks_(local_tracks) {}

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 45 of file ForwardProton.cc.

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

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

                                                                                {
  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;
}

chi2()

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

inline

chi-squared of the fit

Definition at line 72 of file ForwardProton.h.

References chi2_.

{ return chi2_; }

contributingLocalTracks()

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

inline

list of RP tracks that contributed to this global track

Definition at line 131 of file ForwardProton.h.

References contributing_local_tracks_.

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

{ return contributing_local_tracks_; }

error()

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

inline

return the uncertainty on a given component

Definition at line 88 of file ForwardProton.h.

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

Referenced by Page1Parser.Page1Parser::check_for_whole_start_tag(), thetaXError(), thetaYError(), vxError(), vyError(), and xiError().

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

lhcSector()

LHCSector reco::ForwardProton::lhcSector ( ) const

inline

Definition at line 135 of file ForwardProton.h.

References invalid, pz(), sector45, and sector56.

                                {
      if (pz() < 0.)
        return LHCSector::sector56;
      if (pz() > 0.)
        return LHCSector::sector45;
      return LHCSector::invalid;
    }

mass()

static float reco::ForwardProton::mass ( )

inlinestatic

proton mass in GeV

Definition at line 102 of file ForwardProton.h.

References mass_.

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

{ return mass_; }

method()

ReconstructionMethod reco::ForwardProton::method ( ) const

inline

reconstruction method for this track

Definition at line 126 of file ForwardProton.h.

References method_.

Referenced by setMethod().

{ return method_; }

momentum()

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

inline

fitted track direction

Definition at line 59 of file ForwardProton.h.

References momentum_.

{ return momentum_; }

ndof()

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

inline

number of degrees of freedom for the track fit

Definition at line 74 of file ForwardProton.h.

References ndof_.

{ return ndof_; }

normalizedChi2()

float reco::ForwardProton::normalizedChi2 ( ) const

inline

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

Definition at line 76 of file ForwardProton.h.

References chi2_, and ndof_.

{ return (ndof_ != 0) ? chi2_ / ndof_ : chi2_ * 1.e6; }

p()

float reco::ForwardProton::p ( ) const

inline

scalar norm of fitted track momentum

Definition at line 61 of file ForwardProton.h.

References momentum_.

Referenced by Electron.Electron::ptErr(), t(), thetaX(), and thetaY().

{ return momentum_.r(); }

pt()

float reco::ForwardProton::pt ( ) const

inline

scalar fitted track transverse momentum

Definition at line 63 of file ForwardProton.h.

References momentum_.

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(), datamodel.Object::p4(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

{ return momentum_.rho(); }

px()

float reco::ForwardProton::px ( ) const

inline

fitted track momentum horizontal component

Definition at line 65 of file ForwardProton.h.

References momentum_.

Referenced by thetaX().

{ return momentum_.x(); }

py()

float reco::ForwardProton::py ( ) const

inline

fitted track momentum vertical component

Definition at line 67 of file ForwardProton.h.

References momentum_.

Referenced by thetaY().

{ return momentum_.y(); }

pz()

float reco::ForwardProton::pz ( ) const

inline

fitted track momentum longitudinal component

Definition at line 69 of file ForwardProton.h.

References momentum_.

Referenced by lhcSector().

{ return momentum_.z(); }

setContributingLocalTracks()

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

inline

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

Definition at line 129 of file ForwardProton.h.

References contributing_local_tracks_, and findQualityFiles::v.

Referenced by CTPPSProtonProducer::produce().

{ contributing_local_tracks_ = v; }

setMethod()

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

inline

set the reconstruction method for this track

Definition at line 124 of file ForwardProton.h.

References method(), and method_.

{ method_ = method; }

setTime()

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

inline

Definition at line 112 of file ForwardProton.h.

References time(), and time_.

Referenced by CTPPSProtonProducer::produce().

{ time_ = time; }

setTimeError()

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

inline

Definition at line 116 of file ForwardProton.h.

References time_err_.

Referenced by CTPPSProtonProducer::produce().

{ time_err_ = time_err; }

setValidFit()

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

inline

set the flag for the fit validity

Definition at line 119 of file ForwardProton.h.

References validateGeometry_cfg::valid, and valid_fit_.

{ valid_fit_ = valid; }

t()

float ForwardProton::t

(

)

const

four-momentum transfer squared, in GeV^2

Definition at line 52 of file ForwardProton.cc.

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

Referenced by CTPPSProtonReconstructionSimulationValidator::fillPlots().

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

thetaX()

float reco::ForwardProton::thetaX ( ) const

inline

vertical scattering angle, in rad

Definition at line 81 of file ForwardProton.h.

References p(), and px().

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

{ return px() / p(); }

thetaXError()

float reco::ForwardProton::thetaXError ( ) const

inline

uncertainty on fitted momentum horizontal angle opening

Definition at line 93 of file ForwardProton.h.

References error(), and th_x.

{ return error(Index::th_x); }

thetaY()

float reco::ForwardProton::thetaY ( ) const

inline

horizontal scattering angle, in rad

Definition at line 83 of file ForwardProton.h.

References p(), and py().

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

{ return py() / p(); }

thetaYError()

float reco::ForwardProton::thetaYError ( ) const

inline

uncertainty on fitted momentum vertical angle opening

Definition at line 95 of file ForwardProton.h.

References error(), and th_y.

{ return error(Index::th_y); }

time()

float reco::ForwardProton::time ( ) const

inline

time of proton arrival at forward stations

Definition at line 111 of file ForwardProton.h.

References time_.

Referenced by setTime().

{ return time_; }

timeError()

float reco::ForwardProton::timeError ( ) const

inline

uncertainty on time of proton arrival at forward stations

Definition at line 115 of file ForwardProton.h.

References time_err_.

{ return time_err_; }

validFit()

bool reco::ForwardProton::validFit ( ) const

inline

flag for the fit validity

Definition at line 121 of file ForwardProton.h.

References valid_fit_.

Referenced by CTPPSProtonReconstructionPlotter::SingleMultiCorrelationPlots::fill(), and CTPPSProtonReconstructionPlotter::ArmCorrelationPlots::fill().

{ return valid_fit_; }

vertex()

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

inline

fitted vertex position

Definition at line 51 of file ForwardProton.h.

References vertex_.

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

{ return vertex_; }

vx()

float reco::ForwardProton::vx ( ) const

inline

fitted vertex horizontal position

Definition at line 53 of file ForwardProton.h.

References vertex_.

{ return vertex_.x(); }

vxError()

float reco::ForwardProton::vxError ( ) const

inline

uncertainty on fitted vertex horizontal position

Definition at line 97 of file ForwardProton.h.

References error(), and vtx_x.

{ return error(Index::vtx_x); }

vy()

float reco::ForwardProton::vy ( ) const

inline

fitted vertex vertical position

Definition at line 55 of file ForwardProton.h.

References vertex_.

{ return vertex_.y(); }

vyError()

float reco::ForwardProton::vyError ( ) const

inline

uncertainty on fitted vertex vertical position

Definition at line 99 of file ForwardProton.h.

References error(), and vtx_y.

{ return error(Index::vtx_y); }

vz()

float reco::ForwardProton::vz ( ) const

inline

vertex longitudinal position (conventionally set to 0)

Definition at line 57 of file ForwardProton.h.

References vertex_.

{ return vertex_.z(); }

xi()

float reco::ForwardProton::xi ( ) const

inline

longitudinal fractional momentum loss

Definition at line 79 of file ForwardProton.h.

References xi_.

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

{ return xi_; }

xiError()

float reco::ForwardProton::xiError ( ) const

inline

uncertainty on longitudinal fractional momentum loss

Definition at line 91 of file ForwardProton.h.

References error(), and xi.

{ return error(Index::xi); }

Member Data Documentation

chi2_

float reco::ForwardProton::chi2_

private

chi-squared

Definition at line 160 of file ForwardProton.h.

Referenced by chi2(), and normalizedChi2().

contributing_local_tracks_

CTPPSLocalTrackLiteRefVector reco::ForwardProton::contributing_local_tracks_

private

collection of references to tracks contributing to this object definition

Definition at line 168 of file ForwardProton.h.

Referenced by contributingLocalTracks(), and setContributingLocalTracks().

covariance_

CovarianceMatrix reco::ForwardProton::covariance_

private

5x5 covariance matrix

Definition at line 158 of file ForwardProton.h.

Referenced by error().

mass_

constexpr float reco::ForwardProton::mass_ = 0.938272046

staticprivate

proton mass, GeV

Definition at line 144 of file ForwardProton.h.

Referenced by mass().

massSquared_

constexpr float reco::ForwardProton::massSquared_ = 0.88035443

staticprivate

proton mass squared, GeV^2

Definition at line 145 of file ForwardProton.h.

Referenced by calculateT().

method_

ReconstructionMethod reco::ForwardProton::method_

private

type of reconstruction applied

Definition at line 166 of file ForwardProton.h.

Referenced by method(), and setMethod().

momentum_

Vector reco::ForwardProton::momentum_

private

reconstructed momentum vector

Definition at line 150 of file ForwardProton.h.

Referenced by momentum(), p(), pt(), px(), py(), and pz().

ndof_

unsigned int reco::ForwardProton::ndof_

private

number of degrees of freedom

Definition at line 162 of file ForwardProton.h.

Referenced by ndof(), and normalizedChi2().

time_

float reco::ForwardProton::time_

private

reconstructed time at forward detectors

Definition at line 152 of file ForwardProton.h.

Referenced by setTime(), and time().

time_err_

float reco::ForwardProton::time_err_

private

uncertainty on reconstructed time at forward detectors

Definition at line 154 of file ForwardProton.h.

Referenced by setTimeError(), and timeError().

valid_fit_

bool reco::ForwardProton::valid_fit_

private

fit validity flag

Definition at line 164 of file ForwardProton.h.

Referenced by setValidFit(), and validFit().

vertex_

Point reco::ForwardProton::vertex_

private

reconstructed vertex position at z/s = 0

Definition at line 148 of file ForwardProton.h.

Referenced by vertex(), vx(), vy(), and vz().

xi_

float reco::ForwardProton::xi_

private

fractional momentum loss (positive for diffractive protons)

Definition at line 156 of file ForwardProton.h.

Referenced by xi().