CaloSpecificAlgo Class Reference

#include <RecoMET/METAlgorithms/interface/CaloSpecificAlgo.h>

Public Member Functions
reco::CaloMET	addInfo (edm::Handle< edm::View< reco::Candidate > > towers, const CommonMETData &met, bool noHF, double globalThreshold)

Private Types
typedef math::XYZTLorentzVector	LorentzVector
typedef math::XYZPoint	Point

Private Member Functions
void	add_MET_in_HF (SpecificCaloMETData &specific, double sumEtInpHF, double MExInpHF, double MEyInpHF, double sumEtInmHF, double MExInmHF, double MEyInmHF)
DetId	find_DetId_of_ECAL_cell_in_constituent_of (const CaloTower *calotower)
DetId	find_DetId_of_HCAL_cell_in_constituent_of (const CaloTower *calotower)
void	remove_HF_from_MET (CommonMETData &met, double sumEtInpHF, double MExInpHF, double MEyInpHF, double sumEtInmHF, double MExInmHF, double MEyInmHF)
void	update_EmEtInEB_EmEtInEE (float &EmEtInEB, float &EmEtInEE, const CaloTower *calotower)
void	update_HadEtInHB_HadEtInHE_HadEtInHO_HadEtInHF_EmEtInHF (float &HadEtInHB, float &HadEtInHE, float &HadEtInHO, float &HadEtInHF, float &EmEtInHF, const CaloTower *calotower, bool noHF)
void	update_MaxTowerEm_MaxTowerHad (float &MaxTowerEm, float &MaxTowerHad, const CaloTower *calotower, bool noHF)
void	update_sumEtInpHF_MExInpHF_MEyInpHF_sumEtInmHF_MExInmHF_MEyInmHF (double &sumEtInpHF, double &MExInpHF, double &MEyInpHF, double &sumEtInmHF, double &MExInmHF, double &MEyInmHF, const CaloTower *calotower)
void	update_totalEt_totalEm (double &totalEt, double &totalEm, const CaloTower *calotower, bool noHF)

Detailed Description

Description: Adds Calorimeter specific information to MET base class

Implementation: [Notes on implementation]

Definition at line 37 of file CaloSpecificAlgo.h.

Member Typedef Documentation

LorentzVector

typedef math::XYZTLorentzVector CaloSpecificAlgo::LorentzVector

private

Definition at line 45 of file CaloSpecificAlgo.h.

Point

typedef math::XYZPoint CaloSpecificAlgo::Point

private

Definition at line 46 of file CaloSpecificAlgo.h.

Member Function Documentation

add_MET_in_HF()

void CaloSpecificAlgo::add_MET_in_HF ( SpecificCaloMETData &  specific, double  sumEtInpHF, double  MExInpHF, double  MEyInpHF, double  sumEtInmHF, double  MExInmHF, double  MEyInmHF  )

private

Definition at line 202 of file CaloSpecificAlgo.cc.

                                                      {
  LorentzVector METpHF(MExInpHF, MEyInpHF, 0, sqrt(MExInpHF * MExInpHF + MEyInpHF * MEyInpHF));
  LorentzVector METmHF(MExInmHF, MEyInmHF, 0, sqrt(MExInmHF * MExInmHF + MEyInmHF * MEyInmHF));
  specific.CaloMETInpHF = METpHF.pt();
  specific.CaloMETInmHF = METmHF.pt();
  specific.CaloMETPhiInpHF = METpHF.Phi();
  specific.CaloMETPhiInmHF = METmHF.Phi();
  specific.CaloSETInpHF = sumEtInpHF;
  specific.CaloSETInmHF = sumEtInmHF;
}

References timingPdfMaker::specific, and mathSSE::sqrt().

addInfo()

reco::CaloMET CaloSpecificAlgo::addInfo	(	edm::Handle< edm::View< reco::Candidate > >	towers,
		const CommonMETData &	met,
		bool	noHF,
		double	globalThreshold
	)

Definition at line 26 of file CaloSpecificAlgo.cc.

                                                                {
  SpecificCaloMETData specific;
  CommonMETData met = umet;
  double totalEt = 0.0;
  double totalEm = 0.0;
 
  double sumEtInpHF = 0.0;
  double sumEtInmHF = 0.0;
  double MExInpHF = 0.0;
  double MEyInpHF = 0.0;
  double MExInmHF = 0.0;
  double MEyInmHF = 0.0;
 
  for (edm::View<Candidate>::const_iterator towerCand = towers->begin(); towerCand != towers->end(); ++towerCand) {
    const CaloTower *calotower = dynamic_cast<const CaloTower *>(&(*towerCand));
    if (!calotower)
      continue;
    if (calotower->et() < globalThreshold)
      continue;
    update_totalEt_totalEm(totalEt, totalEm, calotower, noHF);
    update_MaxTowerEm_MaxTowerHad(specific.MaxEtInEmTowers, specific.MaxEtInHadTowers, calotower, noHF);
    update_EmEtInEB_EmEtInEE(specific.EmEtInEB, specific.EmEtInEE, calotower);
    update_HadEtInHB_HadEtInHE_HadEtInHO_HadEtInHF_EmEtInHF(specific.HadEtInHB,
                                                            specific.HadEtInHE,
                                                            specific.HadEtInHO,
                                                            specific.HadEtInHF,
                                                            specific.EmEtInHF,
                                                            calotower,
                                                            noHF);
    update_sumEtInpHF_MExInpHF_MEyInpHF_sumEtInmHF_MExInmHF_MEyInmHF(
        sumEtInpHF, MExInpHF, MEyInpHF, sumEtInmHF, MExInmHF, MEyInmHF, calotower);
  }
 
  double totalHad = totalEt - totalEm;
 
  if (noHF)
    remove_HF_from_MET(met, sumEtInpHF, MExInpHF, MEyInpHF, sumEtInmHF, MExInmHF, MEyInmHF);
 
  if (!noHF)
    add_MET_in_HF(specific, sumEtInpHF, MExInpHF, MEyInpHF, sumEtInmHF, MExInmHF, MEyInmHF);
 
  specific.EtFractionHadronic = (totalEt == 0.0) ? 0.0 : totalHad / totalEt;
  specific.EtFractionEm = (totalEt == 0.0) ? 0.0 : totalEm / totalEt;
 
  const LorentzVector p4(met.mex, met.mey, 0.0, met.met);
  const Point vtx(0.0, 0.0, 0.0);
  CaloMET caloMET(specific, met.sumet, p4, vtx);
  return caloMET;
}

References slimmedMETs_cfi::caloMET, ewkMuLumiMonitorDQM_cfi::calotower, HLT_FULL_cff::globalThreshold, BTaggingMonitor_cfi::met, HLT_FULL_cff::noHF, p4, timingPdfMaker::specific, HLT_FULL_cff::towers, and extraflags_cff::vtx.

Referenced by cms::CaloMETProducer::produce().

find_DetId_of_ECAL_cell_in_constituent_of()

DetId CaloSpecificAlgo::find_DetId_of_ECAL_cell_in_constituent_of ( const CaloTower *  calotower )

private

Definition at line 233 of file CaloSpecificAlgo.cc.

                                                                                            {
  DetId ret;
  for (int cell = calotower->constituentsSize() - 1; cell >= 0; --cell) {
    DetId id = calotower->constituent(cell);
    if (id.det() == DetId::Ecal) {
      ret = id;
      break;
    }
  }
  return ret;
}

References ewkMuLumiMonitorDQM_cfi::calotower, DetId::Ecal, triggerObjects_cff::id, and runTheMatrix::ret.

find_DetId_of_HCAL_cell_in_constituent_of()

DetId CaloSpecificAlgo::find_DetId_of_HCAL_cell_in_constituent_of ( const CaloTower *  calotower )

private

Definition at line 220 of file CaloSpecificAlgo.cc.

                                                                                            {
  DetId ret;
  for (int cell = calotower->constituentsSize() - 1; cell >= 0; --cell) {
    DetId id = calotower->constituent(cell);
    if (id.det() == DetId::Hcal) {
      ret = id;
      break;
    }
  }
  return ret;
}

References ewkMuLumiMonitorDQM_cfi::calotower, DetId::Hcal, triggerObjects_cff::id, and runTheMatrix::ret.

remove_HF_from_MET()

void CaloSpecificAlgo::remove_HF_from_MET ( CommonMETData &  met, double  sumEtInpHF, double  MExInpHF, double  MEyInpHF, double  sumEtInmHF, double  MExInmHF, double  MEyInmHF  )

private

Definition at line 188 of file CaloSpecificAlgo.cc.

                                                           {
  met.mex -= (MExInmHF + MExInpHF);
  met.mey -= (MEyInmHF + MEyInpHF);
  met.sumet -= (sumEtInpHF + sumEtInmHF);
  met.met = sqrt(met.mex * met.mex + met.mey * met.mey);
}

References BTaggingMonitor_cfi::met, and mathSSE::sqrt().

update_EmEtInEB_EmEtInEE()

void CaloSpecificAlgo::update_EmEtInEB_EmEtInEE ( float &  EmEtInEB, float &  EmEtInEE, const CaloTower *  calotower  )

private

Definition at line 119 of file CaloSpecificAlgo.cc.

                                                                                                            {
  DetId detIdEcal = find_DetId_of_ECAL_cell_in_constituent_of(calotower);
  if (detIdEcal.null())
    return;
 
  EcalSubdetector subdet = EcalSubdetector(detIdEcal.subdetId());
  if (subdet == EcalBarrel) {
    EmEtInEB += calotower->emEt();
  } else if (subdet == EcalEndcap) {
    EmEtInEE += calotower->emEt();
  }
}

References ewkMuLumiMonitorDQM_cfi::calotower, EcalBarrel, EcalEndcap, DetId::null(), and DetId::subdetId().

