#include <PFClient.h>

Inheritance diagram for PFClient:

Public Member Functions
	PFClient (const edm::ParameterSet &parameterSet)

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginJob ()

void	createEfficiencyPlots (std::string &folder, std::string &name)

void	createProfilePlots (std::string &folder, std::string &name)

void	createProjectionPlots (std::string &folder, std::string &name)

void	createResolutionPlots (std::string &folder, std::string &name)

void	doEfficiency ()

void	doProfiles ()

void	doProjection ()

void	doSummaries ()

void	endJob ()

void	endRun (edm::Run const &run, edm::EventSetup const &eSetup)

void	getHistogramParameters (MonitorElement *me_slice, double &avarage, double &rms, double &mean, double &sigma)

Private Attributes
DQMStore *	dqmStore_

std::vector< std::string >	effHistogramNames_

bool	efficiencyFlag_

std::vector< std::string >	folderNames_

std::vector< std::string >	histogramNames_

bool	profileFlag_

std::vector< std::string >	profileHistogramNames_

std::vector< std::string >	projectionHistogramNames_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 12 of file PFClient.h.

Constructor & Destructor Documentation

PFClient::PFClient ( const edm::ParameterSet & parameterSet )

Definition at line 16 of file PFClient.cc.

References effHistogramNames_, efficiencyFlag_, folderNames_, edm::ParameterSet::getParameter(), histogramNames_, profileFlag_, profileHistogramNames_, and projectionHistogramNames_.

 {
   folderNames_ = parameterSet.getParameter< std::vector<std::string> >( "FolderNames" );
   histogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNames" );
   efficiencyFlag_ =  parameterSet.getParameter< bool> ("CreateEfficiencyPlots" );
   effHistogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNamesForEfficiencyPlots" );
   projectionHistogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNamesForProjectionPlots" );
   profileFlag_ =  parameterSet.getParameter< bool> ("CreateProfilePlots" );
   profileHistogramNames_  = parameterSet.getParameter< std::vector<std::string> >( "HistogramNamesForProfilePlots" );
 }

Member Function Documentation

void PFClient::analyze	(	edm::Event const &	,
		edm::EventSetup const &
	)

inlineprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 19 of file PFClient.h.

19 {;}

void PFClient::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 29 of file PFClient.cc.

References dqmStore_, and cppFunctionSkipper::operator.

                         {
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 }

void PFClient::createEfficiencyPlots	(	std::string &	folder,
		std::string &	name
	)

private

Definition at line 303 of file PFClient.cc.

References DQMStore::book1D(), MonitorElement::DQM_KIND_TH1F, dqmStore_, DQMStore::get(), MonitorElement::getAxisTitle(), MonitorElement::getTH1F(), MonitorElement::kind(), MonitorElement::Reset(), DQMStore::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripMonitorClusterAlca_cfi::xmax, and SiStripMonitorClusterAlca_cfi::xmin.

Referenced by doEfficiency().

                                                                        {     
   MonitorElement* me1 = dqmStore_->get(folder+"/"+name+"ref_");
   MonitorElement* me2 = dqmStore_->get(folder+"/"+name+"gen_");
   if (!me1 || !me2) return;
 
   TH1F* me1_forEff = (TH1F*)me1->getTH1F()->Rebin(2,"me1_forEff");
   TH1F* me2_forEff = (TH1F*)me2->getTH1F()->Rebin(2,"me2_forEff");
     
   MonitorElement* me_eff;
   if ( (me1->kind() == MonitorElement::DQM_KIND_TH1F) &&
        (me1->kind() == MonitorElement::DQM_KIND_TH1F) ) {
     //TH1* th1 = me1->getTH1F();
     //size_t nbinx = me1->getNbinsX();
     size_t nbinx = me1_forEff->GetNbinsX();
 
     //float xmin = th1->GetXaxis()->GetXmin();
     //float xmax = th1->GetXaxis()->GetXmax();
     float xmin = me1_forEff->GetXaxis()->GetXmin();
     float xmax = me1_forEff->GetXaxis()->GetXmax();
     std::string xtit = me1->getAxisTitle(1);
     std::string tit_new;
     tit_new = ";"+xtit+";"+xtit+" efficiency"; 
 
     dqmStore_->setCurrentFolder(folder);
     me_eff = dqmStore_->book1D("efficiency_"+name,tit_new, nbinx, xmin, xmax); 
         
     me_eff->Reset(); me_eff->setEfficiencyFlag();
     /*
     double  efficiency;
     for (size_t ix = 1; ix < nbinx+1; ++ix) {
       float val1 = me1->getBinContent(ix);
       float val2 = me2->getBinContent(ix);
       if (val2 > 0.0) efficiency = val1/val2;
       else efficiency = 0;   
       me_eff->setBinContent(ix,efficiency);
     }
     */
     // Binomial errors "B" asked by Florian
     /*me1_forEff->Sumw2(); me2_forEff->Sumw2();*/ me_eff->getTH1F()->Sumw2();
     me_eff->getTH1F()->Divide( me1_forEff, me2_forEff, 1, 1, "B");
   }
 }

