HCALRecHitAnalyzer Class Reference

#include <HCALRecHitAnalyzer.h>

Inheritance diagram for HCALRecHitAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endJob ()
	HCALRecHitAnalyzer (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	BookHistos ()
void	FillGeometry (const edm::EventSetup &)

Private Attributes
DQMStore *	dbe_
bool	debug_
bool	finebinning_
std::string	FolderName_
edm::InputTag	hBHERecHitsLabel_
edm::InputTag	hFRecHitsLabel_
MonitorElement *	hHCAL_D1_energy_ieta_iphi
MonitorElement *	hHCAL_D1_energyvsieta
MonitorElement *	hHCAL_D1_Maxenergy_ieta_iphi
MonitorElement *	hHCAL_D1_Maxenergyvsieta
MonitorElement *	hHCAL_D1_METPhivsieta
MonitorElement *	hHCAL_D1_METvsieta
MonitorElement *	hHCAL_D1_MExvsieta
MonitorElement *	hHCAL_D1_MEyvsieta
MonitorElement *	hHCAL_D1_Minenergy_ieta_iphi
MonitorElement *	hHCAL_D1_Minenergyvsieta
MonitorElement *	hHCAL_D1_Occ_ieta_iphi
MonitorElement *	hHCAL_D1_Occvsieta
MonitorElement *	hHCAL_D1_SETvsieta
MonitorElement *	hHCAL_D2_energy_ieta_iphi
MonitorElement *	hHCAL_D2_energyvsieta
MonitorElement *	hHCAL_D2_Maxenergy_ieta_iphi
MonitorElement *	hHCAL_D2_Maxenergyvsieta
MonitorElement *	hHCAL_D2_METPhivsieta
MonitorElement *	hHCAL_D2_METvsieta
MonitorElement *	hHCAL_D2_MExvsieta
MonitorElement *	hHCAL_D2_MEyvsieta
MonitorElement *	hHCAL_D2_Minenergy_ieta_iphi
MonitorElement *	hHCAL_D2_Minenergyvsieta
MonitorElement *	hHCAL_D2_Occ_ieta_iphi
MonitorElement *	hHCAL_D2_Occvsieta
MonitorElement *	hHCAL_D2_SETvsieta
MonitorElement *	hHCAL_D3_energy_ieta_iphi
MonitorElement *	hHCAL_D3_energyvsieta
MonitorElement *	hHCAL_D3_Maxenergy_ieta_iphi
MonitorElement *	hHCAL_D3_Maxenergyvsieta
MonitorElement *	hHCAL_D3_METPhivsieta
MonitorElement *	hHCAL_D3_METvsieta
MonitorElement *	hHCAL_D3_MExvsieta
MonitorElement *	hHCAL_D3_MEyvsieta
MonitorElement *	hHCAL_D3_Minenergy_ieta_iphi
MonitorElement *	hHCAL_D3_Minenergyvsieta
MonitorElement *	hHCAL_D3_Occ_ieta_iphi
MonitorElement *	hHCAL_D3_Occvsieta
MonitorElement *	hHCAL_D3_SETvsieta
MonitorElement *	hHCAL_D4_energy_ieta_iphi
MonitorElement *	hHCAL_D4_energyvsieta
MonitorElement *	hHCAL_D4_Maxenergy_ieta_iphi
MonitorElement *	hHCAL_D4_Maxenergyvsieta
MonitorElement *	hHCAL_D4_METPhivsieta
MonitorElement *	hHCAL_D4_METvsieta
MonitorElement *	hHCAL_D4_MExvsieta
MonitorElement *	hHCAL_D4_MEyvsieta
MonitorElement *	hHCAL_D4_Minenergy_ieta_iphi
MonitorElement *	hHCAL_D4_Minenergyvsieta
MonitorElement *	hHCAL_D4_Occ_ieta_iphi
MonitorElement *	hHCAL_D4_Occvsieta
MonitorElement *	hHCAL_D4_SETvsieta
MonitorElement *	hHCAL_ieta_detaMap
MonitorElement *	hHCAL_ieta_dphiMap
MonitorElement *	hHCAL_ieta_iphi_etaMap
MonitorElement *	hHCAL_ieta_iphi_HBMap
MonitorElement *	hHCAL_ieta_iphi_HEMap
MonitorElement *	hHCAL_ieta_iphi_HFMap
MonitorElement *	hHCAL_ieta_iphi_HOMap
MonitorElement *	hHCAL_ieta_iphi_phiMap
MonitorElement *	hHCAL_Nevents
edm::InputTag	hORecHitsLabel_
int	Nevents

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 21 of file HCALRecHitAnalyzer.h.

Constructor & Destructor Documentation

HCALRecHitAnalyzer::HCALRecHitAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 50 of file HCALRecHitAnalyzer.cc.

References debug_, finebinning_, FolderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hBHERecHitsLabel_, hFRecHitsLabel_, and hORecHitsLabel_.

{

  // Retrieve Information from the Configuration File
  hBHERecHitsLabel_  = iConfig.getParameter<edm::InputTag>("HBHERecHitsLabel");
  hORecHitsLabel_    = iConfig.getParameter<edm::InputTag>("HORecHitsLabel");
  hFRecHitsLabel_    = iConfig.getParameter<edm::InputTag>("HFRecHitsLabel");

  debug_             = iConfig.getParameter<bool>("Debug");
  finebinning_ = iConfig.getUntrackedParameter<bool>("FineBinning");
  FolderName_ = iConfig.getUntrackedParameter<std::string>("FolderName");

}

Member Function Documentation

void HCALRecHitAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Int_t depth = det.depth(); //always 4

Implements edm::EDAnalyzer.

Definition at line 466 of file HCALRecHitAnalyzer.cc.

References edm::SortedCollection< T, SORT >::begin(), DEBUG, HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), ET, eta(), create_public_lumi_plots::exp, MonitorElement::Fill(), MonitorElement::getBinContent(), edm::Event::getByLabel(), egHLT::errCodes::HBHERecHits, hBHERecHitsLabel_, HcalBarrel, egHLT::errCodes::HFRecHits, hFRecHitsLabel_, hHCAL_D1_energy_ieta_iphi, hHCAL_D1_energyvsieta, hHCAL_D1_Maxenergy_ieta_iphi, hHCAL_D1_Maxenergyvsieta, hHCAL_D1_METPhivsieta, hHCAL_D1_METvsieta, hHCAL_D1_MExvsieta, hHCAL_D1_MEyvsieta, hHCAL_D1_Minenergy_ieta_iphi, hHCAL_D1_Minenergyvsieta, hHCAL_D1_Occ_ieta_iphi, hHCAL_D1_Occvsieta, hHCAL_D1_SETvsieta, hHCAL_D2_energy_ieta_iphi, hHCAL_D2_energyvsieta, hHCAL_D2_Maxenergy_ieta_iphi, hHCAL_D2_Maxenergyvsieta, hHCAL_D2_METPhivsieta, hHCAL_D2_METvsieta, hHCAL_D2_MExvsieta, hHCAL_D2_MEyvsieta, hHCAL_D2_Minenergy_ieta_iphi, hHCAL_D2_Minenergyvsieta, hHCAL_D2_Occ_ieta_iphi, hHCAL_D2_Occvsieta, hHCAL_D2_SETvsieta, hHCAL_D3_energy_ieta_iphi, hHCAL_D3_energyvsieta, hHCAL_D3_Maxenergy_ieta_iphi, hHCAL_D3_Maxenergyvsieta, hHCAL_D3_METPhivsieta, hHCAL_D3_METvsieta, hHCAL_D3_MExvsieta, hHCAL_D3_MEyvsieta, hHCAL_D3_Minenergy_ieta_iphi, hHCAL_D3_Minenergyvsieta, hHCAL_D3_Occ_ieta_iphi, hHCAL_D3_Occvsieta, hHCAL_D3_SETvsieta, hHCAL_D4_energy_ieta_iphi, hHCAL_D4_energyvsieta, hHCAL_D4_Maxenergy_ieta_iphi, hHCAL_D4_Maxenergyvsieta, hHCAL_D4_METPhivsieta, hHCAL_D4_METvsieta, hHCAL_D4_MExvsieta, hHCAL_D4_MEyvsieta, hHCAL_D4_Minenergy_ieta_iphi, hHCAL_D4_Minenergyvsieta, hHCAL_D4_Occ_ieta_iphi, hHCAL_D4_Occvsieta, hHCAL_D4_SETvsieta, hHCAL_ieta_iphi_etaMap, hHCAL_ieta_iphi_phiMap, hHCAL_Nevents, hORecHitsLabel_, i, HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), j, Nevents, phi, colinearityKinematic::Phi, edm::Handle< T >::product(), reco::tau::disc::Pt(), MonitorElement::setBinContent(), HcalDetId::subdet(), and theta().

