MVAJetPuId Class Reference

#include <MVAJetPuId.h>

Public Types
typedef std::map< std::string, std::pair< float *, float > >	variables_list_t
enum	version_t { USER = -1, CATEv0 = 0 }

Public Member Functions
int	computeCutIDflag (float betaStarClassic, float dR2Mean, float nvtx, float jetPt, float jetEta)
int	computeIDflag (float mva, float jetPt, float jetEta)
int	computeIDflag (float mva, int ptId, int etaId)
PileupJetIdentifier	computeIdVariables (const reco::Jet *jet, float jec, const reco::Vertex *, const reco::VertexCollection &, double rho, bool calculateMva=false)
PileupJetIdentifier	computeMva ()
std::string	dumpVariables () const
std::pair< int, int >	getJetIdKey (float jetPt, float jetEta)
const variables_list_t &	getVariables () const
const std::string	method () const
	MVAJetPuId (int version=CATEv0, const std::string &tmvaWeight="", const std::string &tmvaMethod="", Float_t impactParTkThreshod_=1., const std::vector< std::string > &tmvaVariables=std::vector< std::string >())
	MVAJetPuId (const edm::ParameterSet &ps)
void	set (const PileupJetIdentifier &)
	~MVAJetPuId ()

Protected Member Functions
void	bookReader ()
void	initVariables ()
void	resetVariables ()
void	runMva ()
void	setup ()

Protected Attributes
Float_t	betaStarCut_ [NWPs][NEtas][NPts]
bool	cutBased_
Float_t	impactParTkThreshod_
PileupJetIdentifier	internalId_
Float_t	mvacut_ [NWPs][NEtas][NPts]
TMVA::Reader *	reader_
Float_t	rmsCut_ [NWPs][NEtas][NPts]
std::string	tmvaMethod_
std::map< std::string, std::string >	tmvaNames_
std::vector< std::string >	tmvaSpectators_
std::vector< std::string >	tmvaVariables_
std::string	tmvaWeights_
variables_list_t	variables_
Int_t	version_

Static Private Attributes
static constexpr int	NEtas = 4
static constexpr int	NPts = 4
static constexpr int	NWPs = 3

Detailed Description

Definition at line 16 of file MVAJetPuId.h.

Member Typedef Documentation

variables_list_t

typedef std::map<std::string, std::pair<float *, float> > MVAJetPuId::variables_list_t

Definition at line 45 of file MVAJetPuId.h.

Member Enumeration Documentation

version_t

enum MVAJetPuId::version_t

Enumerator
USER
CATEv0

Definition at line 22 of file MVAJetPuId.h.

{ USER = -1, CATEv0 = 0 };

Constructor & Destructor Documentation

MVAJetPuId() [1/2]

MVAJetPuId::MVAJetPuId	(	int	version = CATEv0,
		const std::string &	tmvaWeight = "",
		const std::string &	tmvaMethod = "",
		Float_t	impactParTkThreshod_ = 1.,
		const std::vector< std::string > &	tmvaVariables = std::vector<std::string>()
	)

Definition at line 50 of file MVAJetPuId.cc.

References impactParTkThreshod_, reader_, setup(), hltPUIdAlgo_cff::tmvaMethod, tmvaMethod_, hltPUIdAlgo_cff::tmvaVariables, tmvaVariables_, hltPUIdAlgo_cff::tmvaWeights, tmvaWeights_, BeamSplash_cfg::version, and version_.

                                                                  {
  impactParTkThreshod_ = impactParTkThreshod;
  tmvaWeights_ = tmvaWeights;
  tmvaMethod_ = tmvaMethod;
  tmvaVariables_ = tmvaVariables;
  version_ = version;

  reader_ = nullptr;

  setup();
}

MVAJetPuId() [2/2]

MVAJetPuId::MVAJetPuId

(

const edm::ParameterSet &

ps

)

Definition at line 16 of file MVAJetPuId.cc.

References edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, impactParTkThreshod_, PileupJetIdentifier::kLoose, PileupJetIdentifier::kMedium, mvacut_, NPts, NWPs, reader_, setup(), AlCaHLTBitMon_QueryRunRegistry::string, tmvaMethod_, tmvaSpectators_, tmvaVariables_, tmvaWeights_, and version_.

                                                {
  impactParTkThreshod_ = 1.;

  tmvaWeights_ = edm::FileInPath(ps.getParameter<std::string>("tmvaWeights")).fullPath();
  tmvaMethod_ = ps.getParameter<std::string>("tmvaMethod");
  tmvaVariables_ = ps.getParameter<std::vector<std::string> >("tmvaVariables");
  tmvaSpectators_ = ps.getParameter<std::vector<std::string> >("tmvaSpectators");
  version_ = ps.getParameter<int>("version");
  reader_ = nullptr;
  edm::ParameterSet jetConfig = ps.getParameter<edm::ParameterSet>("JetIdParams");
  for (int i0 = 0; i0 < NWPs; i0++) {
    std::string lCutType = "Tight";
    if (i0 == PileupJetIdentifier::kMedium)
      lCutType = "Medium";
    if (i0 == PileupJetIdentifier::kLoose)
      lCutType = "Loose";
    for (int i1 = 0; i1 < 1; i1++) {
      std::vector<double> pt010 = jetConfig.getParameter<std::vector<double> >(("Pt010_" + lCutType).c_str());
      std::vector<double> pt1020 = jetConfig.getParameter<std::vector<double> >(("Pt1020_" + lCutType).c_str());
      std::vector<double> pt2030 = jetConfig.getParameter<std::vector<double> >(("Pt2030_" + lCutType).c_str());
      std::vector<double> pt3050 = jetConfig.getParameter<std::vector<double> >(("Pt3050_" + lCutType).c_str());
      for (int i2 = 0; i2 < NPts; i2++)
        mvacut_[i0][0][i2] = pt010[i2];
      for (int i2 = 0; i2 < NPts; i2++)
        mvacut_[i0][1][i2] = pt1020[i2];
      for (int i2 = 0; i2 < NPts; i2++)
        mvacut_[i0][2][i2] = pt2030[i2];
      for (int i2 = 0; i2 < NPts; i2++)
        mvacut_[i0][3][i2] = pt3050[i2];
    }
  }
  setup();
}

~MVAJetPuId()

MVAJetPuId::~MVAJetPuId

(

)

Definition at line 105 of file MVAJetPuId.cc.

References reader_.

                        {
  if (reader_) {
    delete reader_;
  }
}

Member Function Documentation

bookReader()

void MVAJetPuId::bookReader ( )

protected

Definition at line 125 of file MVAJetPuId.cc.

References cms::cuda::assert(), relativeConstraints::empty, dqmdumpme::first, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, reco::details::loadTMVAWeights(), reader_, tmvaMethod_, tmvaNames_, tmvaSpectators_, tmvaVariables_, tmvaWeights_, and variables_.

Referenced by runMva().

                            {
  reader_ = new TMVA::Reader("!Color:Silent");
  assert(!tmvaMethod_.empty() && !tmvaWeights_.empty());
  for (std::vector<std::string>::iterator it = tmvaVariables_.begin(); it != tmvaVariables_.end(); ++it) {
    if (tmvaNames_[*it].empty()) {
      tmvaNames_[*it] = *it;
    }
    reader_->AddVariable(*it, variables_[tmvaNames_[*it]].first);
  }
  for (std::vector<std::string>::iterator it = tmvaSpectators_.begin(); it != tmvaSpectators_.end(); ++it) {
    if (tmvaNames_[*it].empty()) {
      tmvaNames_[*it] = *it;
    }
    reader_->AddSpectator(*it, variables_[tmvaNames_[*it]].first);
  }
  reco::details::loadTMVAWeights(reader_, tmvaMethod_, tmvaWeights_);
}

computeCutIDflag()

int MVAJetPuId::computeCutIDflag	(	float	betaStarClassic,
		float	dR2Mean,
		float	nvtx,
		float	jetPt,
		float	jetEta
	)

computeIDflag() [1/2]

int MVAJetPuId::computeIDflag	(	float	mva,
		float	jetPt,
		float	jetEta
	)

Definition at line 175 of file MVAJetPuId.cc.

References getJetIdKey(), reco::btau::jetEta, reco::btau::jetPt, and beam_dqm_sourceclient-live_cfg::mva.

Referenced by runMva().

                                                                  {
  std::pair<int, int> jetIdKey = getJetIdKey(jetPt, jetEta);
  return computeIDflag(mva, jetIdKey.first, jetIdKey.second);
}

computeIDflag() [2/2]

int MVAJetPuId::computeIDflag	(	float	mva,
		int	ptId,
		int	etaId
	)

Definition at line 180 of file MVAJetPuId.cc.

References PileupJetIdentifier::kLoose, PileupJetIdentifier::kMedium, PileupJetIdentifier::kTight, beam_dqm_sourceclient-live_cfg::mva, and mvacut_.

                                                            {
  int idFlag(0);
  if (mva > mvacut_[PileupJetIdentifier::kTight][ptId][etaId])
    idFlag += 1 << PileupJetIdentifier::kTight;
  if (mva > mvacut_[PileupJetIdentifier::kMedium][ptId][etaId])
    idFlag += 1 << PileupJetIdentifier::kMedium;
  if (mva > mvacut_[PileupJetIdentifier::kLoose][ptId][etaId])
    idFlag += 1 << PileupJetIdentifier::kLoose;
  return idFlag;
}

computeIdVariables()

PileupJetIdentifier MVAJetPuId::computeIdVariables	(	const reco::Jet *	jet,
		float	jec,
		const reco::Vertex *	vtx,
		const reco::VertexCollection &	allvtx,
		double	rho,
		bool	calculateMva = false
	)

Definition at line 196 of file MVAJetPuId.cc.

