#include <HLTElectronPixelMatchFilter.h>

Inheritance diagram for HLTElectronPixelMatchFilter:

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

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

	~HLTElectronPixelMatchFilter ()

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 Member Functions
float	calDPhi1Sq (reco::ElectronSeedCollection::const_iterator seed, int charge) const

float	calDPhi2Sq (reco::ElectronSeedCollection::const_iterator seed, int charge) const

float	calDZ2Sq (reco::ElectronSeedCollection::const_iterator seed, int charge) const

int	getNrOfMatches (edm::Handle< reco::ElectronSeedCollection > &eleSeeds, reco::SuperClusterRef &candSCRef) const

Private Attributes
edm::InputTag	candTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	candToken_

bool	doIsolated_

bool	isPixelVeto_

edm::InputTag	L1IsoCollTag_

edm::InputTag	L1IsoPixelSeedsTag_

edm::EDGetTokenT < reco::ElectronSeedCollection >	L1IsoPixelSeedsToken_

edm::InputTag	L1NonIsoCollTag_

edm::InputTag	L1NonIsoPixelSeedsTag_

edm::EDGetTokenT < reco::ElectronSeedCollection >	L1NonIsoPixelSeedsToken_

int	ncandcut_

double	npixelmatchcut_

float	s2BarrelThres_

float	s2ForwardThres_

float	s2InterThres_

float	sPhi1B_

float	sPhi1F_

float	sPhi1I_

float	sPhi2B_

float	sPhi2F_

float	sPhi2I_

float	sR2F_

float	sR2I_

float	sZ2B_

bool	useS_

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: Aidan Randle-Condle (ULB)

Id:: HLTElectronPixelMatchFilter.cc,v 1.15 2012/03/06 10:13:59 sharper Exp

Author: Aidan Randle-Conde (ULB)

Definition at line 18 of file HLTElectronPixelMatchFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 35 of file HLTElectronPixelMatchFilter.cc.

References candTag_, candToken_, doIsolated_, edm::ParameterSet::getParameter(), isPixelVeto_, L1IsoCollTag_, L1IsoPixelSeedsTag_, L1IsoPixelSeedsToken_, L1NonIsoCollTag_, L1NonIsoPixelSeedsTag_, L1NonIsoPixelSeedsToken_, ncandcut_, npixelmatchcut_, funct::pow(), s2BarrelThres_, s2ForwardThres_, s2InterThres_, sPhi1B_, sPhi1F_, sPhi1I_, sPhi2B_, sPhi2F_, sPhi2I_, sR2F_, sR2I_, sZ2B_, and useS_.

                                                                                        : HLTFilter(iConfig) {
   candTag_                = iConfig.getParameter< edm::InputTag > ("candTag");
   L1IsoPixelSeedsTag_     = iConfig.getParameter< edm::InputTag > ("L1IsoPixelSeedsTag");
   L1NonIsoPixelSeedsTag_  = iConfig.getParameter< edm::InputTag > ("L1NonIsoPixelSeedsTag");
   npixelmatchcut_         = iConfig.getParameter< double >        ("npixelmatchcut");
   ncandcut_               = iConfig.getParameter< int >           ("ncandcut");
   doIsolated_             = iConfig.getParameter< bool >          ("doIsolated");
   L1IsoCollTag_           = iConfig.getParameter< edm::InputTag > ("L1IsoCand");
   L1NonIsoCollTag_        = iConfig.getParameter< edm::InputTag > ("L1NonIsoCand");
   
   candToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(candTag_);
   L1IsoPixelSeedsToken_ = consumes<reco::ElectronSeedCollection>(L1IsoPixelSeedsTag_);
   L1NonIsoPixelSeedsToken_= consumes<reco::ElectronSeedCollection>(L1NonIsoPixelSeedsTag_);
   
   sPhi1B_ = iConfig.getParameter< double >("s_a_phi1B") ;
   sPhi1I_ = iConfig.getParameter< double >("s_a_phi1I") ;
   sPhi1F_ = iConfig.getParameter< double >("s_a_phi1F") ;
   sPhi2B_ = iConfig.getParameter< double >("s_a_phi2B") ;
   sPhi2I_ = iConfig.getParameter< double >("s_a_phi2I") ;
   sPhi2F_ = iConfig.getParameter< double >("s_a_phi2F") ;
   sZ2B_    = iConfig.getParameter< double >("s_a_zB"   ) ;
   sR2I_    = iConfig.getParameter< double >("s_a_rI"   ) ;
   sR2F_    = iConfig.getParameter< double >("s_a_rF"   ) ;
   s2BarrelThres_ = std::pow(std::atanh(iConfig.getParameter< double >("tanhSO10BarrelThres"))*10.,2);
   s2InterThres_ = std::pow(std::atanh(iConfig.getParameter< double >("tanhSO10InterThres"))*10.,2);
   s2ForwardThres_ = std::pow(std::atanh(iConfig.getParameter< double >("tanhSO10ForwardThres"))*10.,2);
   
   isPixelVeto_ = iConfig.getParameter< bool >("pixelVeto" );
   useS_        = iConfig.getParameter< bool >("useS" );
 
 }

