CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
VertexTableProducer Class Reference

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

Inheritance diagram for VertexTableProducer:
edm::stream::EDProducer<>

Public Member Functions

 VertexTableProducer (const edm::ParameterSet &)
 
 ~VertexTableProducer () 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 Member Functions

void beginStream (edm::StreamID) override
 
void endStream () override
 
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 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<>
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
 

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 48 of file VertexTableProducer.cc.

Constructor & Destructor Documentation

◆ VertexTableProducer()

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

Definition at line 82 of file VertexTableProducer.cc.

83  : pvs_(consumes<std::vector<reco::Vertex>>(params.getParameter<edm::InputTag>("pvSrc"))),
84  pvsScore_(consumes<edm::ValueMap<float>>(params.getParameter<edm::InputTag>("pvSrc"))),
85  svs_(consumes<edm::View<reco::VertexCompositePtrCandidate>>(params.getParameter<edm::InputTag>("svSrc"))),
86  svCut_(params.getParameter<std::string>("svCut"), true),
87  goodPvCut_(params.getParameter<std::string>("goodPvCut"), true),
88  goodPvCutString_(params.getParameter<std::string>("goodPvCut")),
89  pvName_(params.getParameter<std::string>("pvName")),
90  svName_(params.getParameter<std::string>("svName")),
91  svDoc_(params.getParameter<std::string>("svDoc")),
92  dlenMin_(params.getParameter<double>("dlenMin")),
93  dlenSigMin_(params.getParameter<double>("dlenSigMin"))
94 
95 {
96  produces<nanoaod::FlatTable>("pv");
97  produces<nanoaod::FlatTable>("otherPVs");
98  produces<nanoaod::FlatTable>("svs");
99  produces<edm::PtrVector<reco::Candidate>>();
100 }

◆ ~VertexTableProducer()

VertexTableProducer::~VertexTableProducer ( )
override

Definition at line 102 of file VertexTableProducer.cc.

102  {
103  // do anything here that needs to be done at destruction time
104  // (e.g. close files, deallocate resources etc.)
105 }

Member Function Documentation

◆ beginStream()

void VertexTableProducer::beginStream ( edm::StreamID  )
overrideprivate

Definition at line 187 of file VertexTableProducer.cc.

187 {}

◆ endStream()

void VertexTableProducer::endStream ( )
overrideprivate

Definition at line 190 of file VertexTableProducer.cc.

190 {}

◆ fillDescriptions()

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

Definition at line 193 of file VertexTableProducer.cc.

193  {
194  //The following says we do not know what parameters are allowed so do no validation
195  // Please change this to state exactly what you do use, even if it is no parameters
196  edm::ParameterSetDescription desc;
197  desc.setUnknown();
198  descriptions.addDefault(desc);
199 }

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

◆ produce()

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

Definition at line 113 of file VertexTableProducer.cc.

