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

Inheritance diagram for VertexTableProducer:

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

	~VertexTableProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 48 of file VertexTableProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 85 of file VertexTableProducer.cc.

                                                                      :
     pvs_(consumes<std::vector<reco::Vertex>>( params.getParameter<edm::InputTag>("pvSrc") )),
     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::Candidate> >();
  
 }

VertexTableProducer::~VertexTableProducer ( )

override

Definition at line 107 of file VertexTableProducer.cc.

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

Member Function Documentation

void VertexTableProducer::beginStream ( edm::StreamID )

overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 189 of file VertexTableProducer.cc.

190 {

191 }

void VertexTableProducer::endStream ( )

overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 195 of file VertexTableProducer.cc.

195 {

196 }

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

static

Definition at line 200 of file VertexTableProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

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

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

overrideprivate

Definition at line 124 of file VertexTableProducer.cc.

References EnergyCorrector::c, RecoVertex::convertError(), RecoVertex::convertPos(), VertexDistance3D::distance(), dlenMin_, dlenSigMin_, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, nanoaod::FlatTable::FloatColumn, edm::Event::getByToken(), goodPvCut_, goodPvCutString_, AlcaRecoSelection_cff::goodPVs, mps_fire::i, nanoaod::FlatTable::IntColumn, eostools::move(), edm::Event::put(), MetAnalyzer::pv(), pvName_, pvs_, pvsScore_, Measurement1D::significance(), svCut_, svName_, svs_, and Measurement1D::value().

 {
     using namespace edm;
     edm::Handle<edm::ValueMap<float>> pvsScoreIn;
     edm::Handle<std::vector<reco::Vertex>> pvsIn;
     iEvent.getByToken(pvs_, pvsIn);
     iEvent.getByToken(pvsScore_, pvsScoreIn);
     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",nanoaod::FlatTable::FloatColumn,8);
     pvTable->addColumnValue<float>("x",(*pvsIn)[0].position().x(),"main primary vertex position x coordinate",nanoaod::FlatTable::FloatColumn,10);
     pvTable->addColumnValue<float>("y",(*pvsIn)[0].position().y(),"main primary vertex position y coordinate",nanoaod::FlatTable::FloatColumn,10);
     pvTable->addColumnValue<float>("z",(*pvsIn)[0].position().z(),"main primary vertex position z coordinate",nanoaod::FlatTable::FloatColumn,16);
     pvTable->addColumnValue<float>("chi2",(*pvsIn)[0].normalizedChi2(),"main primary vertex reduced chi2",nanoaod::FlatTable::FloatColumn,8);
     int goodPVs=0;
     for (const auto & pv : *pvsIn) if (goodPvCut_(pv)) goodPVs++;
     pvTable->addColumnValue<int>("npvs",(*pvsIn).size(),"total number of reconstructed primary vertices",nanoaod::FlatTable::IntColumn);
     pvTable->addColumnValue<int>("npvsGood",goodPVs,"number of good reconstructed primary vertices. selection:"+goodPvCutString_,nanoaod::FlatTable::IntColumn);
     pvTable->addColumnValue<float>("score",(*pvsScoreIn).get(pvsIn.id(),0),"main primary vertex score, i.e. sum pt2 of clustered objects",nanoaod::FlatTable::FloatColumn,8);
 
     auto otherPVsTable = std::make_unique<nanoaod::FlatTable>((*pvsIn).size() >4?3:(*pvsIn).size()-1,"Other"+pvName_,false);
     std::vector<float> pvsz;
     for(size_t i=1;i < (*pvsIn).size() && i < 4; i++) pvsz.push_back((*pvsIn)[i-1].position().z());
     otherPVsTable->addColumn<float>("z",pvsz,"Z position of other primary vertices, excluding the main PV",nanoaod::FlatTable::FloatColumn,8);
 
 
     edm::Handle<edm::View<reco::VertexCompositePtrCandidate> > svsIn;
     iEvent.getByToken(svs_, svsIn);
     auto selCandSv = std::make_unique<PtrVector<reco::Candidate>>();
     std::vector<float> dlen,dlenSig,pAngle;
     VertexDistance3D vdist;
 
     size_t i=0;
     const auto & PV0 = pvsIn->front();
     for (const auto & sv : *svsIn) {
        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());
                edm::Ptr<reco::Candidate> c =  svsIn->ptrAt(i);
                selCandSv->push_back(c);
                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();
                pAngle.push_back(std::acos(pdotv));
            }
        }
        i++;
     }
 
 
     auto svsTable = std::make_unique<nanoaod::FlatTable>(selCandSv->size(),svName_,false);
     // For SV we fill from here only stuff that cannot be created with the SimpleFlatTableProducer 
     svsTable->addColumn<float>("dlen",dlen,"decay length in cm",nanoaod::FlatTable::FloatColumn,10);
     svsTable->addColumn<float>("dlenSig",dlenSig,"decay length significance",nanoaod::FlatTable::FloatColumn, 10);
     svsTable->addColumn<float>("pAngle",pAngle,"pointing angle, i.e. acos(p_SV * (SV - PV)) ",nanoaod::FlatTable::FloatColumn,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

const double VertexTableProducer::dlenMin_

private

Definition at line 76 of file VertexTableProducer.cc.

Referenced by produce().

const double VertexTableProducer::dlenSigMin_

private

Definition at line 76 of file VertexTableProducer.cc.

Referenced by produce().

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

private

Definition at line 71 of file VertexTableProducer.cc.

Referenced by produce().

const std::string VertexTableProducer::goodPvCutString_

private

Definition at line 72 of file VertexTableProducer.cc.

Referenced by produce().

const std::string VertexTableProducer::pvName_

private

Definition at line 73 of file VertexTableProducer.cc.

Referenced by produce().

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

private

Definition at line 67 of file VertexTableProducer.cc.

Referenced by produce().

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

private

Definition at line 68 of file VertexTableProducer.cc.

Referenced by produce().

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

private

Definition at line 70 of file VertexTableProducer.cc.

Referenced by produce().

const std::string VertexTableProducer::svDoc_

private

Definition at line 75 of file VertexTableProducer.cc.

const std::string VertexTableProducer::svName_

private

Definition at line 74 of file VertexTableProducer.cc.

Referenced by produce().

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

private

Definition at line 69 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

Member Function Documentation

Member Data Documentation