#include <CaloTowersClient.h>

Inheritance diagram for CaloTowersClient:

Public Member Functions
virtual void	beginJob (void)

	CaloTowersClient (const edm::ParameterSet &)

int	CaloTowersEndjob (const std::vector< MonitorElement * > &hcalMEs)

virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

virtual	~CaloTowersClient ()

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 Attributes
edm::ParameterSet	conf_

bool	debug_

std::string	dirName_

std::string	dirNameJet_

std::string	dirNameMET_

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 37 of file CaloTowersClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file CaloTowersClient.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");
 }

CaloTowersClient::~CaloTowersClient ( )

virtual

Definition at line 22 of file CaloTowersClient.cc.

 { 
   
 }

Member Function Documentation

void CaloTowersClient::beginJob ( void )

virtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 27 of file CaloTowersClient.cc.

 {
  
 }

int CaloTowersClient::CaloTowersEndjob ( const std::vector< MonitorElement * > & hcalMEs )

Definition at line 69 of file CaloTowersClient.cc.

References generateEDF::cont, gather_cfg::cout, alignCSCRings::e, MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), i, j, nevent, and verbose_.

Referenced by dqmEndJob().

                                                                                {
 
    int useAllHistos = 0;
    MonitorElement* Ntowers_vs_ieta =0;
    MonitorElement* mapEnergy_N =0, *mapEnergy_E =0, *mapEnergy_H =0, *mapEnergy_EH =0;
    MonitorElement* occupancy_map =0, *occupancy_vs_ieta =0;
    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 fev = float(nevent);
 
    for (int i = 1; i <= nx; i++) {
       cont = Ntowers_vs_ieta -> getBinContent(i) / fev ;
       Ntowers_vs_ieta -> setBinContent(i,cont);
    }
 
    // mean energies & occupancies evaluation
 
    nx = mapEnergy_N->getNbinsX();
    int ny = mapEnergy_N->getNbinsY();
    float cnorm;
    float phi_factor;
 
    for (int i = 1; i <= nx; i++) {
       float sumphi = 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 ;
             mapEnergy_E -> setBinContent(i,j,cont);
 
             cont = mapEnergy_H -> getBinContent(i,j) / cnorm ;
             mapEnergy_H -> setBinContent(i,j,cont);
 
             cont = mapEnergy_EH -> getBinContent(i,j) / cnorm ;
             mapEnergy_EH -> setBinContent(i,j,cont);
          }
 
          // Occupancy (needed for occupancy vs ieta)
          cnorm   = occupancy_map -> getBinContent(i,j) / fev;
          if(cnorm > 1.e-30) occupancy_map -> setBinContent(i,j,cnorm);
 
          sumphi += cnorm;
 
       } // 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, 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.;
        }
        if(ieta >= 0) ieta -= 1; // -41 -1, 0 40  - to bring back to histo num
 
        cnorm = sumphi / phi_factor;
        occupancy_vs_ieta->Fill(double(ieta), cnorm);
 
    } // end of ieta cycle (i)
    
    return 1;
 
 }

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

overridevirtual

Implements DQMEDHarvester.

Definition at line 34 of file CaloTowersClient.cc.

References CaloTowersEndjob(), gather_cfg::cout, dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), 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., 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;
       }
 
     }    
 
   }
 
 }

Member Data Documentation

edm::ParameterSet CaloTowersClient::conf_

private

Definition at line 42 of file CaloTowersClient.h.

bool CaloTowersClient::debug_

private

Definition at line 45 of file CaloTowersClient.h.

Referenced by CaloTowersClient().

std::string CaloTowersClient::dirName_

private

Definition at line 47 of file CaloTowersClient.h.

Referenced by CaloTowersClient(), and dqmEndJob().

std::string CaloTowersClient::dirNameJet_

private

Definition at line 48 of file CaloTowersClient.h.

std::string CaloTowersClient::dirNameMET_

private

Definition at line 49 of file CaloTowersClient.h.

std::string CaloTowersClient::outputFile_

private

Definition at line 40 of file CaloTowersClient.h.

Referenced by CaloTowersClient().

bool CaloTowersClient::verbose_

private

Definition at line 44 of file CaloTowersClient.h.

Referenced by CaloTowersClient(), CaloTowersEndjob(), and dqmEndJob().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation