AnaL1CaloCleaner.cc File Reference

#include <memory>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "TTree.h"
#include "TFile.h"
#include "Math/VectorUtil.h"

Go to the source code of this file.

Classes
class	AnaL1CaloCleaner

Functions
std::pair< double, double >	compMEtProjU (const reco::Candidate::LorentzVector &zP4, double metPx, double metPy, int &errorFlag)
	DEFINE_FWK_MODULE (AnaL1CaloCleaner)
std::string	getKey (const DetId &det)

Function Documentation

std::pair<double, double> compMEtProjU	(	const reco::Candidate::LorentzVector &	zP4,
		double	metPx,
		double	metPy,
		int &	errorFlag
	)

Definition at line 107 of file AnaL1CaloCleaner.cc.

Referenced by AnaL1CaloCleaner::analyze(), and MCEmbeddingValidationAnalyzer::plotEntryTypeL1ETM::fillHistograms().

{
  if ( zP4.pt() == 0. ) {
    edm::LogWarning ("compMEtProjU")
      << " Failed to compute projection, because Z0 candidate has zero Pt --> returning dummy solution !!";
    errorFlag = 1;
    return std::pair<double, double>(0., 0.);
  }
  
  double qX = zP4.px();
  double qY = zP4.py();
  double qT = TMath::Sqrt(qX*qX + qY*qY);
  
  double uX = -metPx;//(qX + metPx);
  double uY = -metPy;//(qY + metPy);
  
  double u1 = (uX*qX + uY*qY)/qT;
  double u2 = (uX*qY - uY*qX)/qT;
  
  return std::pair<double, double>(u1, u2);
}

DEFINE_FWK_MODULE

(

AnaL1CaloCleaner

)

std::string getKey

(

const DetId &

det

)

Definition at line 82 of file AnaL1CaloCleaner.cc.

References DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, EcalLaserPnDiode, EcalPreshower, EcalTriggerTower, DetId::Hcal, HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalOther, HcalOuter, HcalTriggerTower, and DetId::subdetId().

Referenced by DTCalibrationMap::getConsts().

{
  std::map<int, std::string > detMap_;
  std::map<int, std::map<int, std::string> > subDetMap_;

  detMap_[DetId::Hcal]="Hcal";
  detMap_[DetId::Ecal]="Ecal";
  
  subDetMap_[DetId::Ecal][EcalBarrel]="EcalBarrel";
  subDetMap_[DetId::Ecal][EcalEndcap]="EcalEndcap";
  subDetMap_[DetId::Ecal][EcalPreshower ]="EcalPreshower";
  subDetMap_[DetId::Ecal][EcalTriggerTower]="EcalTriggerTower";
  subDetMap_[DetId::Ecal][EcalLaserPnDiode]="EcalLaserPnDiode";
  
  subDetMap_[DetId::Hcal][HcalEmpty]="HcalEmpty";
  subDetMap_[DetId::Hcal][HcalBarrel]="HcalBarrel";
  subDetMap_[DetId::Hcal][HcalEndcap]="HcalEndcap";
  subDetMap_[DetId::Hcal][HcalOuter]="HcalOuter";
  subDetMap_[DetId::Hcal][HcalForward]="HcalForward";
  subDetMap_[DetId::Hcal][HcalTriggerTower]="HcalTriggerTower";
  subDetMap_[DetId::Hcal][HcalOther]="HcalOther";

  return "H_"+detMap_[det.det()]+"_"+subDetMap_[det.det()][det.subdetId()];
}