JetMatchingTools Class Reference

#include <JetMatchingTools.h>

Classes
struct	JetConstituent

Public Member Functions
int	generatorId (unsigned fSimTrackId)
	Generator ID. More...
std::vector< const CaloTower * >	getCaloTowers (int fGeneratorId)
	CaloTowers. More...
std::vector< JetConstituent >	getConstituentHits (const CaloTower &fTower)
	get CaloRecHits contributing to the tower More...
std::vector< DetId >	getConstituentIds (const CaloTower &fTower)
	get cells contributing to the tower More...
std::vector< const CaloTower * >	getConstituents (const reco::CaloJet &fJet)
	get towers contributing to CaloJet More...
const EBRecHitCollection *	getEBRecHitCollection ()
const edm::PCaloHitContainer *	getEBSimHitCollection ()
const EERecHitCollection *	getEERecHitCollection ()
const edm::PCaloHitContainer *	getEESimHitCollection ()
const reco::GenParticle *	getGenParticle (int fGeneratorId)
	GenParticle. More...
std::vector< const reco::GenParticle * >	getGenParticles (const reco::CaloJet &fJet, bool fVerbose=true)
	GenParticles for CaloJet. More...
std::vector< const reco::GenParticle * >	getGenParticles (const reco::GenJet &fJet)
	GenParticles for GenJet. More...
const reco::CandidateCollection *	getGenParticlesCollection ()
const HBHERecHitCollection *	getHBHERecHitCollection ()
const edm::PCaloHitContainer *	getHcalSimHitCollection ()
const HFRecHitCollection *	getHFRecHitCollection ()
const HORecHitCollection *	getHORecHitCollection ()
std::vector< const PCaloHit * >	getPCaloHits (DetId fId)
	get PCaloHits contributing to the detId More...
std::vector< const PCaloHit * >	getPCaloHits (int fGeneratorId)
	CaloSimHits. More...
const edm::SimTrackContainer *	getSimTrackCollection ()
const edm::SimVertexContainer *	getSimVertexCollection ()
const SimTrack *	getTrack (unsigned fSimTrackId)
	convert trackId to SimTrack More...
int	getTrackId (const PCaloHit &fHit)
	GEANT track ID. More...
	JetMatchingTools (const edm::Event &fEvent, edm::ConsumesCollector &&iC)
double	lostEnergyFraction (const reco::CaloJet &fJet)
	energy in broken links More...
double	overlapEnergyFraction (const std::vector< const reco::GenParticle *> &fObject, const std::vector< const reco::GenParticle *> &fReference) const
	energy overlap More...
	~JetMatchingTools ()

Private Attributes
edm::EDGetTokenT< reco::CandidateCollection >	input_cands_token_
edm::EDGetTokenT< EBRecHitCollection >	input_ebrechits_token_
edm::EDGetTokenT< EERecHitCollection >	input_eerechits_token_
edm::EDGetTokenT< HBHERecHitCollection >	input_hbherechits_token_
edm::EDGetTokenT< HFRecHitCollection >	input_hfrechits_token_
edm::EDGetTokenT< HORecHitCollection >	input_horechits_token_
edm::EDGetTokenT< edm::PCaloHitContainer >	input_pcalohits_ebcal_token_
edm::EDGetTokenT< edm::PCaloHitContainer >	input_pcalohits_eecal_token_
edm::EDGetTokenT< edm::PCaloHitContainer >	input_pcalohits_hcal_token_
edm::EDGetTokenT< edm::SimTrackContainer >	input_simtrack_token_
edm::EDGetTokenT< edm::SimVertexContainer >	input_simvertex_token_
const EBRecHitCollection *	mEBRecHitCollection
const edm::PCaloHitContainer *	mEBSimHitCollection
const EERecHitCollection *	mEERecHitCollection
const edm::PCaloHitContainer *	mEESimHitCollection
const edm::Event *	mEvent
const reco::CandidateCollection *	mGenParticleCollection
const HBHERecHitCollection *	mHBHERecHitCollection
const edm::PCaloHitContainer *	mHcalSimHitCollection
const HFRecHitCollection *	mHFRecHitCollection
const HORecHitCollection *	mHORecHitCollection
const edm::SimTrackContainer *	mSimTrackCollection
const edm::SimVertexContainer *	mSimVertexCollection

Detailed Description

Definition at line 28 of file JetMatchingTools.h.

Constructor & Destructor Documentation

