tauImpactParameter::ParticleBuilder Class Reference

#include <ParticleBuilder.h>

Public Types
enum	CMSSWPerigee {   aCurv = 0, aTheta, aPhi, aTip,   aLip }

Public Member Functions
	ParticleBuilder ()
	~ParticleBuilder ()

Static Public Member Functions
static LorentzVectorParticle	createLorentzVectorParticle (const reco::TransientTrack &transTrk, const reco::Vertex &V, bool fromPerigee, bool useTrackHelixPropagation)
static TrackParticle	createTrackParticle (const reco::TransientTrack &transTrk, const GlobalPoint &p, bool fromPerigee=true, bool useTrackHelixPropogation=true)
static reco::Vertex	getVertex (const LorentzVectorParticle &p)

Static Private Member Functions
static TVectorT< double >	convertCMSSWTrackParToSFTrackPar (const TVectorT< double > &inpar)
static TVectorT< double >	convertCMSSWTrackPerigeeToSFTrackPar (const TVectorT< double > &inpar)

Detailed Description

Definition at line 25 of file ParticleBuilder.h.

Member Enumeration Documentation

CMSSWPerigee

enum tauImpactParameter::ParticleBuilder::CMSSWPerigee

Enumerator
aCurv
aTheta
aPhi
aTip
aLip

Definition at line 27 of file ParticleBuilder.h.

{ aCurv = 0, aTheta, aPhi, aTip, aLip };

Constructor & Destructor Documentation

ParticleBuilder()

tauImpactParameter::ParticleBuilder::ParticleBuilder ( )

inline

Definition at line 29 of file ParticleBuilder.h.

{};

~ParticleBuilder()

tauImpactParameter::ParticleBuilder::~ParticleBuilder ( )

inline

Definition at line 30 of file ParticleBuilder.h.

{};

Member Function Documentation

convertCMSSWTrackParToSFTrackPar()

TVectorT< double > ParticleBuilder::convertCMSSWTrackParToSFTrackPar ( const TVectorT< double > &  inpar )

staticprivate

Definition at line 136 of file ParticleBuilder.cc.

References funct::cos(), tauImpactParameter::TrackParticle::dxy, tauImpactParameter::TrackParticle::dz, reco::TrackBase::i_dsz, reco::TrackBase::i_dxy, reco::TrackBase::i_lambda, reco::TrackBase::i_phi, reco::TrackBase::i_qoverp, tauImpactParameter::TrackParticle::kappa, tauImpactParameter::TrackParticle::lambda, tauImpactParameter::TrackParticle::NHelixPar, and tauImpactParameter::TrackParticle::phi.

Referenced by createTrackParticle().

                                                                                                {
  TVectorT<double> par(TrackParticle::NHelixPar);
  par(TrackParticle::kappa) = -1.0 * inpar(TrackParticle::NHelixPar) * inpar(reco::TrackBase::i_qoverp) /
                              fabs(cos(inpar(reco::TrackBase::i_lambda)));
  par(TrackParticle::lambda) = inpar(reco::TrackBase::i_lambda);
  par(TrackParticle::phi) = inpar(reco::TrackBase::i_phi);
  par(TrackParticle::dz) = inpar(reco::TrackBase::i_dsz) / fabs(cos(inpar(reco::TrackBase::i_lambda)));
  par(TrackParticle::dxy) = inpar(reco::TrackBase::i_dxy);
  return par;
}

convertCMSSWTrackPerigeeToSFTrackPar()

TVectorT< double > ParticleBuilder::convertCMSSWTrackPerigeeToSFTrackPar ( const TVectorT< double > &  inpar )

staticprivate

Definition at line 147 of file ParticleBuilder.cc.

References aCurv, aLip, aPhi, aTheta, aTip, tauImpactParameter::TrackParticle::dxy, tauImpactParameter::TrackParticle::dz, tauImpactParameter::TrackParticle::kappa, tauImpactParameter::TrackParticle::lambda, tauImpactParameter::TrackParticle::NHelixPar, tauImpactParameter::TrackParticle::phi, and Pi.

