#include <PhysicsTools/VertexTableProducer/plugins/VertexTableProducer.cc>

Inheritance diagram for VertexTableProducer:

Public Member Functions
	VertexTableProducer (const edm::ParameterSet &)

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 Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const double	dlenMin_

const double	dlenSigMin_

const StringCutObjectSelector< reco::Vertex >	goodPvCut_

const std::string	goodPvCutString_

const edm::EDGetTokenT< pat::PackedCandidateCollection >	pfc_

const std::string	pvName_

const edm::EDGetTokenT< std::vector< reco::Vertex > >	pvs_

const edm::EDGetTokenT< edm::ValueMap< float > >	pvsScore_

const StringCutObjectSelector< reco::Candidate >	svCut_

const std::string	svDoc_

const std::string	svName_

const edm::EDGetTokenT< edm::View< reco::VertexCompositePtrCandidate > >	svs_

Additional Inherited Members
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

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 51 of file VertexTableProducer.cc.

Constructor & Destructor Documentation

◆ VertexTableProducer()

VertexTableProducer::VertexTableProducer ( const edm::ParameterSet & params )

explicit

Definition at line 78 of file VertexTableProducer.cc.

     : pvs_(consumes<std::vector<reco::Vertex>>(params.getParameter<edm::InputTag>("pvSrc"))),
       pfc_(consumes<pat::PackedCandidateCollection>(params.getParameter<edm::InputTag>("pfcSrc"))),
       pvsScore_(consumes<edm::ValueMap<float>>(params.getParameter<edm::InputTag>("pvSrc"))),
       svs_(consumes<edm::View<reco::VertexCompositePtrCandidate>>(params.getParameter<edm::InputTag>("svSrc"))),
       svCut_(params.getParameter<std::string>("svCut"), true),
       goodPvCut_(params.getParameter<std::string>("goodPvCut"), true),
       goodPvCutString_(params.getParameter<std::string>("goodPvCut")),
       pvName_(params.getParameter<std::string>("pvName")),
       svName_(params.getParameter<std::string>("svName")),
       svDoc_(params.getParameter<std::string>("svDoc")),
       dlenMin_(params.getParameter<double>("dlenMin")),
       dlenSigMin_(params.getParameter<double>("dlenSigMin"))
 
 {
   produces<nanoaod::FlatTable>("pv");
   produces<nanoaod::FlatTable>("otherPVs");
   produces<nanoaod::FlatTable>("svs");
   produces<edm::PtrVector<reco::VertexCompositePtrCandidate>>();
 }

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 222 of file VertexTableProducer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                      {
   edm::ParameterSetDescription desc;
 
   desc.add<edm::InputTag>("pvSrc")->setComment(
       "std::vector<reco::Vertex> and ValueMap<float> primary vertex input collections");
   desc.add<edm::InputTag>("pfcSrc")->setComment("packedPFCandidates input collections");
   desc.add<std::string>("goodPvCut")->setComment("selection on the primary vertex");
   desc.add<edm::InputTag>("svSrc")->setComment(
       "reco::VertexCompositePtrCandidate compatible secondary vertex input collection");
   desc.add<std::string>("svCut")->setComment("selection on the secondary vertex");
 
   desc.add<double>("dlenMin")->setComment("minimum value of dl to select secondary vertex");
   desc.add<double>("dlenSigMin")->setComment("minimum value of dl significance to select secondary vertex");
 
   desc.add<std::string>("pvName")->setComment("name of the flat table ouput");
   desc.add<std::string>("svName")->setComment("name of the flat table ouput");
   desc.add<std::string>("svDoc")->setComment("a few words of documentation");
 
   descriptions.addWithDefaultLabel(desc);
 }

◆ produce()

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

overrideprivate