void PFClient::createProfilePlots	(	std::string &	folder,
		std::string &	name
	)

private

Definition at line 223 of file PFClient.cc.

References DQMStore::bookProfile(), MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, dqmStore_, DQMStore::get(), MonitorElement::getNbinsX(), MonitorElement::getTH2F(), MonitorElement::getTProfile(), i, MonitorElement::kind(), MonitorElement::setBinContent(), MonitorElement::setBinEntries(), DQMStore::setCurrentFolder(), MonitorElement::setEntries(), AlCaHLTBitMon_QueryRunRegistry::string, fw3dlego::xbins, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

Referenced by doProfiles().

                                                                     {     
   MonitorElement* me = dqmStore_->get(folder+"/"+name);
   if (!me) return;
 
   if ( (me->kind() == MonitorElement::DQM_KIND_TH2F) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2S) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2D) ) {
     TH2* th = me->getTH2F();
     size_t nbinx = me->getNbinsX();
     
     float ymin = th->GetYaxis()->GetXmin();
     float ymax = th->GetYaxis()->GetXmax();
     std::string xtit = th->GetXaxis()->GetTitle();
     std::string ytit = th->GetYaxis()->GetTitle();
     double* xbins = new double[nbinx+1];
     for (size_t ix = 1; ix < nbinx+1; ++ix) {
        xbins[ix-1] = th->GetXaxis()->GetBinLowEdge(ix);
        if (ix == nbinx) xbins[ix] = th->GetXaxis()->GetBinUpEdge(ix);
     }    
 
     std::string tit_new;
     dqmStore_->setCurrentFolder(folder);
     //TProfiles
     MonitorElement* me_profile[2]; 
     me_profile[0] = dqmStore_->bookProfile("profile_"+name, tit_new,
                        nbinx, xbins, ymin, ymax, "" ); 
     me_profile[1] = dqmStore_->bookProfile("profileRMS_"+name, tit_new,
                        nbinx, xbins, ymin, ymax, "s" ); 
     TProfile* profileX = th->ProfileX();
     //size_t nbiny = me->getNbinsY();
     //TProfile* profileX = th->ProfileX("",0,nbiny+1); add underflow and overflow
     static const Int_t NUM_STAT = 6;
     Double_t stats[NUM_STAT] = {0}; th->GetStats(stats);
 
     for (Int_t i = 0; i<2; i++) {
       if (me_profile[i]) {
     for (size_t ix = 0; ix <= nbinx+1; ++ix) {
       me_profile[i]->setBinContent( ix, profileX->GetBinContent(ix)*profileX->GetBinEntries(ix) ) ;
       //me_profile[i]->Fill( profileX->GetBinCenter(ix), profileX->GetBinContent(ix)*profileX->GetBinEntries(ix) ) ;
       me_profile[i]->setBinEntries( ix, profileX->GetBinEntries(ix) );
       me_profile[i]->getTProfile()->GetSumw2()->fArray[ix] = profileX->GetSumw2()->fArray[ix];
       //me_profile[i]->getTProfile()->GetBinSumw2()->fArray[ix] = profileX->GetBinSumw2()->fArray[ix]; // segmentation violation
     }
       }
       me_profile[i]->getTProfile()->PutStats(stats);
       me_profile[i]->setEntries( profileX->GetEntries() ) ;
     }
 
     delete [] xbins;
   }
 }

void PFClient::createProjectionPlots	(	std::string &	folder,
		std::string &	name
	)

private

Definition at line 177 of file PFClient.cc.

References DQMStore::book1D(), MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, dqmStore_, DQMStore::get(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), MonitorElement::getTH2F(), MonitorElement::kind(), MonitorElement::setBinContent(), DQMStore::setCurrentFolder(), MonitorElement::setEntries(), AlCaHLTBitMon_QueryRunRegistry::string, fw3dlego::xbins, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

