#include <TauJetSelectorForHLTTrackSeeding.h>

Inheritance diagram for TauJetSelectorForHLTTrackSeeding:

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

	~TauJetSelectorForHLTTrackSeeding ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

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

Private Attributes
const double	etaMaxCaloJet_

const double	etaMinCaloJet_

const double	fractionMaxChargedPUInCaloCone_

const double	fractionMinCaloInTauCone_

edm::EDGetTokenT < reco::CaloJetCollection >	inputCaloJetToken_

edm::EDGetTokenT < reco::TrackJetCollection >	inputTrackJetToken_

edm::EDGetTokenT < reco::TrackCollection >	inputTrackToken_

const double	isolationConeSize_

const int	nTrkMaxInCaloCone_

const double	ptMinCaloJet_

const double	ptTrkMaxInCaloCone_

const double	tauConeSize_

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

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

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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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 24 of file TauJetSelectorForHLTTrackSeeding.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file TauJetSelectorForHLTTrackSeeding.cc.

References edm::ParameterSet::getParameter(), inputCaloJetToken_, inputTrackJetToken_, and inputTrackToken_.

                                                                                                 :
   ptMinCaloJet_(iConfig.getParameter< double > ("ptMinCaloJet")),
   etaMinCaloJet_(iConfig.getParameter< double > ("etaMinCaloJet")),
   etaMaxCaloJet_(iConfig.getParameter< double > ("etaMaxCaloJet")),
   tauConeSize_(iConfig.getParameter< double > ("tauConeSize")),
   isolationConeSize_(iConfig.getParameter< double > ("isolationConeSize")),
   fractionMinCaloInTauCone_(iConfig.getParameter< double > ("fractionMinCaloInTauCone")),
   fractionMaxChargedPUInCaloCone_(iConfig.getParameter< double > ("fractionMaxChargedPUInCaloCone")),
   ptTrkMaxInCaloCone_(iConfig.getParameter< double > ("ptTrkMaxInCaloCone")),
   nTrkMaxInCaloCone_(iConfig.getParameter< int > ("nTrkMaxInCaloCone"))
 {
    //now do what ever initialization is needed
   inputTrackJetToken_ = consumes< reco::TrackJetCollection >(iConfig.getParameter< edm::InputTag > ("inputTrackJetTag"));
   inputCaloJetToken_  = consumes< reco::CaloJetCollection >(iConfig.getParameter< edm::InputTag > ("inputCaloJetTag"));
   inputTrackToken_    = consumes< reco::TrackCollection >(iConfig.getParameter< edm::InputTag > ("inputTrackTag"));
 
   produces<reco::TrackJetCollection>();
 }

TauJetSelectorForHLTTrackSeeding::~TauJetSelectorForHLTTrackSeeding ( )

Definition at line 26 of file TauJetSelectorForHLTTrackSeeding.cc.

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

Member Function Documentation

void TauJetSelectorForHLTTrackSeeding::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 144 of file TauJetSelectorForHLTTrackSeeding.cc.

144 {}

void TauJetSelectorForHLTTrackSeeding::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 148 of file TauJetSelectorForHLTTrackSeeding.cc.