Definition at line 105 of file VertexTableProducer.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, reco::Candidate::charge(), RecoVertex::convertError(), RecoVertex::convertPos(), VertexDistanceXY::distance(), VertexDistance3D::distance(), dlenMin_, dlenSigMin_, PVValHelper::dx, PVValHelper::dxy, PVValHelper::dy, PVValHelper::dz, goodPvCut_, goodPvCutString_, AlcaRecoSelection_cff::goodPVs, mps_fire::i, l1ctLayer2EG_cff::id, iEvent, dqmiolumiharvest::j, eostools::move(), getGTfromDQMFile::obj, pfc_, position, pvName_, pvs_, pvsScore_, Measurement1D::significance(), mathSSE::sqrt(), pfDeepBoostedJetPreprocessParams_cfi::sv, svCut_, svDoc_, svName_, svs_, Measurement1D::value(), and z.

                                                                                {
   using namespace edm;
   const auto& pvsScoreProd = iEvent.get(pvsScore_);
   auto pvsIn = iEvent.getHandle(pvs_);
 
   //pf candidates collection
   auto pfcIn = iEvent.getHandle(pfc_);
 
   auto pvTable = std::make_unique<nanoaod::FlatTable>(1, pvName_, true);
   pvTable->addColumnValue<float>("ndof", (*pvsIn)[0].ndof(), "main primary vertex number of degree of freedom", 8);
   pvTable->addColumnValue<float>("x", (*pvsIn)[0].position().x(), "main primary vertex position x coordinate", 10);
   pvTable->addColumnValue<float>("y", (*pvsIn)[0].position().y(), "main primary vertex position y coordinate", 10);
   pvTable->addColumnValue<float>("z", (*pvsIn)[0].position().z(), "main primary vertex position z coordinate", 16);
   pvTable->addColumnValue<float>("chi2", (*pvsIn)[0].normalizedChi2(), "main primary vertex reduced chi2", 8);
   int goodPVs = 0;
   for (const auto& pv : *pvsIn)
     if (goodPvCut_(pv))
       goodPVs++;
   pvTable->addColumnValue<uint8_t>("npvs", pvsIn->size(), "total number of reconstructed primary vertices");
   pvTable->addColumnValue<uint8_t>(
       "npvsGood", goodPVs, "number of good reconstructed primary vertices. selection:" + goodPvCutString_);
   pvTable->addColumnValue<float>(
       "score", pvsScoreProd.get(pvsIn.id(), 0), "main primary vertex score, i.e. sum pt2 of clustered objects", 8);
 
   float pv_sumpt2 = 0.0, pv_sumpx = 0.0, pv_sumpy = 0.0;
   for (const auto& obj : *pfcIn) {
     if (obj.charge() == 0) {
       continue;
     }  // skip neutrals
     double dz = fabs(obj.dz((*pvsIn)[0].position()));
     bool include_pfc = false;
     if (dz < 0.2) {
       include_pfc = true;
       for (size_t j = 1; j < (*pvsIn).size(); j++) {
         double newdz = fabs(obj.dz((*pvsIn)[j].position()));
         if (newdz < dz) {
           include_pfc = false;
           break;
         }
       }  // this pf candidate belongs to other PV
     }
     if (include_pfc) {
       float pfc_pt = obj.pt();
       pv_sumpt2 += pfc_pt * pfc_pt;
       pv_sumpx += obj.px();
       pv_sumpy += obj.py();
     }
   }
   pvTable->addColumnValue<float>(
       "sumpt2", pv_sumpt2, "sum pt2 of pf charged candidates for the main primary vertex", 10);
   pvTable->addColumnValue<float>("sumpx", pv_sumpx, "sum px of pf charged candidates for the main primary vertex", 10);
   pvTable->addColumnValue<float>("sumpy", pv_sumpy, "sum py of pf charged candidates for the main primary vertex", 10);
 
   auto otherPVsTable =
       std::make_unique<nanoaod::FlatTable>((*pvsIn).size() > 4 ? 3 : (*pvsIn).size() - 1, "Other" + pvName_, false);
   std::vector<float> pvsz;
   std::vector<float> pvscores;
   for (size_t i = 1; i < (*pvsIn).size() && i < 4; i++) {
     pvsz.push_back((*pvsIn)[i].position().z());
     pvscores.push_back(pvsScoreProd.get(pvsIn.id(), i));
   }
   otherPVsTable->addColumn<float>("z", pvsz, "Z position of other primary vertices, excluding the main PV", 8);
   otherPVsTable->addColumn<float>("score", pvscores, "scores of other primary vertices, excluding the main PV", 8);
 
   const auto& svsProd = iEvent.get(svs_);
   auto selCandSv = std::make_unique<PtrVector<reco::VertexCompositePtrCandidate>>();
   std::vector<float> dlen, dlenSig, pAngle, dxy, dxySig;
   std::vector<int16_t> charge;
   VertexDistance3D vdist;
   VertexDistanceXY vdistXY;
 
   size_t i = 0;
   const auto& PV0 = pvsIn->front();
   for (const auto& sv : svsProd) {
     if (svCut_(sv)) {
       Measurement1D dl =
           vdist.distance(PV0, VertexState(RecoVertex::convertPos(sv.position()), RecoVertex::convertError(sv.error())));
       if (dl.value() > dlenMin_ and dl.significance() > dlenSigMin_) {
         dlen.push_back(dl.value());
         dlenSig.push_back(dl.significance());
         selCandSv->push_back(svsProd.ptrAt(i));
         double dx = (PV0.x() - sv.vx()), dy = (PV0.y() - sv.vy()), dz = (PV0.z() - sv.vz());
         double pdotv = (dx * sv.px() + dy * sv.py() + dz * sv.pz()) / sv.p() / sqrt(dx * dx + dy * dy + dz * dz);
         pAngle.push_back(std::acos(pdotv));
         Measurement1D d2d = vdistXY.distance(
             PV0, VertexState(RecoVertex::convertPos(sv.position()), RecoVertex::convertError(sv.error())));
         dxy.push_back(d2d.value());
         dxySig.push_back(d2d.significance());
 
         int sum_charge = 0;
         for (unsigned int id = 0; id < sv.numberOfDaughters(); ++id) {
           const reco::Candidate* daughter = sv.daughter(id);
           sum_charge += daughter->charge();
         }
         charge.push_back(sum_charge);
       }
     }
     i++;
   }
 
   auto svsTable = std::make_unique<nanoaod::FlatTable>(selCandSv->size(), svName_, false);
   svsTable->setDoc(svDoc_);
   // For SV we fill from here only stuff that cannot be created with the SimpleFlatTableProducer
   svsTable->addColumn<float>("dlen", dlen, "decay length in cm", 10);
   svsTable->addColumn<float>("dlenSig", dlenSig, "decay length significance", 10);
   svsTable->addColumn<float>("dxy", dxy, "2D decay length in cm", 10);
   svsTable->addColumn<float>("dxySig", dxySig, "2D decay length significance", 10);
   svsTable->addColumn<float>("pAngle", pAngle, "pointing angle, i.e. acos(p_SV * (SV - PV)) ", 10);
   svsTable->addColumn<int16_t>("charge", charge, "sum of the charge of the SV tracks", 10);
 
   iEvent.put(std::move(pvTable), "pv");
   iEvent.put(std::move(otherPVsTable), "otherPVs");
   iEvent.put(std::move(svsTable), "svs");
   iEvent.put(std::move(selCandSv));
 }

