tmtt::TP Class Reference

#include <TP.h>

Public Member Functions
const std::vector< const Stub * > &	assocStubs () const
int	charge () const
float	d0 () const
float	dphi (float rad) const
float	eta () const
void	fillNearestJetInfo (const reco::GenJetCollection *genJets)
void	fillTruth (const std::list< Stub > &vStubs)
unsigned int	index () const
bool	inTimeBx () const
float	mass () const
float	nearestJetPt () const
unsigned int	numAssocStubs () const
unsigned int	numLayers () const
bool	operator== (const TP &tpOther) const
int	pdgId () const
float	phi0 () const
bool	physicsCollision () const
float	pt () const
float	qOverPt () const
float	tanLambda () const
float	theta () const
	TP (const TrackingParticlePtr &tpPtr, unsigned int index_in_vTPs, const Settings *settings)
float	tpInJet (float genJetPtCut=30.) const
const TrackingParticlePtr &	trackingParticlePtr () const
float	trkPhiAtR (float rad) const
float	trkPhiAtStub (const Stub *stub) const
float	trkRAtStub (const Stub *stub) const
float	trkZAtR (float rad) const
float	trkZAtStub (const Stub *stub) const
bool	use () const
bool	useForAlgEff () const
bool	useForEff () const
float	vx () const
float	vy () const
float	vz () const
float	z0 () const

Private Member Functions
void	calcNumLayers ()
void	fillUse ()
void	fillUseForAlgEff ()
void	fillUseForEff ()

Private Attributes
std::vector< const Stub * >	assocStubs_
int	charge_
float	d0_
float	eta_
unsigned int	index_in_vTPs_
bool	inTimeBx_
float	mass_
float	nearestJetPt_
unsigned int	nLayersWithStubs_
int	pdgId_
float	phi0_
bool	physicsCollision_
float	pt_
const Settings *	settings_
float	tanLambda_
float	theta_
bool	tpInJet_
TrackingParticlePtr	trackingParticlePtr_
bool	use_
bool	useForAlgEff_
bool	useForEff_
float	vx_
float	vy_
float	vz_
float	z0_

Detailed Description

Definition at line 23 of file TP.h.

Constructor & Destructor Documentation

TP()

tmtt::TP::TP	(	const TrackingParticlePtr &	tpPtr,
		unsigned int	index_in_vTPs,
		const Settings *	settings
	)

Definition at line 16 of file TP.cc.

References runTauDisplay::eid, fillUse(), fillUseForEff(), TrackingParticle::g4Tracks(), inTimeBx_, physicsCollision_, and x.

      : trackingParticlePtr_(tpPtr),
        index_in_vTPs_(index_in_vTPs),
        settings_(settings),
        pdgId_(tpPtr->pdgId()),
        charge_(tpPtr->charge()),
        mass_(tpPtr->mass()),
        pt_(tpPtr->pt()),
        eta_(tpPtr->eta()),
        theta_(tpPtr->theta()),
        tanLambda_(1. / tan(theta_)),
        phi0_(tpPtr->phi()),
        vx_(tpPtr->vertex().x()),
        vy_(tpPtr->vertex().y()),
        vz_(tpPtr->vertex().z()),
        d0_(vx_ * sin(phi0_) - vy_ * cos(phi0_)),                      // Copied from CMSSW class TrackBase::d0().
        z0_(vz_ - (vx_ * cos(phi0_) + vy_ * sin(phi0_)) * tanLambda_)  // Copied from CMSSW class TrackBase::dz().
  {
    const vector<SimTrack>& vst = tpPtr->g4Tracks();
    EncodedEventId eid = vst.at(0).eventId();
    inTimeBx_ = (eid.bunchCrossing() == 0);  // TP from in-time or out-of-time Bx.
    physicsCollision_ = (eid.event() == 0);  // TP from physics collision or from pileup.

    this->fillUse();        // Fill use_ flag, indicating if TP is worth keeping.
    this->fillUseForEff();  // Fill useForEff_ flag, indicating if TP is good for tracking efficiency measurement.
  }

