#include <MuonUpdatorAtVertex.h>

Public Member Functions
	MuonUpdatorAtVertex (const edm::ParameterSet &pset, const MuonServiceProxy *service)
	Constructor. More...

std::pair< bool, FreeTrajectoryState >	propagate (const TrajectoryStateOnSurface &tsos, const reco::BeamSpot &beamSpot) const
	Propagate the state to the 2D-PCA. More...

std::pair< bool, FreeTrajectoryState >	propagate (const TrajectoryStateOnSurface &tsos) const __attribute__((deprecated))
	Propagate the state to the 2D-PCA (nominal CMS axis) - DEPRECATED -. More...

std::pair< bool, FreeTrajectoryState >	propagateToNominalLine (const TrajectoryStateOnSurface &tsos) const
	Propagate the state to the 2D-PCA (nominal CMS axis) More...

std::pair< bool, FreeTrajectoryState >	propagateWithUpdate (const TrajectoryStateOnSurface &tsos, const reco::BeamSpot &beamSpot) const
	Propagate to the 2D-PCA and apply the vertex constraint. More...

std::pair< bool, FreeTrajectoryState >	update (const reco::TransientTrack &track, const reco::BeamSpot &beamSpot) const
	Applies the vertex constraint. More...

std::pair< bool, FreeTrajectoryState >	update (const FreeTrajectoryState &ftsAtVtx, const reco::BeamSpot &beamSpot) const
	Applies the vertex constraint. More...

virtual	~MuonUpdatorAtVertex ()
	Destructor. More...

Private Attributes
double	theChi2Cut

SingleTrackVertexConstraint	theConstrictor

GlobalError	thePositionErrors

std::string	thePropagatorName

const MuonServiceProxy *	theService

TransientTrackFromFTSFactory	theTransientTrackFactory

Detailed Description

This class do the extrapolation of a TrajectoryStateOnSurface to the PCA and can apply, with a different method, the vertex constraint. The vertex constraint is applyed using the Kalman Filter tools used for the vertex reconstruction.

For the time being the propagator is the SteppingHelixPropagator because the method propagate(TSOS,GlobalPoint) it is in its specific interface. Once the interface of the Propagator base class will be updated, then propagator will become generic.

For what concern the beam spot, it is possible set via cff the relevant parameters:

BeamSpotPosition[0] <=> x BeamSpotPosition[1] <=> y BeamSpotPosition[2] <=> z

BeamSpotPositionErrors[0] = sigma(x) BeamSpotPositionErrors[1] = sigma(y) BeamSpotPositionErrors[2] = sigma(z)

Author: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

This class do the extrapolation of a TrajectoryStateOnSurface to the PCA and can apply, with a different method, the vertex constraint. The vertex constraint is applyed using the Kalman Filter tools used for the vertex reconstruction.

Author: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 42 of file MuonUpdatorAtVertex.h.

Constructor & Destructor Documentation

MuonUpdatorAtVertex::MuonUpdatorAtVertex	(	const edm::ParameterSet &	pset,
		const MuonServiceProxy *	service
	)

Constructor.

Definition at line 28 of file MuonUpdatorAtVertex.cc.

References edm::ParameterSet::getParameter(), theChi2Cut, thePositionErrors, and thePropagatorName.

     : theService(service) {
   thePropagatorName = pset.getParameter<string>("Propagator");
 
   // Errors on the Beam spot position
   vector<double> errors = pset.getParameter<vector<double> >("BeamSpotPositionErrors");
   if (errors.size() != 3)
     edm::LogError("Muon|RecoMuon|MuonUpdatorAtVertex")
         << "MuonUpdatorAtVertex::BeamSpotPositionErrors wrong number of parameters!!";
 
   // assume:
   // errors[0] = sigma(x)
   // errors[1] = sigma(y)
   // errors[2] = sigma(z)
 
   AlgebraicSymMatrix33 mat;
   mat(0, 0) = errors[0] * errors[0];
   mat(1, 1) = errors[1] * errors[1];
   mat(2, 2) = errors[2] * errors[2];
   GlobalError glbErrPos(mat);
 
   thePositionErrors = glbErrPos;
 
   // cut on chi^2
   theChi2Cut = pset.getParameter<double>("MaxChi2");
 }

