PFTauTransverseImpactParameters Class Reference

Inheritance diagram for PFTauTransverseImpactParameters:

Public Types
enum	CMSSWPerigee {   aCurv = 0, aTheta, aPhi, aTip,   aLip }
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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_
const edm::ESGetToken < TransientTrackBuilder, TransientTrackRecord >	transTrackBuilderToken_
bool	useFullCalculation_

Detailed Description

Definition at line 53 of file PFTauTransverseImpactParameters.cc.

Member Enumeration Documentation

enum PFTauTransverseImpactParameters::CMSSWPerigee

Enumerator
aCurv
aTheta
aPhi
aTip
aLip

Definition at line 55 of file PFTauTransverseImpactParameters.cc.

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

Constructor & Destructor Documentation

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"))),
      transTrackBuilderToken_(esConsumes(edm::ESInputTag{"", "TransientTrackBuilder"})),
      useFullCalculation_(iConfig.getParameter<bool>("useFullCalculation")) {
  produces<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef>>>();
  produces<PFTauTransverseImpactParameterCollection>("PFTauTIP");
}

PFTauTransverseImpactParameters::~PFTauTransverseImpactParameters ( )

override

Definition at line 81 of file PFTauTransverseImpactParameters.cc.

{}

Member Function Documentation

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

static

Definition at line 193 of file PFTauTransverseImpactParameters.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

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

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

override

Definition at line 102 of file PFTauTransverseImpactParameters.cc.

References RecoVertex::convertPos(), reco::TrackBase::d0Error(), reco::Vertex::dimension, PVValHelper::dxy, reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::TransientTrack::field(), edm::Event::getByToken(), edm::EventSetup::getData(), edm::Event::getRefBeforePut(), getTrack(), mps_fire::i, reco::TransientTrack::impactPointState(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), dqmiolumiharvest::j, edm::Ref< C, T, F >::key(), eostools::move(), PFTauPVAToken_, PFTauSVAToken_, PFTauToken_, edm::Event::put(), IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), HLT_FULL_cff::Tau, HLT_FULL_cff::track, transTrackBuilderToken_, useFullCalculation_, and findQualityFiles::v.

                                                                                           {
  // Obtain Collections
  auto const& transTrackBuilder = iSetup.getData(transTrackBuilderToken_);

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

Member Data Documentation

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

private

Definition at line 63 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

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

private

Definition at line 64 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

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

private

Definition at line 62 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> PFTauTransverseImpactParameters::transTrackBuilderToken_

private

Definition at line 65 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

bool PFTauTransverseImpactParameters::useFullCalculation_

private

Definition at line 66 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().