#include <BeamSpot.h>

Public Types
enum	{ dimension = 7 }
enum	{ dim3 = 3 }
enum	{ resdim = 2 }
enum	BeamType { Unknown = -1, Fake = 0, LHC = 1, Tracker = 2 }
	beam spot flags More...
typedef math::Error< dim3 >::type	Covariance3DMatrix
typedef math::Error< dimension > ::type	CovarianceMatrix
typedef math::XYZPoint	Point
	point in the space
typedef math::Error< resdim >::type	ResCovMatrix
Public Member Functions
	BeamSpot ()
	default constructor
	BeamSpot (const Point &point, double sigmaZ, double dxdz, double dydz, double BeamWidthX, const CovarianceMatrix &error, BeamType type=Unknown)
	constructor from values
double	beamWidthFromBeta (double z, double e) const
double	BeamWidthX () const
	beam width X
double	BeamWidthXError () const
	error on beam width X, assume error in X = Y
double	BeamWidthY () const
	beam width Y
double	BeamWidthYError () const
	error on beam width Y, assume error in X = Y
double	betaStar () const
double	covariance (int i, int j) const
	(i,j)-th element of error matrix
CovarianceMatrix	covariance () const
	return full covariance matrix of dim 7
Covariance3DMatrix	covariance3D () const
	return only 3D position covariance matrix
double	dxdz () const
	dxdz slope
double	dxdzError () const
	error on dxdz
double	dydz () const
	dydz slope
double	dydzError () const
	error on dydz
double	emittanceX () const
	additional information
double	emittanceY () const
const Point	position (const double z) const
	position of the beam spot at a given z value (it takes into account the dxdz and dydz slopes)
const Point &	position () const
	position
void	print (std::stringstream &ss) const
	print information
Covariance3DMatrix	rotatedCovariance3D () const
void	setBeamWidthX (double v)
void	setBeamWidthY (double v)
void	setbetaStar (double v)
void	setEmittanceX (double v)
void	setEmittanceY (double v)
void	setType (BeamType type)
	set beam type
double	sigmaZ () const
	sigma z
double	sigmaZ0Error () const
	error on sigma z
BeamType	type () const
	return beam type
double	x (const double z) const
	x coordinate of the beeam spot position at a given z value (it takes into account the dxdz slope)
double	x0 () const
	x coordinate
double	x0Error () const
	error on x
double	y (const double z) const
	y coordinate of the beeam spot position at a given z value (it takes into account the dydz slope)
double	y0 () const
	y coordinate
double	y0Error () const
	error on y
double	z0 () const
	z coordinate
double	z0Error () const
	error on z
Private Attributes
Double32_t	BeamWidthX_
Double32_t	BeamWidthY_
Double32_t	betaStar_
Double32_t	dxdz_
Double32_t	dydz_
Double32_t	emittanceX_
Double32_t	emittanceY_
CovarianceMatrix	error_
	errors
Point	position_
	position
Double32_t	sigmaZ_
BeamType	type_

Detailed Description

Reconstructed beam spot object which provides position, error, and width of the beam position.

Author:: Francisco Yumiceva, Fermilab (yumiceva@fnal.gov)

Version:

Id:: BeamSpot.h,v 1.12 2011/04/06 20:31:35 venturia Exp

Definition at line 23 of file BeamSpot.h.

Member Typedef Documentation

typedef math::Error<dim3>::type reco::BeamSpot::Covariance3DMatrix

Definition at line 34 of file BeamSpot.h.

typedef math::Error<dimension>::type reco::BeamSpot::CovarianceMatrix

Definition at line 32 of file BeamSpot.h.

typedef math::XYZPoint reco::BeamSpot::Point

point in the space

Definition at line 30 of file BeamSpot.h.

typedef math::Error<resdim>::type reco::BeamSpot::ResCovMatrix

Definition at line 36 of file BeamSpot.h.

Member Enumeration Documentation

anonymous enum

Enumerator:

dimension

Definition at line 31 of file BeamSpot.h.

{ dimension = 7 };

anonymous enum

Enumerator:

dim3

Definition at line 33 of file BeamSpot.h.

{ dim3 = 3 };

anonymous enum

Enumerator:

resdim

Definition at line 35 of file BeamSpot.h.

{ resdim = 2 };

