#include <TrajSeedMatcher.h>

Inheritance diagram for TrajSeedMatcher::MatchingCutsV2:

Public Member Functions
	MatchingCutsV2 (const edm::ParameterSet &pset)

bool	operator() (const SCHitMatch &scHitMatch) const

Public Member Functions inherited from TrajSeedMatcher::MatchingCuts
	MatchingCuts ()

virtual	~MatchingCuts ()

Private Member Functions
size_t	getBinNr (float eta) const

float	getCutValue (float et, float highEt, float highEtThres, float lowEtGrad) const

Private Attributes
std::vector< double >	dPhiHighEt_

std::vector< double >	dPhiHighEtThres_

std::vector< double >	dPhiLowEtGrad_

std::vector< double >	dRZHighEt_

std::vector< double >	dRZHighEtThres_

std::vector< double >	dRZLowEtGrad_

std::vector< double >	etaBins_

Detailed Description

Definition at line 191 of file TrajSeedMatcher.h.

Constructor & Destructor Documentation

TrajSeedMatcher::MatchingCutsV2::MatchingCutsV2 ( const edm::ParameterSet & pset )

explicit

Definition at line 444 of file TrajSeedMatcher.cc.

References dPhiHighEt_, dPhiHighEtThres_, dPhiLowEtGrad_, dRZHighEt_, dRZHighEtThres_, dRZLowEtGrad_, etaBins_, Exception, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                         :
   dPhiHighEt_(pset.getParameter<std::vector<double> >("dPhiMaxHighEt")),
   dPhiHighEtThres_(pset.getParameter<std::vector<double> >("dPhiMaxHighEtThres")),
   dPhiLowEtGrad_(pset.getParameter<std::vector<double> >("dPhiMaxLowEtGrad")),
   dRZHighEt_(pset.getParameter<std::vector<double> >("dRZMaxHighEt")),
   dRZHighEtThres_(pset.getParameter<std::vector<double> >("dRZMaxHighEtThres")),
   dRZLowEtGrad_(pset.getParameter<std::vector<double> >("dRZMaxLowEtGrad")),
   etaBins_(pset.getParameter<std::vector<double> >("etaBins"))
 {
   auto binSizeCheck = [](size_t sizeEtaBins,const std::vector<double>& vec,const std::string& name){
     if(vec.size()!=sizeEtaBins+1){ 
       throw cms::Exception("InvalidConfig")<<" when constructing TrajSeedMatcher::MatchingCutsV2 "<< name<<" has "<<vec.size()<<" bins, it should be equal to #bins of etaBins+1"<<sizeEtaBins+1;
     }
   };
   binSizeCheck(etaBins_.size(),dPhiHighEt_,"dPhiMaxHighEt");
   binSizeCheck(etaBins_.size(),dPhiHighEtThres_,"dPhiMaxHighEtThres");
   binSizeCheck(etaBins_.size(),dPhiLowEtGrad_,"dPhiMaxLowEtGrad");
   binSizeCheck(etaBins_.size(),dRZHighEt_,"dRZMaxHighEt");
   binSizeCheck(etaBins_.size(),dRZHighEtThres_,"dRZMaxHighEtThres");
   binSizeCheck(etaBins_.size(),dRZLowEtGrad_,"dRZMaxLowEtGrad");
 }

Member Function Documentation

size_t TrajSeedMatcher::MatchingCutsV2::getBinNr ( float eta ) const

private

Definition at line 478 of file TrajSeedMatcher.cc.

References funct::abs(), and etaBins_.

Referenced by operator()().

 {
   const float absEta = std::abs(eta);
   for(size_t etaNr=0;etaNr<etaBins_.size();etaNr++){
     if(absEta<etaBins_[etaNr]) return etaNr;
   }
   return etaBins_.size();
 }

float TrajSeedMatcher::MatchingCutsV2::getCutValue	(	float	et,
		float	highEt,
		float	highEtThres,
		float	lowEtGrad
	)		const

inlineprivate

Definition at line 197 of file TrajSeedMatcher.h.

References f, and min().

Referenced by operator()().

                                                                                       {
       return  highEt + std::min(0.f,et-highEtThres)*lowEtGrad;
     }

bool TrajSeedMatcher::MatchingCutsV2::operator() ( const SCHitMatch & scHitMatch ) const

virtual

Implements TrajSeedMatcher::MatchingCuts.

Definition at line 466 of file TrajSeedMatcher.cc.

References funct::abs(), TrajSeedMatcher::SCHitMatch::dPhi(), dPhiHighEt_, dPhiHighEtThres_, dPhiLowEtGrad_, CSCSegmentAlgorithmRU_cfi::dPhiMax, TrajSeedMatcher::SCHitMatch::dRZ(), dRZHighEt_, dRZHighEtThres_, dRZLowEtGrad_, TrajSeedMatcher::SCHitMatch::et(), TrajSeedMatcher::SCHitMatch::eta(), getBinNr(), and getCutValue().

 {
   size_t binNr=getBinNr(scHitMatch.eta());
   float dPhiMax = getCutValue(scHitMatch.et(),dPhiHighEt_[binNr],dPhiHighEtThres_[binNr],dPhiLowEtGrad_[binNr]); 
   if(dPhiMax>=0 && std::abs(scHitMatch.dPhi()) > dPhiMax) return false;  
   float dRZMax = getCutValue(scHitMatch.et(),dRZHighEt_[binNr],dRZHighEtThres_[binNr],dRZLowEtGrad_[binNr]);
   if(dRZMax>=0 && std::abs(scHitMatch.dRZ()) > dRZMax) return false;
   
   return true;
 }