#include <HLTPixlMBFilt.h>

Inheritance diagram for HLTPixlMBFilt:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)

	HLTPixlMBFilt (const edm::ParameterSet &)

	~HLTPixlMBFilt ()

Public Member Functions inherited from HLTFilter
	HLTFilter ()

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

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

virtual	~ProducerBase ()

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 edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

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

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

See header file for documentation

Date:: 2007/04/05 16:43:02

Revision:: 1.2

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_.

                                                            :
     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_;
 
    //register your products
   produces<trigger::TriggerFilterObjectWithRefs>();
 }

HLTPixlMBFilt::~HLTPixlMBFilt ( )

Definition at line 49 of file HLTPixlMBFilt.cc.

50 {

51 }

Member Function Documentation

bool HLTPixlMBFilt::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements HLTFilter.

Definition at line 58 of file HLTPixlMBFilt.cc.

References accept(), edm::Event::getByLabel(), gen::k, LogDebug, min_sep_, min_trks_, module(), path(), pixlTag_, edm::Event::put(), dt_dqm_sourceclient_common_cff::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 filterobject
 
    // The filter object
    auto_ptr<TriggerFilterObjectWithRefs> filterobject (new TriggerFilterObjectWithRefs(path(),module()));
 
    if (accept) {
      for (unsigned int ipos=0; ipos < itstore.size(); ipos++) {
        int iaddr=itstore.at(ipos);
        filterobject->addObject(TriggerTrack,RecoChargedCandidateRef(tracks,iaddr));
      }
    }
    // put filter object into the Event
    iEvent.put(filterobject);
 
 
   LogDebug("") << "Number of pixel-track objects accepted:"
                << " " << npixl_tot;
 
   // return with final filter decision
   return accept;
 
 }