113  {
114  using namespace edm;
115  edm::Handle<edm::ValueMap<float>> pvsScoreIn;
116  edm::Handle<std::vector<reco::Vertex>> pvsIn;
117  iEvent.getByToken(pvs_, pvsIn);
118  iEvent.getByToken(pvsScore_, pvsScoreIn);
119  auto pvTable = std::make_unique<nanoaod::FlatTable>(1, pvName_, true);
120  pvTable->addColumnValue<float>("ndof", (*pvsIn)[0].ndof(), "main primary vertex number of degree of freedom", 8);
121  pvTable->addColumnValue<float>("x", (*pvsIn)[0].position().x(), "main primary vertex position x coordinate", 10);
122  pvTable->addColumnValue<float>("y", (*pvsIn)[0].position().y(), "main primary vertex position y coordinate", 10);
123  pvTable->addColumnValue<float>("z", (*pvsIn)[0].position().z(), "main primary vertex position z coordinate", 16);
124  pvTable->addColumnValue<float>("chi2", (*pvsIn)[0].normalizedChi2(), "main primary vertex reduced chi2", 8);
125  int goodPVs = 0;
126  for (const auto& pv : *pvsIn)
127  if (goodPvCut_(pv))
128  goodPVs++;
129  pvTable->addColumnValue<int>("npvs", pvsIn->size(), "total number of reconstructed primary vertices");
130  pvTable->addColumnValue<int>(
131  "npvsGood", goodPVs, "number of good reconstructed primary vertices. selection:" + goodPvCutString_);
132  pvTable->addColumnValue<float>(
133  "score", pvsScoreIn->get(pvsIn.id(), 0), "main primary vertex score, i.e. sum pt2 of clustered objects", 8);
134 
135  auto otherPVsTable =
136  std::make_unique<nanoaod::FlatTable>((*pvsIn).size() > 4 ? 3 : (*pvsIn).size() - 1, "Other" + pvName_, false);
137  std::vector<float> pvsz;
138  for (size_t i = 1; i < (*pvsIn).size() && i < 4; i++)
139  pvsz.push_back((*pvsIn)[i - 1].position().z());
140  otherPVsTable->addColumn<float>("z", pvsz, "Z position of other primary vertices, excluding the main PV", 8);
141 
142  edm::Handle<edm::View<reco::VertexCompositePtrCandidate>> svsIn;
143  iEvent.getByToken(svs_, svsIn);
144  auto selCandSv = std::make_unique<PtrVector<reco::Candidate>>();
145  std::vector<float> dlen, dlenSig, pAngle, dxy, dxySig;
146  VertexDistance3D vdist;
147  VertexDistanceXY vdistXY;
148 
149  size_t i = 0;
150  const auto& PV0 = pvsIn->front();
151  for (const auto& sv : *svsIn) {
152  if (svCut_(sv)) {
153  Measurement1D dl =
154  vdist.distance(PV0, VertexState(RecoVertex::convertPos(sv.position()), RecoVertex::convertError(sv.error())));
155  if (dl.value() > dlenMin_ and dl.significance() > dlenSigMin_) {
156  dlen.push_back(dl.value());
157  dlenSig.push_back(dl.significance());
158  edm::Ptr<reco::Candidate> c = svsIn->ptrAt(i);
159  selCandSv->push_back(c);
160  double dx = (PV0.x() - sv.vx()), dy = (PV0.y() - sv.vy()), dz = (PV0.z() - sv.vz());
161  double pdotv = (dx * sv.px() + dy * sv.py() + dz * sv.pz()) / sv.p() / sqrt(dx * dx + dy * dy + dz * dz);
162  pAngle.push_back(std::acos(pdotv));
163  Measurement1D d2d = vdistXY.distance(
164  PV0, VertexState(RecoVertex::convertPos(sv.position()), RecoVertex::convertError(sv.error())));
165  dxy.push_back(d2d.value());
166  dxySig.push_back(d2d.significance());
167  }
168  }
169  i++;
170  }
171 
172  auto svsTable = std::make_unique<nanoaod::FlatTable>(selCandSv->size(), svName_, false);
173  // For SV we fill from here only stuff that cannot be created with the SimpleFlatTableProducer
174  svsTable->addColumn<float>("dlen", dlen, "decay length in cm", 10);
175  svsTable->addColumn<float>("dlenSig", dlenSig, "decay length significance", 10);
176  svsTable->addColumn<float>("dxy", dxy, "2D decay length in cm", 10);
177  svsTable->addColumn<float>("dxySig", dxySig, "2D decay length significance", 10);
178  svsTable->addColumn<float>("pAngle", pAngle, "pointing angle, i.e. acos(p_SV * (SV - PV)) ", 10);
179 
180  iEvent.put(std::move(pvTable), "pv");
181  iEvent.put(std::move(otherPVsTable), "otherPVs");
182  iEvent.put(std::move(svsTable), "svs");
183  iEvent.put(std::move(selCandSv));
184 }

References HltBtagPostValidation_cff::c, 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, iEvent, eostools::move(), MetAnalyzer::pv(), pvName_, pvs_, pvsScore_, Measurement1D::significance(), mathSSE::sqrt(), pfDeepBoostedJetPreprocessParams_cfi::sv, svCut_, svName_, svs_, and Measurement1D::value().

Member Data Documentation

◆ dlenMin_

const double VertexTableProducer::dlenMin_
private

Definition at line 76 of file VertexTableProducer.cc.

Referenced by produce().

◆ dlenSigMin_

const double VertexTableProducer::dlenSigMin_
private

Definition at line 76 of file VertexTableProducer.cc.

Referenced by produce().

◆ goodPvCut_

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

Definition at line 71 of file VertexTableProducer.cc.

Referenced by produce().

◆ goodPvCutString_

const std::string VertexTableProducer::goodPvCutString_
private

Definition at line 72 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvName_

const std::string VertexTableProducer::pvName_
private