References a, cms::cuda::assert(), Assign(), b, DummyCfis::c, reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMultiplicity(), PileupJetIdentifier::chgEMfrac_, PileupJetIdentifier::chgHadrfrac_, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), PileupJetIdentifier::d0_, dumpMFGeometry_cfg::delta, reco::deltaPhi(), reco::deltaR(), PileupJetIdentifier::dRLead2nd_, PileupJetIdentifier::dRMeanEm_, PileupJetIdentifier::dRMeanNeut_, l1ctLayer1_cff::dZ, spr::find(), PileupJetIdentifier::fourthFrac_, DivergingColor::frac, reco::PFCandidate::gamma, reco::PFJet::getPFConstituent(), reco::PFJet::getPFConstituents(), reco::PFCandidate::h0, mps_fire::i, internalId_, edm::Ptr< T >::isAvailable(), edm::Ref< C, T, F >::isAvailable(), edm::Ptr< T >::isNonnull(), edm::Ref< C, T, F >::isNonnull(), edm::Ptr< T >::isNull(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, jetMETDQMOfflineSource_cff::jec, metsig::jet, PileupJetIdentifier::jetM_, PileupJetIdentifier::jetPhi_, reco::btau::jetPt, PileupJetIdentifier::leadChEta_, PileupJetIdentifier::leadChPhi_, PileupJetIdentifier::leadChPt_, PileupJetIdentifier::leadEmEta_, PileupJetIdentifier::leadEmPhi_, PileupJetIdentifier::leadEmPt_, PileupJetIdentifier::leadEta_, PileupJetIdentifier::leadFrac_, PileupJetIdentifier::leadNeutEta_, PileupJetIdentifier::leadNeutPhi_, PileupJetIdentifier::leadNeutPt_, PileupJetIdentifier::leadPhi_, PileupJetIdentifier::leadPt_, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, SiStripPI::min, PileupJetIdentifier::neuEMfrac_, PileupJetIdentifier::neuHadrfrac_, reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), reco::PFJet::neutralMultiplicity(), DiDispStaMuonMonitor_cfi::pt, PileupJetIdentifier::pt2A_, PileupJetIdentifier::ptMean_, PileupJetIdentifier::ptRMS_, alignCSCRings::r, resetVariables(), rho, runMva(), PileupJetIdentifier::secondEta_, PileupJetIdentifier::secondFrac_, PileupJetIdentifier::secondPhi_, PileupJetIdentifier::secondPt_, SetPtEtaPhi(), jetUpdater_cfi::sort, mathSSE::sqrt(), PileupJetIdentifier::sumChPt_, PileupJetIdentifier::sumNePt_, TtFullHadEvtBuilder_cfi::sumPt, PileupJetIdentifier::sumPt_, edm::swap(), PileupJetIdentifier::thirdFrac_, L1BJetProducer_cff::vtx, and mps_merge::weight.