enum reco::BeamSpot::BeamType

beam spot flags

Enumerator:

Unknown
Fake
LHC
Tracker

Definition at line 27 of file BeamSpot.h.

{ Unknown=-1, Fake=0, LHC=1, Tracker=2 };

Constructor & Destructor Documentation

reco::BeamSpot::BeamSpot ( )

default constructor

Definition at line 25 of file BeamSpot.cc.

References j, gen::k, and Unknown.

                     {
          // initialize
          position_ = Point(0.,0.,0.);
          sigmaZ_ = 0.;
          dxdz_ = 0.;
          dydz_ = 0.;
          BeamWidthX_ = 0.; BeamWidthY_ = 0;
          for (int j=0; j<7; j++) {
                        for (int k=j; k<7; k++) {
                                error_(j,k) = 0.;
                        }
          }
          type_ = Unknown;
          emittanceX_ = 0;
          emittanceY_ = 0;
          betaStar_ = 0;
  }

reco::BeamSpot::BeamSpot	(	const Point &	point,
		double	sigmaZ,
		double	dxdz,
		double	dydz,
		double	BeamWidthX,
		const CovarianceMatrix &	error,
		BeamType	type = `Unknown`
	)		`[inline]`

constructor from values

Definition at line 42 of file BeamSpot.h.

References BeamWidthX(), BeamWidthX_, BeamWidthY_, betaStar_, dxdz(), dxdz_, dydz(), dydz_, emittanceX_, emittanceY_, error, error_, point, position_, sigmaZ(), sigmaZ_, type(), and type_.

                                           { 
      position_ = point;
      sigmaZ_ = sigmaZ;
      dxdz_ = dxdz;
      dydz_ = dydz;
      BeamWidthX_ = BeamWidthX;
      BeamWidthY_ = BeamWidthX;
      error_ = error;
          type_ = type;
          emittanceX_ = emittanceY_ = 0;
          betaStar_ = 0;
    };

Member Function Documentation

double reco::BeamSpot::beamWidthFromBeta	(	double	z,
		double	e
	)		const `[inline]`

Definition at line 140 of file BeamSpot.h.

References betaStar_, position_, funct::pow(), and mathSSE::sqrt().

                                                             {
                return sqrt( e*betaStar_*(1 + pow((z-position_.Z())/betaStar_,2) ) );
        }

double reco::BeamSpot::BeamWidthX ( ) const [inline]

beam width X

Definition at line 87 of file BeamSpot.h.

References BeamWidthX_.

Referenced by PrimaryVertexAnalyzer4PU::analyze(), PrimaryVertexMonitor::analyze(), PrimaryVertexAnalyzer4PU::analyzeVertexCollection(), BeamSpot(), BeamSpotTransientTrackingRecHit::BeamSpotTransientTrackingRecHit(), BeamFitter::dumpBWTxtFile(), BeamFitter::dumpTxtFile(), TwoBodyDecayVirtualMeasurement::extractBeamSpotError(), MuonVPlusJetsIDSelectionFunctor::fall10Cuts(), BeamMonitorBx::FillTables(), BeamMonitorBx::FillTrendHistos(), BSFitter::Fit(), BeamMonitor::FitAndFill(), MuonVPlusJetsIDSelectionFunctor::kitQCDCuts(), KFBasedPixelFitter::MyBeamSpotHit::MyBeamSpotHit(), PrimaryVertexAnalyzer4PU::printPVTrks(), HIBestVertexProducer::produce(), BeamFitter::runFitterNoTxt(), BeamFitter::runPVandTrkFitter(), HIProtoTrackSelector::select(), MuonVPlusJetsIDSelectionFunctor::spring10Cuts(), BeamMonitorBx::weight(), and BeamFitter::write2DB().

{ return BeamWidthX_; }

double reco::BeamSpot::BeamWidthXError ( ) const [inline]

error on beam width X, assume error in X = Y

Definition at line 104 of file BeamSpot.h.

References error_, and mathSSE::sqrt().

Referenced by BeamFitter::dumpBWTxtFile(), BeamMonitorBx::FillTables(), BeamMonitorBx::FillTrendHistos(), BeamMonitor::FitAndFill(), BeamFitter::runFitterNoTxt(), and BeamMonitorBx::weight().

