CaloTowersValidation Class Reference

#include <CaloTowersValidation.h>

Inheritance diagram for CaloTowersValidation:

Public Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &c) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	CaloTowersValidation (edm::ParameterSet const &conf)
	~CaloTowersValidation () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Private Types
typedef math::RhoEtaPhiVector	Vector

Private Member Functions
double	dR (double eta1, double phi1, double eta2, double phi2)

Private Attributes
MonitorElement *	emean_vs_ieta_E
MonitorElement *	emean_vs_ieta_E1
MonitorElement *	emean_vs_ieta_EH
MonitorElement *	emean_vs_ieta_EH1
MonitorElement *	emean_vs_ieta_H
MonitorElement *	emean_vs_ieta_H1
MonitorElement *	emEnergyTiming_HB
MonitorElement *	emEnergyTiming_HE
MonitorElement *	emEnergyTiming_HF
MonitorElement *	emEnergyTiming_High_HB
MonitorElement *	emEnergyTiming_Low_HB
MonitorElement *	emEnergyTiming_Low_HE
MonitorElement *	emEnergyTiming_profile_HB
MonitorElement *	emEnergyTiming_profile_HE
MonitorElement *	emEnergyTiming_profile_HF
MonitorElement *	emEnergyTiming_profile_High_HB
MonitorElement *	emEnergyTiming_profile_Low_HB
MonitorElement *	emEnergyTiming_profile_Low_HE
MonitorElement *	emTiming_HB
MonitorElement *	emTiming_HE
MonitorElement *	emTiming_HF
double	etaMax [3]
double	etaMin [3]
MonitorElement *	hadEnergyTiming_HB
MonitorElement *	hadEnergyTiming_HE
MonitorElement *	hadEnergyTiming_HF
MonitorElement *	hadEnergyTiming_High_HB
MonitorElement *	hadEnergyTiming_Low_HB
MonitorElement *	hadEnergyTiming_Low_HE
MonitorElement *	hadEnergyTiming_Low_HF
MonitorElement *	hadEnergyTiming_profile_HB
MonitorElement *	hadEnergyTiming_profile_HE
MonitorElement *	hadEnergyTiming_profile_HF
MonitorElement *	hadEnergyTiming_profile_High_HB
MonitorElement *	hadEnergyTiming_profile_Low_HB
MonitorElement *	hadEnergyTiming_profile_Low_HE
MonitorElement *	hadEnergyTiming_profile_Low_HF
MonitorElement *	hadTiming_HB
MonitorElement *	hadTiming_HE
MonitorElement *	hadTiming_HF
std::string	hcalselector_
int	imc
int	isub
MonitorElement *	mapEnergy_E
MonitorElement *	mapEnergy_EH
MonitorElement *	mapEnergy_H
MonitorElement *	mapEnergy_HB
MonitorElement *	mapEnergy_HE
MonitorElement *	mapEnergy_HF
MonitorElement *	mapEnergy_N
MonitorElement *	mapEnergyEcal_HB
MonitorElement *	mapEnergyEcal_HE
MonitorElement *	mapEnergyEcal_HF
MonitorElement *	mapEnergyHcal_HB
MonitorElement *	mapEnergyHcal_HE
MonitorElement *	mapEnergyHcal_HF
std::string	mc_
MonitorElement *	meEnergyEcal_HB
MonitorElement *	meEnergyEcal_HE
MonitorElement *	meEnergyEcal_HF
MonitorElement *	meEnergyEcalTower_HB
MonitorElement *	meEnergyEcalTower_HE
MonitorElement *	meEnergyEcalTower_HF
MonitorElement *	meEnergyHcal_HB
MonitorElement *	meEnergyHcal_HE
MonitorElement *	meEnergyHcal_HF
MonitorElement *	meEnergyHcalTower_HB
MonitorElement *	meEnergyHcalTower_HE
MonitorElement *	meEnergyHcalTower_HF
MonitorElement *	meEnergyHcalvsEcal_HB
MonitorElement *	meEnergyHcalvsEcal_HE
MonitorElement *	meEnergyHcalvsEcal_HF
MonitorElement *	meEnergyHO_HB
MonitorElement *	meEnergyHO_HE
MonitorElement *	meEnergyHO_HF
MonitorElement *	meNumFiredTowers_HB
MonitorElement *	meNumFiredTowers_HE
MonitorElement *	meNumFiredTowers_HF
MonitorElement *	MET_HB
MonitorElement *	MET_HE
MonitorElement *	MET_HF
MonitorElement *	meTotEnergy_HB
MonitorElement *	meTotEnergy_HE
MonitorElement *	meTotEnergy_HF
int	nevent
MonitorElement *	Ntowers_vs_ieta
MonitorElement *	numBadCellsEcal_EB
MonitorElement *	numBadCellsEcal_EE
MonitorElement *	numBadCellsHcal_HB
MonitorElement *	numBadCellsHcal_HE
MonitorElement *	numBadCellsHcal_HF
MonitorElement *	numPrbCellsEcal_EB
MonitorElement *	numPrbCellsEcal_EE
MonitorElement *	numPrbCellsHcal_HB
MonitorElement *	numPrbCellsHcal_HE
MonitorElement *	numPrbCellsHcal_HF
MonitorElement *	numRcvCellsEcal_EB
MonitorElement *	numRcvCellsEcal_EE
MonitorElement *	numRcvCellsHcal_HB
MonitorElement *	numRcvCellsHcal_HE
MonitorElement *	numRcvCellsHcal_HF
MonitorElement *	occupancy_map
MonitorElement *	occupancy_vs_ieta
std::string	outputFile_
MonitorElement *	phiMET_HB
MonitorElement *	phiMET_HE
MonitorElement *	phiMET_HF
MonitorElement *	SET_HB
MonitorElement *	SET_HE
MonitorElement *	SET_HF
edm::EDGetTokenT< CaloTowerCollection >	tok_calo_
edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_
bool	useAllHistos_

