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

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

tauImpactParameter::ParticleBuilder::ParticleBuilder ( )

inline

Definition at line 29 of file ParticleBuilder.h.

{};

tauImpactParameter::ParticleBuilder::~ParticleBuilder ( )

inline

Definition at line 30 of file ParticleBuilder.h.

{};

Member Function Documentation

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

staticprivate

Definition at line 123 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;
}

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

staticprivate

Definition at line 133 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;
}

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(), reco::Vertex::covariance(), createTrackParticle(), i, 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, reco::Vertex::position(), tauImpactParameter::ErrorMatrixPropagator::propagateError(), 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());
}

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

static

Definition at line 50 of file ParticleBuilder.cc.

References funct::abs(), trackerHits::c, reco::TransientTrack::charge(), tauImpactParameter::TrackHelixVertexFitter::computedxydz(), convertCMSSWTrackParToSFTrackPar(), convertCMSSWTrackPerigeeToSFTrackPar(), reco::TrackBase::covariance(), PerigeeTrajectoryError::covarianceMatrix(), tauImpactParameter::TrackParticle::dxy, tauImpactParameter::TrackParticle::dz, reco::TransientTrack::field(), i, reco::TrackBase::i_qoverp, reco::TransientTrack::impactPointTSCP(), MagneticField::inInverseGeV(), j, tauImpactParameter::TrackParticle::kappa, kappa, tauImpactParameter::TrackHelixVertexFitter::kappa0, tauImpactParameter::TrackParticle::lambda, tauImpactParameter::TrackHelixVertexFitter::lambda0, LogDebug, tauImpactParameter::TrackHelixVertexFitter::NFreeTrackPar, tauImpactParameter::TrackParticle::NHelixPar, reco::TrackBase::parameter(), TrajectoryStateClosestToPoint::perigeeError(), TrajectoryStateClosestToPoint::perigeeParameters(), phi, 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, detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), tauImpactParameter::TrackHelixVertexFitter::y0, detailsBasic3DVector::z, 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());
}

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

static

Definition at line 112 of file ParticleBuilder.cc.

References i, j, findQualityFiles::v, tauImpactParameter::LorentzVectorParticle::vertex(), and tauImpactParameter::LorentzVectorParticle::vertexCov().

                                                                   {
  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);
}