{ return sqrt ( error_(6,6) );}

double reco::BeamSpot::BeamWidthY ( ) const [inline]

beam width Y

Definition at line 89 of file BeamSpot.h.

References BeamWidthY_.

Referenced by PrimaryVertexAnalyzer4PU::analyze(), PrimaryVertexMonitor::analyze(), PrimaryVertexAnalyzer4PU::analyzeVertexCollection(), BeamSpotTransientTrackingRecHit::BeamSpotTransientTrackingRecHit(), BeamFitter::dumpTxtFile(), TwoBodyDecayVirtualMeasurement::extractBeamSpotError(), MuonVPlusJetsIDSelectionFunctor::fall10Cuts(), BeamMonitorBx::FillTables(), BeamMonitorBx::FillTrendHistos(), BeamMonitor::FitAndFill(), MuonVPlusJetsIDSelectionFunctor::kitQCDCuts(), PrimaryVertexAnalyzer4PU::printPVTrks(), HIBestVertexProducer::produce(), BeamFitter::runFitterNoTxt(), BeamFitter::runPVandTrkFitter(), HIProtoTrackSelector::select(), MuonVPlusJetsIDSelectionFunctor::spring10Cuts(), BeamMonitorBx::weight(), and BeamFitter::write2DB().

{ return BeamWidthY_; }

double reco::BeamSpot::BeamWidthYError ( ) const [inline]

error on beam width Y, assume error in X = Y

Definition at line 106 of file BeamSpot.h.

References error_, and mathSSE::sqrt().

Referenced by BeamMonitorBx::FillTables(), BeamMonitorBx::FillTrendHistos(), BeamMonitor::FitAndFill(), BeamFitter::runFitterNoTxt(), and BeamMonitorBx::weight().

{ return sqrt ( error_(6,6) );}

double reco::BeamSpot::betaStar ( ) const [inline]

Definition at line 139 of file BeamSpot.h.

References betaStar_.

Referenced by BeamFitter::dumpTxtFile().

{ return betaStar_; }

double reco::BeamSpot::covariance	(	int	i,
		int	j
	)		const `[inline]`

(i,j)-th element of error matrix

Definition at line 113 of file BeamSpot.h.

References error_.

Referenced by BeamFitter::dumpTxtFile(), BSFitter::Fit(), BeamFitter::runPVandTrkFitter(), and BeamFitter::write2DB().

                                           {
      return error_(i,j);
    }

CovarianceMatrix reco::BeamSpot::covariance ( void ) const [inline]

return full covariance matrix of dim 7

Definition at line 117 of file BeamSpot.h.

References error_.

{ return error_; }

Covariance3DMatrix reco::BeamSpot::covariance3D ( ) const [inline]

return only 3D position covariance matrix

Definition at line 119 of file BeamSpot.h.

References error_, j, gen::k, and makeMuonMisalignmentScenario::matrix.

Referenced by pat::PATElectronProducer::embedHighLevel().

                                            {

                Covariance3DMatrix matrix;
                for (int j=0; j<3; j++) {
                        for (int k=j; k<3; k++) {
                                matrix(j,k) = error_(j,k);
                        }
                }
      return matrix;
    };

double reco::BeamSpot::dxdz ( ) const [inline]

double reco::BeamSpot::dxdzError ( ) const [inline]

error on dxdz

Definition at line 99 of file BeamSpot.h.

References error_, and mathSSE::sqrt().

Referenced by BeamMonitor::FitAndFill(), BeamFitter::runFitterNoTxt(), and BeamMonitorBx::weight().

{ return sqrt ( error_(4,4) ); }

double reco::BeamSpot::dydz ( ) const [inline]

double reco::BeamSpot::dydzError ( ) const [inline]

error on dydz

Definition at line 101 of file BeamSpot.h.

References error_, and mathSSE::sqrt().

Referenced by BeamMonitor::FitAndFill(), BeamFitter::runFitterNoTxt(), and BeamMonitorBx::weight().

{ return sqrt ( error_(5,5) ); }

double reco::BeamSpot::emittanceX ( ) const [inline]

additional information

Definition at line 137 of file BeamSpot.h.

References emittanceX_.

Referenced by BeamFitter::dumpTxtFile().

{ return emittanceX_; }

double reco::BeamSpot::emittanceY ( ) const [inline]