JetMatchingTools()

JetMatchingTools::JetMatchingTools	(	const edm::Event &	fEvent,
		edm::ConsumesCollector &&	iC
	)

Definition at line 37 of file JetMatchingTools.cc.

References input_cands_token_, input_ebrechits_token_, input_eerechits_token_, input_hbherechits_token_, input_hfrechits_token_, input_horechits_token_, input_pcalohits_ebcal_token_, input_pcalohits_eecal_token_, input_pcalohits_hcal_token_, input_simtrack_token_, input_simvertex_token_, and HLT_2022v12_cff::InputTag.

    : mEvent(&fEvent),
      mEBRecHitCollection(nullptr),
      mEERecHitCollection(nullptr),
      mHBHERecHitCollection(nullptr),
      mHORecHitCollection(nullptr),
      mHFRecHitCollection(nullptr),
      mEBSimHitCollection(nullptr),
      mEESimHitCollection(nullptr),
      mHcalSimHitCollection(nullptr),
      mSimTrackCollection(nullptr),
      mSimVertexCollection(nullptr),
      mGenParticleCollection(nullptr) {
  input_ebrechits_token_ = iC.mayConsume<EBRecHitCollection>(edm::InputTag("ecalRecHit:EcalRecHitsEB"));
  input_eerechits_token_ = iC.mayConsume<EERecHitCollection>(edm::InputTag("ecalRecHit:EcalRecHitsEE"));
  input_hbherechits_token_ = iC.mayConsume<HBHERecHitCollection>(edm::InputTag("hbhereco"));
  input_horechits_token_ = iC.mayConsume<HORecHitCollection>(edm::InputTag("horeco"));
  input_hfrechits_token_ = iC.mayConsume<HFRecHitCollection>(edm::InputTag("hfreco"));
  input_pcalohits_ebcal_token_ = iC.mayConsume<edm::PCaloHitContainer>(edm::InputTag("g4SimHits:EcalHitsEB"));
  input_pcalohits_eecal_token_ = iC.mayConsume<edm::PCaloHitContainer>(edm::InputTag("g4SimHits:EcalHitsEE"));
  input_pcalohits_hcal_token_ = iC.mayConsume<edm::PCaloHitContainer>(edm::InputTag("g4SimHits:HcalHits"));
  input_simtrack_token_ = iC.mayConsume<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
  input_simvertex_token_ = iC.mayConsume<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
  input_cands_token_ = iC.mayConsume<reco::CandidateCollection>(edm::InputTag("genParticleCandidates"));
}

~JetMatchingTools()

JetMatchingTools::~JetMatchingTools

(

)

Definition at line 63 of file JetMatchingTools.cc.

{}

Member Function Documentation

generatorId()

int JetMatchingTools::generatorId

(

unsigned

fSimTrackId

)

Generator ID.

Definition at line 240 of file JetMatchingTools.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, getSimVertexCollection(), getTrack(), HLT_2022v12_cff::track, and bphysicsOniaDQM_cfi::vertex.

Referenced by getGenParticles(), and lostEnergyFraction().

                                                      {
  const SimTrack* track = getTrack(fSimTrackId);
  if (!track)
    return -1;
  while (track->noGenpart()) {
    if (track->noVertex()) {
      std::cerr << "JetMatchingTools::generatorId-> No vertex for track " << *track << std::endl;
      return -1;
    }
    const SimVertex* vertex = &((*getSimVertexCollection())[track->vertIndex()]);
    if (vertex->noParent()) {
      std::cerr << "JetMatchingTools::generatorId-> No track for vertex " << *vertex << std::endl;
      return -1;
    }
    track = getTrack(vertex->parentIndex());
  }
  return track->genpartIndex();
}

getCaloTowers()

std::vector<const CaloTower*> JetMatchingTools::getCaloTowers

(

int

fGeneratorId

)

CaloTowers.

getConstituentHits()

std::vector< JetMatchingTools::JetConstituent > JetMatchingTools::getConstituentHits

(

const CaloTower &

fTower

)

get CaloRecHits contributing to the tower

Definition at line 164 of file JetMatchingTools.cc.