Detailed Description

Definition at line 33 of file CaloTowersValidation.h.

Member Typedef Documentation

typedef math::RhoEtaPhiVector CaloTowersValidation::Vector

private

Definition at line 48 of file CaloTowersValidation.h.

Constructor & Destructor Documentation

CaloTowersValidation::CaloTowersValidation

(

edm::ParameterSet const &

conf

)

Definition at line 6 of file CaloTowersValidation.cc.

References etaMax, etaMin, edm::ParameterSet::getUntrackedParameter(), hcalselector_, imc, isub, mc_, nevent, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, tok_calo_, tok_evt_, and useAllHistos_.

{

  tok_calo_ = consumes<CaloTowerCollection>(conf.getUntrackedParameter<edm::InputTag>("CaloTowerCollectionLabel"));
    tok_evt_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));

  // DQM ROOT output
  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");

  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");

  mc_           = conf.getUntrackedParameter<std::string>("mc", "yes");
  useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);

  etaMin[0] = 0.;
  etaMax[0] = 1.4;
  etaMin[1] = 1.4;
  etaMax[1] = 2.9;
  etaMin[2] = 2.9;
  etaMax[2] = 5.2;

  isub = 0;
  if(hcalselector_ == "HB") isub = 1;
  if(hcalselector_ == "HE") isub = 2;
  if(hcalselector_ == "HF") isub = 3;

  imc = 1;
  if(mc_ == "no") imc = 0;
  
  if ( !outputFile_.empty() ) {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
  }

  nevent = 0;
    // const char * sub = hcalselector_.c_str();

}

CaloTowersValidation::~CaloTowersValidation ( )

override

Definition at line 47 of file CaloTowersValidation.cc.

                                            {

}

Member Function Documentation

void CaloTowersValidation::analyze ( edm::Event const &  e, edm::EventSetup const &  c  )

override