{
    Nevents++;
    hHCAL_Nevents->Fill(0);
  // ==========================================================
  // Retrieve!
  // ==========================================================
                                                                                                                                                             
  const HBHERecHitCollection *HBHERecHits;
  const HORecHitCollection *HORecHits;
  const HFRecHitCollection *HFRecHits;

  edm::Handle<HBHERecHitCollection> HBHERecHitsHandle;
  iEvent.getByLabel(hBHERecHitsLabel_,HBHERecHitsHandle);
  if (!HBHERecHitsHandle.isValid()) {
    edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
                   << "HCALRecHitAnalyzer::analyze!\n"; return;
  } else {
    HBHERecHits = HBHERecHitsHandle.product();
  }
  edm::Handle<HORecHitCollection> HORecHitsHandle;
  iEvent.getByLabel(hORecHitsLabel_,HORecHitsHandle);
  if (!HORecHitsHandle.isValid()) {
    edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
                   << "HCALRecHitAnalyzer::analyze!\n"; return;
  } else {
    HORecHits = HORecHitsHandle.product();
  }
  edm::Handle<HFRecHitCollection> HFRecHitsHandle;
  iEvent.getByLabel(hFRecHitsLabel_,HFRecHitsHandle);
  if (!HFRecHitsHandle.isValid()) {
    edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
                   << "HCALRecHitAnalyzer::analyze!\n"; return;
  } else {
    HFRecHits = HFRecHitsHandle.product();
  }

  // ==========================================================
  // Fill Histograms!
  // ==========================================================
                                                   
  TLorentzVector vHBHEMET_EtaRing[83][4];
  int HBHEActiveRing[83][4];
  int HBHENActiveCells[83][4];
  double HBHESET_EtaRing[83][4];
  double HBHEMinEnergy_EtaRing[83][4];
  double HBHEMaxEnergy_EtaRing[83][4];

  for (int i = 0; i < 83; i++) {
    for (int j = 0; j < 4; j++) {

      HBHEActiveRing[i][j] = 0;
      HBHENActiveCells[i][j] = 0;
      HBHESET_EtaRing[i][j] = 0; 
      HBHEMinEnergy_EtaRing[i][j] = 14E3; 
      HBHEMaxEnergy_EtaRing[i][j] = -999; 
    }
  }
  
  // Loop over HBHERecHit's
  HBHERecHitCollection::const_iterator hbherechit;
  int nHBrechit = 0, nHErechit = 0;
  
  for (hbherechit = HBHERecHits->begin(); hbherechit != HBHERecHits->end(); hbherechit++) {
    
    HcalDetId det = hbherechit->id();
    double Energy = hbherechit->energy();
    Int_t depth = det.depth();
    Int_t ieta = det.ieta();
    Int_t iphi = det.iphi();
    int EtaRing = 41 + ieta; // this counts from 0    
    double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
    double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
    double theta = 2*TMath::ATan(exp(-1*eta));
    double ET = Energy*TMath::Sin(theta);
    HcalSubdetector HcalNum = det.subdet();
    TLorentzVector v_;
    
    
    
    if (Energy>0) // zero suppress
      {
    HBHEActiveRing[EtaRing][depth-1] = 1;
    HBHENActiveCells[EtaRing][depth-1]++;
    HBHESET_EtaRing[EtaRing][depth-1]+=ET;
    v_.SetPtEtaPhiE(ET, 0, phi, ET);
    vHBHEMET_EtaRing[EtaRing][depth-1]-=v_;
    
    
    DEBUG( EtaRing << " " << Energy << " " << ET << " " << theta << " " << eta << " " << phi << " : " << vHBHEMET_EtaRing[EtaRing][depth-1].Phi() << " " << vHBHEMET_EtaRing[EtaRing][depth-1].Pt() );  
    
    switch (depth) {
    case 1:
      hHCAL_D1_Occ_ieta_iphi->Fill(ieta,iphi);
      break;
    case 2:
      hHCAL_D2_Occ_ieta_iphi->Fill(ieta,iphi);
      break;
    case 3:
      hHCAL_D3_Occ_ieta_iphi->Fill(ieta,iphi);
      break;
    } // end switch
      }
    
    if (Energy > HBHEMaxEnergy_EtaRing[EtaRing][depth-1])
      HBHEMaxEnergy_EtaRing[EtaRing][depth-1] = Energy;
    if (Energy < HBHEMinEnergy_EtaRing[EtaRing][depth-1])
      HBHEMinEnergy_EtaRing[EtaRing][depth-1] = Energy;
    
    switch (depth) {
    case 1:
      hHCAL_D1_energy_ieta_iphi->Fill(ieta,iphi,Energy);           
      hHCAL_D1_energyvsieta->Fill(ieta,Energy);    
      if (Energy>hHCAL_D1_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D1_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D1_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      
      break;
    case 2:
      hHCAL_D2_energy_ieta_iphi->Fill(ieta,iphi,Energy);
      
      hHCAL_D2_energyvsieta->Fill(ieta,Energy);
      if (Energy>hHCAL_D2_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D2_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D2_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      break;
    case 3:
      hHCAL_D3_energy_ieta_iphi->Fill(ieta,iphi,Energy);
      
      hHCAL_D3_energyvsieta->Fill(ieta,Energy);
      if (Energy>hHCAL_D3_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D3_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D3_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D3_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      break;
    } // end switch
    
    
    if (HcalNum == HcalBarrel) {
      nHBrechit++;
    } else { // HcalEndcap
      nHErechit++;
    }
  } // end loop over HBHERecHit's
  
  // Fill eta-ring MET quantities
  for (int iEtaRing=0; iEtaRing < 83; iEtaRing++) {
    for (int jDepth=0; jDepth < 3; jDepth++) {
      
      switch (jDepth+1) {
      case 1:
    hHCAL_D1_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
    hHCAL_D1_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
    break;
      case 2:
        hHCAL_D2_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
        hHCAL_D2_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
        break;
      case 3:
        hHCAL_D3_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
        hHCAL_D3_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
        break;
      }
      
      if (HBHEActiveRing[iEtaRing][jDepth]) {
    
    switch (jDepth+1) {
    case 1:
      hHCAL_D1_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
      hHCAL_D1_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
      hHCAL_D1_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
      hHCAL_D1_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
      hHCAL_D1_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
      hHCAL_D1_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
      break;
    case 2:
      hHCAL_D2_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
      hHCAL_D2_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
      hHCAL_D2_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
      hHCAL_D2_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
      hHCAL_D2_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
      hHCAL_D2_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
      break;
    case 3:
      hHCAL_D3_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
      hHCAL_D3_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
      hHCAL_D3_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
      hHCAL_D3_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
      hHCAL_D3_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
      hHCAL_D3_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
      break;
    } //switch
      } //activering
    }//depth
  } //etaring    

  
  //-------------------------------------------------HO
  TLorentzVector vHOMET_EtaRing[83];
  int HOActiveRing[83];
  int HONActiveCells[83];
  double HOSET_EtaRing[83];
  double HOMinEnergy_EtaRing[83];
  double HOMaxEnergy_EtaRing[83];
  
  for (int i=0;i<83; i++) {

    HOActiveRing[i] = 0;
    HONActiveCells[i] = 0;
    HOSET_EtaRing[i] = 0; 
    HOMinEnergy_EtaRing[i] = 14E3;
    HOMaxEnergy_EtaRing[i] = -999;
  }

  // Loop over HORecHit's
  HORecHitCollection::const_iterator horechit;
  int nHOrechit = 0;

  for (horechit = HORecHits->begin(); horechit != HORecHits->end(); horechit++) {
                                                                                                                                                             
    nHOrechit++;
                                                                                                                                                             
    HcalDetId det = horechit->id();
    double Energy = horechit->energy();
    Int_t ieta = det.ieta();
    Int_t iphi = det.iphi();
    int EtaRing = 41+ieta; // this counts from 0    
    double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
    double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
    double theta = 2*TMath::ATan(exp(-1*eta));
    double ET = Energy*TMath::Sin(theta);
    TLorentzVector v_; 


    if (Energy>0) // zero suppress
      {
    HOActiveRing[EtaRing] = 1;
    HONActiveCells[EtaRing]++;
    HOSET_EtaRing[EtaRing]+=ET;
    v_.SetPtEtaPhiE(ET, 0, phi, ET);
    vHOMET_EtaRing[EtaRing]-=v_;
      
    hHCAL_D4_Occ_ieta_iphi->Fill(ieta,iphi);
      }

    if (Energy > HOMaxEnergy_EtaRing[EtaRing])
      HOMaxEnergy_EtaRing[EtaRing] = Energy;
    if (Energy < HOMinEnergy_EtaRing[EtaRing])
      HOMinEnergy_EtaRing[EtaRing] = Energy;
    
    hHCAL_D4_energy_ieta_iphi->Fill(ieta,iphi,Energy);
    
    hHCAL_D4_energyvsieta->Fill(ieta,Energy);     
    if (Energy>hHCAL_D4_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D4_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D4_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D4_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
                                                                                                                                                   
  } // end loop over HORecHit's

  // Fill eta-ring MET quantities
  for (int iEtaRing=0; iEtaRing<83; iEtaRing++) {
    hHCAL_D4_Maxenergyvsieta->Fill(iEtaRing-41, HOMaxEnergy_EtaRing[iEtaRing]);
    hHCAL_D4_Minenergyvsieta->Fill(iEtaRing-41, HOMinEnergy_EtaRing[iEtaRing]);

    if (HOActiveRing[iEtaRing]) {

      hHCAL_D4_METPhivsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Phi());
      hHCAL_D4_MExvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Px());
      hHCAL_D4_MEyvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Py());
      hHCAL_D4_METvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Pt());
      hHCAL_D4_SETvsieta->Fill(iEtaRing-41, HOSET_EtaRing[iEtaRing]);
      hHCAL_D4_Occvsieta->Fill(iEtaRing-41, HONActiveCells[iEtaRing]);

    }

  }

  //------------------------------------------------HF

  TLorentzVector vHFMET_EtaRing[83][2];
  int HFActiveRing[83][2];
  int HFNActiveCells[83][2];
  double HFSET_EtaRing[83][2];
  double HFMinEnergy_EtaRing[83][2];
  double HFMaxEnergy_EtaRing[83][2];

  for (int i=0;i<83; i++) {
    for (int j=0;j<2; j++) {

      HFActiveRing[i][j] = 0;
      HFNActiveCells[i][j] = 0;
      HFSET_EtaRing[i][j] = 0; 
      HFMinEnergy_EtaRing[i][j] = 14E3; 
      HFMaxEnergy_EtaRing[i][j] = -999; 
      
    }
  }

  // Loop over HFRecHit's
  HFRecHitCollection::const_iterator hfrechit;
  int nHFrechit = 0;

  for (hfrechit = HFRecHits->begin(); hfrechit != HFRecHits->end(); hfrechit++) {
                                                                                                                                                             
    nHFrechit++;
                                                                                                                                                             
    HcalDetId det = hfrechit->id();
    double Energy = hfrechit->energy();
    Int_t depth = det.depth();
    Int_t ieta = det.ieta();
    Int_t iphi = det.iphi();
    int EtaRing = 41+ieta; // this counts from 0    
    double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
    double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
    double theta = 2*TMath::ATan(exp(-1*eta));
    double ET = Energy*TMath::Sin(theta);

    TLorentzVector v_;
    if (Energy>0) // zero suppress
      {
    HFActiveRing[EtaRing][depth-1] = 1;
    HFNActiveCells[EtaRing][depth-1]++;
    HFSET_EtaRing[EtaRing][depth-1]+=ET;
    v_.SetPtEtaPhiE(ET, 0, phi, ET);
    vHFMET_EtaRing[EtaRing][depth-1]-=v_;                                                                         
    
    switch (depth) {
    case 1:
      hHCAL_D1_Occ_ieta_iphi->Fill(ieta,iphi);
      break;
    case 2:
      hHCAL_D2_Occ_ieta_iphi->Fill(ieta,iphi);
      break;
      
    }
      }
    
    if (Energy > HFMaxEnergy_EtaRing[EtaRing][depth-1])
      HFMaxEnergy_EtaRing[EtaRing][depth-1] = Energy;
    if (Energy < HFMinEnergy_EtaRing[EtaRing][depth-1])
      HFMinEnergy_EtaRing[EtaRing][depth-1] = Energy;
    

    switch (depth) {
    case 1:
      hHCAL_D1_energy_ieta_iphi->Fill(ieta,iphi,Energy);
      
      hHCAL_D1_energyvsieta->Fill(ieta,Energy);  
      if (Energy>hHCAL_D1_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D1_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D1_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy); 
      break;
    case 2:
      hHCAL_D2_energy_ieta_iphi->Fill(ieta,iphi,Energy);
      
      hHCAL_D2_energyvsieta->Fill(ieta,Energy);  
      if (Energy>hHCAL_D2_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      if (Energy<hHCAL_D2_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi)) 
    hHCAL_D2_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
      break;
    }

  } // end loop over HFRecHit's

  // Fill eta-ring MET quantities
  for (int iEtaRing=0; iEtaRing<83; iEtaRing++) {
    for (int jDepth=0; jDepth<2; jDepth++) {
      switch (jDepth+1) {
      case 1:
        hHCAL_D1_Maxenergyvsieta->Fill(iEtaRing-41, HFMaxEnergy_EtaRing[iEtaRing][jDepth]);
        hHCAL_D1_Minenergyvsieta->Fill(iEtaRing-41, HFMinEnergy_EtaRing[iEtaRing][jDepth]);
        break;
      case 2:
        hHCAL_D2_Maxenergyvsieta->Fill(iEtaRing-41, HFMaxEnergy_EtaRing[iEtaRing][jDepth]);
        hHCAL_D2_Minenergyvsieta->Fill(iEtaRing-41, HFMinEnergy_EtaRing[iEtaRing][jDepth]);
        break;
      }


      if (HFActiveRing[iEtaRing][jDepth]) {

    switch (jDepth+1) {
    case 1:

      hHCAL_D1_METPhivsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Phi());
      hHCAL_D1_MExvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Px());
      hHCAL_D1_MEyvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Py());
      hHCAL_D1_METvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Pt());
      hHCAL_D1_SETvsieta->Fill(iEtaRing-41, HFSET_EtaRing[iEtaRing][jDepth]);
      hHCAL_D1_Occvsieta->Fill(iEtaRing-41, HFNActiveCells[iEtaRing][jDepth]);
      break;

    case 2:

      hHCAL_D2_METPhivsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Phi());
      hHCAL_D2_MExvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Px());
      hHCAL_D2_MEyvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Py());
      hHCAL_D2_METvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Pt());
      hHCAL_D2_SETvsieta->Fill(iEtaRing-41, HFSET_EtaRing[iEtaRing][jDepth]);
      hHCAL_D2_Occvsieta->Fill(iEtaRing-41, HFNActiveCells[iEtaRing][jDepth]);
      break;

    }
      }
    }
  }



}