Member Function Documentation

assocStubs()

const std::vector<const Stub*>& tmtt::TP::assocStubs ( ) const

inline

Definition at line 79 of file TP.h.

References assocStubs_.

                                                     {
      return assocStubs_;
    }  // associated stubs. (Includes those failing tightened front-end electronics cuts supplied by user). (Which stubs are returned is affected by "StubMatchStrict" config param.)

calcNumLayers()

void tmtt::TP::calcNumLayers ( )

inlineprivate

Definition at line 103 of file TP.h.

References assocStubs_, tmtt::Utility::countLayers(), nLayersWithStubs_, and settings_.

Referenced by fillTruth().

{ nLayersWithStubs_ = Utility::countLayers(settings_, assocStubs_, false); }

charge()

int tmtt::TP::charge ( void  ) const

inline

Definition at line 46 of file TP.h.

References charge_.

{ return charge_; }

d0()

float tmtt::TP::d0 ( ) const

inline

Definition at line 63 of file TP.h.

References d0_.

{ return d0_; }

dphi()

float tmtt::TP::dphi ( float  rad ) const

inline

Definition at line 66 of file TP.h.

References charge_, tmtt::Settings::invPtToDphi(), pt_, and settings_.

Referenced by trkPhiAtR(), and trkPhiAtStub().

{ return asin(settings_->invPtToDphi() * rad * charge_ / pt_); }

eta()

float tmtt::TP::eta ( void  ) const

inline

Definition at line 54 of file TP.h.