Definition at line 405 of file CaloTowersValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, DetId::det(), dR(), DetId::Ecal, EcalBarrel, EcalEndcap, emean_vs_ieta_E, emean_vs_ieta_E1, emean_vs_ieta_EH, emean_vs_ieta_EH1, emean_vs_ieta_H, emean_vs_ieta_H1, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, edm::SortedCollection< T, SORT >::end(), etaMax, etaMin, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::HepMCProduct::GetEvent(), hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, CaloTowerDetId::ieta(), imc, isub, edm::HandleBase::isValid(), mapEnergy_E, mapEnergy_EH, mapEnergy_H, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, MuonErrorMatrixAnalyzer_cfi::maxPt, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, RazorAnalyzer::met, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, nevent, npart, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, AlCaHLTBitMon_ParallelJobs::p, phiMET_HB, phiMET_HE, phiMET_HF, EnergyCorrector::pt, alignCSCRings::r, Gflash::Rmin, SET_HB, SET_HE, SET_HF, mathSSE::sqrt(), HcalDetId::subdet(), DetId::subdetId(), tok_calo_, tok_evt_, and useAllHistos_.

                                                                                {

    nevent++;

    // bool     MC = false; // UNUSED
    double   phi_MC = 9999.;
    double   eta_MC = 9999.;

    if (imc != 0){
        edm::Handle<edm::HepMCProduct> evtMC;
        event.getByToken(tok_evt_,evtMC);  // generator in late 310_preX
        if (!evtMC.isValid()) {
            std::cout << "no HepMCProduct found" << std::endl;    
        } else {
            // MC=true; // UNUSED
            //    std::cout << "*** source HepMCProduct found"<< std::endl;
        }  

        // MC particle with highest pt is taken as a direction reference  
        double maxPt = -99999.;
        int npart    = 0;
        const HepMC::GenEvent * myGenEvent = evtMC->GetEvent();
        for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
                p != myGenEvent->particles_end(); ++p ) {
            double phip = (*p)->momentum().phi();
            double etap = (*p)->momentum().eta();
            //    phi_MC = phip;
            //    eta_MC = etap;
            double pt  = (*p)->momentum().perp();
            if(pt > maxPt) { npart++; maxPt = pt; phi_MC = phip; eta_MC = etap; }
        }
        //  std::cout << "*** Max pT = " << maxPt <<  std::endl;  
    }    

    edm::Handle<CaloTowerCollection> towers;
    event.getByToken(tok_calo_,towers);
    CaloTowerCollection::const_iterator cal;

    double met;
    double phimet;

    // ieta scan 
    double partR  = 0.3;
    double Rmin   = 9999.;
    double Econe  = 0.;
    double Hcone  = 0.;
    double Ee1    = 0.;
    double Eh1    = 0.;
    double ieta_MC = 9999;
    double iphi_MC = 9999;
    //  double  etaM   = 9999.;


    // HB   
    double sumEnergyHcal_HB = 0.;
    double sumEnergyEcal_HB = 0.;
    double sumEnergyHO_HB   = 0.;
    Int_t numFiredTowers_HB = 0;
    double metx_HB   =  0.;
    double mety_HB   =  0.;
    double metz_HB   =  0.;
    double sEt_HB    =  0.;
    // HE   
    double sumEnergyHcal_HE = 0.;
    double sumEnergyEcal_HE = 0.;
    double sumEnergyHO_HE   = 0.;
    Int_t numFiredTowers_HE = 0;
    double metx_HE   =  0.;
    double mety_HE   =  0.;
    double metz_HE   =  0.;
    double sEt_HE    =  0.;
    // HF   
    double sumEnergyHcal_HF = 0.;
    double sumEnergyEcal_HF = 0.;
    double sumEnergyHO_HF   = 0.;
    Int_t numFiredTowers_HF = 0;
    double metx_HF   =  0.;
    double mety_HF   =  0.;
    double metz_HF   =  0.;
    double sEt_HF    =  0.;

    for ( cal = towers->begin(); cal != towers->end(); ++cal ) {

        double eE     = cal->emEnergy();
        double eH     = cal->hadEnergy();
        double eHO    = cal->outerEnergy();
        double etaT   = cal->eta();
        double phiT   = cal->phi();
        double en     = cal->energy();
        double etT    = cal->et();
        double had_tm = cal->hcalTime();
        double em_tm  = cal->ecalTime();

        int numBadEcalCells = cal->numBadEcalCells();
        int numRcvEcalCells = cal->numRecoveredEcalCells();
        int numPrbEcalCells = cal->numProblematicEcalCells();

        int numBadHcalCells = cal->numBadHcalCells();
        int numRcvHcalCells = cal->numRecoveredHcalCells();
        int numPrbHcalCells = cal->numProblematicHcalCells(); 

        math::RhoEtaPhiVector mom(cal->et(), cal->eta(), cal->phi());
        //  Vector mom  = cal->momentum(); 

        // cell properties    
        CaloTowerDetId idT = cal->id();
        int ieta = idT.ieta();
        int iphi = idT.iphi();

        // ecal:  0 EcalBarrel  1 EcalEndcap
        // hcal:  0 hcalBarrel  1 HcalEndcap  2 HcalForward
        std::vector<int> inEcals(2), inHcals(3);
        unsigned int constitSize = cal->constituentsSize();
        for(unsigned int ic=0; ic<constitSize; ic++){
            DetId detId = cal->constituent(ic);
            if(detId.det() == DetId::Ecal){
                if( detId.subdetId()==EcalBarrel ) inEcals[0] =1;
                else if( detId.subdetId()==EcalEndcap ) inEcals[1] =1;
            }
            if(detId.det() == DetId::Hcal){
                if( HcalDetId(detId).subdet()==HcalBarrel ) inHcals[0] =1;
                else if( HcalDetId(detId).subdet()==HcalEndcap ) inHcals[1] =1;
                else if( HcalDetId(detId).subdet()==HcalForward ) inHcals[2] =1;
            }
        }
        //All cell histos are used
        if( inEcals[0] ) {
            numBadCellsEcal_EB->Fill(numBadEcalCells);
            numRcvCellsEcal_EB->Fill(numRcvEcalCells);
            numPrbCellsEcal_EB->Fill(numPrbEcalCells);
        }
        if( inEcals[1] ) {
            numBadCellsEcal_EE->Fill(numBadEcalCells);
            numRcvCellsEcal_EE->Fill(numRcvEcalCells);
            numPrbCellsEcal_EE->Fill(numPrbEcalCells);
        } 

        if (imc != 0){
            double r    = dR(eta_MC, phi_MC, etaT, phiT);

            if( r < partR ){
                Econe += eE; 
                Hcone += eH;

                // closest to MC
                if(r < Rmin) { 
                    if( fabs(eta_MC) < 3.0 && (ieta > 29 || ieta < -29)) {;}
                    else {    
                        Rmin = r;
                        ieta_MC = ieta;
                        iphi_MC = iphi; 
                        Ee1     = eE;
                        Eh1     = eH;
                    }
                }
            }
        }

        //Ntowers is used in EndJob, occupancy_map is used for occupancy vs ieta
        Ntowers_vs_ieta -> Fill(double(ieta),1.);
        occupancy_map -> Fill(double(ieta),double(iphi));


        if((isub == 0 || isub == 1) 
                && (fabs(etaT) <  etaMax[0] && fabs(etaT) >= etaMin[0] )) {

            //All cell histos are used
            numBadCellsHcal_HB->Fill(numBadHcalCells);
            numRcvCellsHcal_HB->Fill(numRcvHcalCells);
            numPrbCellsHcal_HB->Fill(numPrbHcalCells);

            //Map energy histos are not used
            if (useAllHistos_){
                mapEnergy_HB     -> Fill(double(ieta), double(iphi), en); 
                mapEnergyHcal_HB -> Fill(double(ieta), double(iphi), eH); 
                mapEnergyEcal_HB -> Fill(double(ieta), double(iphi), eE); 
            }
            //      std::cout << " e_ecal = " << eE << std::endl;

            //  simple sums
            sumEnergyHcal_HB += eH;
            sumEnergyEcal_HB += eE;
            sumEnergyHO_HB   += eHO;

            numFiredTowers_HB++;

            //Not used
            if (useAllHistos_){
                meEnergyEcalTower_HB->Fill(eE);
                meEnergyHcalTower_HB->Fill(eH);    
            }

            // MET, SET & phimet
            //  double  etT = cal->et();
            metx_HB += mom.x();   
            mety_HB += mom.y();  //etT * sin(phiT);          
            sEt_HB  += etT;

            //Timing (all histos are used)
            emTiming_HB->Fill(em_tm);
            hadTiming_HB->Fill(had_tm);

            emEnergyTiming_Low_HB->Fill(eE, em_tm);
            emEnergyTiming_HB->Fill(eE, em_tm);
            emEnergyTiming_High_HB->Fill(eE, em_tm);
            emEnergyTiming_profile_Low_HB->Fill(eE, em_tm);
            emEnergyTiming_profile_HB->Fill(eE, em_tm);
            emEnergyTiming_profile_High_HB->Fill(eE, em_tm);

            hadEnergyTiming_Low_HB->Fill(eH, had_tm); 
            hadEnergyTiming_HB->Fill(eH, had_tm); 
            hadEnergyTiming_High_HB->Fill(eH, had_tm); 
            hadEnergyTiming_profile_Low_HB->Fill(eH, had_tm);     
            hadEnergyTiming_profile_HB->Fill(eH, had_tm);     
            hadEnergyTiming_profile_High_HB->Fill(eH, had_tm);     
        }

        if((isub == 0 || isub == 2) 
                && (fabs(etaT) <  etaMax[1] && fabs(etaT) >= etaMin[1] )) {

            //All cell histos are used
            numBadCellsHcal_HE->Fill(numBadHcalCells);
            numRcvCellsHcal_HE->Fill(numRcvHcalCells);
            numPrbCellsHcal_HE->Fill(numPrbHcalCells);

            //Map energy histos are not used
            if (useAllHistos_){
                mapEnergy_HE     -> Fill(double(ieta), double(iphi), en); 
                mapEnergyHcal_HE -> Fill(double(ieta), double(iphi), eH); 
                mapEnergyEcal_HE -> Fill(double(ieta), double(iphi), eE); 
            }
            //      std::cout << " e_ecal = " << eE << std::endl;

            //  simple sums
            sumEnergyHcal_HE += eH;
            sumEnergyEcal_HE += eE;
            sumEnergyHO_HE   += eHO;

            numFiredTowers_HE++;

            //Not used
            if (useAllHistos_){
                meEnergyEcalTower_HE->Fill(eE);
                meEnergyHcalTower_HE->Fill(eH);    
            }
            // MET, SET & phimet
            //  double  etT = cal->et();
            metx_HE += mom.x();   
            mety_HE += mom.y();  //etT * sin(phiT);          
            sEt_HE  += etT;    

            //Timing (all histos are used)
            emTiming_HE->Fill(em_tm);
            hadTiming_HE->Fill(had_tm);

            emEnergyTiming_Low_HE->Fill(eE, em_tm);
            emEnergyTiming_HE->Fill(eE, em_tm);
            emEnergyTiming_profile_Low_HE->Fill(eE, em_tm);
            emEnergyTiming_profile_HE->Fill(eE, em_tm);

            hadEnergyTiming_Low_HE->Fill(eH, had_tm);     
            hadEnergyTiming_HE->Fill(eH, had_tm);   
            hadEnergyTiming_profile_Low_HE->Fill(eH, had_tm);     
            hadEnergyTiming_profile_HE->Fill(eH, had_tm);     
        }

        if((isub == 0 || isub == 3) 
                && (fabs(etaT) <  etaMax[2] && fabs(etaT) >= etaMin[2] )) {

            //All cell histos are used
            numBadCellsHcal_HF->Fill(numBadHcalCells);
            numRcvCellsHcal_HF->Fill(numRcvHcalCells);
            numPrbCellsHcal_HF->Fill(numPrbHcalCells);

            //Map energy histos are not used
            if (useAllHistos_){
                mapEnergy_HF     -> Fill(double(ieta), double(iphi), en); 
                mapEnergyHcal_HF -> Fill(double(ieta), double(iphi), eH); 
                mapEnergyEcal_HF -> Fill(double(ieta), double(iphi), eE); 
            }
            //      std::cout << " e_ecal = " << eE << std::endl;

            //  simple sums
            sumEnergyHcal_HF += eH;
            sumEnergyEcal_HF += eE;
            sumEnergyHO_HF   += eHO;

            numFiredTowers_HF++;

            //Not used
            if (useAllHistos_){
                meEnergyEcalTower_HF->Fill(eE);
                meEnergyHcalTower_HF->Fill(eH);    
            }
            // MET, SET & phimet
            //  double  etT = cal->et();
            metx_HF += mom.x();   
            mety_HF += mom.y();  //etT * sin(phiT);          
            sEt_HF  += etT;    

            //Timing (all histos are used)
            emTiming_HF->Fill(em_tm);
            hadTiming_HF->Fill(had_tm);
            emEnergyTiming_HF->Fill(eE, em_tm);
            emEnergyTiming_profile_HF->Fill(eE, em_tm);

            hadEnergyTiming_Low_HF->Fill(eH, had_tm);     
            hadEnergyTiming_HF->Fill(eH, had_tm);     
            hadEnergyTiming_profile_Low_HF->Fill(eH, had_tm);     
            hadEnergyTiming_profile_HF->Fill(eH, had_tm);     
        }



    } // end of Towers cycle 

    //These are the six single pion histos
    emean_vs_ieta_E  -> Fill(double(ieta_MC), Econe); 
    emean_vs_ieta_H  -> Fill(double(ieta_MC), Hcone); 
    emean_vs_ieta_EH -> Fill(double(ieta_MC), Econe+Hcone); 
    emean_vs_ieta_E1  -> Fill(double(ieta_MC), Ee1); 
    emean_vs_ieta_H1  -> Fill(double(ieta_MC), Eh1); 
    emean_vs_ieta_EH1 -> Fill(double(ieta_MC), Ee1+Eh1); 

    //Map histos are not used except the last one in EndJob
    if (useAllHistos_){
        mapEnergy_E -> Fill(double(ieta_MC), double(iphi_MC), Ee1);
        mapEnergy_H -> Fill(double(ieta_MC), double(iphi_MC), Eh1); 
        mapEnergy_EH -> Fill(double(ieta_MC), double(iphi_MC), Ee1+Eh1);
    } 
    mapEnergy_N  -> Fill(double(ieta_MC), double(iphi_MC), 1.); 


    if(isub == 0 || isub == 1) {
        met    = sqrt(metx_HB*metx_HB + mety_HB*mety_HB);
        Vector metv(metx_HB,mety_HB,metz_HB);
        phimet = metv.phi();

        //Five oldest drawn histos first; the rest are not used
        meEnergyHcal_HB->      Fill(sumEnergyHcal_HB);
        meEnergyEcal_HB->      Fill(sumEnergyEcal_HB);
        meNumFiredTowers_HB->  Fill(numFiredTowers_HB);
        MET_HB    -> Fill (met); 
        SET_HB    -> Fill (sEt_HB); 

        if (useAllHistos_){
            meEnergyHcalvsEcal_HB->Fill(sumEnergyEcal_HB, sumEnergyHcal_HB);
            meEnergyHO_HB->        Fill(sumEnergyHO_HB);
            meTotEnergy_HB->       Fill(sumEnergyHcal_HB+sumEnergyEcal_HB
                    +sumEnergyHO_HB);    
            phiMET_HB -> Fill (phimet);
        } 
    }    


    if(isub == 0 || isub == 2) {
        met    = sqrt(metx_HE*metx_HE + mety_HE*mety_HE);
        Vector metv(metx_HE,mety_HE,metz_HE);
        phimet = metv.phi();

        //Five oldest drawn histos first; the rest are not used
        meEnergyHcal_HE->      Fill(sumEnergyHcal_HE);
        meEnergyEcal_HE->      Fill(sumEnergyEcal_HE);
        meNumFiredTowers_HE->  Fill(numFiredTowers_HE);
        MET_HE    -> Fill (met); 
        SET_HE    -> Fill (sEt_HE); 

        if (useAllHistos_){
            meEnergyHcalvsEcal_HE->Fill(sumEnergyEcal_HE, sumEnergyHcal_HE);
            meEnergyHO_HE->        Fill(sumEnergyHO_HE);
            meTotEnergy_HE->       Fill(sumEnergyHcal_HE+sumEnergyEcal_HE
                    +sumEnergyHO_HE);    
            phiMET_HE -> Fill (phimet); 
        }
    }

    if(isub == 0 || isub == 3) {
        met    = sqrt(metx_HF*metx_HF + mety_HF*mety_HF);
        Vector metv(metx_HF,mety_HF,metz_HF);
        phimet = metv.phi();

        //Five oldest drawn histos first; the rest are not used
        meEnergyHcal_HF->      Fill(sumEnergyHcal_HF);
        meEnergyEcal_HF->      Fill(sumEnergyEcal_HF);
        meNumFiredTowers_HF->  Fill(numFiredTowers_HF);
        MET_HF    -> Fill (met); 
        SET_HF    -> Fill (sEt_HF); 

        if (useAllHistos_){
            meEnergyHcalvsEcal_HF->Fill(sumEnergyEcal_HF, sumEnergyHcal_HF);
            meEnergyHO_HF->        Fill(sumEnergyHO_HF);
            meTotEnergy_HF->       Fill(sumEnergyHcal_HF+sumEnergyEcal_HF
                    +sumEnergyHO_HF);    
            phiMET_HF -> Fill (phimet);
        }
    }

}

void CaloTowersValidation::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  irun, edm::EventSetup const &  isetup  )

override

Definition at line 52 of file CaloTowersValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), emean_vs_ieta_E, emean_vs_ieta_E1, emean_vs_ieta_EH, emean_vs_ieta_EH1, emean_vs_ieta_H, emean_vs_ieta_H1, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, trackerHits::histo, isub, mapEnergy_E, mapEnergy_EH, mapEnergy_H, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, occupancy_vs_ieta, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, DQMStore::IBooker::setCurrentFolder(), and useAllHistos_.

{
    Char_t histo[100];

    ibooker.setCurrentFolder("CaloTowersV/CaloTowersTask");

    //These two histos are not drawn by our macros, but they are used
    //in the EndJob for norms and such so I am leaving them alone for now
    //-------------------------------------------------------------------------------------------
    sprintf  (histo, "Ntowers_per_event_vs_ieta" );
    Ntowers_vs_ieta = ibooker.book1D(histo, histo, 83, -41.5, 41.5);

    sprintf  (histo, "CaloTowersTask_map_Nentries" );
    mapEnergy_N = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
    //-------------------------------------------------------------------------------------------

    //These the single pion scan histos
    //-------------------------------------------------------------------------------------------
    sprintf  (histo, "emean_vs_ieta_E" );
    emean_vs_ieta_E = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    sprintf  (histo, "emean_vs_ieta_H" );
    emean_vs_ieta_H = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    sprintf  (histo, "emean_vs_ieta_EH" );
    emean_vs_ieta_EH = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    //These are drawn
    sprintf  (histo, "emean_vs_ieta_E1" );
    emean_vs_ieta_E1 = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    sprintf  (histo, "emean_vs_ieta_H1" );
    emean_vs_ieta_H1 = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    sprintf  (histo, "emean_vs_ieta_EH1" );
    emean_vs_ieta_EH1 = ibooker.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
    //-------------------------------------------------------------------------------------------

    //Map energy histos are not drawn
    if (useAllHistos_){
        sprintf  (histo, "CaloTowersTask_map_energy_E" );
        mapEnergy_E = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
        sprintf  (histo, "CaloTowersTask_map_energy_H");
        mapEnergy_H = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
        sprintf  (histo, "CaloTowersTask_map_energy_EH" );
        mapEnergy_EH = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
    }

    //All ECAL cell histos are used
    // XXX: ECAL 0-25 [0-26, 26 bins]   HCAL 0-4 [0-5, 5 bins]
    sprintf  (histo, "number_of_bad_cells_Ecal_EB");
    numBadCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
    sprintf  (histo, "number_of_bad_cells_Ecal_EE");
    numBadCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
    sprintf  (histo, "number_of_recovered_cells_Ecal_EB");
    numRcvCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
    sprintf  (histo, "number_of_recovered_cells_Ecal_EE");
    numRcvCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
    sprintf  (histo, "number_of_problematic_cells_Ecal_EB");
    numPrbCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
    sprintf  (histo, "number_of_problematic_cells_Ecal_EE");
    numPrbCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26); 

    //Occupancy vs. ieta is drawn, occupancy map is needed to draw it 
    sprintf  (histo, "CaloTowersTask_map_occupancy" );
    occupancy_map = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);

    sprintf  (histo, "CaloTowersTask_occupancy_vs_ieta" );
    occupancy_vs_ieta = ibooker.book1D(histo, histo, 83, -41.5, 41.5);

    if( isub == 1 || isub == 0) {
        //All cell histos are used
        sprintf  (histo, "number_of_bad_cells_Hcal_HB");
        numBadCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_recovered_cells_Hcal_HB");
        numRcvCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_problematic_cells_Hcal_HB");
        numPrbCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5); 

        //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
        //-------------------------------------------------------------------------------------------
        sprintf (histo, "CaloTowersTask_energy_HCAL_HB" ) ;
        meEnergyHcal_HB    = ibooker.book1D(histo, histo, 4100, -200, 8000);  

        sprintf (histo, "CaloTowersTask_energy_ECAL_HB" ) ;
        meEnergyEcal_HB    = ibooker.book1D(histo, histo, 3100, -200, 6000); 

        sprintf (histo, "CaloTowersTask_number_of_fired_towers_HB" ) ;
        meNumFiredTowers_HB = ibooker.book1D(histo, histo, 1000, 0, 2000); 

        sprintf  (histo, "CaloTowersTask_MET_HB" ) ;
        MET_HB = ibooker.book1D(histo, histo, 3000, 0. , 3000. ) ;

        sprintf  (histo, "CaloTowersTask_SET_HB" ) ;
        SET_HB = ibooker.book1D(histo, histo, 8000, 0. , 8000. ) ;
        //-------------------------------------------------------------------------------------------

        //Timing histos and profiles -- all six are necessary
        //-------------------------------------------------------------------------------------------
        sprintf  (histo, "CaloTowersTask_EM_Timing_HB" ) ;
        emTiming_HB = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Timing_HB" ) ;
        hadTiming_HB = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

        //Energy-Timing histos are divided into low, medium and high to reduce memory usage
        //EM
        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_Low_HB" ) ;
        emEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HB" ) ;
        emEnergyTiming_HB = ibooker.book2D(histo, histo, 200, 0. , 400., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_High_HB" ) ;
        emEnergyTiming_High_HB = ibooker.book2D(histo, histo, 200, 0. , 3000., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HB" ) ;
        emEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HB" ) ;
        emEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 200, 0. , 400., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_High_HB" ) ;
        emEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 200, 0. , 3000., 110, -120., 100.  ) ;

        //HAD
        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HB" ) ;
        hadEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HB" ) ;
        hadEnergyTiming_HB = ibooker.book2D(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_High_HB" ) ;
        hadEnergyTiming_High_HB = ibooker.book2D(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HB" ) ;
        hadEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HB" ) ;
        hadEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_High_HB" ) ;
        hadEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;
        //-------------------------------------------------------------------------------------------

        //Everything else is not drawn
        if (useAllHistos_){
            sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HB") ;
            meEnergyHcalvsEcal_HB    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

            sprintf (histo, "CaloTowersTask_energy_OUTER_HB" ) ;
            meEnergyHO_HB    = ibooker.book1D(histo, histo, 1640, -200, 8000);   

            sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HB" ) ;
            meEnergyEcalTower_HB = ibooker.book1D(histo, histo, 440, -200, 2000); 

            sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HB" ) ;
            meEnergyHcalTower_HB = ibooker.book1D(histo, histo, 440 , -200 , 2000); 

            sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HB" ) ;
            meTotEnergy_HB = ibooker.book1D(histo, histo,400, 0., 2000.) ;

            sprintf  (histo, "CaloTowersTask_map_energy_HB" );
            mapEnergy_HB = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HB");
            mapEnergyHcal_HB = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HB" );
            mapEnergyEcal_HB = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);

            sprintf  (histo, "CaloTowersTask_phi_MET_HB" ) ;
            phiMET_HB = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
        }
    } 

    if( isub == 2 || isub == 0) {
        //All cell histos are used
        sprintf  (histo, "number_of_bad_cells_Hcal_HE");
        numBadCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_recovered_cells_Hcal_HE");
        numRcvCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_problematic_cells_Hcal_HE");
        numPrbCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5); 

        //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
        //-------------------------------------------------------------------------------------------
        sprintf (histo, "CaloTowersTask_energy_HCAL_HE" ) ;
        meEnergyHcal_HE    = ibooker.book1D(histo, histo, 1240, -200, 6000); 

        sprintf (histo, "CaloTowersTask_energy_ECAL_HE" ) ;
        meEnergyEcal_HE    = ibooker.book1D(histo, histo, 1240, -200, 6000);

        sprintf (histo, "CaloTowersTask_number_of_fired_towers_HE" ) ;
        meNumFiredTowers_HE = ibooker.book1D(histo, histo, 1000, 0, 2000); 

        sprintf  (histo, "CaloTowersTask_MET_HE" ) ;
        MET_HE = ibooker.book1D(histo, histo, 1000, 0. , 1000. ) ;

        sprintf  (histo, "CaloTowersTask_SET_HE" ) ;
        SET_HE = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
        //-------------------------------------------------------------------------------------------

        //Timing histos and profiles -- all six are necessary
        //-------------------------------------------------------------------------------------------
        sprintf  (histo, "CaloTowersTask_EM_Timing_HE" ) ;
        emTiming_HE = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Timing_HE" ) ;
        hadTiming_HE = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

        //Energy-Timing histos are divided into low and normal to reduce memory usage
        //EM
        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_Low_HE" ) ;
        emEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HE" ) ;
        emEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HE" ) ;
        emEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HE" ) ;
        emEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 110, -120., 100.  ) ;

        //HAD
        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HE" ) ;
        hadEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HE" ) ;
        hadEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HE" ) ;
        hadEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HE" ) ;
        hadEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;
        //-------------------------------------------------------------------------------------------

        //Everything else is not drawn
        if (useAllHistos_){
            sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HE") ;
            meEnergyHcalvsEcal_HE    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

            sprintf (histo, "CaloTowersTask_energy_OUTER_HE" ) ;
            meEnergyHO_HE    = ibooker.book1D(histo, histo, 440, -200, 2000);   

            sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HE" ) ;
            meEnergyEcalTower_HE = ibooker.book1D(histo, histo, 1100, -200, 2000); 

            sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HE" ) ;
            meEnergyHcalTower_HE = ibooker.book1D(histo, histo, 1100 , -200 , 2000); 

            sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HE" ) ;
            meTotEnergy_HE = ibooker.book1D(histo, histo,400, 0., 2000.) ;

            sprintf  (histo, "CaloTowersTask_map_energy_HE" );
            mapEnergy_HE = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HE");
            mapEnergyHcal_HE = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HE" );
            mapEnergyEcal_HE = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);

            sprintf  (histo, "CaloTowersTask_phi_MET_HE" ) ;
            phiMET_HE = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
        }
    }


    if( isub == 3 || isub == 0) {
        //All cell histos are used
        sprintf  (histo, "number_of_bad_cells_Hcal_HF");
        numBadCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_recovered_cells_Hcal_HF");
        numRcvCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
        sprintf  (histo, "number_of_problematic_cells_Hcal_HF");
        numPrbCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5); 

        //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
        //-------------------------------------------------------------------------------------------
        sprintf (histo, "CaloTowersTask_energy_HCAL_HF" ) ;
        meEnergyHcal_HF    = ibooker.book1D(histo, histo, 4040, -200, 20000);

        sprintf (histo, "CaloTowersTask_energy_ECAL_HF" ) ;
        meEnergyEcal_HF    = ibooker.book1D(histo, histo, 2440, -200, 12000);

        sprintf (histo, "CaloTowersTask_number_of_fired_towers_HF" ) ;
        meNumFiredTowers_HF = ibooker.book1D(histo, histo, 1000, 0, 2000); 

        sprintf  (histo, "CaloTowersTask_MET_HF" ) ;
        MET_HF = ibooker.book1D(histo, histo, 500, 0. , 500. ) ;

        sprintf  (histo, "CaloTowersTask_SET_HF" ) ;
        SET_HF = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
        //-------------------------------------------------------------------------------------------

        //Timing histos and profiles -- all six are necessary
        //-------------------------------------------------------------------------------------------
        sprintf  (histo, "CaloTowersTask_EM_Timing_HF" ) ;
        emTiming_HF = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Timing_HF" ) ;
        hadTiming_HF = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

        //EM
        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HF" ) ;
        emEnergyTiming_HF = ibooker.book2D(histo, histo, 150, 0. , 300., 110, -120., 100.  ) ;

        sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HF" ) ;
        emEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 150, 0. , 300., 110, -120., 100.  ) ;

        //HAD (requires two different sets of histograms to lower RAM usage)
        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HF" ) ;
        hadEnergyTiming_Low_HF = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HF" ) ;
        hadEnergyTiming_HF = ibooker.book2D(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HF" ) ;
        hadEnergyTiming_profile_Low_HF = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

        sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HF" ) ;
        hadEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;
        //-------------------------------------------------------------------------------------------

        //Everything else is not drawn
        if (useAllHistos_){
            sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HF") ;
            meEnergyHcalvsEcal_HF    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

            sprintf (histo, "CaloTowersTask_energy_OUTER_HF" ) ;
            meEnergyHO_HF    = ibooker.book1D(histo, histo, 440, -200, 2000);   

            sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HF" ) ;
            meEnergyEcalTower_HF = ibooker.book1D(histo, histo, 440, -200, 2000); 

            sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HF" ) ;
            meEnergyHcalTower_HF = ibooker.book1D(histo, histo, 440 , -200 , 2000); 

            sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HF" ) ;
            meTotEnergy_HF = ibooker.book1D(histo, histo, 400, 0., 2000.) ;

            sprintf  (histo, "CaloTowersTask_map_energy_HF" );
            mapEnergy_HF = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HF");
            mapEnergyHcal_HF = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);
            sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HF" );
            mapEnergyEcal_HF = ibooker.book2D(histo, histo, 83, -41.5, 41.5, 72, 0., 72.);

            sprintf  (histo, "CaloTowersTask_phi_MET_HF" ) ;
            phiMET_HF = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
        }
    }



}