Referenced by MVAJetPuIdProducer::produce().

                                                                      {
  typedef std::vector<reco::PFCandidatePtr> constituents_type;
  typedef std::vector<reco::PFCandidatePtr>::iterator constituents_iterator;

  resetVariables();

  const reco::PFJet *pfjet = dynamic_cast<const reco::PFJet *>(jet);

  if (jec < 0.) {
    jec = 1.;
  }

  constituents_type constituents = pfjet->getPFConstituents();

  reco::PFCandidatePtr lLead, lSecond, lLeadNeut, lLeadEm, lLeadCh, lTrail;
  std::vector<float> frac, fracCh, fracEm, fracNeut;
  constexpr int ncones = 4;
  std::array<float, ncones> cones{{0.1, 0.2, 0.3, 0.4}};
  std::array<float *, ncones> coneFracs{
      {&internalId_.frac01_, &internalId_.frac02_, &internalId_.frac03_, &internalId_.frac04_}};
  TMatrixDSym covMatrix(2);
  covMatrix = 0.;

  reco::TrackRef impactTrack;
  float jetPt = jet->pt() / jec;  // use uncorrected pt for shape variables
  float sumPt = 0., sumPt2 = 0., sumPtCh = 0, sumPtNe = 0;
  float sum_deta = 0;
  float sum_dphi = 0;
  float sum_deta2 = 0;
  float sum_detadphi = 0;
  float sum_dphi2 = 0;
  SetPtEtaPhi(
      *jet, internalId_.jetPt_, internalId_.jetEta_, internalId_.jetPhi_);  // use corrected pt for jet kinematics
  internalId_.jetM_ = jet->mass();
  internalId_.rho_ = rho;  //allvtx.size();
  for (constituents_iterator it = constituents.begin(); it != constituents.end(); ++it) {
    reco::PFCandidatePtr &icand = *it;
    float candPt = icand->pt();
    float candPtFrac = candPt / jetPt;
    float candDr = reco::deltaR(**it, *jet);
    float candDeta = fabs((*it)->eta() - jet->eta());
    float candDphi = reco::deltaPhi(**it, *jet);
    float candPtDr = candPt * candDr;
    size_t icone = std::lower_bound(&cones[0], &cones[ncones], candDr) - &cones[0];
    float weight2 = candPt * candPt;

    if (lLead.isNull() || candPt > lLead->pt()) {
      lSecond = lLead;
      lLead = icand;
    } else if ((lSecond.isNull() || candPt > lSecond->pt()) && (candPt < lLead->pt())) {
      lSecond = icand;
    }

    //internalId_.dRMean_     += candPtDr;
    internalId_.dR2Mean_ += candPtDr * candPtDr;

    internalId_.ptD_ += candPt * candPt;
    sumPt += candPt;
    sumPt2 += candPt * candPt;
    sum_deta += candDeta * weight2;
    sum_dphi += candDphi * weight2;
    sum_deta2 += candDeta * candDeta * weight2;
    sum_detadphi += candDeta * candDphi * weight2;
    sum_dphi2 += candDphi * candDphi * weight2;
    //Teta         += candPt * candDR * candDeta;
    //Tphi         += candPt * candDR * candDphi;

    frac.push_back(candPtFrac);
    if (icone < ncones) {
      *coneFracs[icone] += candPt;
    }

    if (icand->particleId() == reco::PFCandidate::h0) {
      if (lLeadNeut.isNull() || candPt > lLeadNeut->pt()) {
        lLeadNeut = icand;
      }
      internalId_.dRMeanNeut_ += candPtDr;
      fracNeut.push_back(candPtFrac);
      sumPtNe += candPt;
    }

    if (icand->particleId() == reco::PFCandidate::gamma) {
      if (lLeadEm.isNull() || candPt > lLeadEm->pt()) {
        lLeadEm = icand;
      }
      internalId_.dRMeanEm_ += candPtDr;
      fracEm.push_back(candPtFrac);
      sumPtNe += candPt;
    }

    if (icand->trackRef().isNonnull() && icand->trackRef().isAvailable()) {
      if (lLeadCh.isNull() || candPt > lLeadCh->pt()) {
        lLeadCh = icand;
      }
      //internalId_.jetRchg_  += candPtDr;
      fracCh.push_back(candPtFrac);
      sumPtCh += candPt;
    }

    if (icand->trackRef().isNonnull() && icand->trackRef().isAvailable()) {
      float tkpt = icand->trackRef()->pt();
      bool inVtx0 =
          find(vtx->tracks_begin(), vtx->tracks_end(), reco::TrackBaseRef(icand->trackRef())) != vtx->tracks_end();
      bool inAnyOther = false;

      double dZ0 = fabs(icand->trackRef()->dz(vtx->position()));
      double dZ = dZ0;
      for (reco::VertexCollection::const_iterator vi = allvtx.begin(); vi != allvtx.end(); ++vi) {
        const reco::Vertex &iv = *vi;
        if (iv.isFake() || iv.ndof() < 4) {
          continue;
        }

        bool isVtx0 = (iv.position() - vtx->position()).r() < 0.02;

        if (!isVtx0 && !inAnyOther) {
          inAnyOther =
              find(iv.tracks_begin(), iv.tracks_end(), reco::TrackBaseRef(icand->trackRef())) != iv.tracks_end();
        }

        dZ = std::min(dZ, fabs(icand->trackRef()->dz(iv.position())));
      }
      if (inVtx0 && !inAnyOther) {
        internalId_.betaClassic_ += tkpt;
      } else if (!inVtx0 && inAnyOther) {
        internalId_.betaStarClassic_ += tkpt;
      }

      if (dZ0 < 0.2) {
        internalId_.beta_ += tkpt;
      } else if (dZ < 0.2) {
        internalId_.betaStar_ += tkpt;
      }
    }

    if (lTrail.isNull() || candPt < lTrail->pt()) {
      lTrail = icand;
    }
  }

  assert(lLead.isNonnull());
  if (lSecond.isNull()) {
    lSecond = lTrail;
  }
  if (lLeadNeut.isNull()) {
    lLeadNeut = lTrail;
  }
  if (lLeadEm.isNull()) {
    lLeadEm = lTrail;
  }
  if (lLeadCh.isNull()) {
    lLeadCh = lTrail;
  }
  impactTrack = lLeadCh->trackRef();

  internalId_.nCharged_ = pfjet->chargedMultiplicity();
  internalId_.nNeutrals_ = pfjet->neutralMultiplicity();
  internalId_.chgEMfrac_ = pfjet->chargedEmEnergy() / jet->energy();
  internalId_.neuEMfrac_ = pfjet->neutralEmEnergy() / jet->energy();
  internalId_.chgHadrfrac_ = pfjet->chargedHadronEnergy() / jet->energy();
  internalId_.neuHadrfrac_ = pfjet->neutralHadronEnergy() / jet->energy();

  if (impactTrack.isNonnull() && impactTrack.isAvailable()) {
    internalId_.d0_ = fabs(impactTrack->dxy(vtx->position()));
    internalId_.dZ_ = fabs(impactTrack->dz(vtx->position()));
  } else {
    internalId_.nParticles_ = constituents.size();
    SetPtEtaPhi(*lLead, internalId_.leadPt_, internalId_.leadEta_, internalId_.leadPhi_);
    SetPtEtaPhi(*lSecond, internalId_.secondPt_, internalId_.secondEta_, internalId_.secondPhi_);
    SetPtEtaPhi(*lLeadNeut, internalId_.leadNeutPt_, internalId_.leadNeutEta_, internalId_.leadNeutPhi_);
    SetPtEtaPhi(*lLeadEm, internalId_.leadEmPt_, internalId_.leadEmEta_, internalId_.leadEmPhi_);
    SetPtEtaPhi(*lLeadCh, internalId_.leadChPt_, internalId_.leadChEta_, internalId_.leadChPhi_);
    std::sort(frac.begin(), frac.end(), std::greater<float>());
    std::sort(fracCh.begin(), fracCh.end(), std::greater<float>());
    std::sort(fracEm.begin(), fracEm.end(), std::greater<float>());
    std::sort(fracNeut.begin(), fracNeut.end(), std::greater<float>());
    Assign(frac, internalId_.leadFrac_, internalId_.secondFrac_, internalId_.thirdFrac_, internalId_.fourthFrac_);

    //covMatrix(0,0) /= sumPt2;
    //covMatrix(0,1) /= sumPt2;
    //covMatrix(1,1) /= sumPt2;
    //covMatrix(1,0)  = covMatrix(0,1);
    //internalId_.etaW_ = sqrt(covMatrix(0,0));
    //internalId_.phiW_ = sqrt(covMatrix(1,1));
    //internalId_.jetW_ = 0.5*(internalId_.etaW_+internalId_.phiW_);
    //TVectorD eigVals(2); eigVals = TMatrixDSymEigen(covMatrix).GetEigenValues();
    //
    if (internalId_.majW_ < internalId_.minW_) {
      std::swap(internalId_.majW_, internalId_.minW_);
    }

    //internalId_.dRLeadCent_ = reco::deltaR(*jet,*lLead);
    if (lSecond.isNonnull()) {
      internalId_.dRLead2nd_ = reco::deltaR(*jet, *lSecond);
    }
    internalId_.dRMeanNeut_ /= jetPt;
    internalId_.dRMeanEm_ /= jetPt;
    //internalId_.jetRchg_   /= jetPt;
    internalId_.dR2Mean_ /= sumPt2;
    for (size_t ic = 0; ic < ncones; ++ic) {
      *coneFracs[ic] /= jetPt;
    }

    double ptMean = sumPt / internalId_.nParticles_;
    double ptRMS = 0;
    for (unsigned int i0 = 0; i0 < frac.size(); i0++) {
      ptRMS += (frac[i0] - ptMean) * (frac[i0] - ptMean);
    }
    ptRMS /= internalId_.nParticles_;
    ptRMS = sqrt(ptRMS);
    internalId_.jetRchg_ = internalId_.leadChPt_ / sumPt;
    internalId_.jetR_ = internalId_.leadPt_ / sumPt;

    internalId_.ptMean_ = ptMean;
    internalId_.ptRMS_ = ptRMS / jetPt;
    internalId_.pt2A_ = sqrt(internalId_.ptD_ / internalId_.nParticles_) / jetPt;
    internalId_.ptD_ = sqrt(internalId_.ptD_) / sumPt;
    internalId_.sumPt_ = sumPt;
    internalId_.sumChPt_ = sumPtCh;
    internalId_.sumNePt_ = sumPtNe;
    if (sumPt > 0) {
      internalId_.beta_ /= sumPt;
      internalId_.betaStar_ /= sumPt;
    } else {
      assert(internalId_.beta_ == 0. && internalId_.betaStar_ == 0.);
    }
    float ave_deta = sum_deta / sumPt2;
    float ave_dphi = sum_dphi / sumPt2;
    float ave_deta2 = sum_deta2 / sumPt2;
    float ave_dphi2 = sum_dphi2 / sumPt2;
    float a = ave_deta2 - ave_deta * ave_deta;
    float b = ave_dphi2 - ave_dphi * ave_dphi;
    float c = -(sum_detadphi / sumPt2 - ave_deta * ave_dphi);
    float axis1 = 0;
    float axis2 = 0;
    if ((((a - b) * (a - b) + 4 * c * c)) > 0) {
      float delta = sqrt(((a - b) * (a - b) + 4 * c * c));
      if (a + b + delta > 0) {
        axis1 = sqrt(0.5 * (a + b + delta));
      }
      if (a + b - delta > 0) {
        axis2 = sqrt(0.5 * (a + b - delta));
      }
    } else {
      axis1 = -1;
      axis2 = -1;
    }
    internalId_.majW_ = axis1;  //sqrt(fabs(eigVals(0)));
    internalId_.minW_ = axis2;  //sqrt(fabs(eigVals(1)));
    //compute Pull

    float ddetaR_sum(0.0), ddphiR_sum(0.0);
    for (int i = 0; i < internalId_.nParticles_; ++i) {
      reco::PFCandidatePtr part = pfjet->getPFConstituent(i);
      float weight = part->pt() * part->pt();
      float deta = part->eta() - jet->eta();
      float dphi = reco::deltaPhi(*part, *jet);
      float ddeta, ddphi, ddR;
      ddeta = deta - ave_deta;
      ddphi = reco::deltaPhi(dphi, ave_dphi);  //2*atan(tan((dphi - ave_dphi)/2.)) ;
      ddR = sqrt(ddeta * ddeta + ddphi * ddphi);
      ddetaR_sum += ddR * ddeta * weight;
      ddphiR_sum += ddR * ddphi * weight;
    }
    if (sumPt2 > 0) {
      float ddetaR_ave = ddetaR_sum / sumPt2;
      float ddphiR_ave = ddphiR_sum / sumPt2;
      internalId_.dRMean_ = sqrt(ddetaR_ave * ddetaR_ave + ddphiR_ave * ddphiR_ave);
    }
  }

  if (calculateMva) {
    runMva();
  }

  return PileupJetIdentifier(internalId_);
}