Referenced by doProjection().

                                                                        {     
   MonitorElement* me = dqmStore_->get(folder+"/"+name);
   if (!me) return;
 
   MonitorElement* projection = 0;
 
   if ( (me->kind() == MonitorElement::DQM_KIND_TH2F) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2S) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2D) ) {
     TH2* th = me->getTH2F();
     size_t nbinx = me->getNbinsX();
     size_t nbiny = me->getNbinsY();
     
     float ymin = th->GetYaxis()->GetXmin();
     float ymax = th->GetYaxis()->GetXmax();
     std::string xtit = th->GetXaxis()->GetTitle();
     std::string ytit = th->GetYaxis()->GetTitle();
     float* xbins = new float[nbinx+1];
     for (size_t ix = 1; ix < nbinx+1; ++ix) {
        xbins[ix-1] = th->GetXaxis()->GetBinLowEdge(ix);
        if (ix == nbinx) xbins[ix] = th->GetXaxis()->GetBinUpEdge(ix);
     }    
 
     std::string tit_new;
     dqmStore_->setCurrentFolder(folder);
 
     if (folder == "ParticleFlow/PFElectronValidation/CompWithGenElectron") {
       if (name == "delta_et_Over_et_VS_et_") projection = dqmStore_->book1D("delta_et_Over_et","E_{T} resolution;#DeltaE_{T}/E_{T}",nbiny,ymin,ymax) ;
       if (name == "delta_et_VS_et_") projection = dqmStore_->book1D("delta_et_","#DeltaE_{T};#DeltaE_{T}",nbiny,ymin,ymax) ;
       if (name == "delta_eta_VS_et_") projection = dqmStore_->book1D("delta_eta_","#Delta#eta;#Delta#eta",nbiny,ymin,ymax) ;
       if (name == "delta_phi_VS_et_") projection = dqmStore_->book1D("delta_phi_","#Delta#phi;#Delta#phi",nbiny,ymin,ymax) ;
     }
 
     if (projection) {
       for (size_t iy = 1; iy < nbiny+1; ++iy) {
     projection->setBinContent(iy, th->ProjectionY("e")->GetBinContent(iy) ); }
       projection->setEntries( me->getEntries() ) ;
     }
     
     delete [] xbins;
   }
 }

void PFClient::createResolutionPlots	(	std::string &	folder,
		std::string &	name
	)

private

Definition at line 114 of file PFClient.cc.

References PDRates::average, DQMStore::book1D(), MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, dqmStore_, DQMStore::get(), getHistogramParameters(), MonitorElement::getName(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), MonitorElement::getTH2F(), MonitorElement::kind(), timingPdfMaker::mean, DQMStore::removeElement(), MonitorElement::Reset(), plotscripts::rms(), MonitorElement::setBinContent(), DQMStore::setCurrentFolder(), MonitorElement::setEfficiencyFlag(), AlCaHLTBitMon_QueryRunRegistry::string, fw3dlego::xbins, SiStripMonitorClusterAlca_cfi::ymax, and SiStripMonitorClusterAlca_cfi::ymin.

Referenced by doSummaries().

                                                                        {     
   MonitorElement* me = dqmStore_->get(folder+"/"+name);
   if (!me) return;
 
   MonitorElement* me_average; 
   MonitorElement* me_rms; 
   MonitorElement* me_mean; 
   MonitorElement* me_sigma; 
 
   if ( (me->kind() == MonitorElement::DQM_KIND_TH2F) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2S) ||
        (me->kind() == MonitorElement::DQM_KIND_TH2D) ) {
     TH2* th = me->getTH2F();
     size_t nbinx = me->getNbinsX();
     size_t nbiny = me->getNbinsY();
     
     float ymin = th->GetYaxis()->GetXmin();
     float ymax = th->GetYaxis()->GetXmax();
     std::string xtit = th->GetXaxis()->GetTitle();
     std::string ytit = th->GetYaxis()->GetTitle();
     float* xbins = new float[nbinx+1];
     for (size_t ix = 1; ix < nbinx+1; ++ix) {
        xbins[ix-1] = th->GetXaxis()->GetBinLowEdge(ix);
        if (ix == nbinx) xbins[ix] = th->GetXaxis()->GetBinUpEdge(ix);
     }
     std::string tit_new = ";"+xtit+";"+ytit;
     dqmStore_->setCurrentFolder(folder);
     MonitorElement* me_slice = dqmStore_->book1D("PFlowSlice","PFlowSlice",nbiny,ymin,ymax); 
 
     tit_new = ";"+xtit+";Average_"+ytit; 
     me_average = dqmStore_->book1D("average_"+name,tit_new, nbinx, xbins); 
     me_average->setEfficiencyFlag();
     tit_new = ";"+xtit+";RMS_"+ytit; 
     me_rms     = dqmStore_->book1D("rms_"+name,tit_new, nbinx, xbins); 
     me_rms->setEfficiencyFlag();
     tit_new = ";"+xtit+";Mean_"+ytit; 
     me_mean    = dqmStore_->book1D("mean_"+name,tit_new, nbinx, xbins); 
     me_mean->setEfficiencyFlag();
     tit_new = ";"+xtit+";Sigma_"+ytit;               
     me_sigma   = dqmStore_->book1D("sigma_"+name,tit_new, nbinx, xbins); 
     me_sigma->setEfficiencyFlag();
 
     double  average, rms, mean, sigma;
 
     for (size_t ix = 1; ix < nbinx+1; ++ix) {
       me_slice->Reset();
       for (size_t iy = 1; iy < nbiny+1; ++iy) {
     me_slice->setBinContent(iy,th->GetBinContent(ix,iy)); 
       }
       getHistogramParameters(me_slice, average, rms, mean, sigma);
       me_average->setBinContent(ix,average);
       me_rms->setBinContent(ix,rms);
       me_mean->setBinContent(ix,mean);
       me_sigma->setBinContent(ix,sigma);
     }
     if (me_slice) dqmStore_->removeElement(me_slice->getName());
     delete [] xbins;
   }
 }