double CaloTowersValidation::dR ( double  eta1, double  phi1, double  eta2, double  phi2  )

private

Definition at line 804 of file CaloTowersValidation.cc.

References DEFINE_FWK_MODULE, PI, mathSSE::sqrt(), and tmp.

Referenced by analyze().

                                                                                  { 
    double PI = 3.1415926535898;
    double deltaphi= phi1 - phi2;
    if( phi2 > phi1 ) { deltaphi= phi2 - phi1;}
    if(deltaphi > PI) { deltaphi = 2.*PI - deltaphi;}
    double deltaeta = eta2 - eta1;
    double tmp = sqrt(deltaeta* deltaeta + deltaphi*deltaphi);
    return tmp;
}

Member Data Documentation

MonitorElement* CaloTowersValidation::emean_vs_ieta_E

private

Definition at line 63 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emean_vs_ieta_E1

private

Definition at line 67 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emean_vs_ieta_EH

private

Definition at line 65 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emean_vs_ieta_EH1

private

Definition at line 69 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emean_vs_ieta_H

private

Definition at line 64 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emean_vs_ieta_H1

private

Definition at line 68 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_HB

private

Definition at line 121 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_HE

private

Definition at line 156 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_HF

private

Definition at line 185 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_High_HB

private

Definition at line 122 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_Low_HB

