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, 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 &tp, std::vector< SimVertex > const &simVertices) const

Private Attributes
bool	chargedOnly_
bool	intimeOnly_
double	lip_
float	maxRapidity_
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 17 of file CaloParticleSelector.h.

Constructor & Destructor Documentation

CaloParticleSelector::CaloParticleSelector ( )

inline

Definition at line 19 of file CaloParticleSelector.h.

Referenced by reco::modules::ParameterAdapter< CaloParticleSelector >::make().

{}

CaloParticleSelector::CaloParticleSelector ( double  ptMin, double  ptMax, double  minRapidity, double  maxRapidity, double  lip, double  tip, int  minHit, 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 20 of file CaloParticleSelector.h.

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

      : 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),
        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.";
    }
  }

Member Function Documentation

bool CaloParticleSelector::operator() ( const CaloParticle &  tp, std::vector< SimVertex > const &  simVertices  ) const

inline

Definition at line 67 of file CaloParticleSelector.h.

References funct::abs(), EncodedEventId::bunchCrossing(), CaloParticle::charge(), chargedOnly_, SiPixelRawToDigiRegional_cfi::deltaPhi, PVValHelper::eta, EncodedEventId::event(), CaloParticle::eventId(), CaloParticle::genParticle_begin(), CaloParticle::genParticle_end(), intimeOnly_, dqmiolumiharvest::j, maxRapidity_, meanPhi_, minRapidity_, AlCaHLTBitMon_ParallelJobs::p, CaloParticle::p4(), EgammaValidation_cff::pdgid, CaloParticle::pdgId(), pdgId_, ntuple::pdgids, HLT_2018_cff::pt2, ptMax2_, ptMin2_, rangePhi_, signalOnly_, stableOnly_, and CaloParticle::status().

                                                                                         {
    // signal only means no PU particles
    if (signalOnly_ && !(tp.eventId().bunchCrossing() == 0 && tp.eventId().event() == 0))
      return false;
    // intime only means no OOT PU particles
    if (intimeOnly_ && !(tp.eventId().bunchCrossing() == 0))
      return false;

    auto pdgid = tp.pdgId();
    if (!pdgId_.empty()) {
      bool testId = false;
      for (auto id : pdgId_) {
        if (id == pdgid) {
          testId = true;
          break;
        }
      }
      if (!testId)
        return false;
    }

    if (chargedOnly_ && tp.charge() == 0)
      return false;  //select only if charge!=0

    // select only stable particles
    if (stableOnly_) {
      for (CaloParticle::genp_iterator j = tp.genParticle_begin(); j != tp.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 (tp.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 = tp.p4().perp2();
      return pt2 >= ptMin2_ && pt2 <= ptMax2_;
    };

    return (ptOk(tp) && etaOk(tp) && phiOk(tp));
  }

Member Data Documentation

bool CaloParticleSelector::chargedOnly_

private

Definition at line 134 of file CaloParticleSelector.h.

Referenced by operator()().

bool CaloParticleSelector::intimeOnly_

private

Definition at line 133 of file CaloParticleSelector.h.

Referenced by operator()().

double CaloParticleSelector::lip_

private

Definition at line 127 of file CaloParticleSelector.h.

float CaloParticleSelector::maxRapidity_

private

Definition at line 126 of file CaloParticleSelector.h.

Referenced by operator()().

float CaloParticleSelector::meanPhi_

private

Definition at line 129 of file CaloParticleSelector.h.

Referenced by operator()().

int CaloParticleSelector::minHit_

private

Definition at line 131 of file CaloParticleSelector.h.

float CaloParticleSelector::minRapidity_

private

Definition at line 125 of file CaloParticleSelector.h.

Referenced by operator()().

std::vector<int> CaloParticleSelector::pdgId_

private

Definition at line 136 of file CaloParticleSelector.h.

Referenced by operator()().

double CaloParticleSelector::ptMax2_

private

Definition at line 124 of file CaloParticleSelector.h.

Referenced by operator()().

double CaloParticleSelector::ptMin2_

private

Definition at line 123 of file CaloParticleSelector.h.

Referenced by operator()().

float CaloParticleSelector::rangePhi_

private

Definition at line 130 of file CaloParticleSelector.h.

Referenced by operator()().

bool CaloParticleSelector::signalOnly_

private

Definition at line 132 of file CaloParticleSelector.h.

Referenced by operator()().

bool CaloParticleSelector::stableOnly_

private

Definition at line 135 of file CaloParticleSelector.h.

Referenced by operator()().

double CaloParticleSelector::tip2_

private

Definition at line 128 of file CaloParticleSelector.h.