computeMva()

PileupJetIdentifier MVAJetPuId::computeMva

(

)

Definition at line 191 of file MVAJetPuId.cc.

References internalId_, and runMva().

Referenced by MVAJetPuIdProducer::produce().

                                           {
  runMva();
  return PileupJetIdentifier(internalId_);
}

dumpVariables()

std::string MVAJetPuId::dumpVariables

(

)

const

Definition at line 479 of file MVAJetPuId.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, MillePedeFileConverter_cfg::out, and variables_.

                                          {
  std::stringstream out;
  for (variables_list_t::const_iterator it = variables_.begin(); it != variables_.end(); ++it) {
    out << std::setw(15) << it->first << std::setw(3) << "=" << std::setw(5) << *it->second.first << " ("
        << std::setw(5) << it->second.second << ")" << std::endl;
  }
  return out.str();
}

getJetIdKey()

std::pair< int, int > MVAJetPuId::getJetIdKey	(	float	jetPt,
		float	jetEta
	)

Definition at line 156 of file MVAJetPuId.cc.

References reco::btau::jetEta, and reco::btau::jetPt.

Referenced by computeIDflag().

                                                                   {
  int ptId = 0;
  if (jetPt > 10 && jetPt < 20)
    ptId = 1;
  if (jetPt >= 20 && jetPt < 30)
    ptId = 2;
  if (jetPt >= 30)
    ptId = 3;

  int etaId = 0;
  if (fabs(jetEta) > 2.5 && fabs(jetEta) < 2.75)
    etaId = 1;
  if (fabs(jetEta) >= 2.75 && fabs(jetEta) < 3.0)
    etaId = 2;
  if (fabs(jetEta) >= 3.0 && fabs(jetEta) < 5.0)
    etaId = 3;
  return std::pair<int, int>(ptId, etaId);
}

