#include <HLTDisplacedEgammaFilter.h>

Inheritance diagram for HLTDisplacedEgammaFilter:

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

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

	~HLTDisplacedEgammaFilter ()

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

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

virtual	~ProducerBase ()

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Private Attributes
bool	EBOnly

edm::InputTag	inputTag_

edm::InputTag	inputTrk

edm::InputTag	L1IsoCollTag_

edm::InputTag	L1NonIsoCollTag_

int	maxTrkCut

int	ncandcut_

edm::InputTag	rechitsEB

edm::InputTag	rechitsEE

bool	relaxed_

double	seedTimeMax

double	seedTimeMin

double	sMaj_max

double	sMaj_min

double	sMin_max

double	sMin_min

double	trkdRCut

double	trkPtCut

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

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

Author: Monica Vazquez Acosta (CERN)

Id:: HLTDisplacedEgammaFilter.cc,v 1.10 2012/04/11 08:57:07 gruen Exp

Author: Monica Vazquez Acosta (CERN)

Definition at line 18 of file HLTDisplacedEgammaFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 29 of file HLTDisplacedEgammaFilter.cc.

References EBOnly, edm::ParameterSet::getParameter(), inputTag_, inputTrk, L1IsoCollTag_, L1NonIsoCollTag_, maxTrkCut, ncandcut_, rechitsEB, rechitsEE, relaxed_, seedTimeMax, seedTimeMin, sMaj_max, sMaj_min, sMin_max, sMin_min, trkdRCut, and trkPtCut.

                                                                                  : HLTFilter(iConfig) 
 {
   inputTag_    = iConfig.getParameter< edm::InputTag > ("inputTag");
   ncandcut_    = iConfig.getParameter<int> ("ncandcut");
   relaxed_     = iConfig.getParameter<bool> ("relaxed") ;
   L1IsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1IsoCand"); 
   L1NonIsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1NonIsoCand"); 
 
   inputTrk   = iConfig.getParameter< edm::InputTag > ("inputTrack");
   trkPtCut   = iConfig.getParameter<double> ("trackPtCut");
   trkdRCut   = iConfig.getParameter<double> ("trackdRCut");
   maxTrkCut  = iConfig.getParameter<int> ("maxTrackCut");
 
   rechitsEB  = iConfig.getParameter< edm::InputTag > ("RecHitsEB");
   rechitsEE  = iConfig.getParameter< edm::InputTag > ("RecHitsEE");
   
   EBOnly       = iConfig.getParameter<bool> ("EBOnly") ;
   sMin_min     = iConfig.getParameter<double> ("sMin_min");
   sMin_max     = iConfig.getParameter<double> ("sMin_max");
   sMaj_min     = iConfig.getParameter<double> ("sMaj_min");
   sMaj_max     = iConfig.getParameter<double> ("sMaj_max");
   seedTimeMin  = iConfig.getParameter<double> ("seedTimeMin");
   seedTimeMax  = iConfig.getParameter<double> ("seedTimeMax");
 
 }

HLTDisplacedEgammaFilter::~HLTDisplacedEgammaFilter ( )

Definition at line 55 of file HLTDisplacedEgammaFilter.cc.

55 {}

Member Function Documentation

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

static

