#include <GenParticleCustomSelector.h>

Public Member Functions
	GenParticleCustomSelector ()

	GenParticleCustomSelector (double ptMin, double minRapidity, double maxRapidity, double tip, double lip, bool chargedOnly, int status, const std::vector< int > &pdgId=std::vector< int >(), bool invertRapidityCut=false, double minPhi=-3.2, double maxPhi=3.2)

bool	operator() (const reco::GenParticle &tp) const
	Operator() performs the selection: e.g. if (tPSelector(tp)) {...}. More...

Private Attributes
bool	chargedOnly_

bool	invertRapidityCut_

double	lip_

double	maxRapidity_

float	meanPhi_

double	minRapidity_

std::vector< int >	pdgId_

double	ptMin_

float	rangePhi_

int	status_

double	tip_

Detailed Description

Definition at line 11 of file GenParticleCustomSelector.h.

Constructor & Destructor Documentation

◆ GenParticleCustomSelector() [1/2]

GenParticleCustomSelector::GenParticleCustomSelector ( )

inline

Definition at line 13 of file GenParticleCustomSelector.h.

13 {}

◆ GenParticleCustomSelector() [2/2]

GenParticleCustomSelector::GenParticleCustomSelector	(	double	ptMin,
		double	minRapidity,
		double	maxRapidity,
		double	tip,
		double	lip,
		bool	chargedOnly,
		int	status,
		const std::vector< int > &	pdgId = `std::vector<int>()`,
		bool	invertRapidityCut = `false`,
		double	minPhi = `-3.2`,
		double	maxPhi = `3.2`
	)

inline

Definition at line 14 of file GenParticleCustomSelector.h.

       : ptMin_(ptMin),
         minRapidity_(minRapidity),
         maxRapidity_(maxRapidity),
         meanPhi_((minPhi + maxPhi) / 2.),
         rangePhi_((maxPhi - minPhi) / 2.),
         tip_(tip),
         lip_(lip),
         chargedOnly_(chargedOnly),
         status_(status),
         pdgId_(pdgId),
         invertRapidityCut_(invertRapidityCut) {
     if (minPhi >= maxPhi) {
       throw cms::Exception("Configuration")
           << "GenParticleCustomSelector: minPhi (" << minPhi << ") must be smaller than maxPhi (" << maxPhi
           << "). The range is constructed from minPhi to maxPhi around their "
              "average.";
     }
     if (minPhi >= M_PI) {
       throw cms::Exception("Configuration") << "GenParticleCustomSelector: minPhi (" << minPhi
                                             << ") must be smaller than PI. The range is constructed from minPhi "
                                                "to maxPhi around their average.";
     }
     if (maxPhi <= -M_PI) {
       throw cms::Exception("Configuration") << "GenParticleCustomSelector: maxPhi (" << maxPhi
                                             << ") must be larger than -PI. The range is constructed from minPhi "
                                                "to maxPhi around their average.";
     }
   }

References Exception, M_PI, HLT_FULL_cff::maxPhi, and HLT_FULL_cff::minPhi.

Member Function Documentation

◆ operator()()

bool GenParticleCustomSelector::operator() ( const reco::GenParticle & tp ) const

inline

Operator() performs the selection: e.g. if (tPSelector(tp)) {...}.

Definition at line 55 of file GenParticleCustomSelector.h.

                                                    {
     if (chargedOnly_ && tp.charge() == 0)
       return false;  //select only if charge!=0
     bool testId = false;
     unsigned int idSize = pdgId_.size();
     if (idSize == 0)
       testId = true;
     else
       for (unsigned int it = 0; it != idSize; ++it) {
         if (tp.pdgId() == pdgId_[it])
           testId = true;
       }
  
     auto etaOk = [&](const reco::GenParticle& p) -> bool {
       float eta = p.eta();
       if (!invertRapidityCut_)
         return (eta >= minRapidity_) && (eta <= maxRapidity_);
       else
         return (eta < minRapidity_ || eta > maxRapidity_);
     };
     auto phiOk = [&](const reco::GenParticle& p) {
       float dphi = deltaPhi(atan2f(p.py(), p.px()), meanPhi_);
       return dphi >= -rangePhi_ && dphi <= rangePhi_;
     };
     auto ptOk = [&](const reco::GenParticle& p) {
       double pt = p.pt();
       return pt >= ptMin_;
     };
  
     return (ptOk(tp) && etaOk(tp) && phiOk(tp) && sqrt(tp.vertex().perp2()) <= tip_ && fabs(tp.vertex().z()) <= lip_ &&
             tp.status() == status_ && testId);
   }