HcalDigisClient Class Reference

#include <Validation/HcalDigisClient/src/HcalDigisClient.cc>

Inheritance diagram for HcalDigisClient:

Classes
struct	HistLim

Public Member Functions
	HcalDigisClient (const edm::ParameterSet &)
	~HcalDigisClient ()
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
virtual void	beginRun (edm::Run const &, edm::EventSetup const &)
	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 ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	book1D (DQMStore::IBooker &ib, std::string name, int n, double min, double max)
void	book1D (DQMStore::IBooker &ib, std::string name, const HistLim &limX)
void	book2D (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
void	booking (DQMStore::IBooker &ib, std::string subdetopt)
void	bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
void	bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY, const char *option)
virtual void	dqmEndJob (DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter)
void	fill1D (std::string name, double X, double weight=1)
void	fill2D (std::string name, double X, double Y, double weight=1)
void	fillPf (std::string name, double X, double Y)
int	HcalDigisEndjob (const std::vector< MonitorElement * > &hcalMEs, std::string subdet_)
double	integralMETH2D (MonitorElement *ME, int i0, int i1, int j0, int j1)
MonitorElement *	monitor (std::string name)
virtual void	runClient (DQMStore::IBooker &ib, DQMStore::IGetter &ig)
void	scaleMETH2D (MonitorElement *ME, double s)
std::string	str (int x)

Private Attributes
std::string	dirName_
std::map< std::string, MonitorElement * > *	msm_
std::string	outputFile_

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 27 of file HcalDigisClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 of file HcalDigisClient.cc.

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

                                                               {
    outputFile_ = iConfig.getUntrackedParameter<std::string > ("outputFile", "HcalDigisClient.root");
    dirName_ = iConfig.getParameter<std::string > ("DQMDirName");
    msm_ = new std::map<std::string, MonitorElement*>();


}

HcalDigisClient::~HcalDigisClient

(

)

Definition at line 32 of file HcalDigisClient.cc.

References msm_.

                                  {
    delete msm_;
}

Member Function Documentation

void HcalDigisClient::book1D ( DQMStore::IBooker &  ib, std::string  name, int  n, double  min, double  max  )

inlineprivate

Definition at line 63 of file HcalDigisClient.h.

References DQMStore::IBooker::book1D(), bookConverter::max, min(), msm_, gen::n, and mergeVDriftHistosByStation::name.

Referenced by booking().

                                                                                    {
        if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), n, min, max);
    }

void HcalDigisClient::book1D ( DQMStore::IBooker &  ib, std::string  name, const HistLim &  limX  )

inlineprivate

Definition at line 67 of file HcalDigisClient.h.

References DQMStore::IBooker::book1D(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and mergeVDriftHistosByStation::name.

                                                                          {
        if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max);
    }

void HcalDigisClient::book2D ( DQMStore::IBooker &  ib, std::string  name, const HistLim &  limX, const HistLim &  limY  )

inlineprivate

Definition at line 75 of file HcalDigisClient.h.

