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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Implements edm::EDProducer.

Definition at line 21 of file DQMTauProducer.cc.

References MuTauSkim_cff::discriminator, edm::Event::getByToken(), i, isolationCone_, HLT_25ns10e33_v2_cff::jetTracks, matchingCone_, nTracks(), ptMin_, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, 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::auto_ptr<reco::HLTTauCollection> selectedTaus(jetCollection);
   
   iEvent.put(selectedTaus);
   
 
 
 
 }