148 {}

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 41 of file TauJetSelectorForHLTTrackSeeding.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, reco::deltaR2(), etaMaxCaloJet_, etaMinCaloJet_, fractionMaxChargedPUInCaloCone_, fractionMinCaloInTauCone_, edm::Event::getByToken(), inputCaloJetToken_, inputTrackJetToken_, inputTrackToken_, isolationConeSize_, nTrkMaxInCaloCone_, p4, phi, ptMinCaloJet_, ptTrkMaxInCaloCone_, edm::Event::put(), tauConeSize_, reco::btau::trackEta, reco::btau::trackPhi, and testEve_cfg::tracks.

 {
    std::auto_ptr< reco::TrackJetCollection > augmentedTrackJets (new reco::TrackJetCollection);
 
    edm::Handle<reco::TrackJetCollection> trackjets;
    iEvent.getByToken(inputTrackJetToken_, trackjets);
 
    for (reco::TrackJetCollection::const_iterator trackjet = trackjets->begin();
       trackjet != trackjets->end(); trackjet++) {
      augmentedTrackJets->push_back(*trackjet);
    }
 
    edm::Handle<reco::TrackCollection> tracks;
    iEvent.getByToken(inputTrackToken_,tracks);
 
    edm::Handle<reco::CaloJetCollection> calojets;
    iEvent.getByToken(inputCaloJetToken_,calojets);
 
    const double tauConeSize2       = tauConeSize_ * tauConeSize_;
    const double isolationConeSize2 = isolationConeSize_ * isolationConeSize_;
 
    for (reco::CaloJetCollection::const_iterator calojet = calojets->begin();
       calojet != calojets->end(); calojet++) {
 
      if ( calojet->pt() < ptMinCaloJet_ ) continue;
      double etaJet = calojet->eta();
      double phiJet = calojet->phi();
      if ( etaJet < etaMinCaloJet_ ) continue;
      if ( etaJet > etaMaxCaloJet_ ) continue;
 
      std::vector <CaloTowerPtr> const & theTowers = calojet->getCaloConstituents();
      double ptIn = 0.;
      double ptOut = 0.;
      for ( unsigned int itwr = 0; itwr < theTowers.size(); ++itwr ) { 
        double etaTwr = theTowers[itwr]->eta() - etaJet;
        double phiTwr = deltaPhi(theTowers[itwr]->phi(), phiJet);
        double deltaR2 = etaTwr*etaTwr + phiTwr*phiTwr;
        //std::cout << "Tower eta/phi/et : " << etaTwr << " " << phiTwr << " " << theTowers[itwr]->pt() << std::endl;
        if ( deltaR2 < tauConeSize2 ) { 
      ptIn += theTowers[itwr]->pt(); 
        } else if ( deltaR2 < isolationConeSize2 ) { 
      ptOut += theTowers[itwr]->pt(); 
        }
      }
      double ptTot = ptIn+ptOut;
      double fracIn = ptIn/ptTot;
 
      // We are looking for isolated tracks
      if ( fracIn < fractionMinCaloInTauCone_) continue;
 
      int ntrk = 0;
      double ptTrk = 0.;
 
      for (reco::TrackJetCollection::const_iterator trackjet = trackjets->begin();
       trackjet != trackjets->end(); trackjet++) {
        for (unsigned itr=0; itr<trackjet->numberOfTracks(); ++itr) { 
      edm::Ptr<reco::Track> track = trackjet->track(itr);
      double trackEta = track->eta() - etaJet;
      double trackPhi = deltaPhi(track->phi(), phiJet);
      double deltaR2 = trackEta*trackEta + trackPhi*trackPhi;
      if ( deltaR2 < isolationConeSize2 ) { 
        ntrk++; 
        ptTrk += track->pt();
      }
        }
      }
      // We are looking for calojets without signal tracks already in
      if ( ntrk > nTrkMaxInCaloCone_ ) continue;
      if ( ptTrk > ptTrkMaxInCaloCone_ ) continue;
 
      int ntrk2 = 0;
      double ptTrk2 = 0.;
 
      for (reco::TrackCollection::const_iterator track = tracks->begin();
       track != tracks->end(); track++) {
        double trackEta = track->eta() - etaJet;
        double trackPhi = deltaPhi(track->phi(), phiJet);
        double deltaR2 = trackEta*trackEta + trackPhi*trackPhi;
        if ( deltaR2 < isolationConeSize2 ) { 
      ntrk2++; 
      ptTrk2 += track->pt();
        }
      }
      // We are looking for signal jets, not PU jets
      double fractionChargedPU = ptTrk2/calojet->pt(); 
      if ( fractionChargedPU > fractionMaxChargedPUInCaloCone_ ) continue;
      /*
      std::cout << "Calo Jet " << calojet->pt() << " " << calojet->eta() 
            << " " << ptIn << " " << ptOut << " " << fracIn
            << " " << ptTrk << " " << ntrk 
            << " " << fractionChargedPU
            << std::endl;
      */
      math::XYZTLorentzVector p4(calojet->p4());
      math::XYZPoint vertex(calojet->vertex());
      augmentedTrackJets->push_back(reco::TrackJet(p4,vertex));
    }
 
    iEvent.put(augmentedTrackJets);
 }