void HCALRecHitAnalyzer::beginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 64 of file HCALRecHitAnalyzer.cc.

References BookHistos(), FillGeometry(), and Nevents.

{
  Nevents = 0;
  BookHistos();
  FillGeometry(iSetup);

}

void HCALRecHitAnalyzer::BookHistos ( )

privatevirtual

Definition at line 72 of file HCALRecHitAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe_, finebinning_, FolderName_, hHCAL_D1_energy_ieta_iphi, hHCAL_D1_energyvsieta, hHCAL_D1_Maxenergy_ieta_iphi, hHCAL_D1_Maxenergyvsieta, hHCAL_D1_METPhivsieta, hHCAL_D1_METvsieta, hHCAL_D1_MExvsieta, hHCAL_D1_MEyvsieta, hHCAL_D1_Minenergy_ieta_iphi, hHCAL_D1_Minenergyvsieta, hHCAL_D1_Occ_ieta_iphi, hHCAL_D1_Occvsieta, hHCAL_D1_SETvsieta, hHCAL_D2_energy_ieta_iphi, hHCAL_D2_energyvsieta, hHCAL_D2_Maxenergy_ieta_iphi, hHCAL_D2_Maxenergyvsieta, hHCAL_D2_METPhivsieta, hHCAL_D2_METvsieta, hHCAL_D2_MExvsieta, hHCAL_D2_MEyvsieta, hHCAL_D2_Minenergy_ieta_iphi, hHCAL_D2_Minenergyvsieta, hHCAL_D2_Occ_ieta_iphi, hHCAL_D2_Occvsieta, hHCAL_D2_SETvsieta, hHCAL_D3_energy_ieta_iphi, hHCAL_D3_energyvsieta, hHCAL_D3_Maxenergy_ieta_iphi, hHCAL_D3_Maxenergyvsieta, hHCAL_D3_METPhivsieta, hHCAL_D3_METvsieta, hHCAL_D3_MExvsieta, hHCAL_D3_MEyvsieta, hHCAL_D3_Minenergy_ieta_iphi, hHCAL_D3_Minenergyvsieta, hHCAL_D3_Occ_ieta_iphi, hHCAL_D3_Occvsieta, hHCAL_D3_SETvsieta, hHCAL_D4_energy_ieta_iphi, hHCAL_D4_energyvsieta, hHCAL_D4_Maxenergy_ieta_iphi, hHCAL_D4_Maxenergyvsieta, hHCAL_D4_METPhivsieta, hHCAL_D4_METvsieta, hHCAL_D4_MExvsieta, hHCAL_D4_MEyvsieta, hHCAL_D4_Minenergy_ieta_iphi, hHCAL_D4_Minenergyvsieta, hHCAL_D4_Occ_ieta_iphi, hHCAL_D4_Occvsieta, hHCAL_D4_SETvsieta, hHCAL_ieta_detaMap, hHCAL_ieta_dphiMap, hHCAL_ieta_iphi_etaMap, hHCAL_ieta_iphi_HBMap, hHCAL_ieta_iphi_HEMap, hHCAL_ieta_iphi_HFMap, hHCAL_ieta_iphi_HOMap, hHCAL_ieta_iphi_phiMap, hHCAL_Nevents, i, j, cppFunctionSkipper::operator, MonitorElement::setBinContent(), and DQMStore::setCurrentFolder().

