#include <HcalSimHitsClient.h>

Inheritance diagram for HcalSimHitsClient:

Public Member Functions
virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &)

	HcalSimHitsClient (const edm::ParameterSet &)

virtual void	runClient_ (DQMStore::IBooker &, DQMStore::IGetter &)

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

virtual	~HcalSimHitsClient ()

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	outputFile_

bool	verbose_

Static Private Attributes
static const int	nTime = 4

static const int	nType = 25

static const int	nType1 = 4

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 35 of file HcalSimHitsClient.h.

Constructor & Destructor Documentation

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

explicit

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

HcalSimHitsClient::~HcalSimHitsClient ( )

virtual

Definition at line 25 of file HcalSimHitsClient.cc.

25 { }

Member Function Documentation

void HcalSimHitsClient::dqmEndJob	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

virtual

Implements DQMEDHarvester.

Definition at line 34 of file HcalSimHitsClient.cc.

References runClient_().

 {
   runClient_(ib,ig);
 }

void HcalSimHitsClient::runClient_	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

virtual

Definition at line 42 of file HcalSimHitsClient.cc.

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

Referenced by dqmEndJob().

                                                                            {
 
   ig.setCurrentFolder(dirName_);
   
   if (verbose_) std::cout << "\nrunClient" << std::endl; 
 
   std::vector<MonitorElement*> hcalMEs;
   
   std::vector<std::string> fullPathHLTFolders = ig.getSubdirs();
   for (unsigned int i=0;i<fullPathHLTFolders.size();i++) {
     if (verbose_) std::cout <<"\nfullPath: "<< fullPathHLTFolders[i] << std::endl;
     ig.setCurrentFolder(fullPathHLTFolders[i]);
 
     std::vector<std::string> fullSubPathHLTFolders = ig.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(), "HcalHitsV/SimHitsValidationHcal") == 0) {
     hcalMEs = ig.getContents(fullSubPathHLTFolders[j]);
     if (verbose_) std::cout <<"hltMES size : "<<hcalMEs.size()<<std::endl;
     if( !SimHitsEndjob(hcalMEs) ) std::cout<<"\nError in SimhitEndjob!"<<std::endl<<std::endl;
       }
 
     }    
 
   }
   
 }

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

Definition at line 74 of file HcalSimHitsClient.cc.

References generateEDF::cont, gather_cfg::cout, MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), i, j, relval_steps::k, mergeVDriftHistosByStation::name, nevent, nTime, nType, nType1, MonitorElement::setBinContent(), AlCaHLTBitMon_QueryRunRegistry::string, cond::rpcobgas::time, and verbose_.

