HcalRecHitsDQMClient Class Reference

#include <HcalRecHitsDQMClient.h>

Inheritance diagram for HcalRecHitsDQMClient:

Public Member Functions
virtual void	beginJob (void) override
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override
virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
	HcalRecHitsDQMClient (const edm::ParameterSet &)
int	HcalRecHitsEndjob (const std::vector< MonitorElement * > &hcalMEs)
float	phifactor (int ieta)
virtual	~HcalRecHitsDQMClient ()
Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
	DQMEDHarvester (void)
virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	endJob () final
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
virtual void	endRun (edm::Run const &, edm::EventSetup const &)
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::ParameterSet	conf_
bool	debug_
std::string	dirName_
std::string	dirNameJet_
std::string	dirNameMET_
const HcalDDDRecConstants *	hcons
int	maxDepthAll_
int	maxDepthHB_
int	maxDepthHE_
int	maxDepthHF_
int	maxDepthHO_
int	nChannels_ [5]
std::string	outputFile_
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 49 of file HcalRecHitsDQMClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file HcalRecHitsDQMClient.cc.

References debug_, dirName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

                                                                        :conf_(iConfig)
{

  outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile", "myfile.root");
  debug_ = false;
  verbose_ = false;
  dirName_=iConfig.getParameter<std::string>("DQMDirName");
}

HcalRecHitsDQMClient::~HcalRecHitsDQMClient ( )

virtual

Definition at line 21 of file HcalRecHitsDQMClient.cc.

{ 
  
}

Member Function Documentation

void HcalRecHitsDQMClient::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 26 of file HcalRecHitsDQMClient.cc.

{
 

}

void HcalRecHitsDQMClient::beginRun ( edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Reimplemented from DQMEDHarvester.

Definition at line 33 of file HcalRecHitsDQMClient.cc.

References geometry, edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), CaloGeometry::getValidDetIds(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, maxDepthAll_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, and nChannels_.

                                                                             {
   
  edm::ESHandle<HcalDDDRecConstants> pHRNDC;
  es.get<HcalRecNumberingRecord>().get( pHRNDC );
  hcons = &(*pHRNDC);
  maxDepthHB_ = hcons->getMaxDepth(0);
  maxDepthHE_ = hcons->getMaxDepth(1);
  maxDepthHF_ = hcons->getMaxDepth(2);
  maxDepthHO_ = hcons->getMaxDepth(3);
 
  edm::ESHandle<CaloGeometry> geometry;

  es.get<CaloGeometryRecord > ().get(geometry);
 
  const std::vector<DetId>& hbCells = geometry->getValidDetIds(DetId::Hcal, HcalBarrel);
  const std::vector<DetId>& heCells = geometry->getValidDetIds(DetId::Hcal, HcalEndcap);
  const std::vector<DetId>& hoCells = geometry->getValidDetIds(DetId::Hcal, HcalOuter);
  const std::vector<DetId>& hfCells = geometry->getValidDetIds(DetId::Hcal, HcalForward);
 
  nChannels_[1] = hbCells.size(); 
  nChannels_[2] = heCells.size(); 
  nChannels_[3] = hoCells.size(); 
  nChannels_[4] = hfCells.size();
  nChannels_[0] = nChannels_[1] + nChannels_[2] + nChannels_[3] + nChannels_[4];
 
  //std::cout << "Channels HB:" << nChannels_[1] << " HE:" << nChannels_[2] << " HO:" << nChannels_[3] << " HF:" << nChannels_[4] << std::endl;
 
  //We hardcode the HF depths because in the dual readout configuration, rechits are not defined for depths 3&4
  maxDepthHF_ = (maxDepthHF_ > 2 ? 2 : maxDepthHF_); //We reatin the dynamic possibility that HF might have 0 or 1 depths
 
  maxDepthAll_ = ( maxDepthHB_ + maxDepthHO_ > maxDepthHE_ ? maxDepthHB_ + maxDepthHO_ : maxDepthHE_ );
  maxDepthAll_ = ( maxDepthAll_ > maxDepthHF_ ? maxDepthAll_ : maxDepthHF_ );
 
}

void HcalRecHitsDQMClient::dqmEndJob ( DQMStore::IBooker &  ibooker, DQMStore::IGetter &  igetter  )

overridevirtual

Implements DQMEDHarvester.

Definition at line 68 of file HcalRecHitsDQMClient.cc.

References gather_cfg::cout, dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), HcalRecHitsEndjob(), i, j, DQMStore::IGetter::setCurrentFolder(), and verbose_.