Referenced by beginRun().

                                   {
  // get ahold of back-end interface
  dbe_ = edm::Service<DQMStore>().operator->();


  if (dbe_) {

    dbe_->setCurrentFolder(FolderName_+"/geometry");
    hHCAL_ieta_iphi_HBMap = dbe_->book2D("METTask_HCAL_ieta_iphi_HBMap","",83,-41,42,72,1,73); 
    hHCAL_ieta_iphi_HEMap = dbe_->book2D("METTask_HCAL_ieta_iphi_HEMap","",83,-41,42,72,1,73); 
    hHCAL_ieta_iphi_HFMap = dbe_->book2D("METTask_HCAL_ieta_iphi_HFMap","",83,-41,42,72,1,73); 
    hHCAL_ieta_iphi_HOMap = dbe_->book2D("METTask_HCAL_ieta_iphi_HOMap","",83,-41,42,72,1,73); 
    hHCAL_ieta_iphi_etaMap = dbe_->book2D("METTask_HCAL_ieta_iphi_etaMap","",83,-41,42,72,1,73);
    hHCAL_ieta_iphi_phiMap = dbe_->book2D("METTask_HCAL_ieta_iphi_phiMap","",83,-41,42,72,1,73);
    hHCAL_ieta_detaMap = dbe_->book1D("METTask_HCAL_ieta_detaMap","",83,-41,42);
    hHCAL_ieta_dphiMap = dbe_->book1D("METTask_HCAL_ieta_dphiMap","",83,-41,42);  

    // Initialize bins for geometry to -999 because z = 0 is a valid entry 
    for (int i = 1; i <= 83; i++) {

      hHCAL_ieta_detaMap->setBinContent(i,-999);
      hHCAL_ieta_dphiMap->setBinContent(i,-999);

      for (int j = 1; j <= 72; j++) {

    hHCAL_ieta_iphi_HBMap->setBinContent(i,j,0);
    hHCAL_ieta_iphi_HEMap->setBinContent(i,j,0);
    hHCAL_ieta_iphi_HFMap->setBinContent(i,j,0);
    hHCAL_ieta_iphi_HOMap->setBinContent(i,j,0);
        hHCAL_ieta_iphi_etaMap->setBinContent(i,j,-999);
        hHCAL_ieta_iphi_phiMap->setBinContent(i,j,-999);

      }

    }
  
    //    dbe_->setCurrentFolder("RecoMETV/MET_HCAL/data");
    dbe_->setCurrentFolder(FolderName_);
    //--Store number of events used
    hHCAL_Nevents          = dbe_->book1D("METTask_HCAL_Nevents","",1,0,1);  
    //--Data integrated over all events and stored by HCAL(ieta,iphi) 

    hHCAL_D1_energy_ieta_iphi = dbe_->book2D("METTask_HCAL_D1_energy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D2_energy_ieta_iphi = dbe_->book2D("METTask_HCAL_D2_energy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D3_energy_ieta_iphi = dbe_->book2D("METTask_HCAL_D3_energy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D4_energy_ieta_iphi = dbe_->book2D("METTask_HCAL_D4_energy_ieta_iphi","",83,-41,42,72,1,73);  
    
    hHCAL_D1_Minenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D1_Minenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D2_Minenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D2_Minenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D3_Minenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D3_Minenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D4_Minenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D4_Minenergy_ieta_iphi","",83,-41,42,72,1,73); 
    
    hHCAL_D1_Maxenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D1_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D2_Maxenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D2_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D3_Maxenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D3_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D4_Maxenergy_ieta_iphi = dbe_->book2D("METTask_HCAL_D4_Maxenergy_ieta_iphi","",83,-41,42,72,1,73); 
    
    // need to initialize those
    for (int i=1; i<=83; i++)
      for (int j=1; j<=73; j++)
    {
      hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
      hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
      hHCAL_D3_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
      hHCAL_D4_Maxenergy_ieta_iphi->setBinContent(i,j,-999);

      hHCAL_D1_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
      hHCAL_D2_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
      hHCAL_D3_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
      hHCAL_D4_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
    }

    hHCAL_D1_Occ_ieta_iphi = dbe_->book2D("METTask_HCAL_D1_Occ_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D2_Occ_ieta_iphi = dbe_->book2D("METTask_HCAL_D2_Occ_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D3_Occ_ieta_iphi = dbe_->book2D("METTask_HCAL_D3_Occ_ieta_iphi","",83,-41,42,72,1,73);  
    hHCAL_D4_Occ_ieta_iphi = dbe_->book2D("METTask_HCAL_D4_Occ_ieta_iphi","",83,-41,42,72,1,73);  
    //--Data over eta-rings
    
    
    // CaloTower values

    if(finebinning_)
      {
    hHCAL_D1_energyvsieta = dbe_->book2D("METTask_HCAL_D1_energyvsieta","",83,-41,42,20101,-100,2001);  
    hHCAL_D2_energyvsieta = dbe_->book2D("METTask_HCAL_D2_energyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D3_energyvsieta = dbe_->book2D("METTask_HCAL_D3_energyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D4_energyvsieta = dbe_->book2D("METTask_HCAL_D4_energyvsieta","",83,-41,42,20101,-100,2001);
    
    hHCAL_D1_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D1_Minenergyvsieta","",83,-41,42,20101,-100,2001);  
    hHCAL_D2_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D2_Minenergyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D3_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D3_Minenergyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D4_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D4_Minenergyvsieta","",83,-41,42,20101,-100,2001);
    
    hHCAL_D1_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D1_Maxenergyvsieta","",83,-41,42,20101,-100,2001);  
    hHCAL_D2_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D2_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D3_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D3_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
    hHCAL_D4_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D4_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
    
    // Integrated over phi
    hHCAL_D1_Occvsieta = dbe_->book2D("METTask_HCAL_D1_Occvsieta","",83,-41,42,73,0,73);
    hHCAL_D2_Occvsieta = dbe_->book2D("METTask_HCAL_D2_Occvsieta","",83,-41,42,73,0,73);
    hHCAL_D3_Occvsieta = dbe_->book2D("METTask_HCAL_D3_Occvsieta","",83,-41,42,73,0,73);
    hHCAL_D4_Occvsieta = dbe_->book2D("METTask_HCAL_D4_Occvsieta","",83,-41,42,73,0,73);
    
    hHCAL_D1_SETvsieta = dbe_->book2D("METTask_HCAL_D1_SETvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D2_SETvsieta = dbe_->book2D("METTask_HCAL_D2_SETvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D3_SETvsieta = dbe_->book2D("METTask_HCAL_D3_SETvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D4_SETvsieta = dbe_->book2D("METTask_HCAL_D4_SETvsieta","",83,-41,42,20001,0,2001); 
    
    hHCAL_D1_METvsieta = dbe_->book2D("METTask_HCAL_D1_METvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D2_METvsieta = dbe_->book2D("METTask_HCAL_D2_METvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D3_METvsieta = dbe_->book2D("METTask_HCAL_D3_METvsieta","",83,-41,42,20001,0,2001); 
    hHCAL_D4_METvsieta = dbe_->book2D("METTask_HCAL_D4_METvsieta","",83,-41,42,20001,0,2001); 
    
    hHCAL_D1_METPhivsieta      = dbe_->book2D("METTask_HCAL_D1_METPhivsieta","",83,-41,42, 80, -4,4);  
    hHCAL_D2_METPhivsieta      = dbe_->book2D("METTask_HCAL_D2_METPhivsieta","",83,-41,42, 80, -4,4);  
    hHCAL_D3_METPhivsieta      = dbe_->book2D("METTask_HCAL_D3_METPhivsieta","",83,-41,42, 80, -4,4);  
    hHCAL_D4_METPhivsieta      = dbe_->book2D("METTask_HCAL_D4_METPhivsieta","",83,-41,42, 80, -4,4);  
    
    hHCAL_D1_MExvsieta         = dbe_->book2D("METTask_HCAL_D1_MExvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D2_MExvsieta         = dbe_->book2D("METTask_HCAL_D2_MExvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D3_MExvsieta         = dbe_->book2D("METTask_HCAL_D3_MExvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D4_MExvsieta         = dbe_->book2D("METTask_HCAL_D4_MExvsieta","",83,-41,42, 10001,-500,501);  
    
    hHCAL_D1_MEyvsieta         = dbe_->book2D("METTask_HCAL_D1_MEyvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D2_MEyvsieta         = dbe_->book2D("METTask_HCAL_D2_MEyvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D3_MEyvsieta         = dbe_->book2D("METTask_HCAL_D3_MEyvsieta","",83,-41,42, 10001,-500,501);  
    hHCAL_D4_MEyvsieta         = dbe_->book2D("METTask_HCAL_D4_MEyvsieta","",83,-41,42, 10001,-500,501);  
      }
    else
      {
    hHCAL_D1_energyvsieta = dbe_->book2D("METTask_HCAL_D1_energyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D2_energyvsieta = dbe_->book2D("METTask_HCAL_D2_energyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D3_energyvsieta = dbe_->book2D("METTask_HCAL_D3_energyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D4_energyvsieta = dbe_->book2D("METTask_HCAL_D4_energyvsieta","",83,-41,42,1000,-10,1990);

        hHCAL_D1_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D1_Minenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D2_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D2_Minenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D3_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D3_Minenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D4_Minenergyvsieta = dbe_->book2D("METTask_HCAL_D4_Minenergyvsieta","",83,-41,42,1000,-10,1990);

        hHCAL_D1_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D1_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D2_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D2_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D3_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D3_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
        hHCAL_D4_Maxenergyvsieta = dbe_->book2D("METTask_HCAL_D4_Maxenergyvsieta","",83,-41,42,1000,-10,1990);

        // Integrated over phi                                                                                                                                                                                 
        hHCAL_D1_Occvsieta = dbe_->book2D("METTask_HCAL_D1_Occvsieta","",83,-41,42,73,0,73);
        hHCAL_D2_Occvsieta = dbe_->book2D("METTask_HCAL_D2_Occvsieta","",83,-41,42,73,0,73);
        hHCAL_D3_Occvsieta = dbe_->book2D("METTask_HCAL_D3_Occvsieta","",83,-41,42,73,0,73);
        hHCAL_D4_Occvsieta = dbe_->book2D("METTask_HCAL_D4_Occvsieta","",83,-41,42,73,0,73);

        hHCAL_D1_SETvsieta = dbe_->book2D("METTask_HCAL_D1_SETvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D2_SETvsieta = dbe_->book2D("METTask_HCAL_D2_SETvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D3_SETvsieta = dbe_->book2D("METTask_HCAL_D3_SETvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D4_SETvsieta = dbe_->book2D("METTask_HCAL_D4_SETvsieta","",83,-41,42,1000,0,2000);

        hHCAL_D1_METvsieta = dbe_->book2D("METTask_HCAL_D1_METvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D2_METvsieta = dbe_->book2D("METTask_HCAL_D2_METvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D3_METvsieta = dbe_->book2D("METTask_HCAL_D3_METvsieta","",83,-41,42,1000,0,2000);
        hHCAL_D4_METvsieta = dbe_->book2D("METTask_HCAL_D4_METvsieta","",83,-41,42,1000,0,2000);

    hHCAL_D1_METPhivsieta      = dbe_->book2D("METTask_HCAL_D1_METPhivsieta","",83,-41,42, 80, -4,4);
        hHCAL_D2_METPhivsieta      = dbe_->book2D("METTask_HCAL_D2_METPhivsieta","",83,-41,42, 80, -4,4);
        hHCAL_D3_METPhivsieta      = dbe_->book2D("METTask_HCAL_D3_METPhivsieta","",83,-41,42, 80, -4,4);
        hHCAL_D4_METPhivsieta      = dbe_->book2D("METTask_HCAL_D4_METPhivsieta","",83,-41,42, 80, -4,4);

        hHCAL_D1_MExvsieta         = dbe_->book2D("METTask_HCAL_D1_MExvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D2_MExvsieta         = dbe_->book2D("METTask_HCAL_D2_MExvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D3_MExvsieta         = dbe_->book2D("METTask_HCAL_D3_MExvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D4_MExvsieta         = dbe_->book2D("METTask_HCAL_D4_MExvsieta","",83,-41,42, 500,-500,500);

        hHCAL_D1_MEyvsieta         = dbe_->book2D("METTask_HCAL_D1_MEyvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D2_MEyvsieta         = dbe_->book2D("METTask_HCAL_D2_MEyvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D3_MEyvsieta         = dbe_->book2D("METTask_HCAL_D3_MEyvsieta","",83,-41,42, 500,-500,500);
        hHCAL_D4_MEyvsieta         = dbe_->book2D("METTask_HCAL_D4_MEyvsieta","",83,-41,42, 500,-500,500);
    
      }
    
  }
  
  // Inspect Setup for CaloTower Geometry
  //  FillGeometry(iSetup);
  
}

void HCALRecHitAnalyzer::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 884 of file HCALRecHitAnalyzer.cc.

{


} 

void HCALRecHitAnalyzer::FillGeometry ( const edm::EventSetup &  iSetup )

private

Definition at line 258 of file HCALRecHitAnalyzer.cc.

References PV3DBase< T, PVType, FrameType >::eta(), eta(), edm::EventSetup::get(), MonitorElement::getBinContent(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), CaloSubdetectorGeometry::getValidDetIds(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hHCAL_ieta_detaMap, hHCAL_ieta_dphiMap, hHCAL_ieta_iphi_etaMap, hHCAL_ieta_iphi_phiMap, i, HcalDetId::ieta(), edm::ESHandleBase::isValid(), M_PI, phi, PV3DBase< T, PVType, FrameType >::phi(), and MonitorElement::setBinContent().

Referenced by beginRun().

{

  // ==========================================================
  // Retrieve!
  // ==========================================================

  const CaloSubdetectorGeometry* HBgeom;
  const CaloSubdetectorGeometry* HEgeom;
  const CaloSubdetectorGeometry* HOgeom;
  const CaloSubdetectorGeometry* HFgeom;

  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  
  if (!pG.isValid() )
    {
      edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Setup Handle, Aborting Task "
                 << "HCALRecHitAnalyzer::FillGeometry!\n"; 
      return;
    }
  
  const CaloGeometry cG = *pG;
  
  HBgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
  HEgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
  HOgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalOuter);
  HFgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalForward);
    
  
  // ==========================================================
  // Fill Histograms!
  // ==========================================================

  std::vector<DetId>::iterator i;

  int HBmin_ieta = 99, HBmax_ieta = -99;
  int HBmin_iphi = 99, HBmax_iphi = -99;

  // Loop Over all Hcal Barrel DetId's
  int nHBdetid = 0;
  std::vector<DetId> HBids = HBgeom->getValidDetIds(DetId::Hcal,HcalBarrel);

  for (i = HBids.begin(); i != HBids.end(); i++) {

    nHBdetid++;

    const CaloCellGeometry* cell = HBgeom->getGeometry(*i);
    HcalDetId HcalID(i->rawId());
    //GlobalPoint p = cell->getPosition();

    int Calo_ieta = 42 + HcalID.ieta();
    int Calo_iphi = HcalID.iphi();
    double Calo_eta = cell->getPosition().eta();
    double Calo_phi = cell->getPosition().phi();

    if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {

      hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
      hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);

    }

    if (Calo_ieta > HBmax_ieta) HBmax_ieta = Calo_ieta;
    if (Calo_ieta < HBmin_ieta) HBmin_ieta = Calo_ieta;
    if (Calo_iphi > HBmax_iphi) HBmax_iphi = Calo_iphi;
    if (Calo_iphi > HBmax_iphi) HBmin_iphi = Calo_iphi;

  }

  int HEmin_ieta = 99, HEmax_ieta = -99;
  int HEmin_iphi = 99, HEmax_iphi = -99;

  // Loop Over all Hcal Endcap DetId's
  int nHEdetid = 0;
  std::vector<DetId> HEids = HEgeom->getValidDetIds(DetId::Hcal,HcalEndcap);

  for (i = HEids.begin(); i != HEids.end(); i++) {

    nHEdetid++;

    const CaloCellGeometry* cell = HEgeom->getGeometry(*i);
    HcalDetId HcalID(i->rawId());
    //GlobalPoint p = cell->getPosition();

    int Calo_ieta = 42 + HcalID.ieta();
    int Calo_iphi = HcalID.iphi();
    double Calo_eta = cell->getPosition().eta();
    double Calo_phi = cell->getPosition().phi();

    // HCAL to HE eta, phi map comparison
    if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
      hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
      hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
    } 

    if (Calo_ieta > HEmax_ieta) HEmax_ieta = Calo_ieta;
    if (Calo_ieta < HEmin_ieta) HEmin_ieta = Calo_ieta;
    if (Calo_iphi > HEmax_iphi) HEmax_iphi = Calo_iphi;
    if (Calo_iphi > HEmax_iphi) HEmin_iphi = Calo_iphi;

  }

  int HFmin_ieta = 99, HFmax_ieta = -99;
  int HFmin_iphi = 99, HFmax_iphi = -99;

  // Loop Over all Hcal Forward DetId's
  int nHFdetid = 0;
  std::vector<DetId> HFids = HFgeom->getValidDetIds(DetId::Hcal,HcalForward);

  for (i = HFids.begin(); i != HFids.end(); i++) {

    nHFdetid++;

    const CaloCellGeometry* cell = HFgeom->getGeometry(*i);
    HcalDetId HcalID(i->rawId());
    //GlobalPoint p = cell->getPosition();

    int Calo_ieta = 42 + HcalID.ieta();
    int Calo_iphi = HcalID.iphi();
    double Calo_eta = cell->getPosition().eta();
    double Calo_phi = cell->getPosition().phi();

    // HCAL to HF eta, phi map comparison
    if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
      hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
      hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
    } 

    if (Calo_ieta > HFmax_ieta) HFmax_ieta = Calo_ieta;
    if (Calo_ieta < HFmin_ieta) HFmin_ieta = Calo_ieta;
    if (Calo_iphi > HFmax_iphi) HFmax_iphi = Calo_iphi;
    if (Calo_iphi > HFmax_iphi) HFmin_iphi = Calo_iphi;

  }

  int HOmin_ieta = 99, HOmax_ieta = -99;
  int HOmin_iphi = 99, HOmax_iphi = -99;

  // Loop Over all Hcal Outer DetId's
  int nHOdetid = 0;
  std::vector<DetId> HOids = HOgeom->getValidDetIds(DetId::Hcal,HcalOuter);

  for (i = HOids.begin(); i != HOids.end(); i++) {

    nHOdetid++;

    const CaloCellGeometry* cell = HOgeom->getGeometry(*i);
    HcalDetId HcalID(i->rawId());
    //GlobalPoint p = cell->getPosition();

    int Calo_ieta = 42 + HcalID.ieta();
    int Calo_iphi = HcalID.iphi();
    double Calo_eta = cell->getPosition().eta();
    double Calo_phi = cell->getPosition().phi();

    // HCAL to HO eta, phi map comparison
    if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
      hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
      hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
    } 

    if (Calo_ieta > HOmax_ieta) HOmax_ieta = Calo_ieta;
    if (Calo_ieta < HOmin_ieta) HOmin_ieta = Calo_ieta;
    if (Calo_iphi > HOmax_iphi) HOmax_iphi = Calo_iphi;
    if (Calo_iphi > HOmax_iphi) HOmin_iphi = Calo_iphi;

  }

  // Set the Cell Size for each (ieta, iphi) Bin
  double currentLowEdge_eta = 0; //double currentHighEdge_eta = 0;
  for (int ieta = 1; ieta < 42 ; ieta++) {
    
    int ieta_ = 42 + ieta;
    double eta = hHCAL_ieta_iphi_etaMap->getBinContent(ieta_,3);
    double phi = hHCAL_ieta_iphi_phiMap->getBinContent(ieta_,3);
    double deta = 2.0*(eta-currentLowEdge_eta);
    deta = ((float)((int)(1.0E3*deta + 0.5)))/1.0E3;
    double dphi = 2.0*phi;
    if (ieta==40 || ieta==41) dphi = 20;
    if (ieta<=39 && ieta>=21) dphi = 10;
    if (ieta<=20) dphi = 5;
    // BS: This is WRONG...need to correct overlap 
    if (ieta == 28) deta = 0.218;
    if (ieta == 29) deta= 0.096;      
    currentLowEdge_eta += deta;

    // BS: This is WRONG...need to correct overlap 
    if (ieta == 29) currentLowEdge_eta = 2.964;

    hHCAL_ieta_detaMap->setBinContent(ieta_,deta); // positive rings
    hHCAL_ieta_dphiMap->setBinContent(ieta_,dphi); // positive rings
    hHCAL_ieta_detaMap->setBinContent(42-ieta,deta); // negative rings
    hHCAL_ieta_dphiMap->setBinContent(42-ieta,dphi); // negative rings

  } // end loop over ieta

  edm::LogInfo("OutputInfo") << "HB ieta range: " << HBmin_ieta << " " << HBmax_ieta;
  edm::LogInfo("OutputInfo") << "HB iphi range: " << HBmin_iphi << " " << HBmax_iphi;
  edm::LogInfo("OutputInfo") << "HE ieta range: " << HEmin_ieta << " " << HEmax_ieta;
  edm::LogInfo("OutputInfo") << "HE iphi range: " << HEmin_iphi << " " << HEmax_iphi;
  edm::LogInfo("OutputInfo") << "HF ieta range: " << HFmin_ieta << " " << HFmax_ieta;
  edm::LogInfo("OutputInfo") << "HF iphi range: " << HFmin_iphi << " " << HFmax_iphi;
  edm::LogInfo("OutputInfo") << "HO ieta range: " << HOmin_ieta << " " << HOmax_ieta;
  edm::LogInfo("OutputInfo") << "HO iphi range: " << HOmin_iphi << " " << HOmax_iphi;

}

Member Data Documentation

DQMStore* HCALRecHitAnalyzer::dbe_

private

Definition at line 34 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos().

bool HCALRecHitAnalyzer::debug_

private

Definition at line 40 of file HCALRecHitAnalyzer.h.

Referenced by HCALRecHitAnalyzer().

bool HCALRecHitAnalyzer::finebinning_

private

Definition at line 41 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos(), and HCALRecHitAnalyzer().

std::string HCALRecHitAnalyzer::FolderName_

private

Definition at line 42 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos(), and HCALRecHitAnalyzer().

edm::InputTag HCALRecHitAnalyzer::hBHERecHitsLabel_

private

Definition at line 37 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and HCALRecHitAnalyzer().

edm::InputTag HCALRecHitAnalyzer::hFRecHitsLabel_

private

Definition at line 38 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and HCALRecHitAnalyzer().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_energy_ieta_iphi

private

Definition at line 60 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_energyvsieta

private

Definition at line 80 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Maxenergy_ieta_iphi

private

Definition at line 70 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Maxenergyvsieta

private

Definition at line 90 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_METPhivsieta

private

Definition at line 110 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_METvsieta

private

Definition at line 105 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_MExvsieta

private

Definition at line 115 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_MEyvsieta

private

Definition at line 120 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Minenergy_ieta_iphi

private

Definition at line 65 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Minenergyvsieta

private

Definition at line 85 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Occ_ieta_iphi

private

Definition at line 75 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_Occvsieta

private

Definition at line 95 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D1_SETvsieta

private

Definition at line 100 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_energy_ieta_iphi

private

Definition at line 61 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_energyvsieta

private

Definition at line 81 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Maxenergy_ieta_iphi

private

Definition at line 71 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Maxenergyvsieta

private

Definition at line 91 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_METPhivsieta

private

Definition at line 111 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_METvsieta

private

Definition at line 106 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_MExvsieta

private

Definition at line 116 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_MEyvsieta

private

Definition at line 121 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Minenergy_ieta_iphi

private

Definition at line 66 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Minenergyvsieta

private

Definition at line 86 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Occ_ieta_iphi

private

Definition at line 76 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_Occvsieta

private

Definition at line 96 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D2_SETvsieta

private

Definition at line 101 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_energy_ieta_iphi

private

Definition at line 62 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_energyvsieta

private

Definition at line 82 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Maxenergy_ieta_iphi

private

Definition at line 72 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Maxenergyvsieta

private

Definition at line 92 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_METPhivsieta

private

Definition at line 112 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_METvsieta

private

Definition at line 107 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_MExvsieta

private

Definition at line 117 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_MEyvsieta

private

Definition at line 122 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Minenergy_ieta_iphi

private

Definition at line 67 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Minenergyvsieta

private

Definition at line 87 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Occ_ieta_iphi

private

Definition at line 77 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_Occvsieta

private

Definition at line 97 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D3_SETvsieta

private

Definition at line 102 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_energy_ieta_iphi

private

Definition at line 63 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_energyvsieta

private

Definition at line 83 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Maxenergy_ieta_iphi

private

Definition at line 73 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Maxenergyvsieta

private

Definition at line 93 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_METPhivsieta

private

Definition at line 113 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_METvsieta

private

Definition at line 108 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_MExvsieta

private

Definition at line 118 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_MEyvsieta

private

Definition at line 123 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Minenergy_ieta_iphi

private

Definition at line 68 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Minenergyvsieta

private

Definition at line 88 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Occ_ieta_iphi

private

Definition at line 78 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_Occvsieta

private

Definition at line 98 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_D4_SETvsieta

private

Definition at line 103 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_detaMap

private

Definition at line 55 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos(), and FillGeometry().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_dphiMap

private

Definition at line 56 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos(), and FillGeometry().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_etaMap

private

Definition at line 53 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), BookHistos(), and FillGeometry().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_HBMap

private

Definition at line 49 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_HEMap

private

Definition at line 50 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_HFMap

private

Definition at line 51 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_HOMap

private

Definition at line 52 of file HCALRecHitAnalyzer.h.

Referenced by BookHistos().

MonitorElement* HCALRecHitAnalyzer::hHCAL_ieta_iphi_phiMap

private

Definition at line 54 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), BookHistos(), and FillGeometry().

MonitorElement* HCALRecHitAnalyzer::hHCAL_Nevents

private

Definition at line 58 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and BookHistos().

edm::InputTag HCALRecHitAnalyzer::hORecHitsLabel_

private

Definition at line 39 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and HCALRecHitAnalyzer().

int HCALRecHitAnalyzer::Nevents

private

Definition at line 46 of file HCALRecHitAnalyzer.h.

Referenced by analyze(), and beginRun().