Referenced by runClient_().

                                                                               {
   
   MonitorElement *Occupancy_map[nTime][nType];
   MonitorElement *Energy[nType1], *Time_weighteden[nType1];
   MonitorElement *HitEnergyvsieta[nType], *HitTimevsieta[nType];
   std::string divisions[nType]={"HB0","HB1","HE0+z","HE1+z","HE2+z","HE0-z","HE1-z",
                 "HE2-z","HO0","HFL0+z","HFS0+z","HFL1+z","HFS1+z",
                 "HFL2+z","HFS2+z","HFL3+z","HFS3+z","HFL0-z","HFS0-z",
                 "HFL1-z","HFS1-z","HFL2-z","HFS2-z","HFL3-z","HFS3-z"};
   
   std::string time[nTime]={"25","50","100","250"};
   std::string detdivision[nType1]={"HB","HE","HF","HO"};
   char name[40], name1[40], name2[40], name3[40], name4[40];
 
   for (int k=0; k<nType1;k++) {
     Energy[k] = 0;
     Time_weighteden[k] = 0;
     for (unsigned int ih=0; ih<hcalMEs.size(); ih++) {
       sprintf (name1, "Energy_%s", detdivision[k].c_str());
       sprintf (name2, "Time_Enweighted_%s", detdivision[k].c_str());
       if (strcmp(hcalMEs[ih]->getName().c_str(), name1) == 0) {
     Energy[k] = hcalMEs[ih];
       }
       if (strcmp(hcalMEs[ih]->getName().c_str(), name2) == 0) {
     Time_weighteden[k] = hcalMEs[ih];
       }
     }
   }
     
   for (int i=0; i<nTime; i++) {
     for (int j=0; j<nType;j++) {
       Occupancy_map[i][j]= 0;
       for (unsigned int ih=0; ih<hcalMEs.size(); ih++) {
     sprintf (name, "HcalHitE%s%s", time[i].c_str(),divisions[j].c_str());
     if (strcmp(hcalMEs[ih]->getName().c_str(), name) == 0) {
       Occupancy_map[i][j]= hcalMEs[ih];
     }
       }
     }
   }
 
   for (int k=0; k<nType;k++) {
     HitEnergyvsieta[k]= 0;
     HitTimevsieta[k]= 0;
     for (unsigned int ih=0; ih<hcalMEs.size(); ih++) {
       sprintf (name3, "HcalHitEta%s",divisions[k].c_str());
       sprintf (name4, "HcalHitTimeAEta%s",divisions[k].c_str());
       if (strcmp(hcalMEs[ih]->getName().c_str(), name3) == 0) {
     HitEnergyvsieta[k]= hcalMEs[ih];
       }
       if (strcmp(hcalMEs[ih]->getName().c_str(), name4) == 0) {
     HitTimevsieta[k]= hcalMEs[ih];
       }
     }
   }
     
   //mean energy 
  
   double nevent = Energy[0]->getEntries();
   if (verbose_) std::cout<<"nevent : "<<nevent<<std::endl;
 
   float cont[nTime][nType];
   float en[nType1], tme[nType1];
   float hitenergy[nType], hittime[nType];
   float fev = float(nevent);
   
   for (int dettype=0; dettype<nType1; dettype++) {
     int nx1=Energy[dettype]->getNbinsX();
     for (int i=0; i<=nx1; i++) {
       en[dettype]= Energy[dettype]->getBinContent(i)/fev;
       Energy[dettype]->setBinContent(i,en[dettype]);
     }
     int nx2= Time_weighteden[dettype]->getNbinsX();
     for (int i=0; i<=nx2; i++) {
       //std::cout<<" time_eneweighted for bin: "<< i<<"is:"<<Time_weighteden[dettype]->getBinContent(i)<<std::endl;
       tme[dettype]= Time_weighteden[dettype]->getBinContent(i)/fev;
       //std::cout<<" averagetime for bin : "<<i<<"is:"<<tme[dettype]<<std::endl;
       Time_weighteden[dettype]->setBinContent(i,tme[dettype]);
     }
   }
   
   for (int dettype=0; dettype<nType; dettype++)  {
     if (HitEnergyvsieta[dettype] != 0) {
       int nx1=HitEnergyvsieta[dettype]->getNbinsX();
       for (int i=0; i<=nx1; i++) {
     hitenergy[dettype]= HitEnergyvsieta[dettype]->getBinContent(i)/fev;
     HitEnergyvsieta[dettype]->setBinContent(i,hitenergy[dettype]);
       }
       int nx2= HitTimevsieta[dettype]->getNbinsX();
       for (int i=0; i<=nx2; i++) {
     hittime[dettype]= HitTimevsieta[dettype]->getBinContent(i)/fev;
     HitTimevsieta[dettype]->setBinContent(i,hittime[dettype]);
       }
     }
   }
   
   for (int itime=0; itime<nTime; itime++) {
     for (int det=0; det<nType;det++) {
       if (Occupancy_map[itime][det] != 0) {
     int ny= Occupancy_map[itime][det]->getNbinsY();
     int nx= Occupancy_map[itime][det]->getNbinsX(); 
     for (int i=1; i<nx+1; i++) {
       for (int j=1; j<ny+1; j++) {
         cont[itime][det] = Occupancy_map[itime][det]->getBinContent(i,j)/fev;
         Occupancy_map[itime][det]->setBinContent(i,j,cont[itime][det]);
       }
     }
       }
     }
   }
  
   return 1;
 }