References eta_.

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), fillNearestJetInfo(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

{ return eta_; }

fillNearestJetInfo()

void tmtt::TP::fillNearestJetInfo

(

const reco::GenJetCollection *

genJets

)

Definition at line 166 of file TP.cc.

References funct::abs(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), reco::deltaR(), PbPb_ZMuSkimMuonDPG_cff::deltaR, eta(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenJet, l1tCaloJetHTTProducer_cfi::genJets, razorScouting_cff::maxEta, HLT_2023v12_cff::minDR, PV_cfg::minPt, nearestJetPt_, reco::LeafCandidate::phi(), phi0(), reco::LeafCandidate::pt(), and tpInJet_.

                                                                 {
    double minDR = 999.;
    double ptOfNearestJet = -1;

    reco::GenJetCollection::const_iterator iterGenJet;
    for (iterGenJet = genJets->begin(); iterGenJet != genJets->end(); ++iterGenJet) {
      reco::GenJet myJet = reco::GenJet(*iterGenJet);

      // Don't consider GenJets failing these cuts.
      constexpr float minPt = 30.0;
      constexpr float maxEta = 2.5;

      if (myJet.pt() > minPt && std::abs(myJet.eta()) > maxEta) {
        double deltaR = reco::deltaR(this->eta(), this->phi0(), myJet.eta(), myJet.phi());

        if (deltaR < minDR) {
          minDR = deltaR;
          ptOfNearestJet = myJet.pt();
        }
      }
    }

    // Only consider GenJets within this distance of TP.
    constexpr float cutDR = 0.4;
    tpInJet_ = (minDR < cutDR);
    nearestJetPt_ = tpInJet_ ? ptOfNearestJet : -1.;
  }

fillTruth()

void tmtt::TP::fillTruth

(

const std::list< Stub > &

vStubs

)

Definition at line 45 of file TP.cc.

References assocStubs_, calcNumLayers(), fillUseForAlgEff(), index(), and alignCSCRings::s.

                                             {
    for (const Stub& s : vStubs) {
      for (const TP* tp_i : s.assocTPs()) {
        if (tp_i->index() == this->index())
          assocStubs_.push_back(&s);
      }
    }

    this->fillUseForAlgEff();  // Fill useForAlgEff_ flag.

    this->calcNumLayers();  // Calculate number of tracker layers this TP has stubs in.
  }

fillUse()

void tmtt::TP::fillUse ( )

private

Definition at line 61 of file TP.cc.

References funct::abs(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), l1pfProducerBarrel_cfi::etaExtra, ALCARECOTkAlBeamHalo_cff::etaMax, tmtt::Settings::etaRegions(), tmtt::Settings::genMaxAbsEta(), tmtt::Settings::genMaxVertR(), tmtt::Settings::genMaxVertZ(), tmtt::Settings::genMinPt(), tmtt::Settings::houghMinPt(), SiStripPI::max, SiStripPI::min, AlignmentTrackSelector_cfi::ptMax, ptMin, settings_, trackingParticlePtr_, trackingParticleSelector_cfi::trackingParticleSelector, and use_.

Referenced by TP().

                   {
    constexpr bool useOnlyInTimeParticles = false;
    constexpr bool useOnlyTPfromPhysicsCollisionFalse = false;
    // Use looser cuts here those those used for tracking efficiency measurement.
    // Keep only those TP that have a chance (allowing for finite track resolution) of being reconstructed as L1 tracks. L1 tracks not matching these TP will be defined as fake.

    // Include all possible particle types here, as if some are left out, L1 tracks matching one of missing types will be declared fake.
    constexpr std::array<int, 5> genPdgIdsAllUnsigned = {{11, 13, 211, 321, 2212}};
    vector<int> genPdgIdsAll;
    for (const int& iPdg : genPdgIdsAllUnsigned) {
      genPdgIdsAll.push_back(iPdg);
      genPdgIdsAll.push_back(-iPdg);
    }

    // Range big enough to include all TP needed to measure tracking efficiency
    // and big enough to include any TP that might be reconstructed for fake rate measurement.
    constexpr float ptMinScale = 0.7;
    const float ptMin = min(settings_->genMinPt(), ptMinScale * settings_->houghMinPt());
    constexpr double ptMax = 9.9e9;
    const float etaExtra = 0.2;
    const float etaMax = max(settings_->genMaxAbsEta(), etaExtra + std::abs(settings_->etaRegions()[0]));
    constexpr double fixedVertRcut = 10.;
    constexpr double fixedVertZcut = 35.;

    static const TrackingParticleSelector trackingParticleSelector(ptMin,
                                                                   ptMax,
                                                                   -etaMax,
                                                                   etaMax,
                                                                   max(fixedVertRcut, settings_->genMaxVertR()),
                                                                   max(fixedVertZcut, settings_->genMaxVertZ()),
                                                                   0,
                                                                   useOnlyTPfromPhysicsCollisionFalse,
                                                                   useOnlyInTimeParticles,
                                                                   true,
                                                                   false,
                                                                   genPdgIdsAll);

    use_ = trackingParticleSelector(*trackingParticlePtr_);
  }

fillUseForAlgEff()

void tmtt::TP::fillUseForAlgEff ( )

private

Definition at line 134 of file TP.cc.

References assocStubs_, tmtt::Utility::countLayers(), tmtt::Settings::genMinStubLayers(), settings_, useForAlgEff_, and useForEff_.

Referenced by fillTruth().

                            {
    useForAlgEff_ = false;
    if (useForEff_) {
      useForAlgEff_ = (Utility::countLayers(settings_, assocStubs_, true) >= settings_->genMinStubLayers());
    }
  }

fillUseForEff()

void tmtt::TP::fillUseForEff ( )

private

Definition at line 103 of file TP.cc.

References funct::abs(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), d0_, tmtt::Settings::genMaxAbsEta(), tmtt::Settings::genMaxD0(), tmtt::Settings::genMaxVertR(), tmtt::Settings::genMaxVertZ(), tmtt::Settings::genMaxZ0(), tmtt::Settings::genMinPt(), tmtt::Settings::genPdgIds(), AlignmentTrackSelector_cfi::ptMax, settings_, trackingParticlePtr_, trackingParticleSelector_cfi::trackingParticleSelector, use_, useForEff_, and z0_.

Referenced by TP().

                         {
    useForEff_ = false;
    if (use_) {
      constexpr bool useOnlyInTimeParticles = true;
      constexpr bool useOnlyTPfromPhysicsCollision = true;
      constexpr double ptMax = 9.9e9;
      static const TrackingParticleSelector trackingParticleSelector(settings_->genMinPt(),
                                                                     ptMax,
                                                                     -settings_->genMaxAbsEta(),
                                                                     settings_->genMaxAbsEta(),
                                                                     settings_->genMaxVertR(),
                                                                     settings_->genMaxVertZ(),
                                                                     0,
                                                                     useOnlyTPfromPhysicsCollision,
                                                                     useOnlyInTimeParticles,
                                                                     true,
                                                                     false,
                                                                     settings_->genPdgIds());

      useForEff_ = trackingParticleSelector(*trackingParticlePtr_);

      // Add additional cut on particle transverse impact parameter.
      if (std::abs(d0_) > settings_->genMaxD0())
        useForEff_ = false;
      if (std::abs(z0_) > settings_->genMaxZ0())
        useForEff_ = false;
    }
  }

index()

unsigned int tmtt::TP::index ( ) const

inline

Definition at line 39 of file TP.h.

References index_in_vTPs_.

Referenced by tmtt::L1track3D::cheat(), fillTruth(), and operator==().

{ return index_in_vTPs_; }

inTimeBx()

bool tmtt::TP::inTimeBx ( ) const

inline

Definition at line 43 of file TP.h.

References inTimeBx_.

{ return inTimeBx_; }

mass()

float tmtt::TP::mass ( ) const

inline

Definition at line 47 of file TP.h.

References mass_.

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), and datamodel.Object::p4().

{ return mass_; }

nearestJetPt()

float tmtt::TP::nearestJetPt ( ) const

inline

Definition at line 95 of file TP.h.

References nearestJetPt_.

{ return nearestJetPt_; }  // -ve if no nearest jet.

numAssocStubs()

unsigned int tmtt::TP::numAssocStubs ( ) const

inline

Definition at line 82 of file TP.h.

References assocStubs_.

{ return assocStubs_.size(); }

numLayers()

unsigned int tmtt::TP::numLayers ( ) const

inline

Definition at line 83 of file TP.h.

References nLayersWithStubs_.

{ return nLayersWithStubs_; }

operator==()

bool tmtt::TP::operator== ( const TP &  tpOther ) const

inline

Definition at line 31 of file TP.h.

References index().

{ return (this->index() == tpOther.index()); }

pdgId()

int tmtt::TP::pdgId ( ) const

inline

Definition at line 41 of file TP.h.

References pdgId_.

Referenced by Particle.Particle::__str__().

{ return pdgId_; }

phi0()

float tmtt::TP::phi0 ( ) const

inline

Definition at line 57 of file TP.h.

References phi0_.

Referenced by fillNearestJetInfo(), and tmtt::SimpleLR4::fit().

{ return phi0_; }

physicsCollision()

bool tmtt::TP::physicsCollision ( ) const

inline

Definition at line 45 of file TP.h.

References physicsCollision_.

{ return physicsCollision_; }

pt()

float tmtt::TP::pt ( ) const

inline

Definition at line 48 of file TP.h.

References pt_.

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), tmtt::KFParamsComb::isGoodState(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), datamodel.Object::p4(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

{ return pt_; }

qOverPt()

float tmtt::TP::qOverPt ( ) const

inline

Definition at line 50 of file TP.h.

References charge_, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), and pt_.