{
  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., CaloTowersD/CaloTowersTask, HcalRecHitsD/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(), "HcalRecHitsD/HcalRecHitTask") ==0  ){
         hcalMEs = igetter.getContents(fullSubPathHLTFolders[j]);
         if (verbose_) std::cout <<"hltMES size : "<<hcalMEs.size()<<std::endl;
         if( !HcalRecHitsEndjob(hcalMEs) ) std::cout<<"\nError in HcalRecHitsEndjob!"<<std::endl<<std::endl;
      }

    }    

  }

}

int HcalRecHitsDQMClient::HcalRecHitsEndjob

(

const std::vector< MonitorElement * > &

hcalMEs

)

Definition at line 104 of file HcalRecHitsDQMClient.cc.

References ws_sso_content_reader::content, gather_cfg::cout, spr::find(), objects.autophobj::float, MonitorElement::getEntries(), i, j, nChannels_, phifactor(), tablePrinter::prefix, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

Referenced by dqmEndJob().

                                                                                    {

   MonitorElement* Nhf=0;

   //Search for emap histograms, and collect them into this vector
   //All subdtectors are plotted together in these histograms. We only need to look for different depths
   std::vector<MonitorElement*> emap_depths;

   //This vector is filled occupancy_maps identified by both subdetector and depth
   std::vector<MonitorElement*> occupancy_maps;
   std::vector<std::string> occupancyID;

   //This vector is filled with emean_vs_ieta histograms, they are divided by both subdetector and depth
   std::vector<MonitorElement*> emean_vs_ieta;

   //These are the only histograms filled in this module; however, the histograms are created empty in HcalRecHitsAnalyzer
   //occupancy_vs_ieta, divided by both subdetector and depth
   std::vector<MonitorElement*> occupancy_vs_ieta;
   std::vector<std::string> occupancy_vs_ietaID;

   //RecHit_StatusWord & RecHit_Aux_StatusWord
   //Divided by subdectector
   std::vector<MonitorElement*> RecHit_StatusWord;
   std::vector<float>           RecHit_StatusWord_Channels;
   std::vector<MonitorElement*> RecHit_Aux_StatusWord;
   std::vector<float>           RecHit_Aux_StatusWord_Channels;

   for(unsigned int ih=0; ih<hcalMEs.size(); ih++){

      //N_HF is not special, it is just convient to get the total number of events
      //The number of entries in N_HF is equal to the number of events
      if( hcalMEs[ih]->getName() == "N_HF" ){ 
         Nhf= hcalMEs[ih];
         continue;
      }

      // ***********************
      // * We fill the various MonitorElement vectors by searching for a matching substring
      // * The methods that are used are agnostic to the ordering of vectors
      // ***********************     

      if( hcalMEs[ih]->getName().find("emap_depth") != std::string::npos ){
         emap_depths.push_back(hcalMEs[ih]);
         continue;  
      } 

      if( hcalMEs[ih]->getName().find("occupancy_map_H") != std::string::npos ){
         occupancy_maps.push_back(hcalMEs[ih]);

         // Use occupancyID to save the subdetector and depth information
         // This will help preserve both indifference to vector ordering and specific details of the detector topology
         // The position in occupancyID must correspond to the histogram position in occupancy_maps

         // Save the string after "occupancy_map_"
         
         std::string prefix = "occupancy_map_";

         occupancyID.push_back( hcalMEs[ih]->getName().substr(prefix.size()) );

         continue;
      }

      if( hcalMEs[ih]->getName().find("emean_vs_ieta_H") != std::string::npos ){
         emean_vs_ieta.push_back(hcalMEs[ih]);
         continue;
      }

      if( hcalMEs[ih]->getName().find("occupancy_vs_ieta_H") != std::string::npos ){
         occupancy_vs_ieta.push_back(hcalMEs[ih]);

         // Use occupancy_vs_ietaID to save the subdetector and depth information
         // This will help preserve both indifference to vector ordering and specific details of the detector topology
         // The position in occupancyID must correspond to the histogram position in occupancy_vs_ieta

         // Save the string after "occupancy_vs_ieta_"
         
         std::string prefix = "occupancy_vs_ieta_";

         occupancy_vs_ietaID.push_back( hcalMEs[ih]->getName().substr(prefix.size()) );

         continue;
      }

      if( hcalMEs[ih]->getName().find("HcalRecHitTask_RecHit_StatusWord_H") != std::string::npos ){
         RecHit_StatusWord.push_back(hcalMEs[ih]);

         if(hcalMEs[ih]->getName().find("HB") != std::string::npos ){
            RecHit_StatusWord_Channels.push_back((float)nChannels_[1]);
         }else if(hcalMEs[ih]->getName().find("HE") != std::string::npos ){
            RecHit_StatusWord_Channels.push_back((float)nChannels_[2]);
         }else if(hcalMEs[ih]->getName().find("H0") != std::string::npos ){
            RecHit_StatusWord_Channels.push_back((float)nChannels_[3]);
         }else if(hcalMEs[ih]->getName().find("HF") != std::string::npos ){
            RecHit_StatusWord_Channels.push_back((float)nChannels_[4]);
         } else {
            RecHit_StatusWord_Channels.push_back(1.);
         }

         continue;
      }

      if( hcalMEs[ih]->getName().find("HcalRecHitTask_RecHit_Aux_StatusWord_H") != std::string::npos ){
         RecHit_Aux_StatusWord.push_back(hcalMEs[ih]);

         if(hcalMEs[ih]->getName().find("HB") != std::string::npos ){
            RecHit_Aux_StatusWord_Channels.push_back((float)nChannels_[1]);
         }else if(hcalMEs[ih]->getName().find("HE") != std::string::npos ){
            RecHit_Aux_StatusWord_Channels.push_back((float)nChannels_[2]);
         }else if(hcalMEs[ih]->getName().find("H0") != std::string::npos ){
            RecHit_Aux_StatusWord_Channels.push_back((float)nChannels_[3]);
         }else if(hcalMEs[ih]->getName().find("HF") != std::string::npos ){
            RecHit_Aux_StatusWord_Channels.push_back((float)nChannels_[4]);
         } else {
            RecHit_Aux_StatusWord_Channels.push_back(1.);
         }

         continue;
      }

   } 

   // mean energies and occupancies evaluation

   double nevtot = Nhf->getEntries(); // Use the number of entries in the Nhf histogram to give the total number of events

   if(verbose_) std::cout<<"nevtot : "<<nevtot<<std::endl;

   // emap histograms are scaled by the number of events
   float  fev           = float (nevtot);
   double scaleBynevtot = 1 / fev;

   // In this and the following histogram vectors, recognize that the for-loop index
   // does not have to correspond to any particular depth
   for(unsigned int depthIdx = 0; depthIdx < emap_depths.size(); depthIdx++){

      int nx = emap_depths[depthIdx]->getNbinsX();
      int ny = emap_depths[depthIdx]->getNbinsY();

      float cnorm;
      float enorm;

      for (int i = 1; i <= nx; i++) {      
         for (int j = 1; j <= ny; j++) {
        cnorm = emap_depths[depthIdx]->getBinContent(i,j) * scaleBynevtot;
        enorm = emap_depths[depthIdx]->getBinError(i,j) * scaleBynevtot;
            emap_depths[depthIdx]->setBinContent(i,j,cnorm);
            emap_depths[depthIdx]->setBinError(i,j,enorm);

         }
      }
   }

   // occupancy_maps & matched occupancy_vs_ieta

   bool omatched = false;

   for(unsigned int occupancyIdx = 0; occupancyIdx < occupancy_maps.size(); occupancyIdx++){

      int nx = occupancy_maps[occupancyIdx]->getNbinsX();
      int ny = occupancy_maps[occupancyIdx]->getNbinsY();

      float cnorm;
      float enorm;

      unsigned int vsIetaIdx = occupancy_vs_ieta.size();
      omatched = false;

      for(vsIetaIdx = 0; vsIetaIdx < occupancy_vs_ieta.size(); vsIetaIdx++){
         if(occupancyID[occupancyIdx] == occupancy_vs_ietaID[vsIetaIdx]){
            omatched = true;
            break;
         }
      }// match occupancy_vs_ieta histogram

      for (int i = 1; i <= nx; i++) {      
         for (int j = 1; j <= ny; j++) {
        cnorm = occupancy_maps[occupancyIdx]->getBinContent(i,j) * scaleBynevtot;
        enorm = occupancy_maps[occupancyIdx]->getBinError(i,j) * scaleBynevtot;
            occupancy_maps[occupancyIdx]->setBinContent(i,j,cnorm);
            occupancy_maps[occupancyIdx]->setBinError(i,j,enorm);

         }
      }

      //Fill occupancy_vs_ieta

      if(omatched){

         //We run over all of the ieta values
         for (int ieta = -41; ieta <= 41; ieta++) {
            float phi_factor = 1.;
            float sumphi = 0.;

            if(ieta == 0) continue; //ieta=0 is not defined

            phi_factor = phifactor(ieta);

            //the rechits occupancy map defines iphi as 0..71
            for (int iphi = 0; iphi <= 71; iphi++) {
               int binIeta = occupancy_maps[occupancyIdx]->getTH2F()->GetXaxis()->FindBin(float(ieta));
               int binIphi = occupancy_maps[occupancyIdx]->getTH2F()->GetYaxis()->FindBin(float(iphi));

               float content = occupancy_maps[occupancyIdx]->getBinContent(binIeta,binIphi);

               sumphi += content;
            }//for loop over phi

            double deta = double(ieta);

            // fill occupancies vs ieta
            cnorm = sumphi / phi_factor;
            occupancy_vs_ieta[vsIetaIdx]->Fill(deta, cnorm);

         }//Fill occupancy_vs_ieta
      }//if omatched
   }

   // Status Word
   // Normalized by number of events and by number of channels per subdetector as well

   for(unsigned int StatusWordIdx = 0; StatusWordIdx < RecHit_StatusWord.size(); StatusWordIdx++){
         
      int nx = RecHit_StatusWord[StatusWordIdx]->getNbinsX();

      float cnorm;
      float enorm;

      for (int i = 1; i <= nx; i++) {      
         cnorm = RecHit_StatusWord[StatusWordIdx]->getBinContent(i) * scaleBynevtot / RecHit_StatusWord_Channels[StatusWordIdx];
         enorm = RecHit_StatusWord[StatusWordIdx]->getBinError(i) * scaleBynevtot / RecHit_StatusWord_Channels[StatusWordIdx];
         RecHit_StatusWord[StatusWordIdx]->setBinContent(i,cnorm);
         RecHit_StatusWord[StatusWordIdx]->setBinError(i,enorm);

      }

   }

   for(unsigned int AuxStatusWordIdx = 0; AuxStatusWordIdx < RecHit_Aux_StatusWord.size(); AuxStatusWordIdx++){
         
      int nx = RecHit_Aux_StatusWord[AuxStatusWordIdx]->getNbinsX();

      float cnorm;
      float enorm;

      for (int i = 1; i <= nx; i++) {      
         cnorm = RecHit_Aux_StatusWord[AuxStatusWordIdx]->getBinContent(i) * scaleBynevtot / RecHit_Aux_StatusWord_Channels[AuxStatusWordIdx];
         enorm = RecHit_Aux_StatusWord[AuxStatusWordIdx]->getBinError(i) * scaleBynevtot / RecHit_Aux_StatusWord_Channels[AuxStatusWordIdx];
         RecHit_Aux_StatusWord[AuxStatusWordIdx]->setBinContent(i,cnorm);
         RecHit_Aux_StatusWord[AuxStatusWordIdx]->setBinError(i,enorm);

      }


   }

   return 1;
}