void PFClient::doEfficiency ( )

private

Definition at line 99 of file PFClient.cc.

References createEfficiencyPlots(), effHistogramNames_, folderNames_, cmsHarvester::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                             {
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                     ifolder != folderNames_.end(); ifolder++) {
     std::string path  = "ParticleFlow/"+(*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = effHistogramNames_.begin();
      ihist != effHistogramNames_.end(); ihist++) {
       std::string hname = (*ihist);
       createEfficiencyPlots(path, hname);
     }
   }
 }

void PFClient::doProfiles ( )

private

Definition at line 83 of file PFClient.cc.

References createProfilePlots(), folderNames_, cmsHarvester::path, profileHistogramNames_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                           {
 
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                     ifolder != folderNames_.end(); ifolder++) {
     std::string path  = "ParticleFlow/"+(*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = profileHistogramNames_.begin();
        ihist != profileHistogramNames_.end(); ihist++) {
       std::string hname = (*ihist);
       createProfilePlots(path, hname);
     }
   }
 }

void PFClient::doProjection ( )

private

Definition at line 67 of file PFClient.cc.

References createProjectionPlots(), folderNames_, cmsHarvester::path, projectionHistogramNames_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                             {
 
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                     ifolder != folderNames_.end(); ifolder++) {
     std::string path  = "ParticleFlow/"+(*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = projectionHistogramNames_.begin();
        ihist != projectionHistogramNames_.end(); ihist++) {
       std::string hname = (*ihist);
       createProjectionPlots(path, hname);
     }
   }
 }

void PFClient::doSummaries ( )

private

Definition at line 51 of file PFClient.cc.

References createResolutionPlots(), folderNames_, histogramNames_, cmsHarvester::path, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endRun().

                            {
 
   for (std::vector<std::string>::const_iterator ifolder = folderNames_.begin();
                                     ifolder != folderNames_.end(); ifolder++) {
     std::string path  = "ParticleFlow/"+(*ifolder);
 
     for (std::vector<std::string>::const_iterator ihist = histogramNames_.begin();
        ihist != histogramNames_.end(); ihist++) {
       std::string hname = (*ihist);
       createResolutionPlots(path, hname);
     }
   }
 }

void PFClient::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 45 of file PFClient.cc.

                       {
 
 }

void PFClient::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file PFClient.cc.

References doEfficiency(), doProfiles(), doProjection(), doSummaries(), efficiencyFlag_, and profileFlag_.

                                                                     {
   doSummaries(); 
   doProjection();
   if (efficiencyFlag_) doEfficiency();
   if (profileFlag_) doProfiles();
 }

void PFClient::getHistogramParameters	(	MonitorElement *	me_slice,
		double &	avarage,
		double &	rms,
		double &	mean,
		double &	sigma
	)

private

Definition at line 278 of file PFClient.cc.

References MonitorElement::DQM_KIND_TH1F, MonitorElement::getMean(), MonitorElement::getRMS(), MonitorElement::getTH1F(), and MonitorElement::kind().

Referenced by createResolutionPlots().

                                                                   {
   average = 0.0;
   rms     = 0.0;
   mean    = 0.0;
   sigma   = 0.0;
 
   if (!me_slice) return;
   if  (me_slice->kind() == MonitorElement::DQM_KIND_TH1F) {
     average = me_slice->getMean();
     rms     = me_slice->getRMS();  
     TH1F* th_slice = me_slice->getTH1F();
     if (th_slice && th_slice->GetEntries() > 0) {
       th_slice->Fit( "gaus","Q0");
       TF1* gaus = th_slice->GetFunction( "gaus" );
       if (gaus) {
     sigma = gaus->GetParameter(2);
         mean  = gaus->GetParameter(1);
       }
     }
   }
 }