TSCBLBuilderNoMaterial Class Reference

#include <TSCBLBuilderNoMaterial.h>

Inheritance diagram for TSCBLBuilderNoMaterial:

Public Member Functions
TrajectoryStateClosestToBeamLine	operator() (const FTS &originalFTS, const reco::BeamSpot &beamSpot) const override
	~TSCBLBuilderNoMaterial () override
Public Member Functions inherited from TrajectoryStateClosestToBeamLineBuilder
virtual	~TrajectoryStateClosestToBeamLineBuilder ()=default

Additional Inherited Members
Public Types inherited from TrajectoryStateClosestToBeamLineBuilder
typedef FreeTrajectoryState	FTS

Detailed Description

This class builds a TrajectoryStateClosestToBeamLine given an original FreeTrajectoryState. This new state is then defined at the point of closest approach to the beam line. It is to be used when there is no material between the state and the BeamLine

Definition at line 13 of file TSCBLBuilderNoMaterial.h.

Constructor & Destructor Documentation

~TSCBLBuilderNoMaterial()

TSCBLBuilderNoMaterial::~TSCBLBuilderNoMaterial ( )

inlineoverride

Definition at line 17 of file TSCBLBuilderNoMaterial.h.

{};

Member Function Documentation

operator()()

TrajectoryStateClosestToBeamLine TSCBLBuilderNoMaterial::operator() ( const FTS &  originalFTS, const reco::BeamSpot &  beamSpot  ) const

overridevirtual

Implements TrajectoryStateClosestToBeamLineBuilder.

Definition at line 8 of file TSCBLBuilderNoMaterial.cc.

                                                                                                        {
  TwoTrackMinimumDistance ttmd;
  bool status = ttmd.calculate(
      originalFTS.parameters(),
      GlobalTrajectoryParameters(GlobalPoint(beamSpot.position().x(), beamSpot.position().y(), beamSpot.position().z()),
                                 GlobalVector(beamSpot.dxdz(), beamSpot.dydz(), 1.),
                                 0,
                                 &(originalFTS.parameters().magneticField())));
  if (!status) {
    LogDebug("TrackingTools|PatternTools")
        << "TSCBLBuilderNoMaterial: Failure in TTMD when searching for PCA of track to beamline.\n"
        << "TrajectoryStateClosestToBeamLine is now invalid.";
    return TrajectoryStateClosestToBeamLine();
  }
 
  pair<GlobalPoint, GlobalPoint> points = ttmd.points();
 
  GlobalPoint xTrack = points.first;
  GlobalVector pTrack = GlobalVector(
      GlobalVector::Cylindrical(originalFTS.momentum().perp(), ttmd.firstAngle(), originalFTS.momentum().z()));
 
  double s = ttmd.pathLength().first;
 
  FreeTrajectoryState theFTS;
  if (originalFTS.hasError()) {
    const AlgebraicSymMatrix55& errorMatrix = originalFTS.curvilinearError().matrix();
    AnalyticalCurvilinearJacobian curvilinJacobian(originalFTS.parameters(), xTrack, pTrack, s);
    const AlgebraicMatrix55& jacobian = curvilinJacobian.jacobian();
    CurvilinearTrajectoryError cte(ROOT::Math::Similarity(jacobian, errorMatrix));
 
    theFTS = FreeTrajectoryState(
        GlobalTrajectoryParameters(xTrack, pTrack, originalFTS.charge(), &(originalFTS.parameters().magneticField())),
        cte);
  } else {
    theFTS = FreeTrajectoryState(
        GlobalTrajectoryParameters(xTrack, pTrack, originalFTS.charge(), &(originalFTS.parameters().magneticField())));
  }
  return TrajectoryStateClosestToBeamLine(theFTS, points.second, beamSpot);
}

References pwdgSkimBPark_cfi::beamSpot, TwoTrackMinimumDistance::calculate(), FreeTrajectoryState::charge(), FreeTrajectoryState::curvilinearError(), TwoTrackMinimumDistance::firstAngle(), FreeTrajectoryState::hasError(), AnalyticalCurvilinearJacobian::jacobian(), LogDebug, GlobalTrajectoryParameters::magneticField(), CurvilinearTrajectoryError::matrix(), FreeTrajectoryState::momentum(), FreeTrajectoryState::parameters(), TwoTrackMinimumDistance::pathLength(), PV3DBase< T, PVType, FrameType >::perp(), TwoTrackMinimumDistance::points(), HLT_2018_cff::points, alignCSCRings::s, mps_update::status, and PV3DBase< T, PVType, FrameType >::z().