Referenced by createTrackParticle().

                                                                                                    {
  TVectorT<double> par(TrackParticle::NHelixPar);
  par(TrackParticle::kappa) = inpar(aCurv);
  par(TrackParticle::lambda) = TMath::Pi() / 2 - inpar(aTheta);
  par(TrackParticle::phi) = inpar(aPhi);
  par(TrackParticle::dxy) = -inpar(aTip);
  par(TrackParticle::dz) = inpar(aLip);
  return par;
}

createLorentzVectorParticle()

LorentzVectorParticle ParticleBuilder::createLorentzVectorParticle ( const reco::TransientTrack &  transTrk, const reco::Vertex &  V, bool  fromPerigee, bool  useTrackHelixPropagation  )

static

Definition at line 19 of file ParticleBuilder.cc.

References tauImpactParameter::TrackHelixVertexFitter::BField0, tauImpactParameter::TrackHelixVertexFitter::computeLorentzVectorPar(), createTrackParticle(), mps_fire::i, dqmiolumiharvest::j, tauImpactParameter::TrackParticle::kappa, tauImpactParameter::TrackHelixVertexFitter::kappa0, tauImpactParameter::TrackParticle::lambda, tauImpactParameter::TrackHelixVertexFitter::lambda0, tauImpactParameter::TrackHelixVertexFitter::MassOffSet, N, tauImpactParameter::TrackHelixVertexFitter::NExtraPar, tauImpactParameter::TrackHelixVertexFitter::NFreeTrackPar, tauImpactParameter::TrackHelixVertexFitter::NFreeVertexPar, AlCaHLTBitMon_ParallelJobs::p, tauImpactParameter::TrackParticle::phi, tauImpactParameter::TrackHelixVertexFitter::phi0, tauImpactParameter::ErrorMatrixPropagator::propagateError(), cmsswSequenceInfo::tp, cms::cuda::V, tauImpactParameter::TrackHelixVertexFitter::x0, tauImpactParameter::TrackHelixVertexFitter::y0, and tauImpactParameter::TrackHelixVertexFitter::z0.

                                                                                                  {
  GlobalPoint p(V.position().x(), V.position().y(), V.position().z());
  TrackParticle tp = createTrackParticle(transTrk, p, fromPerigee, useTrackHelixPropagation);

  int N =
      TrackHelixVertexFitter::NFreeTrackPar + TrackHelixVertexFitter::NExtraPar + TrackHelixVertexFitter::MassOffSet;
  TVectorT<double> par(N);
  par(TrackHelixVertexFitter::x0) = V.position().x();
  par(TrackHelixVertexFitter::y0) = V.position().y();
  par(TrackHelixVertexFitter::z0) = V.position().z();
  par(TrackHelixVertexFitter::kappa0) = tp.parameter(TrackParticle::kappa);
  par(TrackHelixVertexFitter::lambda0) = tp.parameter(TrackParticle::lambda);
  par(TrackHelixVertexFitter::phi0) = tp.parameter(TrackParticle::phi);
  par(TrackHelixVertexFitter::NFreeTrackPar + TrackHelixVertexFitter::MassOffSet) = tp.mass();
  par(TrackHelixVertexFitter::NFreeTrackPar + TrackHelixVertexFitter::BField0) = tp.bField();

  TMatrixTSym<double> parcov(N);
  for (int i = 0; i < TrackHelixVertexFitter::NFreeVertexPar; i++) {
    for (int j = 0; j < TrackHelixVertexFitter::NFreeVertexPar; j++) {
      parcov(i, j) = V.covariance(i, j);
    }
  }
  parcov(TrackHelixVertexFitter::kappa0, TrackHelixVertexFitter::kappa0) =
      tp.covariance(TrackParticle::kappa, TrackParticle::kappa);
  parcov(TrackHelixVertexFitter::lambda0, TrackHelixVertexFitter::lambda0) =
      tp.covariance(TrackParticle::lambda, TrackParticle::lambda);
  parcov(TrackHelixVertexFitter::phi0, TrackHelixVertexFitter::phi0) =
      tp.covariance(TrackParticle::phi, TrackParticle::phi);

  TVectorT<double> LVPar = TrackHelixVertexFitter::computeLorentzVectorPar(par);
  TMatrixTSym<double> LVCov =
      ErrorMatrixPropagator::propagateError(&TrackHelixVertexFitter::computeLorentzVectorPar, par, parcov);
  return LorentzVectorParticle(LVPar, LVCov, tp.pdgId(), tp.charge(), tp.bField());
}

