PFTauSecondaryVertexProducer Class Reference

Inheritance diagram for PFTauSecondaryVertexProducer:

Public Types
enum	Alg { useInputPV =0, usePVwithMaxSumPt, useTauPV }
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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

Public Member Functions
	PFTauSecondaryVertexProducer (const edm::ParameterSet &iConfig)
void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
	~PFTauSecondaryVertexProducer () override
Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::global::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)

Private Attributes
const edm::InputTag	PFTauTag_
const edm::EDGetTokenT< std::vector< reco::PFTau > >	PFTauToken_

Additional Inherited Members
Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Definition at line 50 of file PFTauSecondaryVertexProducer.cc.

Member Enumeration Documentation

enum PFTauSecondaryVertexProducer::Alg

Enumerator
useInputPV
usePVwithMaxSumPt
useTauPV

Definition at line 52 of file PFTauSecondaryVertexProducer.cc.

{useInputPV=0, usePVwithMaxSumPt, useTauPV};

Constructor & Destructor Documentation

PFTauSecondaryVertexProducer::PFTauSecondaryVertexProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 62 of file PFTauSecondaryVertexProducer.cc.

                                                                                        :
  PFTauTag_(iConfig.getParameter<edm::InputTag>("PFTauTag")),
  PFTauToken_(consumes<std::vector<reco::PFTau> >(iConfig.getParameter<edm::InputTag>("PFTauTag")))
{
  produces<edm::AssociationVector<PFTauRefProd, std::vector<std::vector<reco::VertexRef> > > >();
  produces<VertexCollection>("PFTauSecondaryVertices");
}

PFTauSecondaryVertexProducer::~PFTauSecondaryVertexProducer ( )

override

Definition at line 70 of file PFTauSecondaryVertexProducer.cc.

                                                           {

}

Member Function Documentation

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

override

Definition at line 74 of file PFTauSecondaryVertexProducer.cc.

References TransientTrackBuilder::build(), DEFINE_FWK_MODULE, edm::EventSetup::get(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), TransientVertex::hasRefittedTracks(), edm::HandleBase::isValid(), eostools::move(), PFTauToken_, edm::Event::put(), TransientVertex::refittedTracks(), nanoDQM_cfi::SV, and KalmanVertexFitter::vertex().

                                                                                                          {
  // Obtain 
  edm::ESHandle<TransientTrackBuilder> transTrackBuilder;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",transTrackBuilder);

  edm::Handle<std::vector<reco::PFTau> > Tau;
  iEvent.getByToken(PFTauToken_,Tau);

  // Set Association Map
  auto AVPFTauSV = std::make_unique<edm::AssociationVector<PFTauRefProd, std::vector<std::vector<reco::VertexRef>>>>(PFTauRefProd(Tau));
  auto VertexCollection_out = std::make_unique<VertexCollection>();
  reco::VertexRefProd VertexRefProd_out = iEvent.getRefBeforePut<reco::VertexCollection>("PFTauSecondaryVertices");

  // For each Tau Run Algorithim
  if(Tau.isValid()) {
    for(reco::PFTauCollection::size_type iPFTau = 0; iPFTau < Tau->size(); iPFTau++) {
      reco::PFTauRef RefPFTau(Tau, iPFTau);
      std::vector<reco::VertexRef> SV;
      if(RefPFTau->decayMode()>=5){
    // Get tracks form PFTau daugthers
    std::vector<reco::TransientTrack> transTrk;
    TransientVertex transVtx;
    const std::vector<edm::Ptr<reco::PFCandidate> > cands = RefPFTau->signalPFChargedHadrCands(); 
    for (std::vector<edm::Ptr<reco::PFCandidate> >::const_iterator iter = cands.begin(); iter!=cands.end(); ++iter) {
      if(iter->get()->trackRef().isNonnull())transTrk.push_back(transTrackBuilder->build(iter->get()->trackRef()));
      else if(iter->get()->gsfTrackRef().isNonnull())transTrk.push_back(transTrackBuilder->build(iter->get()->gsfTrackRef()));
    }
    // Fit the secondary vertex
    bool FitOk(true);
    KalmanVertexFitter kvf(true);
    try{
      transVtx = kvf.vertex(transTrk); //KalmanVertexFitter  
    }catch(...){
      FitOk=false;
    }
    if(!transVtx.hasRefittedTracks()) FitOk=false;
    if(transVtx.refittedTracks().size()!=transTrk.size()) FitOk=false;
    if(FitOk){
      SV.push_back(reco::VertexRef(VertexRefProd_out, VertexCollection_out->size()));
      VertexCollection_out->push_back(transVtx);
    }
      }
      AVPFTauSV->setValue(iPFTau, SV);
    }
  }
  iEvent.put(std::move(VertexCollection_out),"PFTauSecondaryVertices");
  iEvent.put(std::move(AVPFTauSV));
}

Member Data Documentation

const edm::InputTag PFTauSecondaryVertexProducer::PFTauTag_

private

Definition at line 58 of file PFTauSecondaryVertexProducer.cc.

const edm::EDGetTokenT<std::vector<reco::PFTau> > PFTauSecondaryVertexProducer::PFTauToken_

private

Definition at line 59 of file PFTauSecondaryVertexProducer.cc.

Referenced by produce().