getVariables()

const variables_list_t& MVAJetPuId::getVariables ( ) const

inline

Definition at line 52 of file MVAJetPuId.h.

References variables_.

{ return variables_; };

initVariables()

void MVAJetPuId::initVariables ( )

protected

Definition at line 499 of file MVAJetPuId.cc.

References HLT_2024v13_cff::beta, simBeamSpotPI::betaStar, d0, l1ctLayer1_cff::dZ, PileupJetIdentifier::idFlag_, INIT_VARIABLE, internalId_, reco::btau::jetEta, reco::btau::jetPhi, reco::btau::jetPt, LHEJetFilter_cfi::jetR, large_val, beam_dqm_sourceclient-live_cfg::mva, HLT_2024v13_cff::nCharged, jetsAK4_CHS_cff::ptD, rho, and TtFullHadEvtBuilder_cfi::sumPt.

Referenced by setup().

                               {
  internalId_.idFlag_ = 0;
  INIT_VARIABLE(mva, "", -100.);

  INIT_VARIABLE(jetPt, "jetPt", 0.);
  INIT_VARIABLE(jetEta, "jetEta", large_val);
  INIT_VARIABLE(jetPhi, "", large_val);
  INIT_VARIABLE(jetM, "", 0.);
  INIT_VARIABLE(nCharged, "nCharged", 0.);
  INIT_VARIABLE(nNeutrals, "", 0.);

  INIT_VARIABLE(chgEMfrac, "", 0.);
  INIT_VARIABLE(neuEMfrac, "", 0.);
  INIT_VARIABLE(chgHadrfrac, "", 0.);
  INIT_VARIABLE(neuHadrfrac, "", 0.);

  INIT_VARIABLE(d0, "", -1000.);
  INIT_VARIABLE(dZ, "", -1000.);
  INIT_VARIABLE(nParticles, "nParticles", 0.);

  INIT_VARIABLE(leadPt, "", 0.);
  INIT_VARIABLE(leadEta, "", large_val);
  INIT_VARIABLE(leadPhi, "", large_val);
  INIT_VARIABLE(secondPt, "", 0.);
  INIT_VARIABLE(secondEta, "", large_val);
  INIT_VARIABLE(secondPhi, "", large_val);
  INIT_VARIABLE(leadNeutPt, "", 0.);
  INIT_VARIABLE(leadNeutEta, "", large_val);

  INIT_VARIABLE(jetR, "jetR", 0.);
  INIT_VARIABLE(pull, "pull", 0.);
  INIT_VARIABLE(jetRchg, "jetRchg", 0.);
  INIT_VARIABLE(dR2Mean, "dR2Mean", 0.);

  INIT_VARIABLE(ptD, "ptD", 0.);
  INIT_VARIABLE(ptMean, "", 0.);
  INIT_VARIABLE(ptRMS, "", 0.);
  INIT_VARIABLE(pt2A, "", 0.);
  INIT_VARIABLE(ptDCh, "", 0.);
  INIT_VARIABLE(ptDNe, "", 0.);
  INIT_VARIABLE(sumPt, "", 0.);
  INIT_VARIABLE(sumChPt, "", 0.);
  INIT_VARIABLE(sumNePt, "", 0.);
  INIT_VARIABLE(secondFrac, "", 0.);
  INIT_VARIABLE(thirdFrac, "", 0.);
  INIT_VARIABLE(fourthFrac, "", 0.);
  INIT_VARIABLE(leadChFrac, "", 0.);
  INIT_VARIABLE(secondChFrac, "", 0.);
  INIT_VARIABLE(thirdChFrac, "", 0.);
  INIT_VARIABLE(fourthChFrac, "", 0.);
  INIT_VARIABLE(leadNeutFrac, "", 0.);
  INIT_VARIABLE(secondNeutFrac, "", 0.);
  INIT_VARIABLE(thirdNeutFrac, "", 0.);
  INIT_VARIABLE(fourthNeutFrac, "", 0.);
  INIT_VARIABLE(leadEmFrac, "", 0.);
  INIT_VARIABLE(secondEmFrac, "", 0.);
  INIT_VARIABLE(thirdEmFrac, "", 0.);
  INIT_VARIABLE(fourthEmFrac, "", 0.);
  INIT_VARIABLE(jetW, "", 1.);
  INIT_VARIABLE(etaW, "", 1.);
  INIT_VARIABLE(phiW, "", 1.);
  INIT_VARIABLE(majW, "majW", 1.);
  INIT_VARIABLE(minW, "minW", 1.);
  INIT_VARIABLE(frac01, "frac01", 0.);
  INIT_VARIABLE(frac02, "frac02", 0.);
  INIT_VARIABLE(frac03, "frac03", 0.);
  INIT_VARIABLE(frac04, "frac04", 0.);

  INIT_VARIABLE(beta, "beta", 0.);
  INIT_VARIABLE(betaStar, "betaStar", 0.);
  INIT_VARIABLE(betaClassic, "betaClassic", 0.);
  INIT_VARIABLE(betaStarClassic, "betaStarClassic", 0.);
  INIT_VARIABLE(rho, "rho", 0.);
}