Referenced by tmtt::SimpleLR4::fit().

                          {
      constexpr float big = 9.9e9;
      return (pt_ > 0) ? charge_ / pt_ : big;
    }

tanLambda()

float tmtt::TP::tanLambda ( ) const

inline

Definition at line 56 of file TP.h.

References tanLambda_.

{ return tanLambda_; }

theta()

float tmtt::TP::theta ( void  ) const

inline

Definition at line 55 of file TP.h.

References theta_.

Referenced by Tau.Tau::zImpact().

{ return theta_; }

tpInJet()

float tmtt::TP::tpInJet ( float  genJetPtCut = 30. ) const

inline

Definition at line 94 of file TP.h.

References nearestJetPt_, and tpInJet_.

{ return (tpInJet_ && nearestJetPt_ > genJetPtCut); }

trackingParticlePtr()

const TrackingParticlePtr& tmtt::TP::trackingParticlePtr ( ) const

inline

Definition at line 29 of file TP.h.

References trackingParticlePtr_.

{ return trackingParticlePtr_; }

trkPhiAtR()

float tmtt::TP::trkPhiAtR ( float  rad ) const

inline

Definition at line 68 of file TP.h.

References d0_, reco::deltaPhi(), dphi(), and phi0_.

{ return reco::deltaPhi(phi0_ - this->dphi(rad) - d0_ / rad, 0.); }