MuonUpdatorAtVertex::~MuonUpdatorAtVertex ( )

virtual

Destructor.

Definition at line 56 of file MuonUpdatorAtVertex.cc.

56 {}

Member Function Documentation

pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::propagate	(	const TrajectoryStateOnSurface &	tsos,
		const reco::BeamSpot &	beamSpot
	)		const

Propagate the state to the 2D-PCA.

Propagate the state to the PCA in 2D, i.e. to the beam line.

Definition at line 61 of file MuonUpdatorAtVertex.cc.

References beam_dqm_sourceclient-live_cfg::beamSpot, TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), FreeTrajectoryState::hasError(), TrackerBounds::isInside(), TrajectoryStateClosestToBeamLine::isValid(), LogTrace, metname, FreeTrajectoryState::parameters(), mps_fire::result, thePropagatorName, theService, and TrajectoryStateClosestToBeamLine::trackStateAtPCA().

Referenced by AlignmentMonitorTracksFromTrajectories::event(), and propagateWithUpdate().

                                                                                                    {
   const string metname = "Muon|RecoMuon|MuonUpdatorAtVertex";
 
   if (TrackerBounds::isInside(tsos.globalPosition())) {
     LogTrace(metname) << "Trajectory inside the Tracker";
 
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tscbl = tscblBuilder(*(tsos.freeState()), beamSpot);
 
     if (tscbl.isValid())
       return pair<bool, FreeTrajectoryState>(true, tscbl.trackStateAtPCA());
     else
       edm::LogWarning(metname) << "Propagation to the PCA using TSCPBuilderNoMaterial failed!"
                                << " This can cause a severe bug.";
   } else {
     LogTrace(metname) << "Trajectory inside the muon system";
 
     FreeTrajectoryState result = theService->propagator(thePropagatorName)->propagate(*tsos.freeState(), beamSpot);
 
     LogTrace(metname) << "MuonUpdatorAtVertex::propagate, path: " << result << " parameters: " << result.parameters();
 
     if (result.hasError())
       return pair<bool, FreeTrajectoryState>(true, result);
     else
       edm::LogInfo(metname) << "Propagation to the PCA failed!";
   }
   return pair<bool, FreeTrajectoryState>(false, FreeTrajectoryState());
 }

std::pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::propagate ( const TrajectoryStateOnSurface & tsos ) const

Propagate the state to the 2D-PCA (nominal CMS axis) - DEPRECATED -.

Definition at line 177 of file MuonUpdatorAtVertex.cc.

References propagateToNominalLine().

                                                                                                             {
   return propagateToNominalLine(tsos);
 }

std::pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::propagateToNominalLine ( const TrajectoryStateOnSurface & tsos ) const

Propagate the state to the 2D-PCA (nominal CMS axis)

Definition at line 141 of file MuonUpdatorAtVertex.cc.

References TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), TrackerBounds::isInside(), TrajectoryStateClosestToPoint::isValid(), LogTrace, metname, fireworks::p1, fireworks::p2, mps_fire::result, thePropagatorName, theService, and TrajectoryStateClosestToPoint::theState().

