#include <HLTPixlMBFilt.h>

Inheritance diagram for HLTPixlMBFilt:

Public Member Functions
virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct)

	HLTPixlMBFilt (const edm::ParameterSet &)

	~HLTPixlMBFilt ()

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

int	module () const

const std::string *	moduleLabel () const

int	path () const

const std::string *	pathName () const

std::pair< int, int >	pmid () const

bool	saveTags () const

virtual	~HLTFilter ()

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

virtual	~EDFilter ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
double	min_Pt_

float	min_sep_

unsigned int	min_trks_

edm::InputTag	pixlTag_

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

typedef WorkerT< EDFilter >	WorkerType

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

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

See header file for documentation

Date:: 2012/01/21 15:00:22

Revision:: 1.4

Author: Mika Huhtinen

Definition at line 22 of file HLTPixlMBFilt.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file HLTPixlMBFilt.cc.

References edm::InputTag::encode(), LogDebug, min_Pt_, min_sep_, min_trks_, and pixlTag_.

                                                            : HLTFilter(iConfig),
     pixlTag_ (iConfig.getParameter<edm::InputTag>("pixlTag")),
     min_Pt_  (iConfig.getParameter<double>("MinPt")),
     min_trks_  (iConfig.getParameter<unsigned int>("MinTrks")),
     min_sep_  (iConfig.getParameter<double>("MinSep"))
 
 {
   LogDebug("") << "MinPt cut " << min_Pt_   << "pixl: " << pixlTag_.encode();
   LogDebug("") << "Requesting : " << min_trks_ << " tracks from same vertex ";
   LogDebug("") << "Requesting tracks from same vertex eta-phi separation by " << min_sep_;
 }

HLTPixlMBFilt::~HLTPixlMBFilt ( )

Definition at line 46 of file HLTPixlMBFilt.cc.

47 {

48 }

Member Function Documentation

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

virtual

Implements HLTFilter.

Definition at line 55 of file HLTPixlMBFilt.cc.

References accept(), trigger::TriggerRefsCollections::addObject(), edm::Event::getByLabel(), gen::k, LogDebug, min_sep_, min_trks_, pixlTag_, L1Trigger_dataformats::reco, mathSSE::sqrt(), testEve_cfg::tracks, and trigger::TriggerTrack.

 {
    using namespace std;
    using namespace edm;
    using namespace reco;
    using namespace trigger;
 
    // All HLT filters must create and fill an HLT filter object,
    // recording any reconstructed physics objects satisfying (or not)
    // this HLT filter, and place it in the Event.
 
 
 
    // Specific filter code
 
    // get hold of products from Event
 
    Handle<RecoChargedCandidateCollection> tracks;
 
    iEvent.getByLabel(pixlTag_,tracks);
 
    // pixel tracks
    int npixl_tot = 0;
    vector<double> etastore;
    vector<double> phistore;
    vector<int> itstore;
    bool accept;
    RecoChargedCandidateCollection::const_iterator apixl(tracks->begin());
    RecoChargedCandidateCollection::const_iterator epixl(tracks->end());
    RecoChargedCandidateCollection::const_iterator ipixl, jpixl;
    unsigned int nsame_vtx=0;
    int itrk = -1;
    if (tracks->size() >= min_trks_) {
      for (ipixl=apixl; ipixl!=epixl; ipixl++){ 
        itrk++;
        const double& ztrk1 = ipixl->vz();           
        const double& etatrk1 = ipixl->momentum().eta();
        const double& phitrk1 = ipixl->momentum().phi();
        nsame_vtx=1;
        etastore.clear();
        phistore.clear();
        itstore.clear();
        etastore.push_back(etatrk1);
        phistore.push_back(phitrk1);
        itstore.push_back(itrk);
        if (fabs(ztrk1) < 15.0) {
          //  check this track against all others to see if others start from same point
      int jtrk=-1;
          for (jpixl=apixl; jpixl!=epixl; jpixl++) {
        jtrk++;
        if (jpixl==ipixl) continue;
            const double& ztrk2 = jpixl->vz();           
            const double& etatrk2 = jpixl->momentum().eta();
            const double& phitrk2 = jpixl->momentum().phi();
            double eta_dist=etatrk2-etatrk1;
            double phi_dist=phitrk2-phitrk1;
            double etaphi_dist=sqrt(eta_dist*eta_dist + phi_dist*phi_dist);
            if (fabs(ztrk2-ztrk1) < 1.0 && etaphi_dist > min_sep_) {
           if (min_trks_ <= 2 || itstore.size() <= 1) {
             etastore.push_back(etatrk2);
             phistore.push_back(phitrk2);
         itstore.push_back(jtrk);
                 nsame_vtx++;
               } else {
                 // check also separation to already found 'second' tracks
         LogDebug("") << "HLTPixlMBFilt: with mintrks=2 we should not be here...";
         bool isok = true;
         for (unsigned int k=1; k < itstore.size(); k++) {
                   eta_dist=etatrk2-etastore.at(k);
                   phi_dist=phitrk2-phistore.at(k);
                   etaphi_dist=sqrt(eta_dist*eta_dist + phi_dist*phi_dist);
           if (etaphi_dist < min_sep_) {
             isok=false;
             break;
                   }
         }
         if (isok) {
               etastore.push_back(etatrk2);
               phistore.push_back(phitrk2);
                   itstore.push_back(jtrk);
                   nsame_vtx++;
                 }
           }
        }
            if (nsame_vtx >= min_trks_) break;
          }
        }
        npixl_tot++;
 
        if (nsame_vtx >= min_trks_) break;
      }
 
      //   final filter decision:
      //   request at least min_trks_ tracks compatible with vertex-region
      accept = (nsame_vtx >= min_trks_ ) ;
 
    } else {
      accept = false;
    }
 
    // At this point we have the indices of the accepted tracks stored in itstore
    // we now move them to the filterproduct
 
    if (accept) {
      for (unsigned int ipos=0; ipos < itstore.size(); ipos++) {
        int iaddr=itstore.at(ipos);
        filterproduct.addObject(TriggerTrack,RecoChargedCandidateRef(tracks,iaddr));
      }
    }
 
   LogDebug("") << "Number of pixel-track objects accepted:"
                << " " << npixl_tot;
 
   // return with final filter decision
   return accept;
 
 }