PFTauTransverseImpactParameters Class Reference

Inheritance diagram for PFTauTransverseImpactParameters:

Public Types
enum	CMSSWPerigee {   aCurv = 0, aTheta, aPhi, aTip,   aLip }
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
	PFTauTransverseImpactParameters (const edm::ParameterSet &iConfig)
void	produce (edm::Event &, const edm::EventSetup &) override
	~PFTauTransverseImpactParameters () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< reco::VertexRef > > >	PFTauPVAToken_
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< std::vector< reco::VertexRef > > > >	PFTauSVAToken_
edm::EDGetTokenT< std::vector< reco::PFTau > >	PFTauToken_
bool	useFullCalculation_

Detailed Description

Definition at line 54 of file PFTauTransverseImpactParameters.cc.

Member Enumeration Documentation

CMSSWPerigee

enum PFTauTransverseImpactParameters::CMSSWPerigee

Enumerator
aCurv
aTheta
aPhi
aTip
aLip

Definition at line 56 of file PFTauTransverseImpactParameters.cc.

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

Constructor & Destructor Documentation

PFTauTransverseImpactParameters()

PFTauTransverseImpactParameters::PFTauTransverseImpactParameters ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 69 of file PFTauTransverseImpactParameters.cc.

    : PFTauToken_(consumes<std::vector<reco::PFTau>>(iConfig.getParameter<edm::InputTag>("PFTauTag"))),
      PFTauPVAToken_(consumes<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef>>>(
          iConfig.getParameter<edm::InputTag>("PFTauPVATag"))),
      PFTauSVAToken_(consumes<edm::AssociationVector<PFTauRefProd, std::vector<std::vector<reco::VertexRef>>>>(
          iConfig.getParameter<edm::InputTag>("PFTauSVATag"))),
      useFullCalculation_(iConfig.getParameter<bool>("useFullCalculation")) {
  produces<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef>>>();
  produces<PFTauTransverseImpactParameterCollection>("PFTauTIP");
}

~PFTauTransverseImpactParameters()

PFTauTransverseImpactParameters::~PFTauTransverseImpactParameters ( )

override

Definition at line 80 of file PFTauTransverseImpactParameters.cc.

{}

Member Function Documentation

fillDescriptions()

void PFTauTransverseImpactParameters::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 193 of file PFTauTransverseImpactParameters.cc.

                                                                                                 {
  // PFTauTransverseImpactParameters
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("PFTauPVATag", edm::InputTag("PFTauPrimaryVertexProducer"));
  desc.add<bool>("useFullCalculation", false);
  desc.add<edm::InputTag>("PFTauTag", edm::InputTag("hpsPFTauProducer"));
  desc.add<edm::InputTag>("PFTauSVATag", edm::InputTag("PFTauSecondaryVertexProducer"));
  descriptions.add("PFTauTransverseImpactParameters", desc);
}

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_2018_cff::InputTag.

produce()