trkPhiAtStub()

float tmtt::TP::trkPhiAtStub

(

const Stub *

stub

)

const

Definition at line 143 of file TP.cc.

References dphi(), phi0_, run3scouting_cff::trkPhi, and trkRAtStub().

                                               {
    float trkPhi = phi0_ - this->dphi(this->trkRAtStub(stub));
    return trkPhi;
  }

trkRAtStub()

float tmtt::TP::trkRAtStub

(

const Stub *

stub

)

const

Definition at line 152 of file TP.cc.

References tmtt::Stub::barrel(), tmtt::Stub::r(), tanLambda_, tmtt::Stub::z(), and z0_.

Referenced by trkPhiAtStub().

                                             {
    float rTrk = (stub->barrel()) ? stub->r() : (stub->z() - z0_) / tanLambda_;
    return rTrk;
  }

trkZAtR()

float tmtt::TP::trkZAtR ( float  rad ) const

inline

Definition at line 70 of file TP.h.

References tanLambda_, and vz_.

{ return (vz_ + rad * tanLambda_); }

trkZAtStub()

float tmtt::TP::trkZAtStub

(

const Stub *

stub

)

const

Definition at line 161 of file TP.cc.

References tmtt::Stub::barrel(), tmtt::Stub::r(), tanLambda_, tmtt::Stub::z(), and z0_.

                                             {
    float zTrk = (stub->barrel()) ? z0_ + tanLambda_ * stub->r() : stub->z();
    return zTrk;
  }

use()

bool tmtt::TP::use ( ) const

inline

Definition at line 85 of file TP.h.

References use_.

{ return use_; }

useForAlgEff()

bool tmtt::TP::useForAlgEff ( ) const

inline

Definition at line 89 of file TP.h.

References useForAlgEff_.

Referenced by tmtt::L1track3D::cheat().

{ return useForAlgEff_; }

useForEff()

bool tmtt::TP::useForEff ( ) const

inline

Definition at line 87 of file TP.h.

References useForEff_.

{ return useForEff_; }

vx()

float tmtt::TP::vx ( ) const

inline

Definition at line 59 of file TP.h.

References vx_.

{ return vx_; }

vy()

float tmtt::TP::vy ( ) const

inline

Definition at line 60 of file TP.h.

References vy_.

{ return vy_; }

vz()

float tmtt::TP::vz ( ) const

inline

Definition at line 61 of file TP.h.

References vz_.

{ return vz_; }

z0()

float tmtt::TP::z0 ( ) const

inline

Definition at line 64 of file TP.h.

References z0_.

{ return z0_; }

Member Data Documentation

assocStubs_

std::vector<const Stub*> tmtt::TP::assocStubs_

private

Definition at line 128 of file TP.h.

