HiGammaJetSignalDef Class Reference

#include <HiPhotonType.h>

Public Member Functions
double	getDeltaPhi (const reco::Candidate &track1, const reco::Candidate &track2)
double	getDeltaR (const reco::Candidate &track1, const reco::Candidate &track2)
	HiGammaJetSignalDef ()
	HiGammaJetSignalDef (const reco::GenParticleCollection *sigPartic)
bool	IsIsolated (const reco::GenParticle &pp)
bool	IsIsolatedJP (const reco::GenParticle &pp)
bool	IsIsolatedPP (const reco::GenParticle &pp)

Public Attributes
double	PI

Private Attributes
const reco::GenParticleCollection *	fSigParticles

Detailed Description

Definition at line 23 of file HiPhotonType.h.

Constructor & Destructor Documentation

HiGammaJetSignalDef::HiGammaJetSignalDef

(

)

Definition at line 50 of file HiPhotonType.cc.

References PI.

                                         : fSigParticles(nullptr) {
  PI = 3.141592653589793238462643383279502884197169399375105820974945;
}

HiGammaJetSignalDef::HiGammaJetSignalDef

(

const reco::GenParticleCollection *

sigPartic

)

Definition at line 53 of file HiPhotonType.cc.

References fSigParticles, and PI.

                                                                                      : fSigParticles(nullptr) {
  using namespace std;

  fSigParticles = sigParticles;
  PI = 3.141592653589793238462643383279502884197169399375105820974945;
}

Member Function Documentation

double HiGammaJetSignalDef::getDeltaPhi	(	const reco::Candidate &	track1,
		const reco::Candidate &	track2
	)

Definition at line 201 of file HiPhotonType.cc.

References reco::Candidate::phi(), and PI.

                                                                                                {
  double dPhi = track1.phi() - track2.phi();

  while (dPhi >= PI)
    dPhi -= (2.0 * PI);
  while (dPhi < (-1.0 * PI))
    dPhi += (2.0 * PI);

  return dPhi;
}

double HiGammaJetSignalDef::getDeltaR	(	const reco::Candidate &	track1,
		const reco::Candidate &	track2
	)

Definition at line 189 of file HiPhotonType.cc.

References reco::Candidate::eta(), reco::Candidate::phi(), PI, and mathSSE::sqrt().

Referenced by IsIsolated(), IsIsolatedJP(), and IsIsolatedPP().

                                                                                              {
  double dEta = track1.eta() - track2.eta();
  double dPhi = track1.phi() - track2.phi();

  while (dPhi >= PI)
    dPhi -= (2.0 * PI);
  while (dPhi < (-1.0 * PI))
    dPhi += (2.0 * PI);

  return sqrt(dEta * dEta + dPhi * dPhi);
}

bool HiGammaJetSignalDef::IsIsolated

(

const reco::GenParticle &

pp

)

Definition at line 60 of file HiPhotonType.cc.

References funct::abs(), reco::GenParticle::collisionId(), reco::LeafCandidate::et(), fSigParticles, getDeltaR(), mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::pdgId(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), dt_dqm_sourceclient_common_cff::reco, and reco::LeafCandidate::status().

                                                              {
  using namespace std;
  using namespace edm;
  using namespace reco;
  // Check if a given particle is isolated.

  double EtCone = 0;
  double EtPhoton = 0;
  double PtMaxHadron = 0;
  double cone = 0.5;
  const int maxindex = (int)fSigParticles->size();
  for (int i = 0; i < maxindex; ++i) {
    const GenParticle &p = (*fSigParticles)[i];

    if (p.status() != 1)
      continue;
    if (p.collisionId() != pp.collisionId())
      continue;

    int apid = abs(p.pdgId());
    if (apid > 11 && apid < 20)
      continue;  //get rid of muons and neutrinos

    if (getDeltaR(p, pp) > cone)
      continue;

    double et = p.et();
    double pt = p.pt();
    EtCone += et;
    bool isHadron = false;
    if (fabs(pp.pdgId()) > 100 && fabs(pp.pdgId()) != 310)
      isHadron = true;

    if (apid > 100 && apid != 310 && pt > PtMaxHadron) {
      if ((isHadron == true) && (fabs(pp.pt() - pt) < 0.001) && (pp.pdgId() == p.pdgId()))
        continue;
      else
        PtMaxHadron = pt;
    }
  }
  EtPhoton = pp.et();

  // isolation cuts
  if (EtCone - EtPhoton > 5 + EtPhoton / 20)
    return kFALSE;
  if (PtMaxHadron > 4.5 + EtPhoton / 40)
    return kFALSE;

  return kTRUE;
}

bool HiGammaJetSignalDef::IsIsolatedJP

(

const reco::GenParticle &

pp

)

Definition at line 153 of file HiPhotonType.cc.

References reco::GenParticle::collisionId(), reco::LeafCandidate::et(), fSigParticles, getDeltaR(), mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, dt_dqm_sourceclient_common_cff::reco, and reco::LeafCandidate::status().

                                                                {
  using namespace std;
  using namespace edm;
  using namespace reco;

  // Check if a given particle is isolated.

  double EtCone = 0;
  double EtPhoton = 0;
  double cone = 0.4;

  const int maxindex = (int)fSigParticles->size();
  for (int i = 0; i < maxindex; ++i) {
    const GenParticle &p = (*fSigParticles)[i];

    if (p.status() != 1)
      continue;
    if (p.collisionId() != pp.collisionId())
      continue;

    if (getDeltaR(p, pp) > cone)
      continue;

    double et = p.et();
    //  double pt = p.pt();
    EtCone += et;
  }
  EtPhoton = pp.et();

  // isolation cuts

  if (EtCone - EtPhoton > 5)
    return kFALSE;
  return kTRUE;
}

bool HiGammaJetSignalDef::IsIsolatedPP

(

const reco::GenParticle &

pp

)

Definition at line 111 of file HiPhotonType.cc.

References reco::GenParticle::collisionId(), reco::LeafCandidate::et(), fSigParticles, getDeltaR(), mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, dt_dqm_sourceclient_common_cff::reco, and reco::LeafCandidate::status().

                                                                {
  using namespace std;
  using namespace edm;
  using namespace reco;

  // Check if a given particle is isolated.

  double EtCone = 0;
  double EtPhoton = 0;
  double cone = 0.4;

  const int maxindex = (int)fSigParticles->size();
  for (int i = 0; i < maxindex; ++i) {
    const GenParticle &p = (*fSigParticles)[i];

    if (p.status() != 1)
      continue;
    if (p.collisionId() != pp.collisionId())
      continue;

    //int apid= abs(p.pdgId());
    //  if(apid>11 &&  apid<20)
    //  continue; //get rid of muons and neutrinos

    if (getDeltaR(p, pp) > cone)
      continue;
    double et = p.et();
    //  double pt = p.pt();
    EtCone += et;
  }

  EtPhoton = pp.et();

  // isolation cuts

  if (EtCone - EtPhoton > 5)
    return kFALSE;
  //if(PtMaxHadron > 4.5+EtPhoton/40)
  //  return kFALSE;

  return kTRUE;
}

Member Data Documentation

const reco::GenParticleCollection* HiGammaJetSignalDef::fSigParticles

private

Definition at line 38 of file HiPhotonType.h.

Referenced by HiGammaJetSignalDef(), IsIsolated(), IsIsolatedJP(), and IsIsolatedPP().

double HiGammaJetSignalDef::PI

Definition at line 35 of file HiPhotonType.h.

Referenced by getDeltaPhi(), getDeltaR(), and HiGammaJetSignalDef().