Member Data Documentation

◆ dlenMin_

const double VertexTableProducer::dlenMin_

private

Definition at line 72 of file VertexTableProducer.cc.

Referenced by produce().

◆ dlenSigMin_

const double VertexTableProducer::dlenSigMin_

private

Definition at line 72 of file VertexTableProducer.cc.

Referenced by produce().

◆ goodPvCut_

const StringCutObjectSelector<reco::Vertex> VertexTableProducer::goodPvCut_

private

Definition at line 67 of file VertexTableProducer.cc.

Referenced by produce().

◆ goodPvCutString_

const std::string VertexTableProducer::goodPvCutString_

private

Definition at line 68 of file VertexTableProducer.cc.

Referenced by produce().

◆ pfc_

const edm::EDGetTokenT<pat::PackedCandidateCollection> VertexTableProducer::pfc_

private

Definition at line 63 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvName_

const std::string VertexTableProducer::pvName_

private

Definition at line 69 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvs_

const edm::EDGetTokenT<std::vector<reco::Vertex> > VertexTableProducer::pvs_

private

Definition at line 62 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvsScore_

const edm::EDGetTokenT<edm::ValueMap<float> > VertexTableProducer::pvsScore_

private

Definition at line 64 of file VertexTableProducer.cc.

Referenced by produce().

◆ svCut_

const StringCutObjectSelector<reco::Candidate> VertexTableProducer::svCut_

private

Definition at line 66 of file VertexTableProducer.cc.

Referenced by produce().

◆ svDoc_

const std::string VertexTableProducer::svDoc_

private

Definition at line 71 of file VertexTableProducer.cc.

Referenced by produce().

◆ svName_

const std::string VertexTableProducer::svName_

private

Definition at line 70 of file VertexTableProducer.cc.

Referenced by produce().

◆ svs_

const edm::EDGetTokenT<edm::View<reco::VertexCompositePtrCandidate> > VertexTableProducer::svs_

private

Definition at line 65 of file VertexTableProducer.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ VertexTableProducer()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ dlenMin_

◆ dlenSigMin_

◆ goodPvCut_

◆ goodPvCutString_

◆ pfc_

◆ pvName_

◆ pvs_

◆ pvsScore_

◆ svCut_

◆ svDoc_

◆ svName_

◆ svs_