#include <CaloTowersClient.h>

Inheritance diagram for CaloTowersClient:

Public Member Functions
void	beginJob (void) override

	CaloTowersClient (const edm::ParameterSet &)

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

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

	~CaloTowersClient () override

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

	DQMEDHarvester ()

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

void	endJob () final

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) final

void	endRun (edm::Run const &, edm::EventSetup const &) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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 &	itemsToGetFrom (BranchType iType) 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) 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_

std::string	outputFile_

bool	verbose_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, ESInputTag const &tag)

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)

Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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 ( )

override

Definition at line 22 of file CaloTowersClient.cc.

 { 
   
 }

Member Function Documentation

void CaloTowersClient::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

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 gather_cfg::cout, DEFINE_FWK_MODULE, MillePedeFileConverter_cfg::e, objects.autophobj::float, MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), mps_fire::i, nevent, mathSSE::sqrt(), and verbose_.

Referenced by dqmEndJob().

                                                                                {
 
    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;
 
 }

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(), mps_fire::i, 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