References DQMStore::IBooker::book2D(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and mergeVDriftHistosByStation::name.

                                                                                               {
        if (!msm_->count(name)) (*msm_)[name] = ib.book2D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
    }

void HcalDigisClient::booking ( DQMStore::IBooker &  ib, std::string  subdetopt  )

private

Definition at line 37 of file HcalDigisClient.cc.

References book1D(), HLT_25ns14e33_v1_cff::depth, str(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dqmEndJob().

                                                                      {

    std::string strtmp;
    HistLim ietaLim(85, -42.5, 42.5);

    for (int depth = 1; depth <= 4; depth++) {
        strtmp = "HcalDigiTask_occupancy_vs_ieta_depth" + str(depth) + "_" + subdetopt;
        book1D(ib,strtmp, ietaLim);
    }
    
}

void HcalDigisClient::bookPf ( DQMStore::IBooker &  ib, std::string  name, const HistLim &  limX, const HistLim &  limY  )

inlineprivate

Definition at line 83 of file HcalDigisClient.h.

References DQMStore::IBooker::bookProfile(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and mergeVDriftHistosByStation::name.

                                                                                               {
        if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
    }

void HcalDigisClient::bookPf ( DQMStore::IBooker &  ib, std::string  name, const HistLim &  limX, const HistLim &  limY, const char *  option  )

inlineprivate

Definition at line 87 of file HcalDigisClient.h.

References DQMStore::IBooker::bookProfile(), HcalDigisClient::HistLim::max, HcalDigisClient::HistLim::min, msm_, HcalDigisClient::HistLim::n, and mergeVDriftHistosByStation::name.

                                                                                                                   {
        if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max, option);
    }

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

inlineprivatevirtual

Implements DQMEDHarvester.

Definition at line 35 of file HcalDigisClient.h.

References booking(), dirName_, runClient(), and DQMStore::IGetter::setCurrentFolder().

                                                                                 {
      igetter.setCurrentFolder("HcalDigisV/HcalDigiTask"); // moved this line from constructor

      // the following booking clas were moved from the constructor
      booking(ibooker, "HB");
      booking(ibooker, "HE");
      booking(ibooker, "HO");
      booking(ibooker, "HF");

      igetter.setCurrentFolder(dirName_); // This sets the DQMStore (should apply to ibooker as well
      runClient(ibooker, igetter);
    }

void HcalDigisClient::fill1D ( std::string  name, double  X, double  weight = 1  )

inlineprivate

Definition at line 71 of file HcalDigisClient.h.

References msm_, and puppiForMET_cff::weight.

Referenced by HcalDigisEndjob().

                                                             {
        msm_->find(name)->second->Fill(X, weight);
    }

void HcalDigisClient::fill2D ( std::string  name, double  X, double  Y, double  weight = 1  )

inlineprivate

Definition at line 79 of file HcalDigisClient.h.

References msm_, and puppiForMET_cff::weight.

                                                                       {
        msm_->find(name)->second->Fill(X, Y, weight);
    }

void HcalDigisClient::fillPf ( std::string  name, double  X, double  Y  )

inlineprivate

Definition at line 91 of file HcalDigisClient.h.

References msm_.

                                                    {
        msm_->find(name)->second->Fill(X, Y);
    }

int HcalDigisClient::HcalDigisEndjob ( const std::vector< MonitorElement * > &  hcalMEs, std::string  subdet_  )

private

Definition at line 73 of file HcalDigisClient.cc.

References ev, fill1D(), MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), MonitorElement::getTH1(), i, j, and MonitorElement::setBinContent().

Referenced by runClient().

                                                                                                 {

    using namespace std;
    string strtmp;


    MonitorElement * nevtot(0);
    MonitorElement * ieta_iphi_occupancy_map1(0);
    MonitorElement * ieta_iphi_occupancy_map2(0);
    MonitorElement * ieta_iphi_occupancy_map3(0);
    MonitorElement * ieta_iphi_occupancy_map4(0);


    // std::cout << " Number of histos " <<     hcalMEs.size() << std::endl;

    for (unsigned int ih = 0; ih < hcalMEs.size(); ih++) {
      
         if (hcalMEs[ih]->getName() == "nevtot") nevtot = hcalMEs[ih];

         strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth1_" + subdet_;
         if (hcalMEs[ih]->getName() == strtmp) ieta_iphi_occupancy_map1 = hcalMEs[ih];
         strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth2_" + subdet_;
         if (hcalMEs[ih]->getName() == strtmp) ieta_iphi_occupancy_map2 = hcalMEs[ih];
         strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth3_" + subdet_;
         if (hcalMEs[ih]->getName() == strtmp) ieta_iphi_occupancy_map3 = hcalMEs[ih];
         strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth4_" + subdet_;
         if (hcalMEs[ih]->getName() == strtmp) ieta_iphi_occupancy_map4 = hcalMEs[ih];

    }

    if (nevtot                   == 0 ||
    ieta_iphi_occupancy_map1 == 0 ||
    ieta_iphi_occupancy_map2 == 0 ||
    ieta_iphi_occupancy_map3 == 0 ||
    ieta_iphi_occupancy_map4 == 0   
    ) {
      edm::LogError("HcalDigisClient") << "No nevtot or maps histo found..."; 
      return 0;
    }

    int ev = nevtot->getEntries();
    if(ev <= 0) {
      edm::LogError("HcalDigisClient") << "normalization factor <= 0!"; 
      return 0;
    }

    float fev = (float) nevtot->getEntries();

    float sumphi_1, sumphi_2, sumphi_3, sumphi_4;
    float phi_factor;
    float cnorm;

    int nx = ieta_iphi_occupancy_map1->getNbinsX();
    int ny = ieta_iphi_occupancy_map1->getNbinsY();
    
    for (int i = 1; i <= nx; i++) {
        for (int j = 1; j <= ny; j++) {

            // occupancies

            strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth1_" + subdet_;
            cnorm = ieta_iphi_occupancy_map1->getBinContent(i, j) / fev;
            ieta_iphi_occupancy_map1->setBinContent(i, j, cnorm);

            strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth2_" + subdet_;
            cnorm = ieta_iphi_occupancy_map2->getBinContent(i, j) / fev;
            ieta_iphi_occupancy_map2->setBinContent(i, j, cnorm);

            strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth3_" + subdet_;
            cnorm = ieta_iphi_occupancy_map3->getBinContent(i, j) / fev; 
            ieta_iphi_occupancy_map3->setBinContent(i, j, cnorm);

            strtmp = "HcalDigiTask_ieta_iphi_occupancy_map_depth4_" + subdet_;
            cnorm = ieta_iphi_occupancy_map4->getBinContent(i, j) / fev; 
            ieta_iphi_occupancy_map4->setBinContent(i, j, cnorm);

        }
    }       

    for (int i = 1; i <= 82; i++) {

        int ieta = i - 42; // -41 -1, 0 40
        if (ieta >= 0) ieta += 1; // -41 -1, 1 41  - to make it detector-like

        if (ieta >= -20 && ieta <= 20) {
          phi_factor = 72.; 
        } else {
          if (ieta >= 40 || ieta <= -40) 
              phi_factor = 18.; 
           else 
              phi_factor = 36.; 
        }    
    
        sumphi_1 = 0.;
        sumphi_2 = 0.;
        sumphi_3 = 0.;
        sumphi_4 = 0.;

        for (int iphi = 1; iphi <= 72; iphi++) {
            sumphi_1 += ieta_iphi_occupancy_map1->getTH1()->GetBinContent(ieta_iphi_occupancy_map1->getTH1()->FindFixBin(double(ieta),double(iphi)));
            sumphi_2 += ieta_iphi_occupancy_map2->getTH1()->GetBinContent(ieta_iphi_occupancy_map2->getTH1()->FindFixBin(double(ieta),double(iphi)));
            sumphi_3 += ieta_iphi_occupancy_map3->getTH1()->GetBinContent(ieta_iphi_occupancy_map3->getTH1()->FindFixBin(double(ieta),double(iphi)));
            sumphi_4 += ieta_iphi_occupancy_map4->getTH1()->GetBinContent(ieta_iphi_occupancy_map4->getTH1()->FindFixBin(double(ieta),double(iphi)));
        }
 
        //REMOVED (JRD) if (ieta >= 0) ieta -= 1; // -41 -1, 0 40  - to bring back to strtmp num !!! 
        double deta = double(ieta);

        // occupancies vs ieta
        cnorm = sumphi_1 / phi_factor;
        strtmp = "HcalDigiTask_occupancy_vs_ieta_depth1_" + subdet_;
        fill1D(strtmp, deta, cnorm);

        cnorm = sumphi_2 / phi_factor;
        strtmp = "HcalDigiTask_occupancy_vs_ieta_depth2_" + subdet_;
        fill1D(strtmp, deta, cnorm);

        cnorm = sumphi_3 / phi_factor;
        strtmp = "HcalDigiTask_occupancy_vs_ieta_depth3_" + subdet_;
        fill1D(strtmp, deta, cnorm);

        cnorm = sumphi_4 / phi_factor;
        strtmp = "HcalDigiTask_occupancy_vs_ieta_depth4_" + subdet_;
        fill1D(strtmp, deta, cnorm);

    } // end of i-loop

  return 1;
}

double HcalDigisClient::integralMETH2D ( MonitorElement *  ME, int  i0, int  i1, int  j0, int  j1  )

private

Definition at line 214 of file HcalDigisClient.cc.

References MonitorElement::getBinContent(), i, and j.

                                                                                         {
    double sum(0);
    for (int i = i0; i <= i1; i++) {
        for (int j = j0; j <= j1; j++) {
            sum += ME->getBinContent(i, j);
        }
    }

    return sum;
}

MonitorElement * HcalDigisClient::monitor ( std::string  name )

private

Definition at line 203 of file HcalDigisClient.cc.

References msm_, and NULL.

                                                       {
    if (!msm_->count(name)) return NULL;
    else return msm_->find(name)->second;
}

void HcalDigisClient::runClient ( DQMStore::IBooker &  ib, DQMStore::IGetter &  ig  )

privatevirtual

Definition at line 49 of file HcalDigisClient.cc.

References dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), HcalDigisEndjob(), i, j, and DQMStore::IGetter::setCurrentFolder().

Referenced by dqmEndJob().

                                                                        {
    ig.setCurrentFolder(dirName_);
    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 = ig.getSubdirs();
    for (unsigned int i = 0; i < fullPathHLTFolders.size(); i++) {
        ig.setCurrentFolder(fullPathHLTFolders[i]);
        std::vector<std::string> fullSubPathHLTFolders = ig.getSubdirs();
        for (unsigned int j = 0; j < fullSubPathHLTFolders.size(); j++) {
            if (strcmp(fullSubPathHLTFolders[j].c_str(), "HcalDigisV/HcalDigiTask") == 0) {
                hcalMEs = ig.getContents(fullSubPathHLTFolders[j]);
                if (!HcalDigisEndjob(hcalMEs, "HB")) 
          edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HB"; 
                if (!HcalDigisEndjob(hcalMEs, "HE")) 
          edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HE"; 
                if (!HcalDigisEndjob(hcalMEs, "HO")) 
          edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HO"; 
                if (!HcalDigisEndjob(hcalMEs, "HF")) 
          edm::LogError("HcalDigisClient") << "Error in HcalDigisEndjob! HF";             }
        }
    }
}

void HcalDigisClient::scaleMETH2D ( MonitorElement *  ME, double  s  )

private

Definition at line 225 of file HcalDigisClient.cc.

References ws_sso_content_reader::content, MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), i, j, alignCSCRings::s, and MonitorElement::setBinContent().

                                                              {
    int nx = ME->getNbinsX();
    int ny = ME->getNbinsY();

    double content(0);
    for (int i = 1; i <= nx; i++) {
        for (int j = 1; j <= ny; j++) {
            content = ME->getBinContent(i, j);
            content *= s;
            ME->setBinContent(i, j, content);
        }
    }
}

std::string HcalDigisClient::str ( int  x )

private

Definition at line 208 of file HcalDigisClient.cc.

References dbtoconf::out, and x.

Referenced by booking().

                                    {
    std::stringstream out;
    out << x;
    return out.str();
}

Member Data Documentation

std::string HcalDigisClient::dirName_

private

Definition at line 103 of file HcalDigisClient.h.

Referenced by dqmEndJob(), HcalDigisClient(), and runClient().

std::map<std::string, MonitorElement*>* HcalDigisClient::msm_

private

Definition at line 101 of file HcalDigisClient.h.

Referenced by book1D(), book2D(), bookPf(), fill1D(), fill2D(), fillPf(), HcalDigisClient(), monitor(), and ~HcalDigisClient().

std::string HcalDigisClient::outputFile_

private

Definition at line 102 of file HcalDigisClient.h.

Referenced by HcalDigisClient().