Definition at line 138 of file BeamSpot.h.

References emittanceY_.

Referenced by BeamFitter::dumpTxtFile().

{ return emittanceY_; }

const Point& reco::BeamSpot::position ( ) const [inline]

position

Definition at line 63 of file BeamSpot.h.

References position_.

{ return position_; }

const BeamSpot::Point reco::BeamSpot::position ( const double z ) const

position of the beam spot at a given z value (it takes into account the dxdz and dydz slopes)

Definition at line 43 of file BeamSpot.cc.

References pos, x, and detailsBasic3DVector::y.

                                                             {

    Point pos(x(z),y(z),z);
    return pos;

  }

void reco::BeamSpot::print ( std::stringstream & ss ) const

print information

Definition at line 51 of file BeamSpot.cc.

Referenced by reco::operator<<().

                                                {

    ss << "-----------------------------------------------------\n"
       << "              Beam Spot Data\n\n"
           << " Beam type    = " << type() << "\n"
       << "       X0     = " << x0() << " +/- " << x0Error() << " [cm]\n"
       << "       Y0     = " << y0() << " +/- " << y0Error() << " [cm]\n"
       << "       Z0     = " << z0() << " +/- " << z0Error() << " [cm]\n"
       << " Sigma Z0     = " << sigmaZ() << " +/- " << sigmaZ0Error() << " [cm]\n"
       << " dxdz         = " << dxdz() << " +/- " << dxdzError() << " [radians]\n"
       << " dydz         = " << dydz() << " +/- " << dydzError() << " [radians]\n"
       << " Beam width X = " << BeamWidthX() << " +/- " << BeamWidthXError() << " [cm]\n"
           << " Beam width Y = " << BeamWidthY() << " +/- " << BeamWidthYError() << " [cm]\n"
           << " EmittanceX   = " << emittanceX() << " [cm]\n"
           << " EmittanceY   = " << emittanceY() << " [cm]\n"
           << " beta-star    = " << betaStar() << " [cm]\n"
       << "-----------------------------------------------------\n\n";
        
  }

BeamSpot::Covariance3DMatrix reco::BeamSpot::rotatedCovariance3D ( ) const

Definition at line 79 of file BeamSpot.cc.

References a, asHepVector(), makeMuonMisalignmentScenario::matrix, funct::pow(), idealTransformation::rotation, TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), TkRotation< T >::yx(), TkRotation< T >::yy(), TkRotation< T >::yz(), TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

Referenced by GsfConstraintAtVertex::constrainAtBeamSpot(), pat::PATMuonProducer::embedHighLevel(), PixelVertexProducer::produce(), PrimaryVertexProducer::produce(), and TrajectoryStateClosestToBeamLine::transverseImpactParameter().

  {
      AlgebraicVector3 newZ(dxdz(), dydz(), 1.);
      AlgebraicVector3 globalZ(0.,0.,1.);
      AlgebraicVector3 rotationAxis = ROOT::Math::Cross(globalZ.Unit(), newZ.Unit());
      float rotationAngle = -acos( ROOT::Math::Dot(globalZ.Unit(),newZ.Unit()));
      AlgebraicVector a = asHepVector(rotationAxis);
      Basic3DVector<float> aa(a[0], a[1], a[2]);
      TkRotation<float> rotation(aa ,rotationAngle);
      AlgebraicMatrix33 rotationMatrix;
      rotationMatrix(0,0) = rotation.xx();
      rotationMatrix(0,1) = rotation.xy();
      rotationMatrix(0,2) = rotation.xz();
      rotationMatrix(1,0) = rotation.yx();
      rotationMatrix(1,1) = rotation.yy();
      rotationMatrix(1,2) = rotation.yz();
      rotationMatrix(2,0) = rotation.zx();
      rotationMatrix(2,1) = rotation.zy();
      rotationMatrix(2,2) = rotation.zz();

      AlgebraicSymMatrix33 diagError ;
      diagError(0,0) = pow(BeamWidthX(),2);
      diagError(1,1) = pow(BeamWidthY(),2);
      diagError(2,2) = pow(sigmaZ(),2);

      Covariance3DMatrix matrix;
      matrix = ROOT::Math::Similarity(rotationMatrix, diagError) + covariance3D();
      return matrix;
  }