|
|
|
Make Jets from protoobjects. More...
Functions | |
| HcalSubdetector | hcalSubdetector (int fEta) |
| converts eta to the corresponding HCAL subdetector. | |
| bool | makeSpecific (const JetReco::InputCollection &fConstituents, reco::PFJet::Specific *fJetSpecific) |
| Make PFlowJet specifics. Assumes ProtoJet is made from ParticleFlowCandidates. | |
| bool | makeSpecific (const JetReco::InputCollection &fConstituents, reco::GenJet::Specific *fJetSpecific) |
| Make GenJet specifics. Assumes ProtoJet is made from HepMCCandidate. | |
| bool | makeSpecific (const JetReco::InputCollection &fConstituents, const CaloSubdetectorGeometry &fTowerGeometry, reco::CaloJet::Specific *fJetSpecific) |
| Make CaloJet specifics. Assumes ProtoJet is made from CaloTowerCandidates. | |
Make Jets from protoobjects.
| HcalSubdetector JetMaker::hcalSubdetector | ( | int | fEta | ) |
converts eta to the corresponding HCAL subdetector.
Referenced by makeSpecific().
| bool JetMaker::makeSpecific | ( | const JetReco::InputCollection & | fConstituents, |
| reco::PFJet::Specific * | fJetSpecific | ||
| ) |
Make PFlowJet specifics. Assumes ProtoJet is made from ParticleFlowCandidates.
@@ PFJET *************************
Definition at line 127 of file JetMaker.cc.
References benchmark_cfg::cerr, reco::LeafCandidate::energy(), reco::Candidate::get(), h, reco::PFJet::Specific::mChargedEmEnergy, reco::PFJet::Specific::mChargedHadronEnergy, reco::PFJet::Specific::mChargedMuEnergy, reco::PFJet::Specific::mChargedMultiplicity, reco::PFJet::Specific::mMuonMultiplicity, reco::PFJet::Specific::mNeutralEmEnergy, reco::PFJet::Specific::mNeutralHadronEnergy, reco::PFJet::Specific::mNeutralMultiplicity, and reco::PFCandidate::particleId().
{
if (!fJetSpecific) return false;
// 1.- Loop over PFCandidates,
// 2.- Get the corresponding PFCandidate
// 3.- Calculate the different PFJet specific quantities
float chargedHadronEnergy=0.;
float neutralHadronEnergy=0.;
float chargedEmEnergy=0.;
float neutralEmEnergy=0.;
float chargedMuEnergy=0.;
int chargedMultiplicity=0;
int neutralMultiplicity=0;
int muonMultiplicity=0;
JetReco::InputCollection::const_iterator constituent = fPFCandidates.begin();
for (; constituent != fPFCandidates.end(); ++constituent) {
const Candidate* candidate = constituent->get ();
if (candidate) {
const PFCandidate* pfCand = dynamic_cast<const PFCandidate*> (candidate);
if (pfCand) {
switch ( PFCandidate::ParticleType (pfCand->particleId())) {
case PFCandidate::h: // charged hadron
chargedHadronEnergy += pfCand->energy();
chargedMultiplicity++;
break;
case PFCandidate::e: // electron
chargedEmEnergy += pfCand->energy();
chargedMultiplicity++;
break;
case PFCandidate::mu: // muon
chargedMuEnergy += pfCand->energy();
chargedMultiplicity++;
muonMultiplicity++;
break;
case PFCandidate::gamma: // photon
case PFCandidate::egamma_HF : // electromagnetic in HF
neutralEmEnergy += pfCand->energy();
neutralMultiplicity++;
break;
case PFCandidate::h0 : // neutral hadron
case PFCandidate::h_HF : // hadron in HF
neutralHadronEnergy += pfCand->energy();
neutralMultiplicity++;
break;
default:
std::cerr << "JetMaker::makeSpecific (PFJetJet)-> Unknown PFCandidate::ParticleType: " << pfCand->particleId() << " is ignored" << std::endl;
break;
}
}
else {
std::cerr << "JetMaker::makeSpecific (PFJetJet)-> Referred constituent is not PFCandidate" << std::endl;
}
}
else {
std::cerr << "JetMaker::makeSpecific (PFJetJet)-> Referred constituent is not available in the event" << std::endl;
}
}
fJetSpecific->mChargedHadronEnergy=chargedHadronEnergy;
fJetSpecific->mNeutralHadronEnergy= neutralHadronEnergy;
fJetSpecific->mChargedEmEnergy=chargedEmEnergy;
fJetSpecific->mChargedMuEnergy=chargedMuEnergy;
fJetSpecific->mNeutralEmEnergy=neutralEmEnergy;
fJetSpecific->mChargedMultiplicity=chargedMultiplicity;
fJetSpecific->mNeutralMultiplicity=neutralMultiplicity;
fJetSpecific->mMuonMultiplicity=muonMultiplicity;
return true;
}
| bool JetMaker::makeSpecific | ( | const JetReco::InputCollection & | fConstituents, |
| reco::GenJet::Specific * | fJetSpecific | ||
| ) |
Make GenJet specifics. Assumes ProtoJet is made from HepMCCandidate.
Definition at line 204 of file JetMaker.cc.
References abs, benchmark_cfg::cerr, reco::LeafCandidate::energy(), edm::RefToBase< T >::get(), reco::Candidate::get(), reco::Candidate::hasMasterClone(), reco::GenJet::Specific::m_AuxiliaryEnergy, reco::GenJet::Specific::m_EmEnergy, reco::GenJet::Specific::m_HadEnergy, reco::GenJet::Specific::m_InvisibleEnergy, reco::Candidate::masterClone(), and reco::LeafCandidate::pdgId().
{
for (JetReco::InputCollection::const_iterator genCand = fMcParticles.begin(); genCand != fMcParticles.end(); ++genCand) {
const Candidate* candidate = genCand->get ();
if (candidate->hasMasterClone ()) candidate = candidate->masterClone().get ();
if (candidate) {
const GenParticle* genParticle = GenJet::genParticle (candidate);
if (genParticle) {
double e = genParticle->energy();
switch (std::abs (genParticle->pdgId ())) {
case 22: // photon
case 11: // e
fJetSpecific->m_EmEnergy += e;
break;
case 211: // pi
case 321: // K
case 130: // KL
case 2212: // p
case 2112: // n
fJetSpecific->m_HadEnergy += e;
break;
case 13: // muon
case 12: // nu_e
case 14: // nu_mu
case 16: // nu_tau
fJetSpecific->m_InvisibleEnergy += e;
break;
default:
fJetSpecific->m_AuxiliaryEnergy += e;
}
}
else {
std::cerr << "JetMaker::makeSpecific (GenJet)-> Referred GenParticleCandidate is not available in the event" << std::endl;
}
}
else {
std::cerr << "JetMaker::makeSpecific (GenJet)-> Referred constituent is not available in the event" << std::endl;
}
}
return true;
}
| bool JetMaker::makeSpecific | ( | const JetReco::InputCollection & | fConstituents, |
| const CaloSubdetectorGeometry & | fTowerGeometry, | ||
| reco::CaloJet::Specific * | fJetSpecific | ||
| ) |
Make CaloJet specifics. Assumes ProtoJet is made from CaloTowerCandidates.
Definition at line 21 of file JetMaker.cc.
References benchmark_cfg::cerr, dPhi(), CaloTower::emEnergy(), reco::LeafCandidate::energy(), geometry, reco::Candidate::get(), CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloTower::hadEnergy(), HcalBarrel, HcalEndcap, HcalForward, hcalSubdetector(), CaloTower::id(), CaloTowerDetId::ieta(), reco::CaloJet::Specific::mEmEnergyInEB, reco::CaloJet::Specific::mEmEnergyInEE, reco::CaloJet::Specific::mEmEnergyInHF, reco::CaloJet::Specific::mEnergyFractionEm, reco::CaloJet::Specific::mEnergyFractionHadronic, reco::CaloJet::Specific::mHadEnergyInHB, reco::CaloJet::Specific::mHadEnergyInHE, reco::CaloJet::Specific::mHadEnergyInHF, reco::CaloJet::Specific::mHadEnergyInHO, reco::CaloJet::Specific::mMaxEInEmTowers, reco::CaloJet::Specific::mMaxEInHadTowers, reco::CaloJet::Specific::mTowersArea, CaloTower::outerEnergy(), and python::multivaluedict::sort().
{
if (!fJetSpecific) return false;
// 1.- Loop over the tower Ids,
// 2.- Get the corresponding CaloTower
// 3.- Calculate the different CaloJet specific quantities
vector<double> eECal_i;
vector<double> eHCal_i;
double eInHad = 0.;
double eInEm = 0.;
double eInHO = 0.;
double eInHB = 0.;
double eInHE = 0.;
double eHadInHF = 0.;
double eEmInHF = 0.;
double eInEB = 0.;
double eInEE = 0.;
double jetArea = 0.;
for (JetReco::InputCollection::const_iterator towerCand = fTowers.begin(); towerCand != fTowers.end(); ++towerCand) {
const Candidate* candidate = towerCand->get ();
if (candidate) {
const CaloTower* tower = dynamic_cast<const CaloTower*> (candidate);
if (tower) {
//Array of energy in EM Towers:
eECal_i.push_back(tower->emEnergy());
eInEm += tower->emEnergy();
//Array of energy in HCAL Towers:
eHCal_i.push_back(tower->hadEnergy());
eInHad += tower->hadEnergy();
// figure out contributions
switch (JetMaker::hcalSubdetector (tower->id().ieta())) {
case HcalBarrel:
eInHB += tower->hadEnergy();
eInHO += tower->outerEnergy();
eInEB += tower->emEnergy();
break;
case HcalEndcap:
eInHE += tower->hadEnergy();
eInEE += tower->emEnergy();
break;
case HcalForward:
eHadInHF += tower->hadEnergy();
eEmInHF += tower->emEnergy();
break;
default:
break;
}
// get area of the tower (++ minus --)
if ( tower->energy() > 0 ) {
const CaloCellGeometry* geometry = fTowerGeometry.getGeometry(tower->id());
if (geometry) {
float dEta = fabs (geometry->getCorners() [0].eta() - geometry->getCorners() [2].eta());
float dPhi = fabs (geometry->getCorners() [0].phi() - geometry->getCorners() [2].phi());
jetArea += dEta * dPhi;
}
}
else {
std::cerr << "JetMaker::makeSpecific (CaloJet)-> Geometry for cell " << tower->id() << " can not be found. Ignoring cell" << std::endl;
}
}
else {
std::cerr << "JetMaker::makeSpecific (CaloJet)-> Constituent is not of CaloTower type" << std::endl;
}
}
else {
std::cerr << "JetMaker::makeSpecific (CaloJet)-> Referred constituent is not available in the event" << std::endl;
}
}
double towerEnergy = eInHad + eInEm;
fJetSpecific->mHadEnergyInHO = eInHO;
fJetSpecific->mHadEnergyInHB = eInHB;
fJetSpecific->mHadEnergyInHE = eInHE;
fJetSpecific->mHadEnergyInHF = eHadInHF;
fJetSpecific->mEmEnergyInHF = eEmInHF;
fJetSpecific->mEmEnergyInEB = eInEB;
fJetSpecific->mEmEnergyInEE = eInEE;
if (towerEnergy > 0) {
fJetSpecific->mEnergyFractionHadronic = eInHad / towerEnergy;
fJetSpecific->mEnergyFractionEm = eInEm / towerEnergy;
}
else { // HO only jet
fJetSpecific->mEnergyFractionHadronic = 1.;
fJetSpecific->mEnergyFractionEm = 0.;
}
fJetSpecific->mTowersArea = jetArea;
fJetSpecific->mMaxEInEmTowers = 0;
fJetSpecific->mMaxEInHadTowers = 0;
//Sort the arrays
sort(eECal_i.begin(), eECal_i.end(), greater<double>());
sort(eHCal_i.begin(), eHCal_i.end(), greater<double>());
if (!fTowers.empty ()) {
//Highest value in the array is the first element of the array
fJetSpecific->mMaxEInEmTowers = eECal_i.front();
fJetSpecific->mMaxEInHadTowers = eHCal_i.front();
}
return true;
}
1.7.3