createTrackParticle()

TrackParticle ParticleBuilder::createTrackParticle ( const reco::TransientTrack &  transTrk, const GlobalPoint &  p, bool  fromPerigee = true, bool  useTrackHelixPropogation = true  )

static

Definition at line 57 of file ParticleBuilder.cc.

References funct::abs(), c, reco::TransientTrack::charge(), tauImpactParameter::TrackHelixVertexFitter::computedxydz(), convertCMSSWTrackParToSFTrackPar(), convertCMSSWTrackPerigeeToSFTrackPar(), reco::TrackBase::covariance(), PerigeeTrajectoryError::covarianceMatrix(), tauImpactParameter::TrackParticle::dxy, PVValHelper::dxy, tauImpactParameter::TrackParticle::dz, PVValHelper::dz, reco::TransientTrack::field(), mps_fire::i, reco::TrackBase::i_qoverp, reco::TransientTrack::impactPointTSCP(), MagneticField::inInverseGeV(), dqmiolumiharvest::j, tauImpactParameter::TrackParticle::kappa, hgcalLayerClusters_cfi::kappa, tauImpactParameter::TrackHelixVertexFitter::kappa0, tauImpactParameter::TrackParticle::lambda, tauImpactParameter::TrackHelixVertexFitter::lambda0, LogDebug, tauImpactParameter::TrackHelixVertexFitter::NFreeTrackPar, tauImpactParameter::TrackParticle::NHelixPar, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::parameter(), TrajectoryStateClosestToPoint::perigeeError(), TrajectoryStateClosestToPoint::perigeeParameters(), tauImpactParameter::TrackParticle::phi, tauImpactParameter::TrackHelixVertexFitter::phi0, tauImpactParameter::PDGInfo::pi_mass(), PdtPdgMini::pi_plus, TrajectoryStateClosestToPoint::position(), tauImpactParameter::ErrorMatrixPropagator::propagateError(), alignCSCRings::s, reco::TransientTrack::track(), reco::TransientTrack::trajectoryStateClosestToPoint(), PerigeeTrajectoryParameters::vector(), x, PV3DBase< T, PVType, FrameType >::x(), tauImpactParameter::TrackHelixVertexFitter::x0, PV3DBase< T, PVType, FrameType >::y(), tauImpactParameter::TrackHelixVertexFitter::y0, PV3DBase< T, PVType, FrameType >::z(), and tauImpactParameter::TrackHelixVertexFitter::z0.