method()

const std::string MVAJetPuId::method ( ) const

inline

Definition at line 41 of file MVAJetPuId.h.

References tmvaMethod_.

{ return tmvaMethod_; }

resetVariables()

void MVAJetPuId::resetVariables ( )

protected

Definition at line 488 of file MVAJetPuId.cc.

References PileupJetIdentifier::idFlag_, internalId_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and variables_.

Referenced by computeIdVariables().

                                {
  internalId_.idFlag_ = 0;
  for (variables_list_t::iterator it = variables_.begin(); it != variables_.end(); ++it) {
    *it->second.first = it->second.second;
  }
}

runMva()

void MVAJetPuId::runMva ( )

protected

Definition at line 145 of file MVAJetPuId.cc.

References bookReader(), computeIDflag(), PileupJetIdentifier::idFlag_, internalId_, PileupJetIdentifier::mva_, reader_, and tmvaMethod_.

Referenced by computeIdVariables(), and computeMva().

                        {
  if (!reader_) {
    bookReader();
  }
  if (fabs(internalId_.jetEta_) < 5.0)
    internalId_.mva_ = reader_->EvaluateMVA(tmvaMethod_.c_str());
  if (fabs(internalId_.jetEta_) >= 5.0)
    internalId_.mva_ = -2.;
  internalId_.idFlag_ = computeIDflag(internalId_.mva_, internalId_.jetPt_, internalId_.jetEta_);
}

set()

void MVAJetPuId::set

(

const PileupJetIdentifier &

id

)

Definition at line 143 of file MVAJetPuId.cc.

References l1ctLayer2EG_cff::id, and internalId_.

Referenced by MVAJetPuIdProducer::produce().

{ internalId_ = id; }

setup()

void MVAJetPuId::setup ( )

protected

Definition at line 66 of file MVAJetPuId.cc.

References initVariables(), tmvaNames_, and tmvaVariables_.

