#include <TrajSeedMatcher.h>

Inheritance diagram for TrajSeedMatcher::MatchingCutsV2:

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

bool	operator() (const SCHitMatch &scHitMatch) const override

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 106 of file TrajSeedMatcher.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 411 of file TrajSeedMatcher.cc.

References dPhiHighEt_, dPhiHighEtThres_, dPhiLowEtGrad_, dRZHighEt_, dRZHighEtThres_, dRZLowEtGrad_, etaBins_, Exception, mergeVDriftHistosByStation::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 447 of file TrajSeedMatcher.cc.

References funct::abs().

                                                               {
   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 113 of file TrajSeedMatcher.h.

References validate-o2o-wbm::f, and min().

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

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

overridevirtual

Implements TrajSeedMatcher::MatchingCuts.

Definition at line 434 of file TrajSeedMatcher.cc.

References funct::abs(), TrajSeedMatcher::SCHitMatch::dPhi, TrajSeedMatcher::SCHitMatch::dRZ, TrajSeedMatcher::SCHitMatch::et, and TrajSeedMatcher::SCHitMatch::eta.

                                                                                                 {
   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;
 }