Referenced by createLorentzVectorParticle().

                                                                                  {
  // Configured for CMSSW Tracks only
  TVectorT<double> par(TrackParticle::NHelixPar + 1);
  TMatrixTSym<double> cov(TrackParticle::NHelixPar + 1);
  TVectorT<double> SFpar(TrackParticle::NHelixPar);
  TMatrixTSym<double> SFcov(TrackParticle::NHelixPar);
  if (!fromPerigee) {
    for (int i = 0; i < TrackParticle::NHelixPar; i++) {
      par(i) = transTrk.track().parameter(i);
      for (int j = 0; j < TrackParticle::NHelixPar; j++) {
        cov(i, j) = transTrk.track().covariance(i, j);
      }
    }
    par(TrackParticle::NHelixPar) = transTrk.field()->inInverseGeV(p).z();
    SFpar = convertCMSSWTrackParToSFTrackPar(par);
    SFcov = ErrorMatrixPropagator::propagateError(&ParticleBuilder::convertCMSSWTrackParToSFTrackPar, par, cov);
  } else {
    GlobalPoint TrackIPPos = transTrk.impactPointTSCP().position();
    //GlobalPoint TrackIPOrigin=transTrk.impactPointTSCP().referencePoint();
    GlobalPoint origin(0.0, 0.0, 0.0);
    for (int i = 0; i < TrackParticle::NHelixPar; i++) {
      par(i) = transTrk.trajectoryStateClosestToPoint(origin).perigeeParameters().vector()(i);
      for (int j = 0; j < TrackParticle::NHelixPar; j++) {
        cov(i, j) = transTrk.trajectoryStateClosestToPoint(origin).perigeeError().covarianceMatrix()(i, j);
      }
    }
    par(TrackParticle::NHelixPar) = transTrk.field()->inInverseGeV(p).z();
    SFpar = convertCMSSWTrackPerigeeToSFTrackPar(par);
    SFcov = ErrorMatrixPropagator::propagateError(&ParticleBuilder::convertCMSSWTrackPerigeeToSFTrackPar, par, cov);
    if (useTrackHelixPropagation) {
      // correct dxy dz neglecting material and radiative corrections

      LogDebug("RecoTauImpactParameterParticleBuilder")
          << "Offical CMS dxy - " << par(TrackParticle::dxy) << " dz " << par(TrackParticle::dz) << " kappa "
          << par(reco::TrackBase::i_qoverp);
      LogDebug("RecoTauImpactParameterParticleBuilder")
          << "Offical CMS dxy - SimpleFits Format" << SFpar(TrackParticle::dxy) << " dz " << SFpar(TrackParticle::dz)
          << " kappa " << SFpar(reco::TrackBase::i_qoverp);

      double x, y, z, dxy, dz, s, kappa, lambda, phi;
      TVectorT<double> freehelix(TrackHelixVertexFitter::NFreeTrackPar);
      freehelix(TrackHelixVertexFitter::x0) = TrackIPPos.x();
      freehelix(TrackHelixVertexFitter::y0) = TrackIPPos.y();
      freehelix(TrackHelixVertexFitter::z0) = TrackIPPos.z();
      freehelix(TrackHelixVertexFitter::kappa0) = SFpar(TrackParticle::kappa);
      freehelix(TrackHelixVertexFitter::lambda0) = SFpar(TrackParticle::lambda);
      freehelix(TrackHelixVertexFitter::phi0) = SFpar(TrackParticle::phi);
      TrackHelixVertexFitter::computedxydz(freehelix, 0, kappa, lambda, phi, x, y, z, s, dxy, dz);
      SFpar(TrackParticle::dxy) = dxy;
      SFpar(TrackParticle::dz) = dz;
      LogDebug("RecoTauImpactParameterParticleBuilder") << "Found values dxy " << dxy << " dz " << dz;
      //exit(0);
    }
  }

  PDGInfo pdgInfo;
  double c = transTrk.charge();
  return TrackParticle(
      SFpar, SFcov, abs(PdtPdgMini::pi_plus) * c, pdgInfo.pi_mass(), c, transTrk.field()->inInverseGeV(p).z());
}

getVertex()

reco::Vertex ParticleBuilder::getVertex ( const LorentzVectorParticle &  p )

static

Definition at line 123 of file ParticleBuilder.cc.

References mps_fire::i, dqmiolumiharvest::j, AlCaHLTBitMon_ParallelJobs::p, findQualityFiles::v, and HltBtagValidation_cff::Vertex.

                                                                    {
  TVector3 v = p.vertex();
  TMatrixTSym<double> vcov = p.vertexCov();
  reco::Vertex::Point vp(v.X(), v.Y(), v.Z());
  reco::Vertex::Error ve;
  for (int i = 0; i < vcov.GetNrows(); i++) {
    for (int j = 0; j < vcov.GetNrows(); j++) {
      ve(i, j) = vcov(i, j);
    }
  }
  return reco::Vertex(vp, ve);
}