private

Definition at line 120 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_Low_HE

private

Definition at line 155 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_HB

private

Definition at line 124 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_HE

private

Definition at line 158 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_HF

private

Definition at line 186 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_High_HB

private

Definition at line 125 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_Low_HB

private

Definition at line 123 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emEnergyTiming_profile_Low_HE

private

Definition at line 157 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emTiming_HB

private

Definition at line 117 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emTiming_HE

private

Definition at line 152 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::emTiming_HF

private

Definition at line 183 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersValidation::etaMax[3]

private

Definition at line 59 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

double CaloTowersValidation::etaMin[3]

private

Definition at line 60 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

MonitorElement* CaloTowersValidation::hadEnergyTiming_HB

private

Definition at line 128 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_HE

private

Definition at line 161 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_HF

private

Definition at line 189 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_High_HB

private

Definition at line 129 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_Low_HB

private

Definition at line 127 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_Low_HE

private

Definition at line 160 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_Low_HF

private

Definition at line 188 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_HB

private

Definition at line 131 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_HE

private

Definition at line 163 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_HF

private

Definition at line 191 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_High_HB

private

Definition at line 132 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_Low_HB

private

Definition at line 130 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_Low_HE

private

Definition at line 162 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadEnergyTiming_profile_Low_HF

private

