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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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(), submitPVResolutionJobs::desc, and HLT_FULL_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_FULL_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().