float HcalRecHitsDQMClient::phifactor

(

int

ieta

)

Definition at line 362 of file HcalRecHitsDQMClient.cc.

References DEFINE_FWK_MODULE.

Referenced by HcalRecHitsEndjob().

                                             {

   float phi_factor_;

   if(ieta >= -20 && ieta <= 20 ){
      phi_factor_ = 72.;
   } else {
      if(ieta >= 40 || ieta <= -40 ){
         phi_factor_ = 18.;
      } else {
         phi_factor_ = 36.;
      }
   }

   return phi_factor_;

}

Member Data Documentation

edm::ParameterSet HcalRecHitsDQMClient::conf_

private

Definition at line 53 of file HcalRecHitsDQMClient.h.

bool HcalRecHitsDQMClient::debug_

private

Definition at line 56 of file HcalRecHitsDQMClient.h.

Referenced by HcalRecHitsDQMClient().

std::string HcalRecHitsDQMClient::dirName_

private

Definition at line 58 of file HcalRecHitsDQMClient.h.

Referenced by dqmEndJob(), and HcalRecHitsDQMClient().

std::string HcalRecHitsDQMClient::dirNameJet_

private

Definition at line 59 of file HcalRecHitsDQMClient.h.

std::string HcalRecHitsDQMClient::dirNameMET_

private

Definition at line 60 of file HcalRecHitsDQMClient.h.

const HcalDDDRecConstants* HcalRecHitsDQMClient::hcons

private

Definition at line 62 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::maxDepthAll_

private

Definition at line 63 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::maxDepthHB_

private

Definition at line 63 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::maxDepthHE_

private

Definition at line 63 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::maxDepthHF_

private

Definition at line 63 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::maxDepthHO_

private

Definition at line 63 of file HcalRecHitsDQMClient.h.

Referenced by beginRun().

int HcalRecHitsDQMClient::nChannels_[5]

private

Definition at line 65 of file HcalRecHitsDQMClient.h.

Referenced by beginRun(), and HcalRecHitsEndjob().

std::string HcalRecHitsDQMClient::outputFile_

private

Definition at line 52 of file HcalRecHitsDQMClient.h.

Referenced by HcalRecHitsDQMClient().

bool HcalRecHitsDQMClient::verbose_

private

Definition at line 55 of file HcalRecHitsDQMClient.h.

Referenced by dqmEndJob(), HcalRecHitsDQMClient(), and HcalRecHitsEndjob().