HLTElectronPixelMatchFilter::~HLTElectronPixelMatchFilter ( )

Definition at line 67 of file HLTElectronPixelMatchFilter.cc.

68 {}

Member Function Documentation

float HLTElectronPixelMatchFilter::calDPhi1Sq	(	reco::ElectronSeedCollection::const_iterator	seed,
		int	charge
	)		const

private

Definition at line 70 of file HLTElectronPixelMatchFilter.cc.

References sPhi1B_, sPhi1F_, and sPhi1I_.

Referenced by getNrOfMatches().

 {
   const float dPhi1Const = seed->subDet1()==1 ? seed->subDet2()==1 ? sPhi1B_ : sPhi1I_ : sPhi1F_;
   float dPhi1 = charge<0 ? seed->dPhi1()/dPhi1Const :  seed->dPhi1Pos()/dPhi1Const; 
   return dPhi1*dPhi1;
 }

float HLTElectronPixelMatchFilter::calDPhi2Sq	(	reco::ElectronSeedCollection::const_iterator	seed,
		int	charge
	)		const

private

Definition at line 77 of file HLTElectronPixelMatchFilter.cc.

References sPhi2B_, sPhi2F_, and sPhi2I_.

Referenced by getNrOfMatches().

 {
   const float dPhi2Const = seed->subDet1()==1 ? seed->subDet2()==1 ? sPhi2B_ : sPhi2I_ : sPhi2F_;
   float dPhi2 = charge <0 ? seed->dPhi2()/dPhi2Const :  seed->dPhi2Pos()/dPhi2Const;
   return dPhi2*dPhi2;
 }

float HLTElectronPixelMatchFilter::calDZ2Sq	(	reco::ElectronSeedCollection::const_iterator	seed,
		int	charge
	)		const

private

Definition at line 85 of file HLTElectronPixelMatchFilter.cc.

References sR2F_, sR2I_, and sZ2B_.

Referenced by getNrOfMatches().

 {
   const float dRZ2Const = seed->subDet1()==1 ? seed->subDet2()==1 ? sZ2B_ : sR2I_ : sR2F_;
   float dRZ2 = charge<0 ? seed->dRz2()/dRZ2Const : seed->dRz2Pos()/dRZ2Const;
   return dRZ2*dRZ2;
 }

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

static

