#include <DQMTauProducer.h>

Inheritance diagram for DQMTauProducer:

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

virtual void	produce (edm::Event &, const edm::EventSetup &) override

	~DQMTauProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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 ()

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > 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
double	isolationCone_

double	matchingCone_

double	ptMin_

double	ptMinLeadTk_

double	rmax_

double	rmin_

double	signalCone_

edm::EDGetTokenT< reco::IsolatedTauTagInfoCollection >	trackIsolatedJets_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
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 18 of file DQMTauProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 8 of file DQMTauProducer.cc.

References edm::ParameterSet::getParameter(), isolationCone_, matchingCone_, ptMin_, signalCone_, and trackIsolatedJets_.

 {
   trackIsolatedJets_ = consumes<reco::IsolatedTauTagInfoCollection>(iConfig.getParameter<edm::InputTag>("TrackIsoJets"));
   matchingCone_ = iConfig.getParameter<double>("MatchingCone");
   signalCone_ = iConfig.getParameter<double>("SignalCone");
   ptMin_ = iConfig.getParameter<double>("MinPtTracks");
   
   isolationCone_ = iConfig.getParameter<double>("IsolationCone");
   produces<reco::HLTTauCollection>();
 }

DQMTauProducer::~DQMTauProducer ( )

Definition at line 19 of file DQMTauProducer.cc.

19 { }

Member Function Documentation

void DQMTauProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iES
	)

overridevirtual

Definition at line 21 of file DQMTauProducer.cc.

References mvaPFMET_leptons_cfi::discriminator, edm::Event::getByToken(), mps_fire::i, createfilelist::int, isolationCone_, impactParameterTagInfos_cfi::jetTracks, matchingCone_, eostools::move(), nTracks(), edm::Handle< T >::product(), ptMin_, edm::Event::put(), objectSelection_cff::selectedTaus, reco::HLTTau::setNL25TrackIsolation(), reco::HLTTau::setNL3TrackIsolation(), signalCone_, metsig::tau, and trackIsolatedJets_.

 {
 
   using namespace reco;
   using namespace edm;
   using namespace std;
   
 
   HLTTauCollection * jetCollection = new HLTTauCollection;
  
 
   edm::Handle<IsolatedTauTagInfoCollection> tauL25Jets;
   iEvent.getByToken(trackIsolatedJets_, tauL25Jets );
   
 
   IsolatedTauTagInfoCollection tau = *(tauL25Jets.product());
 
   
   float eta_, phi_, pt_;
   float ptLeadTk=0.;
   int trackIsolation=1000.;  
   int nTracks=1000.;
 
  for(unsigned int i=0; i<tau.size();i++)
    {
      JetTracksAssociationRef jetTracks = tau[i].jtaRef();
      math::XYZVector jetDir(jetTracks->first->px(),jetTracks->first->py(),jetTracks->first->pz());   
      eta_ = jetDir.eta();
      phi_ = jetDir.phi();
      pt_ = jetTracks->first->pt();
      
      const TrackRef leadTk = tau[i].leadingSignalTrack(jetDir,matchingCone_, 1.);
      if(!leadTk) 
        {}else{
      trackIsolation = (int)tau[i].discriminator(jetDir,matchingCone_, signalCone_, isolationCone_,1.,1.,0);
      ptLeadTk = (*leadTk).pt();
      nTracks = (tau[i].tracksInCone((*leadTk).momentum(), isolationCone_,   ptMin_ )).size() - (tau[i].tracksInCone((*leadTk).momentum(), signalCone_, ptMin_ )).size();
        }
      HLTTau pippo(eta_,phi_,pt_,-1,trackIsolation,ptLeadTk,trackIsolation,ptLeadTk);
      pippo.setNL25TrackIsolation(nTracks);
      pippo.setNL3TrackIsolation(nTracks);
      jetCollection->push_back(pippo);
      
    }
  
    
   std::unique_ptr<reco::HLTTauCollection> selectedTaus(jetCollection);
   
   iEvent.put(std::move(selectedTaus));
   
 
 
 
 }