Referenced by assocStubs(), calcNumLayers(), fillTruth(), fillUseForAlgEff(), and numAssocStubs().

charge_

int tmtt::TP::charge_

private

Definition at line 115 of file TP.h.

Referenced by charge(), dphi(), and qOverPt().

d0_

float tmtt::TP::d0_

private

Definition at line 125 of file TP.h.

Referenced by d0(), fillUseForEff(), and trkPhiAtR().

eta_

float tmtt::TP::eta_

private

Definition at line 118 of file TP.h.

Referenced by eta().

index_in_vTPs_

unsigned int tmtt::TP::index_in_vTPs_

private

Definition at line 108 of file TP.h.

Referenced by index().

inTimeBx_

bool tmtt::TP::inTimeBx_

private

Definition at line 113 of file TP.h.

Referenced by inTimeBx(), and TP().

mass_

float tmtt::TP::mass_

private

Definition at line 116 of file TP.h.

Referenced by mass().

nearestJetPt_

float tmtt::TP::nearestJetPt_

private

Definition at line 136 of file TP.h.

Referenced by fillNearestJetInfo(), nearestJetPt(), and tpInJet().

nLayersWithStubs_

unsigned int tmtt::TP::nLayersWithStubs_

private

Definition at line 129 of file TP.h.

Referenced by calcNumLayers(), and numLayers().

pdgId_

int tmtt::TP::pdgId_

private

Definition at line 112 of file TP.h.

Referenced by pdgId().

phi0_

float tmtt::TP::phi0_

private

Definition at line 121 of file TP.h.

Referenced by phi0(), trkPhiAtR(), and trkPhiAtStub().

physicsCollision_

bool tmtt::TP::physicsCollision_

private

Definition at line 114 of file TP.h.

Referenced by physicsCollision(), and TP().

pt_

float tmtt::TP::pt_

private

Definition at line 117 of file TP.h.

Referenced by dphi(), pt(), and qOverPt().

settings_

const Settings* tmtt::TP::settings_

private

Definition at line 110 of file TP.h.

Referenced by calcNumLayers(), dphi(), fillUse(), fillUseForAlgEff(), and fillUseForEff().

tanLambda_

float tmtt::TP::tanLambda_

private

Definition at line 120 of file TP.h.

Referenced by tanLambda(), trkRAtStub(), trkZAtR(), and trkZAtStub().

theta_

float tmtt::TP::theta_

private

Definition at line 119 of file TP.h.

Referenced by theta().

tpInJet_

bool tmtt::TP::tpInJet_

private

Definition at line 135 of file TP.h.

Referenced by fillNearestJetInfo(), and tpInJet().

trackingParticlePtr_

TrackingParticlePtr tmtt::TP::trackingParticlePtr_

private

Definition at line 106 of file TP.h.

Referenced by fillUse(), fillUseForEff(), and trackingParticlePtr().

use_

bool tmtt::TP::use_

private

Definition at line 131 of file TP.h.

Referenced by fillUse(), fillUseForEff(), and use().

useForAlgEff_

bool tmtt::TP::useForAlgEff_

private

Definition at line 133 of file TP.h.

Referenced by fillUseForAlgEff(), and useForAlgEff().

useForEff_

bool tmtt::TP::useForEff_

private

Definition at line 132 of file TP.h.

Referenced by fillUseForAlgEff(), fillUseForEff(), and useForEff().

vx_

float tmtt::TP::vx_

private

Definition at line 122 of file TP.h.

Referenced by vx().

vy_

float tmtt::TP::vy_

private

Definition at line 123 of file TP.h.

Referenced by vy().

vz_

float tmtt::TP::vz_

private

Definition at line 124 of file TP.h.

Referenced by trkZAtR(), and vz().

z0_

float tmtt::TP::z0_

private

Definition at line 126 of file TP.h.

Referenced by fillUseForEff(), trkRAtStub(), trkZAtStub(), and z0().