Definition at line 73 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvs_

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

Definition at line 67 of file VertexTableProducer.cc.

Referenced by produce().

◆ pvsScore_

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

Definition at line 68 of file VertexTableProducer.cc.

Referenced by produce().

◆ svCut_

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

Definition at line 70 of file VertexTableProducer.cc.

Referenced by produce().

◆ svDoc_

const std::string VertexTableProducer::svDoc_
private

Definition at line 75 of file VertexTableProducer.cc.

◆ svName_

const std::string VertexTableProducer::svName_
private

Definition at line 74 of file VertexTableProducer.cc.

Referenced by produce().

◆ svs_

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

Definition at line 69 of file VertexTableProducer.cc.

Referenced by produce().

VertexTableProducer::goodPvCut_
const StringCutObjectSelector< reco::Vertex > goodPvCut_
Definition: VertexTableProducer.cc:71
Measurement1D
Definition: Measurement1D.h:11
mps_fire.i
i
Definition: mps_fire.py:428
CalibrationSummaryClient_cfi.params
params
Definition: CalibrationSummaryClient_cfi.py:14
edm
HLT enums.
Definition: AlignableModifier.h:19
Measurement1D::value
double value() const
Definition: Measurement1D.h:25
VertexDistance3D::distance
Measurement1D distance(const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const override
Definition: VertexDistance3D.cc:17
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
VertexTableProducer::pvs_
const edm::EDGetTokenT< std::vector< reco::Vertex > > pvs_
Definition: VertexTableProducer.cc:67
VertexTableProducer::svCut_
const StringCutObjectSelector< reco::Candidate > svCut_
Definition: VertexTableProducer.cc:70
VertexTableProducer::svs_
const edm::EDGetTokenT< edm::View< reco::VertexCompositePtrCandidate > > svs_
Definition: VertexTableProducer.cc:69
edm::Handle
Definition: AssociativeIterator.h:50
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
pfDeepBoostedJetPreprocessParams_cfi.sv
sv
Definition: pfDeepBoostedJetPreprocessParams_cfi.py:352
VertexDistance3D
Definition: VertexDistance3D.h:13
VertexTableProducer::dlenSigMin_
const double dlenSigMin_
Definition: VertexTableProducer.cc:76
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
VertexTableProducer::svName_
const std::string svName_
Definition: VertexTableProducer.cc:74
edm::View
Definition: CaloClusterFwd.h:14
Measurement1D::significance
double significance() const
Definition: Measurement1D.h:29
RecoVertex::convertError
reco::Vertex::Error convertError(const GlobalError &ge)
Definition: ConvertError.h:8
VertexTableProducer::pvName_
const std::string pvName_
Definition: VertexTableProducer.cc:73
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
iEvent
int iEvent
Definition: GenABIO.cc:224
PVValHelper::dy
Definition: PVValidationHelpers.h:49
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
VertexDistanceXY
Definition: VertexDistanceXY.h:11
VertexTableProducer::dlenMin_
const double dlenMin_
Definition: VertexTableProducer.cc:76
edm::Ptr< reco::Candidate >
VertexTableProducer::goodPvCutString_
const std::string goodPvCutString_
Definition: VertexTableProducer.cc:72
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
AlcaRecoSelection_cff.goodPVs
goodPVs
Good Primary Vertex Selection.
Definition: AlcaRecoSelection_cff.py:54
eostools.move
def move(src, dest)
Definition: eostools.py:511
VertexDistanceXY::distance
Measurement1D distance(const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const override
Definition: VertexDistanceXY.cc:19
PVValHelper::dxy
Definition: PVValidationHelpers.h:47
PVValHelper::dz
Definition: PVValidationHelpers.h:50
VertexState
Definition: VertexState.h:13
edm::ValueMap< float >
RecoVertex::convertPos
reco::Vertex::Point convertPos(const GlobalPoint &p)
Definition: ConvertToFromReco.h:7
VertexTableProducer::pvsScore_
const edm::EDGetTokenT< edm::ValueMap< float > > pvsScore_
Definition: VertexTableProducer.cc:68
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
edm::InputTag
Definition: InputTag.h:15
PVValHelper::dx
Definition: PVValidationHelpers.h:48
VertexTableProducer::svDoc_
const std::string svDoc_
Definition: VertexTableProducer.cc:75
Generated for CMSSW Reference Manual by doxygen 1.8.16