#include <HLTDeDxFilter.h>

Inheritance diagram for HLTDeDxFilter:

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

virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	~HLTDeDxFilter ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default

Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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 (const std::string &iProcessName, std::vector< const char * > &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 ()

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

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::InputTag	caloTowersTag_

edm::EDGetToken	caloTowersToken_

edm::InputTag	inputdedxTag_

edm::EDGetToken	inputdedxToken_

edm::InputTag	inputTracksTag_

edm::EDGetToken	inputTracksToken_

double	maxAssocCaloE_

double	maxAssocCaloEDeltaRSize_

double	maxETA_

double	maxNHitMissIn_

double	maxNHitMissMid_

double	maxRelTrkIsoDeltaRp3_

double	minDEDx_

double	minNOM_

double	minNumValidHits_

double	minPT_

double	relTrkIsoDeltaRSize_

bool	saveTags_

edm::InputTag	thisModuleTag_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

Public Types inherited from edm::ProducerBase
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

Author: Claude Nuttens; Gheorghe Lungu

Definition at line 22 of file HLTDeDxFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file HLTDeDxFilter.cc.

References caloTowersTag_, caloTowersToken_, edm::ParameterSet::getParameter(), inputdedxTag_, inputdedxToken_, HLT_25ns14e33_v1_cff::InputTag, inputTracksTag_, inputTracksToken_, maxAssocCaloE_, maxAssocCaloEDeltaRSize_, maxETA_, maxNHitMissIn_, maxNHitMissMid_, maxRelTrkIsoDeltaRp3_, minDEDx_, minNOM_, minNumValidHits_, minPT_, relTrkIsoDeltaRSize_, AlCaHLTBitMon_QueryRunRegistry::string, and thisModuleTag_.

                                                            : HLTFilter(iConfig)
 {
   minDEDx_  = iConfig.getParameter<double> ("minDEDx");
   minPT_        = iConfig.getParameter<double> ("minPT");
   minNOM_       = iConfig.getParameter<double> ("minNOM");
   maxETA_       = iConfig.getParameter<double> ("maxETA");
   minNumValidHits_  = iConfig.getParameter<double> ("minNumValidHits");
   maxNHitMissIn_    = iConfig.getParameter<double> ("maxNHitMissIn");
   maxNHitMissMid_   = iConfig.getParameter<double> ("maxNHitMissMid");
   maxRelTrkIsoDeltaRp3_     = iConfig.getParameter<double> ("maxRelTrkIsoDeltaRp3");
   relTrkIsoDeltaRSize_     = iConfig.getParameter<double> ("relTrkIsoDeltaRSize");
   maxAssocCaloE_     = iConfig.getParameter<double> ("maxAssocCaloE");
   maxAssocCaloEDeltaRSize_     = iConfig.getParameter<double> ("maxAssocCaloEDeltaRSize");
   inputTracksTag_ = iConfig.getParameter< edm::InputTag > ("inputTracksTag");
   inputdedxTag_   = iConfig.getParameter< edm::InputTag > ("inputDeDxTag");
   caloTowersTag_ =  iConfig.getParameter<edm::InputTag>("caloTowersTag");
 
   if(maxAssocCaloE_ >= 0) caloTowersToken_ = consumes<CaloTowerCollection> (iConfig.getParameter<edm::InputTag>("caloTowersTag"));
   inputTracksToken_ = consumes<reco::TrackCollection>(iConfig.getParameter< edm::InputTag > ("inputTracksTag"));
   inputdedxToken_   = consumes<edm::ValueMap<reco::DeDxData> >(iConfig.getParameter< edm::InputTag > ("inputDeDxTag"));
 
   thisModuleTag_ = edm::InputTag(iConfig.getParameter<std::string>("@module_label"));
 
   //register your products
   produces<reco::RecoChargedCandidateCollection>();
 }

HLTDeDxFilter::~HLTDeDxFilter ( )

Definition at line 64 of file HLTDeDxFilter.cc.

64 {}

Member Function Documentation

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

static

Definition at line 66 of file HLTDeDxFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_25ns14e33_v1_cff::InputTag.

                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<bool>("saveTags",false);
   desc.add<double>("minDEDx",0.0);
   desc.add<double>("minPT",0.0);
   desc.add<double>("minNOM",0.0);
   desc.add<double>("maxETA",5.5);
   desc.add<double>("minNumValidHits",0);
   desc.add<double>("maxNHitMissIn",99);
   desc.add<double>("maxNHitMissMid",99);
   desc.add<double>("maxRelTrkIsoDeltaRp3", -1);
   desc.add<double>("relTrkIsoDeltaRSize", 0.3);
   desc.add<double>("maxAssocCaloE", -99);
   desc.add<double>("maxAssocCaloEDeltaRSize", 0.5);
   desc.add<edm::InputTag>("caloTowersTag",edm::InputTag("hltTowerMakerForAll"));
   desc.add<edm::InputTag>("inputTracksTag",edm::InputTag("hltL3Mouns"));
   desc.add<edm::InputTag>("inputDeDxTag",edm::InputTag("HLTdedxHarm2"));
   descriptions.add("hltDeDxFilter",desc);
 }

