CaloTowersClient.cc

Go to the documentation of this file.

#include "Validation/CaloTowers/interface/CaloTowersClient.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DQMServices/Core/interface/MonitorElement.h"

CaloTowersClient::CaloTowersClient(const edm::ParameterSet& iConfig):conf_(iConfig)
{

  outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile", "myfile.root");

  debug_ = false;
  verbose_ = false;
  dirName_=iConfig.getParameter<std::string>("DQMDirName");
}


CaloTowersClient::~CaloTowersClient()
{ 
  
}

void CaloTowersClient::beginJob()
{
 

}


void CaloTowersClient::dqmEndJob(DQMStore::IBooker & ibooker, DQMStore::IGetter & igetter)
{
  igetter.setCurrentFolder(dirName_);

  if (verbose_) std::cout << "\nrunClient" << std::endl; 

  std::vector<MonitorElement*> hcalMEs;

  // Since out folders are fixed to three, we can just go over these three folders
  // i.e., CaloTowersV/CaloTowersTask, HcalRecHitsV/HcalRecHitTask, NoiseRatesV/NoiseRatesTask.
  std::vector<std::string> fullPathHLTFolders = igetter.getSubdirs();
  for(unsigned int i=0;i<fullPathHLTFolders.size();i++) {

    if (verbose_) std::cout <<"\nfullPath: "<< fullPathHLTFolders[i] << std::endl;
    igetter.setCurrentFolder(fullPathHLTFolders[i]);

    std::vector<std::string> fullSubPathHLTFolders = igetter.getSubdirs();
    for(unsigned int j=0;j<fullSubPathHLTFolders.size();j++) {

      if (verbose_) std::cout <<"fullSub: "<<fullSubPathHLTFolders[j] << std::endl;

      if( strcmp(fullSubPathHLTFolders[j].c_str(), "CaloTowersV/CaloTowersTask") ==0  ){
         hcalMEs = igetter.getContents(fullSubPathHLTFolders[j]);
         if (verbose_) std::cout <<"hltMES size : "<<hcalMEs.size()<<std::endl;
         if( !CaloTowersEndjob(hcalMEs) ) std::cout<<"\nError in CaloTowersEndjob!"<<std::endl<<std::endl;
      }

    }    

  }

}

// called after entering the CaloTowersV/CaloTowersTask directory
// hcalMEs are within that directory
int CaloTowersClient::CaloTowersEndjob(const std::vector<MonitorElement*> &hcalMEs){

   int useAllHistos = 0;
   MonitorElement* Ntowers_vs_ieta =nullptr;
   MonitorElement* mapEnergy_N =nullptr, *mapEnergy_E =nullptr, *mapEnergy_H =nullptr, *mapEnergy_EH =nullptr;
   MonitorElement* occupancy_map =nullptr, *occupancy_vs_ieta =nullptr;
   for(unsigned int ih=0; ih<hcalMEs.size(); ih++){
      if( strcmp(hcalMEs[ih]->getName().c_str(), "Ntowers_per_event_vs_ieta") ==0  ){
         Ntowers_vs_ieta = hcalMEs[ih];
      } 
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_Nentries") ==0  ){
         mapEnergy_N = hcalMEs[ih];
      }
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_H") ==0  ){
         useAllHistos++; mapEnergy_H = hcalMEs[ih];
      }
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_E") ==0  ){
         useAllHistos++; mapEnergy_E = hcalMEs[ih];
      }
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_EH") ==0  ){
         useAllHistos++; mapEnergy_EH = hcalMEs[ih];
      }
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_occupancy") ==0  ){
         occupancy_map = hcalMEs[ih];
      }
      if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_occupancy_vs_ieta") ==0  ){
         occupancy_vs_ieta = hcalMEs[ih];
      }
  
   } 
   if( useAllHistos !=0 && useAllHistos !=3 ) return 0;

   double nevent = mapEnergy_N->getEntries();
   if(verbose_) std::cout<<"nevent : "<<nevent<<std::endl;

   // mean number of towers per ieta
   int nx = Ntowers_vs_ieta->getNbinsX();
   float cont;
   float econt;
   float fev = float(nevent);

   for (int i = 1; i <= nx; i++) {
      cont = Ntowers_vs_ieta -> getBinContent(i) / fev ;
      econt = Ntowers_vs_ieta -> getBinError(i) / fev ;
      Ntowers_vs_ieta -> setBinContent(i,cont);
      Ntowers_vs_ieta -> setBinError(i,econt);
   }

   // mean energies & occupancies evaluation

   nx = mapEnergy_N->getNbinsX();
   int ny = mapEnergy_N->getNbinsY();
   float cnorm;
   float enorm;
   float phi_factor;

   for (int i = 1; i <= nx; i++) {
      float sumphi = 0.;
      float sumphie = 0.;

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

         // Emean
         cnorm   = mapEnergy_N -> getBinContent(i,j);
         //Phi histos are not used in the macros
         if(cnorm > 0.000001 && useAllHistos) {

            cont = mapEnergy_E -> getBinContent(i,j) / cnorm ;
            econt = mapEnergy_E -> getBinError(i,j) / cnorm;
            mapEnergy_E -> setBinContent(i,j,cont);
            mapEnergy_E -> setBinError(i,j,econt);


            cont = mapEnergy_H -> getBinContent(i,j) / cnorm ;
            econt = mapEnergy_H -> getBinError(i,j) / cnorm;
            mapEnergy_H -> setBinContent(i,j,cont);
            mapEnergy_H -> setBinError(i,j,econt);

            cont = mapEnergy_EH -> getBinContent(i,j) / cnorm ;
            econt = mapEnergy_EH -> getBinError(i,j) / cnorm;
            mapEnergy_EH -> setBinContent(i,j,cont);
            mapEnergy_EH -> setBinError(i,j,econt);
         }


         // Occupancy (needed for occupancy vs ieta)
         cnorm   = occupancy_map -> getBinContent(i,j) / fev;
         enorm   = occupancy_map -> getBinError(i,j) / fev;
         if(cnorm > 1.e-30) occupancy_map -> setBinContent(i,j,cnorm);

         sumphi += cnorm;
         sumphie += enorm*enorm;

      } // end of iphy cycle (j)

      //Occupancy vs ieta histo is drawn
      // phi-factor evaluation for occupancy_vs_ieta calculation
      int ieta = i - 42;        // -41,-1,1,41 
      if(ieta >= -20 && ieta <= 20 )
         {phi_factor = 72.;}
         else {
            if(ieta >= 40 || ieta <= -40 ) {phi_factor = 18.;}
         else
            phi_factor = 36.;
       }

       cnorm = sumphi / phi_factor;
       enorm = sqrt(sumphie) / phi_factor;
       occupancy_vs_ieta->setBinContent(i, cnorm);
       occupancy_vs_ieta->setBinError(i,enorm);

   } // end of ieta cycle (i)
   
   return 1;

}

DEFINE_FWK_MODULE(CaloTowersClient);