Definition at line 59 of file HLTDisplacedEgammaFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLTFilter::makeHLTFilterDescription().

                                                                                           {
    edm::ParameterSetDescription desc;
    makeHLTFilterDescription(desc);
    desc.add<edm::InputTag>("inputTag",edm::InputTag("hltEGRegionalL1SingleEG22"));
    desc.add<edm::InputTag>("L1IsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
    desc.add<edm::InputTag>("L1NonIsoCand",edm::InputTag("hltL1NonIsoRecoEcalCandidate"));
    desc.add<edm::InputTag>("RecHitsEB",edm::InputTag("hltEcalRecHitAll", "EcalRecHitsEB"));
    desc.add<edm::InputTag>("RecHitsEE",edm::InputTag("hltEcalRecHitAll", "EcalRecHitsEE"));
    desc.add<edm::InputTag>("inputTrack",edm::InputTag("hltL1SeededEgammaRegionalCTFFinalFitWithMaterial"));
    desc.add<bool>("relaxed",false);
    desc.add<int>("ncandcut",1);
    desc.add<bool>("EBOnly",false);
    desc.add<double>("sMin_min",0.1);
    desc.add<double>("sMin_max",0.4);
    desc.add<double>("sMaj_min",0.0);
    desc.add<double>("sMaj_max",999.0);
    desc.add<double>("seedTimeMin", -25.0);
    desc.add<double>("seedTimeMax",  25.0);
    desc.add<int>("maxTrackCut", 0);
    desc.add<double>("trackPtCut", 3.0);
    desc.add<double>("trackdRCut", 0.5);
    descriptions.add("hltDisplacedEgammaFilter",desc);
 }

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

virtual

Implements HLTFilter.

Definition at line 84 of file HLTDisplacedEgammaFilter.cc.

References abs, accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), EcalClusterTools::cluster2ndMoments(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, EBOnly, edm::SortedCollection< T, SORT >::find(), edm::Event::getByLabel(), EcalClusterTools::getMaximum(), i, inputTag_, inputTrk, L1IsoCollTag_, L1NonIsoCollTag_, maxTrkCut, n, ncandcut_, HI_PhotonSkim_cff::rechits, rechitsEB, rechitsEE, relaxed_, HLTFilter::saveTags(), seedTimeMax, Cluster2ndMoments::sMaj, sMaj_max, Cluster2ndMoments::sMin, sMin_max, testEve_cfg::tracks, trigger::TriggerCluster, trkdRCut, and trkPtCut.

 {
   using namespace trigger;
 
   // The filter object
   if (saveTags()) {
     filterproduct.addCollectionTag(L1IsoCollTag_);
     if (relaxed_) filterproduct.addCollectionTag(L1NonIsoCollTag_);
   }
 
   // Ref to Candidate object to be recorded in filter object
    edm::Ref<reco::RecoEcalCandidateCollection> ref;
 
   // get hold of filtered candidates
   //edm::Handle<reco::HLTFilterObjectWithRefs> recoecalcands;
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByLabel (inputTag_,PrevFilterOutput);
 
   // get hold of collection of objects
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByLabel( inputTrk , tracks);
 
   // get the EcalRecHit
   edm::Handle<EcalRecHitCollection>      rechitsEB_ ;
   edm::Handle<EcalRecHitCollection>      rechitsEE_ ;
   iEvent.getByLabel( rechitsEB,     rechitsEB_ );
   iEvent.getByLabel( rechitsEE,     rechitsEE_ );
 
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > recoecalcands;   
   PrevFilterOutput->getObjects(TriggerCluster, recoecalcands);
  
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
   for (unsigned int i=0; i<recoecalcands.size(); i++) {
     
     ref = recoecalcands[i] ;
     if ( EBOnly &&  std::abs( ref->eta() ) >= 1.479  ) continue ;
 
 
     // S_Minor Cuts from the seed cluster
     reco::CaloClusterPtr SCseed = ref->superCluster()->seed() ;
     const EcalRecHitCollection* rechits = ( std::abs( ref->eta() ) < 1.479 ) ? rechitsEB_.product() : rechitsEE_.product() ;
 
     Cluster2ndMoments moments = EcalClusterTools::cluster2ndMoments(*SCseed, *rechits);
     float sMin =  moments.sMin  ;
     float sMaj =  moments.sMaj  ;
     if ( sMin < sMin_min || sMin > sMin_max ) continue ;
     if ( sMaj < sMaj_min || sMaj > sMaj_max ) continue ;
 
     // seed Time 
     std::pair<DetId, float> maxRH = EcalClusterTools::getMaximum( *SCseed, rechits );
     DetId seedCrystalId = maxRH.first;
     EcalRecHitCollection::const_iterator seedRH = rechits->find(seedCrystalId);
     float seedTime = (float)seedRH->time();
     if ( seedTime < seedTimeMin || seedTime > seedTimeMax ) continue ;
  
     //Track Veto
     
     int nTrk = 0 ;
     for (reco::TrackCollection::const_iterator it = tracks->begin(); it != tracks->end(); it++ )  {
         if ( it->pt() < trkPtCut ) continue ;
         LorentzVector trkP4( it->px(), it->py(), it->pz(), it->p() ) ;
         double dR =  ROOT::Math::VectorUtil::DeltaR( trkP4 , ref->p4()  ) ;
         if ( dR < trkdRCut )  nTrk++ ;
         if ( nTrk > maxTrkCut ) break ;
     }
     if ( nTrk > maxTrkCut ) continue ;     
     
 
     n++;
     // std::cout << "Passed eta: " << ref->eta() << std::endl;
     filterproduct.addObject(TriggerCluster, ref);
   }
   
   
   // filter decision
   bool accept(n>=ncandcut_);
   
   return accept;
 }