bool HLTDeDxFilter::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

Implements HLTFilter.

Definition at line 90 of file HLTDeDxFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), edm::SortedCollection< T, SORT >::begin(), caloTowersToken_, reco::deltaR2(), edm::SortedCollection< T, SORT >::end(), eta, edm::Event::getByToken(), i, inputdedxTag_, inputdedxToken_, inputTracksTag_, inputTracksToken_, j, maxAssocCaloE_, maxAssocCaloEDeltaRSize_, maxETA_, maxNHitMissIn_, maxNHitMissMid_, maxRelTrkIsoDeltaRp3_, minDEDx_, minNOM_, minNumValidHits_, minPT_, reco::HitPattern::MISSING_INNER_HITS, listHistos::NTracks, p4, phi, funct::pow(), EnergyCorrector::pt, edm::Event::put(), lumiQueryAPI::q, dt_dqm_sourceclient_common_cff::reco, relTrkIsoDeltaRSize_, HLTFilter::saveTags(), reco::RecoChargedCandidate::setTrack(), mathSSE::sqrt(), thisModuleTag_, reco::HitPattern::TRACK_HITS, HLT_25ns14e33_v3_cff::trackCollection, and trigger::TriggerMuon.

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   auto_ptr<RecoChargedCandidateCollection> chargedCandidates( new std::vector<RecoChargedCandidate> );
 
   ModuleDescription moduleDesc_;
 
   if (saveTags()){
      filterproduct.addCollectionTag(thisModuleTag_);
      filterproduct.addCollectionTag(inputTracksTag_);
      filterproduct.addCollectionTag(inputdedxTag_);
   }
 
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByToken(inputTracksToken_,trackCollectionHandle);
   const reco::TrackCollection &trackCollection = *trackCollectionHandle.product();
 
   edm::Handle<edm::ValueMap<reco::DeDxData> > dEdxTrackHandle;
   iEvent.getByToken(inputdedxToken_, dEdxTrackHandle);
   const edm::ValueMap<reco::DeDxData> &dEdxTrack = *dEdxTrackHandle.product();
 
   edm::Handle<CaloTowerCollection> caloTowersHandle;
   if(maxAssocCaloE_ >= 0) iEvent.getByToken(caloTowersToken_, caloTowersHandle);
 
   bool accept=false;
   int  NTracks = 0;
   //fill local arrays for eta, phi, and pt
   float eta[trackCollection.size()], phi[trackCollection.size()], pt[trackCollection.size()]; 
   for(unsigned int i=0; i<trackCollection.size(); i++){
     eta[i] = trackCollection[i].eta();
     phi[i] = trackCollection[i].phi();
     pt[i] = trackCollection[i].pt(); 
   }   
   for(unsigned int i=0; i<trackCollection.size(); i++){
     reco::TrackRef track  = reco::TrackRef( trackCollectionHandle, i );
      if(pt[i]>minPT_ && fabs(eta[i])<maxETA_ && dEdxTrack[track].numberOfMeasurements()>minNOM_ && dEdxTrack[track].dEdx()>minDEDx_){
       NTracks++;
        if(track->numberOfValidHits() < minNumValidHits_) continue;
        if(track->hitPattern().trackerLayersWithoutMeasurement( reco::HitPattern::MISSING_INNER_HITS) > maxNHitMissIn_) continue;
        if(track->hitPattern().trackerLayersWithoutMeasurement( reco::HitPattern::TRACK_HITS) > maxNHitMissMid_) continue;
        if (saveTags()){
      Particle::Charge q = track->charge();
           //SAVE DEDX INFORMATION AS IF IT WAS THE MASS OF THE PARTICLE
           Particle::LorentzVector p4(track->px(), track->py(), track->pz(), sqrt(pow(track->p(),2) + pow(dEdxTrack[track].dEdx(),2)));
           Particle::Point vtx(track->vx(),track->vy(), track->vz());
           //SAVE NOH, NOM, NOS INFORMATION AS IF IT WAS THE PDGID OF THE PARTICLE
           int Hits  = ((dEdxTrack[track].numberOfSaturatedMeasurements()&0xFF)<<16) | ((dEdxTrack[track].numberOfMeasurements()&0xFF)<<8) | (track->found()&0xFF);
           RecoChargedCandidate cand(q, p4, vtx, Hits, 0);
           cand.setTrack(track);
           chargedCandidates->push_back(cand);
        }
        
        //calculate relative trk isolation only if parameter maxRelTrkIsoDeltaRp3 is greater than 0
        if(maxRelTrkIsoDeltaRp3_ >= 0){
      auto ptCone = trackCollection[i].pt();
      for(unsigned int j=0; j<trackCollection.size(); j++){
      if (i==j) continue; // do not compare track to itself
          auto trkDeltaR2 = deltaR2(eta[i], phi[i], eta[j], phi[j]);
          if (trkDeltaR2 < relTrkIsoDeltaRSize_ * relTrkIsoDeltaRSize_){
            ptCone+=pt[j];
      }
        }
        double relTrkIso = (ptCone - pt[i])/(pt[i]);
        if (relTrkIso > maxRelTrkIsoDeltaRp3_) continue;
        }
        
        //calculate the calorimeter energy associated with the track if maxAssocCaloE_ >= 0 
        if(maxAssocCaloE_ >= 0){
        //Access info about Calo Towers                                                                                       
        double caloEMDeltaRp5  = 0;
        double caloHadDeltaRp5 = 0;
        const CaloTowerCollection &caloTower = *caloTowersHandle.product();
        for (CaloTowerCollection::const_iterator j=caloTower.begin(); j!=caloTower.end(); j++) {
      auto caloDeltaR2 = deltaR2(eta[i], phi[i], j->eta(), j->phi());
      double Eem  = j->emEnergy();
      double Ehad = j->hadEnergy();
 
          if (caloDeltaR2 < (maxAssocCaloEDeltaRSize_ * maxAssocCaloEDeltaRSize_) ) {
            caloEMDeltaRp5  += Eem;
            caloHadDeltaRp5 += Ehad;
          }
        }
        if (caloEMDeltaRp5 + caloHadDeltaRp5 > maxAssocCaloE_) continue;
        }
 
        accept=true;
     }
   }
 
   // put filter object into the Event
    if(saveTags()){
      edm::OrphanHandle<RecoChargedCandidateCollection> chargedCandidatesHandle = iEvent.put(chargedCandidates);
      for(int i=0; i<NTracks; i++){
           filterproduct.addObject(TriggerMuon,RecoChargedCandidateRef(chargedCandidatesHandle,i));
      }
    }
 
   return accept;
 }