void PFTauTransverseImpactParameters::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 101 of file PFTauTransverseImpactParameters.cc.

                                                                                           {
  // Obtain Collections
  edm::ESHandle<TransientTrackBuilder> transTrackBuilder;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", transTrackBuilder);
 
  edm::Handle<std::vector<reco::PFTau>> Tau;
  iEvent.getByToken(PFTauToken_, Tau);
 
  edm::Handle<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef>>> PFTauPVA;
  iEvent.getByToken(PFTauPVAToken_, PFTauPVA);
 
  edm::Handle<edm::AssociationVector<PFTauRefProd, std::vector<std::vector<reco::VertexRef>>>> PFTauSVA;
  iEvent.getByToken(PFTauSVAToken_, PFTauSVA);
 
  // Set Association Map
  auto AVPFTauTIP =
      std::make_unique<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef>>>(
          PFTauRefProd(Tau));
  auto TIPCollection_out = std::make_unique<PFTauTransverseImpactParameterCollection>();
  reco::PFTauTransverseImpactParameterRefProd TIPRefProd_out =
      iEvent.getRefBeforePut<reco::PFTauTransverseImpactParameterCollection>("PFTauTIP");
 
  // For each Tau Run Algorithim
  if (Tau.isValid()) {
    for (reco::PFTauCollection::size_type iPFTau = 0; iPFTau < Tau->size(); iPFTau++) {
      reco::PFTauRef RefPFTau(Tau, iPFTau);
      const reco::VertexRef PV = PFTauPVA->value(RefPFTau.key());
      const std::vector<reco::VertexRef> SV = PFTauSVA->value(RefPFTau.key());
      double dxy(-999), dxy_err(-999);
      reco::Vertex::Point poca(0, 0, 0);
      double ip3d(-999), ip3d_err(-999);
      reco::Vertex::Point ip3d_poca(0, 0, 0);
      if (RefPFTau->leadChargedHadrCand().isNonnull()) {
        const reco::Track* track = getTrack(*RefPFTau->leadChargedHadrCand());
        if (track != nullptr) {
          if (useFullCalculation_) {
            reco::TransientTrack transTrk = transTrackBuilder->build(*track);
            GlobalVector direction(
                RefPFTau->p4().px(), RefPFTau->p4().py(), RefPFTau->p4().pz());  //To compute sign of IP
            std::pair<bool, Measurement1D> signed_IP2D =
                IPTools::signedTransverseImpactParameter(transTrk, direction, (*PV));
            dxy = signed_IP2D.second.value();
            dxy_err = signed_IP2D.second.error();
            std::pair<bool, Measurement1D> signed_IP3D = IPTools::signedImpactParameter3D(transTrk, direction, (*PV));
            ip3d = signed_IP3D.second.value();
            ip3d_err = signed_IP3D.second.error();
            TransverseImpactPointExtrapolator extrapolator(transTrk.field());
            GlobalPoint pos =
                extrapolator.extrapolate(transTrk.impactPointState(), RecoVertex::convertPos(PV->position()))
                    .globalPosition();
            poca = reco::Vertex::Point(pos.x(), pos.y(), pos.z());
            AnalyticalImpactPointExtrapolator extrapolator3D(transTrk.field());
            GlobalPoint pos3d =
                extrapolator3D.extrapolate(transTrk.impactPointState(), RecoVertex::convertPos(PV->position()))
                    .globalPosition();
            ip3d_poca = reco::Vertex::Point(pos3d.x(), pos3d.y(), pos3d.z());
          } else {
            dxy_err = track->d0Error();
            dxy = track->dxy(PV->position());
            ip3d_err = track->dzError();       //store dz, ip3d not available
            ip3d = track->dz(PV->position());  //store dz, ip3d not available
          }
        }
      }
      if (!SV.empty()) {
        reco::Vertex::CovarianceMatrix cov;
        reco::Vertex::Point v(SV.at(0)->x() - PV->x(), SV.at(0)->y() - PV->y(), SV.at(0)->z() - PV->z());
        for (int i = 0; i < reco::Vertex::dimension; i++) {
          for (int j = 0; j < reco::Vertex::dimension; j++) {
            cov(i, j) = SV.at(0)->covariance(i, j) + PV->covariance(i, j);
          }
        }
        GlobalVector direction(RefPFTau->px(), RefPFTau->py(), RefPFTau->pz());
        double vSig = SecondaryVertex::computeDist3d(*PV, *SV.at(0), direction, true).significance();
        reco::PFTauTransverseImpactParameter TIPV(poca, dxy, dxy_err, ip3d_poca, ip3d, ip3d_err, PV, v, vSig, SV.at(0));
        reco::PFTauTransverseImpactParameterRef TIPVRef =
            reco::PFTauTransverseImpactParameterRef(TIPRefProd_out, TIPCollection_out->size());
        TIPCollection_out->push_back(TIPV);
        AVPFTauTIP->setValue(iPFTau, TIPVRef);
      } else {
        reco::PFTauTransverseImpactParameter TIPV(poca, dxy, dxy_err, ip3d_poca, ip3d, ip3d_err, PV);
        reco::PFTauTransverseImpactParameterRef TIPVRef =
            reco::PFTauTransverseImpactParameterRef(TIPRefProd_out, TIPCollection_out->size());
        TIPCollection_out->push_back(TIPV);
        AVPFTauTIP->setValue(iPFTau, TIPVRef);
      }
    }
  }
  iEvent.put(std::move(TIPCollection_out), "PFTauTIP");
  iEvent.put(std::move(AVPFTauTIP));
}

References TransientTrackBuilder::build(), RecoVertex::convertPos(), reco::Vertex::dimension, PVValHelper::dxy, reco::TransientTrack::field(), edm::EventSetup::get(), get, getTrack(), mps_fire::i, iEvent, reco::TransientTrack::impactPointState(), electrons_cff::ip3d, edm::Ref< C, T, F >::isNonnull(), dqmiolumiharvest::j, edm::Ref< C, T, F >::key(), eostools::move(), PFTauPVAToken_, PFTauSVAToken_, PFTauToken_, nanoDQM_cfi::PV, IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), met_cff::significance, nanoDQM_cfi::SV, L1TRate_Offline_cfi::Tau, HLT_2018_cff::track, useFullCalculation_, and findQualityFiles::v.

Member Data Documentation

PFTauPVAToken_

edm::EDGetTokenT<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > > PFTauTransverseImpactParameters::PFTauPVAToken_

private

Definition at line 64 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

PFTauSVAToken_

edm::EDGetTokenT<edm::AssociationVector<PFTauRefProd, std::vector<std::vector<reco::VertexRef> > > > PFTauTransverseImpactParameters::PFTauSVAToken_

private

Definition at line 65 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

PFTauToken_

edm::EDGetTokenT<std::vector<reco::PFTau> > PFTauTransverseImpactParameters::PFTauToken_

private

Definition at line 63 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

useFullCalculation_

bool PFTauTransverseImpactParameters::useFullCalculation_

private

Definition at line 66 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().