SignedDecayLength3D Class Reference

#include <SignedDecayLength3D.h>

Public Member Functions
int	id () const
	SignedDecayLength3D ()

Static Public Member Functions
static std::pair< bool, Measurement1D >	apply (const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)

Static Private Member Functions
static TrajectoryStateOnSurface	closestApproachToJet (const FreeTrajectoryState &, const reco::Vertex &, const GlobalVector &, const MagneticField *)

Detailed Description

Threedimensional track decay length (minimum distance of the closest approach to a jet from the primary vertex) signed according to the jet direction

Definition at line 15 of file SignedDecayLength3D.h.

Constructor & Destructor Documentation

SignedDecayLength3D()

SignedDecayLength3D::SignedDecayLength3D ( )

inline

Definition at line 19 of file SignedDecayLength3D.h.

{};

Member Function Documentation

apply()

pair< bool, Measurement1D > SignedDecayLength3D::apply ( const reco::TransientTrack &  track, const GlobalVector &  direction, const reco::Vertex &  vertex  )

static

chech it!!!!!!!!!!!!!!!!!!!!!!!

Definition at line 19 of file SignedDecayLength3D.cc.

                                                                           {
  double theError = 0.;
  bool theIsValid;
 
  //TrajectoryStateOnSurface TSOS = (aRecTrack).impactPointStateOnSurface();
  TrajectoryStateOnSurface TSOS = transientTrack.impactPointState();
  const FreeTrajectoryState* FTS = TSOS.freeTrajectoryState();
 
  TrajectoryStateOnSurface theTSOS = closestApproachToJet(*FTS, vertex, direction, transientTrack.field());
  theIsValid = theTSOS.isValid();
 
  if (theIsValid) {
    GlobalVector J = direction.unit();
    GlobalPoint vertexPosition(vertex.x(), vertex.y(), vertex.z());
 
    double theValue = J.dot(theTSOS.globalPosition() - vertexPosition);
 
    //error calculation
 
    AlgebraicVector3 j;
    j[0] = J.x();
    j[1] = J.y();
    j[2] = J.z();
    AlgebraicVector6 jj;
    jj[0] = J.x();
    jj[1] = J.y();
    jj[2] = J.z();
    jj[3] = 0.;
    jj[4] = 0.;
    jj[5] = 0.;  
    double E1 = ROOT::Math::Similarity(jj, theTSOS.cartesianError().matrix());
    //  double E2 = (aJet.vertex().positionError().matrix()).similarity(j);
    double E2 = ROOT::Math::Similarity(j, vertex.covariance());
 
    theError = sqrt(E1 + E2);
 
    //cout<< "Error ="<< theError<<endl;
    Measurement1D A(theValue, theError);
    return pair<bool, Measurement1D>(theIsValid, A);
  } else {
    return pair<bool, Measurement1D>(theIsValid, Measurement1D(0., 0.));
  }  // endif (isValid)
}  // end constructor declaration

References MaterialEffects_cfi::A, TrajectoryStateOnSurface::cartesianError(), IPTools::closestApproachToJet(), Vector3DBase< T, FrameTag >::dot(), reco::TransientTrack::field(), TrajectoryStateOnSurface::freeTrajectoryState(), TrajectoryStateOnSurface::globalPosition(), reco::TransientTrack::impactPointState(), TrajectoryStateOnSurface::isValid(), dqmiolumiharvest::j, findQualityFiles::jj, CartesianTrajectoryError::matrix(), mathSSE::sqrt(), Vector3DBase< T, FrameTag >::unit(), bphysicsOniaDQM_cfi::vertex, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by heavyIonTools.ConfigureHeavyIons::__call__(), coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), heavyIonTools.ProductionDefaults::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), trigTools.SwitchOnTrigger::__call__(), runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::__call__(), heavyIonTools.SelectionDefaults::__call__(), heavyIonTools.DisbaleMonteCarloDeps::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), tauTools.AddTauCollection::__call__(), trackTools.MakeTrackCandidates::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), trigTools.SwitchOnTriggerMatchingStandAlone::__call__(), trigTools.SwitchOnTriggerMatchEmbedding::__call__(), jetTools.AddJetCollection::__call__(), jetTools.SwitchJetCollection::__call__(), jetTools.UpdateJetCollection::__call__(), jetTools.AddJetID::__call__(), and jetTools.SetTagInfos::__call__().

closestApproachToJet()

TrajectoryStateOnSurface SignedDecayLength3D::closestApproachToJet ( const FreeTrajectoryState &  aFTS, const reco::Vertex &  vertex, const GlobalVector &  aJetDirection, const MagneticField *  field  )

staticprivate

Definition at line 65 of file SignedDecayLength3D.cc.

                                                                                               {
  GlobalVector J = aJetDirection.unit();
 
  Line::PositionType pos(GlobalPoint(vertex.x(), vertex.y(), vertex.z()));
  Line::DirectionType dir(J);
  Line Jet(pos, dir);
 
  AnalyticalTrajectoryExtrapolatorToLine TETL(field);
 
  return TETL.extrapolate(aFTS, Jet);
}

References DeadROC_duringRun::dir, AnalyticalTrajectoryExtrapolatorToLine::extrapolate(), L1TRate_cfi::Jet, Vector3DBase< T, FrameTag >::unit(), and bphysicsOniaDQM_cfi::vertex.

id()

int SignedDecayLength3D::id ( void  ) const

inline

Definition at line 25 of file SignedDecayLength3D.h.

{ return 3; }