Referenced by propagate().

                                                 {
   const string metname = "Muon|RecoMuon|MuonUpdatorAtVertex";
 
   if (TrackerBounds::isInside(tsos.globalPosition())) {
     LogTrace(metname) << "Trajectory inside the Tracker";
 
     TSCPBuilderNoMaterial tscpBuilder;
     TrajectoryStateClosestToPoint tscp = tscpBuilder(*(tsos.freeState()), GlobalPoint(0., 0., 0.));
 
     if (tscp.isValid())
       return pair<bool, FreeTrajectoryState>(true, tscp.theState());
     else
       edm::LogWarning(metname) << "Propagation to the PCA using TSCPBuilderNoMaterial failed!"
                                << " This can cause a severe bug.";
   } else {
     LogTrace(metname) << "Trajectory inside the muon system";
 
     // Define a line using two 3D-points
     GlobalPoint p1(0., 0., -1500);
     GlobalPoint p2(0., 0., 1500);
 
     pair<FreeTrajectoryState, double> result =
         theService->propagator(thePropagatorName)->propagateWithPath(*tsos.freeState(), p1, p2);
 
     LogTrace(metname) << "MuonUpdatorAtVertex::propagate, path: " << result.second
                       << " parameters: " << result.first.parameters();
 
     if (result.first.hasError())
       return pair<bool, FreeTrajectoryState>(true, result.first);
     else
       edm::LogInfo(metname) << "Propagation to the PCA failed! Path: " << result.second;
   }
   return pair<bool, FreeTrajectoryState>(false, FreeTrajectoryState());
 }

pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::propagateWithUpdate	(	const TrajectoryStateOnSurface &	tsos,
		const reco::BeamSpot &	beamSpot
	)		const

Propagate to the 2D-PCA and apply the vertex constraint.

Definition at line 129 of file MuonUpdatorAtVertex.cc.

References propagate(), and update().

Referenced by AlignmentMonitorTracksFromTrajectories::event().

                                                                                                              {
   pair<bool, FreeTrajectoryState> propagationResult = propagate(tsos, beamSpot);
 
   if (propagationResult.first) {
     return update(propagationResult.second, beamSpot);
   } else {
     edm::LogInfo("Muon|RecoMuon|MuonUpdatorAtVertex") << "Constraint at vertex failed";
     return pair<bool, FreeTrajectoryState>(false, FreeTrajectoryState());
   }
 }

pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::update	(	const reco::TransientTrack &	track,
		const reco::BeamSpot &	beamSpot
	)		const

Applies the vertex constraint.

Definition at line 91 of file MuonUpdatorAtVertex.cc.

References SingleTrackVertexConstraint::constrain(), alignCSCRings::e, cms::Exception::explainSelf(), LogTrace, metname, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, theChi2Cut, theConstrictor, thePositionErrors, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), MatrixUtil.Steps::overwrite(), propagateWithUpdate(), and update().

                                                                                                 {
   const std::string metname = "Muon|RecoMuon|MuonUpdatorAtVertex";
 
   pair<bool, FreeTrajectoryState> result(false, FreeTrajectoryState());
 
   GlobalPoint spotPos(beamSpot.x0(), beamSpot.y0(), beamSpot.z0());
 
   SingleTrackVertexConstraint::BTFtuple constrainedTransientTrack;
 
   // Temporary hack here. A fix in the SingleTrackVertexConstraint code
   // is needed. Once available try&catch will be removed
   try {
     constrainedTransientTrack = theConstrictor.constrain(track, spotPos, thePositionErrors);
   } catch (cms::Exception &e) {
     edm::LogWarning(metname) << "cms::Exception caught in MuonUpdatorAtVertex::update\n"
                              << "Exception from SingleTrackVertexConstraint\n"
                              << e.explainSelf();
     return result;
   }
 
   if (std::get<0>(constrainedTransientTrack))
     if (std::get<2>(constrainedTransientTrack) <= theChi2Cut) {
       result.first = true;
       result.second = *std::get<1>(constrainedTransientTrack).impactPointState().freeState();
     } else
       LogTrace(metname) << "Constraint at vertex failed: too large chi2";
   else
     LogTrace(metname) << "Constraint at vertex failed";
 
   return result;
 }

pair< bool, FreeTrajectoryState > MuonUpdatorAtVertex::update	(	const FreeTrajectoryState &	ftsAtVtx,
		const reco::BeamSpot &	beamSpot
	)		const

Applies the vertex constraint.

Definition at line 124 of file MuonUpdatorAtVertex.cc.

References beam_dqm_sourceclient-live_cfg::beamSpot, TransientTrackFromFTSFactory::build(), theTransientTrackFactory, and update().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                                                                 {
   return update(theTransientTrackFactory.build(ftsAtVtx), beamSpot);
 }