References cms::cuda::assert(), EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, CaloTower::constituent(), CaloTower::constituentsSize(), DetId::Ecal, EcalBarrel, EcalEndcap, getEBRecHitCollection(), getEERecHitCollection(), getHBHERecHitCollection(), getHFRecHitCollection(), getHORecHitCollection(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, and mps_fire::result.

Referenced by lostEnergyFraction().

                                                                                                      {
  std::vector<JetConstituent> result;

  for (unsigned i = 0; i < fTower.constituentsSize(); ++i) {
    DetId id = fTower.constituent(i);

    if (id.det() == DetId::Ecal) {
      const EcalRecHit* hit = nullptr;

      if ((EcalSubdetector)id.subdetId() == EcalBarrel) {
        hit = getHit(*getEBRecHitCollection(), id);
      } else if ((EcalSubdetector)id.subdetId() == EcalEndcap) {
        hit = getHit(*getEERecHitCollection(), id);
      }

      assert(hit != nullptr);
      if (hit)
        result.push_back(JetConstituent(*hit));
      else
        std::cerr << "Can not find rechit for id " << id.rawId() << std::endl;
    } else if (id.det() == DetId::Hcal) {
      const CaloRecHit* hit = nullptr;

      if ((HcalSubdetector)id.subdetId() == HcalBarrel || (HcalSubdetector)id.subdetId() == HcalEndcap) {
        hit = getHit(*getHBHERecHitCollection(), id);
      } else if ((HcalSubdetector)id.subdetId() == HcalOuter) {
        hit = getHit(*getHORecHitCollection(), id);
      }
      if ((HcalSubdetector)id.subdetId() == HcalForward) {
        hit = getHit(*getHFRecHitCollection(), id);
      }

      if (hit)
        result.push_back(JetConstituent(*hit));
      else
        std::cerr << "Can not find rechit for id " << id.rawId() << std::endl;
    }
  }

  return result;
}

getConstituentIds()

std::vector< DetId > JetMatchingTools::getConstituentIds

(

const CaloTower &

fTower

)

get cells contributing to the tower

Definition at line 207 of file JetMatchingTools.cc.

References CaloTower::constituent(), CaloTower::constituentsSize(), mps_fire::i, and mps_fire::result.

Referenced by getGenParticles().

                                                                            {
  std::vector<DetId> result;
  for (unsigned i = 0; i < fTower.constituentsSize(); ++i) {
    DetId id = fTower.constituent(i);
    result.push_back(id);
  }
  return result;
}

getConstituents()

std::vector< const CaloTower * > JetMatchingTools::getConstituents

(

const reco::CaloJet &

fJet

)

get towers contributing to CaloJet

Definition at line 155 of file JetMatchingTools.cc.

References reco::CaloJet::getCaloConstituents(), mps_fire::i, and mps_fire::result.

Referenced by getGenParticles(), and lostEnergyFraction().

                                                                                     {
  std::vector<const CaloTower*> result;
  std::vector<CaloTowerPtr> constituents = fJet.getCaloConstituents();
  for (unsigned i = 0; i < constituents.size(); ++i)
    result.push_back(&*(constituents[i]));
  return result;
}

getEBRecHitCollection()

const EBRecHitCollection * JetMatchingTools::getEBRecHitCollection

(

)

Definition at line 65 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_ebrechits_token_, mEBRecHitCollection, mEvent, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by getConstituentHits().

                                                                  {
  if (!mEBRecHitCollection) {
    edm::Handle<EBRecHitCollection> recHits;
    mEvent->getByToken(input_ebrechits_token_, recHits);
    mEBRecHitCollection = &*recHits;
  }
  return mEBRecHitCollection;
}

getEBSimHitCollection()

const PCaloHitContainer * JetMatchingTools::getEBSimHitCollection

(

)

Definition at line 105 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_pcalohits_ebcal_token_, mEBSimHitCollection, mEvent, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by getPCaloHits().

                                                                 {
  if (!mEBSimHitCollection) {
    edm::Handle<PCaloHitContainer> simHits;
    mEvent->getByToken(input_pcalohits_ebcal_token_, simHits);
    mEBSimHitCollection = &*simHits;
  }
  return mEBSimHitCollection;
}

getEERecHitCollection()

const EERecHitCollection * JetMatchingTools::getEERecHitCollection

(

)

Definition at line 73 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_eerechits_token_, mEERecHitCollection, mEvent, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by getConstituentHits().

                                                                  {
  if (!mEERecHitCollection) {
    edm::Handle<EERecHitCollection> recHits;
    mEvent->getByToken(input_eerechits_token_, recHits);
    mEERecHitCollection = &*recHits;
  }
  return mEERecHitCollection;
}

getEESimHitCollection()

const PCaloHitContainer * JetMatchingTools::getEESimHitCollection

(

)

Definition at line 113 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_pcalohits_eecal_token_, mEESimHitCollection, mEvent, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by getPCaloHits().

                                                                 {
  if (!mEESimHitCollection) {
    edm::Handle<PCaloHitContainer> simHits;
    mEvent->getByToken(input_pcalohits_eecal_token_, simHits);
    mEESimHitCollection = &*simHits;
  }
  return mEESimHitCollection;
}

getGenParticle()

const reco::GenParticle * JetMatchingTools::getGenParticle

(

int

fGeneratorId

)

GenParticle.

Definition at line 260 of file JetMatchingTools.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, reco::GenJet::genParticle(), getGenParticlesCollection(), edm::OwnVector< T, P >::size(), and findQualityFiles::size.

Referenced by getGenParticles().

                                                                        {
  if (fGeneratorId > int(getGenParticlesCollection()->size())) {
    std::cerr << "JetMatchingTools::getGenParticle-> requested index " << fGeneratorId
              << " is grater then container size " << getGenParticlesCollection()->size() << std::endl;
    return nullptr;
  }
  return reco::GenJet::genParticle(
      &(*getGenParticlesCollection())[fGeneratorId - 1]);  // knowhow: index is shifted by 1
}

getGenParticles() [1/2]

std::vector< const reco::GenParticle * > JetMatchingTools::getGenParticles	(	const reco::CaloJet &	fJet,
		bool	fVerbose = true
	)

GenParticles for CaloJet.

Definition at line 271 of file JetMatchingTools.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, generatorId(), getConstituentIds(), getConstituents(), getGenParticle(), getPCaloHits(), getTrackId(), phits, mps_fire::result, and HLT_2022v12_cff::towers.

                                                                                                          {
  std::set<const reco::GenParticle*> result;
  // follow the chain
  std::vector<const CaloTower*> towers = getConstituents(fJet);
  for (unsigned itower = 0; itower < towers.size(); ++itower) {
    std::vector<DetId> detids = getConstituentIds(*(towers[itower]));
    for (unsigned iid = 0; iid < detids.size(); ++iid) {
      std::vector<const PCaloHit*> phits = getPCaloHits(detids[iid]);
      for (unsigned iphit = 0; iphit < phits.size(); ++iphit) {
        int trackId = getTrackId(*(phits[iphit]));
        if (trackId >= 0) {
          int genId = generatorId(trackId);
          if (genId >= 0) {
            const reco::GenParticle* genPart = getGenParticle(genId);
            if (genPart) {
              result.insert(genPart);
            } else if (fVerbose) {
              std::cerr << "JetMatchingTools::getGenParticles-> Can not convert genId " << genId << " to GenParticle"
                        << std::endl;
            }
          } else if (fVerbose) {
            std::cerr << "JetMatchingTools::getGenParticles-> Can not convert trackId " << trackId << " to genId"
                      << std::endl;
          }
        } else if (fVerbose) {
          std::cerr << "JetMatchingTools::getGenParticles-> Unknown trackId for PCaloHit " << *(phits[iphit])
                    << std::endl;
        }
      }
    }
  }
  return std::vector<const reco::GenParticle*>(result.begin(), result.end());
}

getGenParticles() [2/2]

std::vector< const reco::GenParticle * > JetMatchingTools::getGenParticles

(

const reco::GenJet &

fJet

)

GenParticles for GenJet.

Definition at line 306 of file JetMatchingTools.cc.

References reco::GenJet::getGenConstituents().

                                                                                          {
  return fJet.getGenConstituents();
}

getGenParticlesCollection()

const reco::CandidateCollection * JetMatchingTools::getGenParticlesCollection

(

)

Definition at line 145 of file JetMatchingTools.cc.

References edm::Event::getByToken(), patZpeak::handle, input_cands_token_, mEvent, and mGenParticleCollection.

Referenced by getGenParticle().

                                                                           {
  if (!mGenParticleCollection) {
    edm::Handle<reco::CandidateCollection> handle;
    mEvent->getByToken(input_cands_token_, handle);
    mGenParticleCollection = &*handle;
  }
  return mGenParticleCollection;
}

getHBHERecHitCollection()

const HBHERecHitCollection * JetMatchingTools::getHBHERecHitCollection

(

)

Definition at line 81 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_hbherechits_token_, mEvent, mHBHERecHitCollection, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by getConstituentHits().

                                                                      {
  if (!mHBHERecHitCollection) {
    edm::Handle<HBHERecHitCollection> recHits;
    mEvent->getByToken(input_hbherechits_token_, recHits);
    mHBHERecHitCollection = &*recHits;
  }
  return mHBHERecHitCollection;
}

getHcalSimHitCollection()

const PCaloHitContainer * JetMatchingTools::getHcalSimHitCollection

(

)

Definition at line 121 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_pcalohits_hcal_token_, mEvent, mHcalSimHitCollection, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by getPCaloHits().

                                                                   {
  if (!mHcalSimHitCollection) {
    edm::Handle<PCaloHitContainer> simHits;
    mEvent->getByToken(input_pcalohits_hcal_token_, simHits);
    mHcalSimHitCollection = &*simHits;
  }
  return mHcalSimHitCollection;
}

getHFRecHitCollection()

const HFRecHitCollection * JetMatchingTools::getHFRecHitCollection

(

)

Definition at line 97 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_hfrechits_token_, mEvent, mHFRecHitCollection, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by getConstituentHits().

                                                                  {
  if (!mHFRecHitCollection) {
    edm::Handle<HFRecHitCollection> recHits;
    mEvent->getByToken(input_hfrechits_token_, recHits);
    mHFRecHitCollection = &*recHits;
  }
  return mHFRecHitCollection;
}

getHORecHitCollection()

const HORecHitCollection * JetMatchingTools::getHORecHitCollection

(

)

Definition at line 89 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_horechits_token_, mEvent, mHORecHitCollection, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by getConstituentHits().

                                                                  {
  if (!mHORecHitCollection) {
    edm::Handle<HORecHitCollection> recHits;
    mEvent->getByToken(input_horechits_token_, recHits);
    mHORecHitCollection = &*recHits;
  }
  return mHORecHitCollection;
}

getPCaloHits() [1/2]

std::vector< const PCaloHit * > JetMatchingTools::getPCaloHits

(

DetId

fId

)

get PCaloHits contributing to the detId

Definition at line 216 of file JetMatchingTools.cc.

References DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, getEBSimHitCollection(), getEESimHitCollection(), getHcalSimHitCollection(), DetId::Hcal, mps_fire::result, and DetId::subdetId().

Referenced by getGenParticles(), and lostEnergyFraction().

                                                                   {
  std::vector<const PCaloHit*> result;
  if (fId.det() == DetId::Ecal) {
    if ((EcalSubdetector)fId.subdetId() == EcalBarrel) {
      result = getSimHits(*getEBSimHitCollection(), fId);
    } else if ((EcalSubdetector)fId.subdetId() == EcalEndcap) {
      result = getSimHits(*getEESimHitCollection(), fId);
    }
  } else if (fId.det() == DetId::Hcal) {
    result = getSimHits(*getHcalSimHitCollection(), fId);
  }
  return result;
}

getPCaloHits() [2/2]

std::vector<const PCaloHit*> JetMatchingTools::getPCaloHits

(

int

fGeneratorId

)

CaloSimHits.

getSimTrackCollection()

const SimTrackContainer * JetMatchingTools::getSimTrackCollection

(

)

Definition at line 129 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_simtrack_token_, mEvent, mSimTrackCollection, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by getTrack().

                                                                 {
  if (!mSimTrackCollection) {
    edm::Handle<SimTrackContainer> simHits;
    mEvent->getByToken(input_simtrack_token_, simHits);
    mSimTrackCollection = &*simHits;
  }
  return mSimTrackCollection;
}

getSimVertexCollection()

const SimVertexContainer * JetMatchingTools::getSimVertexCollection

(

)

Definition at line 137 of file JetMatchingTools.cc.

References edm::Event::getByToken(), input_simvertex_token_, mEvent, mSimVertexCollection, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by generatorId().

                                                                   {
  if (!mSimVertexCollection) {
    edm::Handle<SimVertexContainer> simHits;
    mEvent->getByToken(input_simvertex_token_, simHits);
    mSimVertexCollection = &*simHits;
  }
  return mSimVertexCollection;
}

getTrack()

const SimTrack * JetMatchingTools::getTrack

(

unsigned

fSimTrackId

)

convert trackId to SimTrack

Definition at line 232 of file JetMatchingTools.cc.

References getSimTrackCollection(), and mps_fire::i.

Referenced by generatorId().

                                                               {
  for (unsigned i = 0; i < getSimTrackCollection()->size(); ++i) {
    if ((*getSimTrackCollection())[i].trackId() == fSimTrackId)
      return &(*getSimTrackCollection())[i];
  }
  return nullptr;
}

getTrackId()

int JetMatchingTools::getTrackId

(

const PCaloHit &

fHit

)

GEANT track ID.

Definition at line 230 of file JetMatchingTools.cc.

References PCaloHit::geantTrackId().

Referenced by getGenParticles(), and lostEnergyFraction().

{ return fHit.geantTrackId(); }

lostEnergyFraction()

double JetMatchingTools::lostEnergyFraction

(

const reco::CaloJet &

fJet

)

energy in broken links

Definition at line 311 of file JetMatchingTools.cc.

References generatorId(), getConstituentHits(), getConstituents(), getPCaloHits(), getTrackId(), phits, FastTrackerRecHitMaskProducer_cfi::recHits, and HLT_2022v12_cff::towers.

                                                                   {
  double totalEnergy = 0;
  double lostEnergy = 0;
  // follow the chain
  std::vector<const CaloTower*> towers = getConstituents(fJet);
  for (unsigned itower = 0; itower < towers.size(); ++itower) {
    std::vector<JetConstituent> recHits = getConstituentHits(*(towers[itower]));
    for (unsigned ihit = 0; ihit < recHits.size(); ++ihit) {
      double foundSimEnergy = 0;
      double lostSimEnergy = 0;
      std::vector<const PCaloHit*> phits = getPCaloHits(recHits[ihit].id);
      for (unsigned iphit = 0; iphit < phits.size(); ++iphit) {
        double simEnergy = phits[iphit]->energy();
        int trackId = getTrackId(*(phits[iphit]));
        if (trackId < 0 || generatorId(trackId) < 0)
          lostSimEnergy += simEnergy;
        else
          foundSimEnergy += simEnergy;
      }
      if (foundSimEnergy > 0 || lostSimEnergy > 0) {
        totalEnergy += recHits[ihit].energy;
        lostEnergy += recHits[ihit].energy * lostSimEnergy / (foundSimEnergy + lostSimEnergy);
      }
    }
  }
  return lostEnergy / totalEnergy;
}

overlapEnergyFraction()

double JetMatchingTools::overlapEnergyFraction	(	const std::vector< const reco::GenParticle *> &	fObject,
		const std::vector< const reco::GenParticle *> &	fReference
	)		const

energy overlap

Definition at line 340 of file JetMatchingTools.cc.

References HCALHighEnergyHPDFilter_cfi::energy, spr::find(), and mps_fire::i.

                                                                                                          {
  if (fObject.empty())
    return 0;
  double totalEnergy = 0;
  double overlapEnergy = 0;
  for (unsigned i = 0; i < fObject.size(); ++i) {
    totalEnergy += fObject[i]->energy();
    if (find(fReference.begin(), fReference.end(), fObject[i]) != fReference.end())
      overlapEnergy += fObject[i]->energy();
  }
  return overlapEnergy / totalEnergy;
}

Member Data Documentation

input_cands_token_

edm::EDGetTokenT<reco::CandidateCollection> JetMatchingTools::input_cands_token_

private

Definition at line 114 of file JetMatchingTools.h.

Referenced by getGenParticlesCollection(), and JetMatchingTools().

input_ebrechits_token_

edm::EDGetTokenT<EBRecHitCollection> JetMatchingTools::input_ebrechits_token_

private

Definition at line 104 of file JetMatchingTools.h.

Referenced by getEBRecHitCollection(), and JetMatchingTools().

input_eerechits_token_

edm::EDGetTokenT<EERecHitCollection> JetMatchingTools::input_eerechits_token_

private

Definition at line 105 of file JetMatchingTools.h.

Referenced by getEERecHitCollection(), and JetMatchingTools().

input_hbherechits_token_

edm::EDGetTokenT<HBHERecHitCollection> JetMatchingTools::input_hbherechits_token_

private

Definition at line 106 of file JetMatchingTools.h.

Referenced by getHBHERecHitCollection(), and JetMatchingTools().

input_hfrechits_token_

edm::EDGetTokenT<HFRecHitCollection> JetMatchingTools::input_hfrechits_token_

private

Definition at line 108 of file JetMatchingTools.h.

Referenced by getHFRecHitCollection(), and JetMatchingTools().

input_horechits_token_

edm::EDGetTokenT<HORecHitCollection> JetMatchingTools::input_horechits_token_

private

Definition at line 107 of file JetMatchingTools.h.

Referenced by getHORecHitCollection(), and JetMatchingTools().

input_pcalohits_ebcal_token_

edm::EDGetTokenT<edm::PCaloHitContainer> JetMatchingTools::input_pcalohits_ebcal_token_

private

Definition at line 110 of file JetMatchingTools.h.

Referenced by getEBSimHitCollection(), and JetMatchingTools().

input_pcalohits_eecal_token_

edm::EDGetTokenT<edm::PCaloHitContainer> JetMatchingTools::input_pcalohits_eecal_token_

private

Definition at line 109 of file JetMatchingTools.h.

Referenced by getEESimHitCollection(), and JetMatchingTools().

input_pcalohits_hcal_token_

edm::EDGetTokenT<edm::PCaloHitContainer> JetMatchingTools::input_pcalohits_hcal_token_

private

Definition at line 111 of file JetMatchingTools.h.

Referenced by getHcalSimHitCollection(), and JetMatchingTools().

input_simtrack_token_

edm::EDGetTokenT<edm::SimTrackContainer> JetMatchingTools::input_simtrack_token_

private

Definition at line 112 of file JetMatchingTools.h.

Referenced by getSimTrackCollection(), and JetMatchingTools().

input_simvertex_token_

edm::EDGetTokenT<edm::SimVertexContainer> JetMatchingTools::input_simvertex_token_

private

Definition at line 113 of file JetMatchingTools.h.

Referenced by getSimVertexCollection(), and JetMatchingTools().

mEBRecHitCollection

const EBRecHitCollection* JetMatchingTools::mEBRecHitCollection

private

Definition at line 92 of file JetMatchingTools.h.

Referenced by getEBRecHitCollection().

mEBSimHitCollection

const edm::PCaloHitContainer* JetMatchingTools::mEBSimHitCollection

private

Definition at line 97 of file JetMatchingTools.h.

Referenced by getEBSimHitCollection().

mEERecHitCollection

const EERecHitCollection* JetMatchingTools::mEERecHitCollection

private

Definition at line 93 of file JetMatchingTools.h.

Referenced by getEERecHitCollection().

mEESimHitCollection

const edm::PCaloHitContainer* JetMatchingTools::mEESimHitCollection

private

Definition at line 98 of file JetMatchingTools.h.

Referenced by getEESimHitCollection().

mEvent

const edm::Event* JetMatchingTools::mEvent

private

Definition at line 91 of file JetMatchingTools.h.

Referenced by getEBRecHitCollection(), getEBSimHitCollection(), getEERecHitCollection(), getEESimHitCollection(), getGenParticlesCollection(), getHBHERecHitCollection(), getHcalSimHitCollection(), getHFRecHitCollection(), getHORecHitCollection(), getSimTrackCollection(), and getSimVertexCollection().

mGenParticleCollection

const reco::CandidateCollection* JetMatchingTools::mGenParticleCollection

private

Definition at line 102 of file JetMatchingTools.h.

Referenced by getGenParticlesCollection().

mHBHERecHitCollection

const HBHERecHitCollection* JetMatchingTools::mHBHERecHitCollection

private

Definition at line 94 of file JetMatchingTools.h.

Referenced by getHBHERecHitCollection().

mHcalSimHitCollection

const edm::PCaloHitContainer* JetMatchingTools::mHcalSimHitCollection

private

Definition at line 99 of file JetMatchingTools.h.

Referenced by getHcalSimHitCollection().

mHFRecHitCollection

const HFRecHitCollection* JetMatchingTools::mHFRecHitCollection

private

Definition at line 96 of file JetMatchingTools.h.

Referenced by getHFRecHitCollection().

mHORecHitCollection

const HORecHitCollection* JetMatchingTools::mHORecHitCollection

private

Definition at line 95 of file JetMatchingTools.h.

Referenced by getHORecHitCollection().

mSimTrackCollection

const edm::SimTrackContainer* JetMatchingTools::mSimTrackCollection

private

Definition at line 100 of file JetMatchingTools.h.

Referenced by getSimTrackCollection().

mSimVertexCollection

const edm::SimVertexContainer* JetMatchingTools::mSimVertexCollection

private

Definition at line 101 of file JetMatchingTools.h.

Referenced by getSimVertexCollection().