update_HadEtInHB_HadEtInHE_HadEtInHO_HadEtInHF_EmEtInHF()

void CaloSpecificAlgo::update_HadEtInHB_HadEtInHE_HadEtInHO_HadEtInHF_EmEtInHF ( float &  HadEtInHB, float &  HadEtInHE, float &  HadEtInHO, float &  HadEtInHF, float &  EmEtInHF, const CaloTower *  calotower, bool  noHF  )

private

Definition at line 133 of file CaloSpecificAlgo.cc.

                                                                                          {
  DetId detIdHcal = find_DetId_of_HCAL_cell_in_constituent_of(calotower);
  if (detIdHcal.null())
    return;
 
  HcalSubdetector subdet = HcalDetId(detIdHcal).subdet();
  if (subdet == HcalBarrel || subdet == HcalOuter) {
    HadEtInHB += calotower->hadEt();
    HadEtInHO += calotower->outerEt();
  }
 
  if (subdet == HcalEndcap) {
    HadEtInHE += calotower->hadEt();
  }
 
  if (subdet == HcalForward && !noHF) {
    HadEtInHF += calotower->hadEt();
    EmEtInHF += calotower->emEt();
  }
}

References ewkMuLumiMonitorDQM_cfi::calotower, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HLT_FULL_cff::noHF, DetId::null(), and HcalDetId::subdet().

update_MaxTowerEm_MaxTowerHad()

void CaloSpecificAlgo::update_MaxTowerEm_MaxTowerHad ( float &  MaxTowerEm, float &  MaxTowerHad, const CaloTower *  calotower, bool  noHF  )

private

Definition at line 95 of file CaloSpecificAlgo.cc.

                                                                {
  DetId detIdHcal = find_DetId_of_HCAL_cell_in_constituent_of(calotower);
  DetId detIdEcal = find_DetId_of_ECAL_cell_in_constituent_of(calotower);
 
  if (!detIdHcal.null()) {
    HcalSubdetector subdet = HcalDetId(detIdHcal).subdet();
    if (subdet == HcalBarrel || subdet == HcalOuter || subdet == HcalEndcap || (!noHF && subdet == HcalForward)) {
      if (calotower->hadEt() > MaxTowerHad)
        MaxTowerHad = calotower->hadEt();
      if (calotower->emEt() > MaxTowerEm)
        MaxTowerEm = calotower->emEt();
    }
  }
 
  if (!detIdEcal.null()) {
    if (calotower->emEt() > MaxTowerEm)
      MaxTowerEm = calotower->emEt();
  }
}

References ewkMuLumiMonitorDQM_cfi::calotower, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HLT_FULL_cff::noHF, DetId::null(), and HcalDetId::subdet().

update_sumEtInpHF_MExInpHF_MEyInpHF_sumEtInmHF_MExInmHF_MEyInmHF()

void CaloSpecificAlgo::update_sumEtInpHF_MExInpHF_MEyInpHF_sumEtInmHF_MExInmHF_MEyInmHF ( double &  sumEtInpHF, double &  MExInpHF, double &  MEyInpHF, double &  sumEtInmHF, double &  MExInmHF, double &  MEyInmHF, const CaloTower *  calotower  )

private

Definition at line 161 of file CaloSpecificAlgo.cc.

                                                                                                                    {
  DetId detIdHcal = find_DetId_of_HCAL_cell_in_constituent_of(calotower);
  if (detIdHcal.null())
    return;
 
  HcalSubdetector subdet = HcalDetId(detIdHcal).subdet();
  if (!(subdet == HcalForward))
    return;
 
  if (calotower->eta() >= 0) {
    sumEtInpHF += calotower->et();
    MExInpHF -= (calotower->et() * cos(calotower->phi()));
    MEyInpHF -= (calotower->et() * sin(calotower->phi()));
  } else {
    sumEtInmHF += calotower->et();
    MExInmHF -= (calotower->et() * cos(calotower->phi()));
    MEyInmHF -= (calotower->et() * sin(calotower->phi()));
  }
}

References ewkMuLumiMonitorDQM_cfi::calotower, funct::cos(), HcalForward, DetId::null(), funct::sin(), and HcalDetId::subdet().

update_totalEt_totalEm()

void CaloSpecificAlgo::update_totalEt_totalEm ( double &  totalEt, double &  totalEm, const CaloTower *  calotower, bool  noHF  )

private

Definition at line 80 of file CaloSpecificAlgo.cc.

                                                                                                                     {
  if (noHF) {
    DetId detIdHcal = find_DetId_of_HCAL_cell_in_constituent_of(calotower);
    if (!detIdHcal.null()) {
      HcalSubdetector subdet = HcalDetId(detIdHcal).subdet();
      if (subdet == HcalForward)
        return;
    }
  }
 
  totalEt += calotower->et();
  totalEm += calotower->emEt();
}

References ewkMuLumiMonitorDQM_cfi::calotower, HcalForward, HLT_FULL_cff::noHF, DetId::null(), and HcalDetId::subdet().