Definition at line 190 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadTiming_HB

private

Definition at line 118 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadTiming_HE

private

Definition at line 153 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::hadTiming_HF

private

Definition at line 184 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersValidation::hcalselector_

private

Definition at line 44 of file CaloTowersValidation.h.

Referenced by CaloTowersValidation().

int CaloTowersValidation::imc

private

Definition at line 56 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

int CaloTowersValidation::isub

private

Definition at line 53 of file CaloTowersValidation.h.

Referenced by analyze(), bookHistograms(), and CaloTowersValidation().

MonitorElement* CaloTowersValidation::mapEnergy_E

private

Definition at line 76 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_EH

private

Definition at line 78 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_H

private

Definition at line 77 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_HB

private

Definition at line 110 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_HE

private

Definition at line 145 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_HF

private

Definition at line 176 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergy_N

private

Definition at line 79 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyEcal_HB

private

Definition at line 111 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyEcal_HE

private

Definition at line 146 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyEcal_HF

private

Definition at line 177 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyHcal_HB

private

Definition at line 112 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyHcal_HE

private

Definition at line 147 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::mapEnergyHcal_HF

private

Definition at line 178 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersValidation::mc_

private

Definition at line 45 of file CaloTowersValidation.h.

Referenced by CaloTowersValidation().