Definition at line 92 of file HLTElectronPixelMatchFilter.cc.

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

                                                                                              {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("candTag",edm::InputTag("hltEgammaHcalIsolFilter"));
   desc.add<edm::InputTag>("L1IsoPixelSeedsTag",edm::InputTag("electronPixelSeeds"));
   desc.add<edm::InputTag>("L1NonIsoPixelSeedsTag",edm::InputTag("electronPixelSeeds"));
   desc.add<double>("npixelmatchcut",1.0);
   desc.add<int>("ncandcut",1);
   desc.add<bool>("doIsolated",true);
   desc.add<edm::InputTag>("L1IsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
   desc.add<edm::InputTag>("L1NonIsoCand",edm::InputTag("hltL1NonIsoRecoEcalCandidate"));
   desc.add<double>("s_a_phi1B",    0.0069) ;
   desc.add<double>("s_a_phi1I",    0.0088) ;
   desc.add<double>("s_a_phi1F",    0.0076) ;
   desc.add<double>("s_a_phi2B",    0.00037) ;
   desc.add<double>("s_a_phi2I",    0.00070) ;
   desc.add<double>("s_a_phi2F",    0.00906) ;
   desc.add<double>("s_a_zB"   ,    0.012) ;
   desc.add<double>("s_a_rI"   ,    0.027) ;
   desc.add<double>("s_a_rF"   ,    0.040) ;
   desc.add<double>("s2_threshold", 0);
   desc.add<double>("tanhSO10BarrelThres",0.35);
   desc.add<double>("tanhSO10InterThres",1);
   desc.add<double>("tanhSO10ForwardThres",1);
   desc.add<bool>  ("useS"     , false);
   desc.add<bool>  ("pixelVeto", false);
 
   descriptions.add("hltElectronPixelMatchFilter",desc);
 }

int HLTElectronPixelMatchFilter::getNrOfMatches	(	edm::Handle< reco::ElectronSeedCollection > &	eleSeeds,
		reco::SuperClusterRef &	candSCRef
	)		const

private

Definition at line 178 of file HLTElectronPixelMatchFilter.cc.

References calDPhi1Sq(), calDPhi2Sq(), calDZ2Sq(), edm::RefToBase< T >::castTo(), s2BarrelThres_, s2ForwardThres_, s2InterThres_, and useS_.

Referenced by hltFilter().

 {
   int nrMatch=0;
   for(reco::ElectronSeedCollection::const_iterator seedIt = eleSeeds->begin(); seedIt != eleSeeds->end(); seedIt++){
     edm::RefToBase<reco::CaloCluster> caloCluster = seedIt->caloCluster() ;
     reco::SuperClusterRef scRef = caloCluster.castTo<reco::SuperClusterRef>() ;
     if(&(*candSCRef) ==  &(*scRef)){
       if(useS_){
     float s2Neg = calDPhi1Sq(seedIt,-1) + calDPhi2Sq(seedIt,-1) + calDZ2Sq(seedIt,-1);
     float s2Pos = calDPhi1Sq(seedIt,1) + calDPhi2Sq(seedIt,1) + calDZ2Sq(seedIt,1);
     
     const float s2Thres = seedIt->subDet1()==1 ? seedIt->subDet2()==1 ? s2BarrelThres_ : s2InterThres_ : s2ForwardThres_; 
     if(s2Neg<s2Thres || s2Pos<s2Thres) nrMatch++ ;
       }
       else nrMatch++;
     }//end sc ref match
   }//end loop over ele seeds
   return nrMatch;
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 122 of file HLTElectronPixelMatchFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), candToken_, doIsolated_, edm::Event::getByToken(), getNrOfMatches(), i, isPixelVeto_, L1IsoCollTag_, L1IsoPixelSeedsToken_, L1NonIsoCollTag_, L1NonIsoPixelSeedsToken_, gen::n, ncandcut_, npixelmatchcut_, HLTFilter::saveTags(), trigger::TriggerCluster, and trigger::TriggerPhoton.

                                                                                                                                                    {
   // The filter object
   using namespace trigger;
   if (saveTags()) {
     filterproduct.addCollectionTag(L1IsoCollTag_);
     if (not doIsolated_) filterproduct.addCollectionTag(L1NonIsoCollTag_);
   }
   
   // Ref to Candidate object to be recorded in filter object
   edm::Ref<reco::RecoEcalCandidateCollection> ref;
   
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken(candToken_,PrevFilterOutput);
   
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > recoecalcands;
   PrevFilterOutput->getObjects(TriggerCluster, recoecalcands);
   if(recoecalcands.empty()) PrevFilterOutput->getObjects(TriggerPhoton,recoecalcands);  //we dont know if its type trigger cluster or trigger photon
   
   //get hold of the pixel seed - supercluster association map
   edm::Handle<reco::ElectronSeedCollection> L1IsoSeeds;
   iEvent.getByToken(L1IsoPixelSeedsToken_,L1IsoSeeds);
 
   edm::Handle<reco::ElectronSeedCollection> L1NonIsoSeeds;
   if(!doIsolated_){
     iEvent.getByToken(L1NonIsoPixelSeedsToken_,L1NonIsoSeeds);
   }
   
   // look at all egammas,  check cuts and add to filter object
   int n = 0;
   for (unsigned int i=0; i<recoecalcands.size(); i++) {
 
     ref = recoecalcands[i];
     reco::SuperClusterRef recr2 = ref->superCluster();
     
     int nmatch = getNrOfMatches(L1IsoSeeds,recr2);
     if(!doIsolated_) nmatch+=getNrOfMatches(L1NonIsoSeeds,recr2);
 
     if (!isPixelVeto_) {
       if ( nmatch >= npixelmatchcut_) {
     n++;
     filterproduct.addObject(TriggerCluster, ref);
       }
     } else {
       if ( nmatch == 0) {
     n++;
     filterproduct.addObject(TriggerCluster, ref);
       }
     }
  
   }//end of loop over candidates
   
   // filter decision
   const bool accept(n>=ncandcut_);
   return accept;
 }