Referenced by MVAJetPuId().

                       {
  initVariables();

  tmvaVariables_.clear();
  tmvaVariables_.push_back("rho");
  tmvaVariables_.push_back("nParticles");
  tmvaVariables_.push_back("nCharged");
  tmvaVariables_.push_back("majW");
  tmvaVariables_.push_back("minW");
  tmvaVariables_.push_back("frac01");
  tmvaVariables_.push_back("frac02");
  tmvaVariables_.push_back("frac03");
  tmvaVariables_.push_back("frac04");
  tmvaVariables_.push_back("ptD");
  tmvaVariables_.push_back("beta");
  tmvaVariables_.push_back("betaStar");
  tmvaVariables_.push_back("dR2Mean");
  tmvaVariables_.push_back("pull");
  tmvaVariables_.push_back("jetR");
  tmvaVariables_.push_back("jetRchg");

  tmvaNames_["rho"] = "rho";
  tmvaNames_["nParticles"] = "nParticles";
  tmvaNames_["nCharged"] = "nCharged";
  tmvaNames_["majW"] = "majW";
  tmvaNames_["minW"] = "minW";
  tmvaNames_["frac01"] = "frac01";
  tmvaNames_["frac02"] = "frac02";
  tmvaNames_["frac03"] = "frac03";
  tmvaNames_["frac04"] = "frac04";
  tmvaNames_["ptD"] = "ptD";
  tmvaNames_["beta"] = "beta";
  tmvaNames_["betaStar"] = "betaStar";
  tmvaNames_["dR2Mean"] = "dR2Mean";
  tmvaNames_["pull"] = "pull";
  tmvaNames_["jetR"] = "jetR";
  tmvaNames_["jetRchg"] = "jetRchg";
}

Member Data Documentation

betaStarCut_

Float_t MVAJetPuId::betaStarCut_[NWPs][NEtas][NPts]

protected

Definition at line 75 of file MVAJetPuId.h.

cutBased_

bool MVAJetPuId::cutBased_

protected

Definition at line 72 of file MVAJetPuId.h.

impactParTkThreshod_

Float_t MVAJetPuId::impactParTkThreshod_

protected

Definition at line 71 of file MVAJetPuId.h.

Referenced by MVAJetPuId().

internalId_

PileupJetIdentifier MVAJetPuId::internalId_

protected

Definition at line 61 of file MVAJetPuId.h.

Referenced by computeIdVariables(), computeMva(), initVariables(), resetVariables(), runMva(), and set().

mvacut_

Float_t MVAJetPuId::mvacut_[NWPs][NEtas][NPts]

protected

Definition at line 73 of file MVAJetPuId.h.

Referenced by computeIDflag(), and MVAJetPuId().

NEtas

constexpr int MVAJetPuId::NEtas = 4

staticprivate

Definition at line 19 of file MVAJetPuId.h.

NPts

constexpr int MVAJetPuId::NPts = 4

staticprivate

Definition at line 18 of file MVAJetPuId.h.

Referenced by MVAJetPuId().

NWPs

constexpr int MVAJetPuId::NWPs = 3

staticprivate

Definition at line 17 of file MVAJetPuId.h.

Referenced by MVAJetPuId().

reader_

TMVA::Reader* MVAJetPuId::reader_

protected

Definition at line 64 of file MVAJetPuId.h.

Referenced by bookReader(), MVAJetPuId(), runMva(), and ~MVAJetPuId().

rmsCut_

Float_t MVAJetPuId::rmsCut_[NWPs][NEtas][NPts]

protected

Definition at line 74 of file MVAJetPuId.h.

tmvaMethod_

std::string MVAJetPuId::tmvaMethod_

protected

Definition at line 65 of file MVAJetPuId.h.

Referenced by bookReader(), method(), MVAJetPuId(), and runMva().

tmvaNames_

std::map<std::string, std::string> MVAJetPuId::tmvaNames_

protected

Definition at line 68 of file MVAJetPuId.h.

Referenced by bookReader(), and setup().

tmvaSpectators_

std::vector<std::string> MVAJetPuId::tmvaSpectators_

protected

Definition at line 67 of file MVAJetPuId.h.

Referenced by bookReader(), and MVAJetPuId().

tmvaVariables_

std::vector<std::string> MVAJetPuId::tmvaVariables_

protected

Definition at line 66 of file MVAJetPuId.h.

Referenced by bookReader(), MVAJetPuId(), and setup().

tmvaWeights_

std::string MVAJetPuId::tmvaWeights_

protected

Definition at line 65 of file MVAJetPuId.h.

Referenced by bookReader(), and MVAJetPuId().

variables_

variables_list_t MVAJetPuId::variables_

protected

Definition at line 62 of file MVAJetPuId.h.

Referenced by bookReader(), dumpVariables(), getVariables(), and resetVariables().

version_

Int_t MVAJetPuId::version_

protected

Definition at line 70 of file MVAJetPuId.h.

Referenced by MVAJetPuId().