MonitorElement* CaloTowersValidation::meEnergyEcal_HB

private

Definition at line 102 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyEcal_HE

private

Definition at line 137 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyEcal_HF

private

Definition at line 168 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyEcalTower_HB

private

Definition at line 106 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyEcalTower_HE

private

Definition at line 141 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyEcalTower_HF

private

Definition at line 172 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcal_HB

private

Definition at line 103 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcal_HE

private

Definition at line 138 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcal_HF

private

Definition at line 169 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalTower_HB

private

Definition at line 107 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalTower_HE

private

Definition at line 142 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalTower_HF

private

Definition at line 173 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalvsEcal_HB

private

Definition at line 100 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalvsEcal_HE

private

Definition at line 135 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHcalvsEcal_HF

private

Definition at line 166 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHO_HB

private

Definition at line 101 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHO_HE

private

Definition at line 136 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meEnergyHO_HF

private

Definition at line 167 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meNumFiredTowers_HB

private

Definition at line 104 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meNumFiredTowers_HE

private

Definition at line 139 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meNumFiredTowers_HF

private

Definition at line 170 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::MET_HB

private

Definition at line 113 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::MET_HE

private

Definition at line 148 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::MET_HF

private

Definition at line 179 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meTotEnergy_HB

private

Definition at line 108 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meTotEnergy_HE

private

Definition at line 143 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::meTotEnergy_HF

private

Definition at line 174 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

int CaloTowersValidation::nevent

private

Definition at line 54 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

MonitorElement* CaloTowersValidation::Ntowers_vs_ieta

private

Definition at line 71 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numBadCellsEcal_EB

private

Definition at line 82 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numBadCellsEcal_EE

private

Definition at line 83 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numBadCellsHcal_HB

private

Definition at line 89 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numBadCellsHcal_HE

private

Definition at line 90 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numBadCellsHcal_HF

private

Definition at line 91 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numPrbCellsEcal_EB

private

Definition at line 86 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numPrbCellsEcal_EE

private

Definition at line 87 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numPrbCellsHcal_HB

private

Definition at line 95 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numPrbCellsHcal_HE

private

Definition at line 96 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numPrbCellsHcal_HF

private

Definition at line 97 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numRcvCellsEcal_EB

private

Definition at line 84 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numRcvCellsEcal_EE

private

Definition at line 85 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numRcvCellsHcal_HB

private

Definition at line 92 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numRcvCellsHcal_HE

private

Definition at line 93 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::numRcvCellsHcal_HF

private

Definition at line 94 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::occupancy_map

private

Definition at line 72 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::occupancy_vs_ieta

private

Definition at line 73 of file CaloTowersValidation.h.

Referenced by bookHistograms().

std::string CaloTowersValidation::outputFile_

private

Definition at line 43 of file CaloTowersValidation.h.

Referenced by CaloTowersValidation().

MonitorElement* CaloTowersValidation::phiMET_HB

private

Definition at line 115 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::phiMET_HE

private

Definition at line 150 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::phiMET_HF

private

Definition at line 181 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::SET_HB

private

Definition at line 114 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::SET_HE

private

Definition at line 149 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersValidation::SET_HF

private

Definition at line 180 of file CaloTowersValidation.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<CaloTowerCollection> CaloTowersValidation::tok_calo_

private

Definition at line 50 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

edm::EDGetTokenT<edm::HepMCProduct> CaloTowersValidation::tok_evt_

private

Definition at line 51 of file CaloTowersValidation.h.

Referenced by analyze(), and CaloTowersValidation().

bool CaloTowersValidation::useAllHistos_

private

Definition at line 46 of file CaloTowersValidation.h.

Referenced by analyze(), bookHistograms(), and CaloTowersValidation().