CaloParticleSelector Class Reference

#include <CaloParticleSelector.h>

Public Member Functions
	CaloParticleSelector ()
	CaloParticleSelector (double ptMin, double ptMax, double minRapidity, double maxRapidity, double lip, double tip, int minHit, unsigned int maxSimClusters, bool signalOnly, bool intimeOnly, bool chargedOnly, bool stableOnly, const std::vector< int > &pdgId=std::vector< int >(), double minPhi=-3.2, double maxPhi=3.2)
bool	operator() (const CaloParticle &cp, std::vector< SimVertex > const &simVertices) const

Private Attributes
bool	chargedOnly_
bool	intimeOnly_
double	lip_
float	maxRapidity_
unsigned int	maxSimClusters_
float	meanPhi_
int	minHit_
float	minRapidity_
std::vector< int >	pdgId_
double	ptMax2_
double	ptMin2_
float	rangePhi_
bool	signalOnly_
bool	stableOnly_
double	tip2_

Detailed Description

Definition at line 10 of file CaloParticleSelector.h.

Constructor & Destructor Documentation

CaloParticleSelector() [1/2]

CaloParticleSelector::CaloParticleSelector ( )

inline

Definition at line 12 of file CaloParticleSelector.h.

{}

CaloParticleSelector() [2/2]

CaloParticleSelector::CaloParticleSelector ( double  ptMin, double  ptMax, double  minRapidity, double  maxRapidity, double  lip, double  tip, int  minHit, unsigned int  maxSimClusters, bool  signalOnly, bool  intimeOnly, bool  chargedOnly, bool  stableOnly, const std::vector< int > &  pdgId = std::vector<int>(), double  minPhi = -3.2, double  maxPhi = 3.2  )

inline

Definition at line 13 of file CaloParticleSelector.h.

      : ptMin2_(ptMin * ptMin),
        ptMax2_(ptMax * ptMax),
        minRapidity_(minRapidity),
        maxRapidity_(maxRapidity),
        lip_(lip),
        tip2_(tip * tip),
        meanPhi_((minPhi + maxPhi) / 2.),
        rangePhi_((maxPhi - minPhi) / 2.),
        minHit_(minHit),
        maxSimClusters_(maxSimClusters),
        signalOnly_(signalOnly),
        intimeOnly_(intimeOnly),
        chargedOnly_(chargedOnly),
        stableOnly_(stableOnly),
        pdgId_(pdgId) {
    if (minPhi >= maxPhi) {
      throw cms::Exception("Configuration")
          << "CaloParticleSelector: 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")
          << "CaloParticleSelector: 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")
          << "CaloParticleSelector: 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 CaloParticleSelector::operator() ( const CaloParticle &  cp, std::vector< SimVertex > const &  simVertices  ) const

inline

Definition at line 62 of file CaloParticleSelector.h.

                                                                                         {
    // signal only means no PU particles
    if (signalOnly_ && !(cp.eventId().bunchCrossing() == 0 && cp.eventId().event() == 0))
      return false;
    // intime only means no OOT PU particles
    if (intimeOnly_ && !(cp.eventId().bunchCrossing() == 0))
      return false;
 
    if (cp.simClusters().size() > maxSimClusters_)
      return false;
 
    auto pdgid = cp.pdgId();
    if (!pdgId_.empty()) {
      bool testId = false;
      for (auto id : pdgId_) {
        if (id == pdgid) {
          testId = true;
          break;
        }
      }
      if (!testId)
        return false;
    }
 
    if (chargedOnly_ && cp.charge() == 0)
      return false;  //select only if charge!=0
 
    // select only stable particles
    if (stableOnly_) {
      for (CaloParticle::genp_iterator j = cp.genParticle_begin(); j != cp.genParticle_end(); ++j) {
        if (j->get() == nullptr || j->get()->status() != 1) {
          return false;
        }
      }
 
      // test for remaining unstabled due to lack of genparticle pointer
      std::vector<int> pdgids{11, 13, 211, 321, 2212, 3112, 3222, 3312, 3334};
      if (cp.status() == -99 && (!std::binary_search(pdgids.begin(), pdgids.end(), std::abs(pdgid)))) {
        return false;
      }
    }
 
    auto etaOk = [&](const CaloParticle& p) -> bool {
      float eta = etaFromXYZ(p.px(), p.py(), p.pz());
      return (eta >= minRapidity_) & (eta <= maxRapidity_);
    };
    auto phiOk = [&](const CaloParticle& p) {
      float dphi = deltaPhi(atan2f(p.py(), p.px()), meanPhi_);
      return dphi >= -rangePhi_ && dphi <= rangePhi_;
    };
    auto ptOk = [&](const CaloParticle& p) {
      double pt2 = cp.p4().perp2();
      return pt2 >= ptMin2_ && pt2 <= ptMax2_;
    };
 
    return (ptOk(cp) && etaOk(cp) && phiOk(cp));
  }

References funct::abs(), chargedOnly_, CommonMethods::cp(), SiPixelRawToDigiRegional_cfi::deltaPhi, PVValHelper::eta, intimeOnly_, dqmiolumiharvest::j, maxRapidity_, maxSimClusters_, meanPhi_, minRapidity_, AlCaHLTBitMon_ParallelJobs::p, EgammaValidation_cff::pdgid, pdgId_, ntuple::pdgids, HLT_FULL_cff::pt2, ptMax2_, ptMin2_, rangePhi_, signalOnly_, and stableOnly_.

Member Data Documentation

chargedOnly_

bool CaloParticleSelector::chargedOnly_

private

Definition at line 133 of file CaloParticleSelector.h.

Referenced by operator()().

intimeOnly_

bool CaloParticleSelector::intimeOnly_

private

Definition at line 132 of file CaloParticleSelector.h.

Referenced by operator()().

lip_

double CaloParticleSelector::lip_

private

Definition at line 125 of file CaloParticleSelector.h.

maxRapidity_

float CaloParticleSelector::maxRapidity_

private

Definition at line 124 of file CaloParticleSelector.h.

Referenced by operator()().

maxSimClusters_

unsigned int CaloParticleSelector::maxSimClusters_

private

Definition at line 130 of file CaloParticleSelector.h.

Referenced by operator()().

meanPhi_

float CaloParticleSelector::meanPhi_

private

Definition at line 127 of file CaloParticleSelector.h.

Referenced by operator()().

minHit_

int CaloParticleSelector::minHit_

private

Definition at line 129 of file CaloParticleSelector.h.

minRapidity_

float CaloParticleSelector::minRapidity_

private

Definition at line 123 of file CaloParticleSelector.h.

Referenced by operator()().

pdgId_

std::vector<int> CaloParticleSelector::pdgId_

private

Definition at line 135 of file CaloParticleSelector.h.

Referenced by operator()().

ptMax2_

double CaloParticleSelector::ptMax2_

private

Definition at line 122 of file CaloParticleSelector.h.

Referenced by operator()().

ptMin2_

double CaloParticleSelector::ptMin2_

private

Definition at line 121 of file CaloParticleSelector.h.

Referenced by operator()().

rangePhi_

float CaloParticleSelector::rangePhi_

private

Definition at line 128 of file CaloParticleSelector.h.

Referenced by operator()().

signalOnly_

bool CaloParticleSelector::signalOnly_

private

Definition at line 131 of file CaloParticleSelector.h.

Referenced by operator()().

stableOnly_

bool CaloParticleSelector::stableOnly_

private

Definition at line 134 of file CaloParticleSelector.h.

Referenced by operator()().

tip2_

double CaloParticleSelector::tip2_

private

Definition at line 126 of file CaloParticleSelector.h.