#include <PhotonOfflineClient.h>

Inheritance diagram for PhotonOfflineClient:

Public Member Functions
void	book2DHistoVector (DQMStore::IBooker &iBooker, std::vector< std::vector< MonitorElement * > > &vecOfHist, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)

void	book3DHistoVector (DQMStore::IBooker &iBooker, std::vector< std::vector< std::vector< MonitorElement * > > > &vecOfHist, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)

MonitorElement *	bookHisto (DQMStore::IBooker &iBooker, std::string histoName, std::string title, int bin, double min, double max)

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

	PhotonOfflineClient (const edm::ParameterSet &pset)

MonitorElement *	retrieveHisto (DQMStore::IGetter &iGetter, std::string dir, std::string name)

virtual void	runClient (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter)

	~PhotonOfflineClient () 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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	dividePlots (MonitorElement dividend, MonitorElement numerator, MonitorElement *denominator)

void	dividePlots (MonitorElement dividend, MonitorElement numerator, double denominator)

Private Attributes
std::string	analyzerName_

bool	batch_

std::stringstream	currentFolder_

double	cutStep_

int	etaBin

double	etaMax

double	etaMin

int	etBin

double	etMax

double	etMin

bool	excludeBkgHistos_

int	histo_index_conversions_

int	histo_index_efficiency_

int	histo_index_invMass_

int	histo_index_photons_

std::string	inputFileName_

bool	minimalSetOfHistos_

int	numberOfSteps_

std::string	outputFileName_

std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsEt_

std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsEta_

std::vector< std::vector< MonitorElement * > >	p_badChannelsFractionVsPhi_

std::vector< std::vector< MonitorElement * > >	p_convFractionVsEt_

std::vector< std::vector< MonitorElement * > >	p_convFractionVsEta_

MonitorElement *	p_convFractionVsEtaLoose_

MonitorElement *	p_convFractionVsEtaTight_

MonitorElement *	p_convFractionVsEtLoose_

MonitorElement *	p_convFractionVsEtTight_

std::vector< std::vector< std::vector< MonitorElement * > > >	p_convFractionVsPhi_

MonitorElement *	p_efficiencyVsEtaHLT_

MonitorElement *	p_efficiencyVsEtaLoose_

MonitorElement *	p_efficiencyVsEtaTight_

MonitorElement *	p_efficiencyVsEtHLT_

MonitorElement *	p_efficiencyVsEtLoose_

MonitorElement *	p_efficiencyVsEtTight_

MonitorElement *	p_vertexReconstructionEfficiencyVsEta_

edm::ParameterSet	parameters_

std::vector< std::string >	parts_

int	phiBin

double	phiMax

double	phiMin

bool	standAlone_

std::vector< std::string >	types_

int	verbosity_

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

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::ProducerBase
ProducesCollector	producesCollector ()

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

$Id: PhotonOfflineClient authors: Nancy Marinelli, U. of Notre Dame, US Jamie Antonelli, U. of Notre Dame, US

Definition at line 60 of file PhotonOfflineClient.h.

Constructor & Destructor Documentation

PhotonOfflineClient::PhotonOfflineClient ( const edm::ParameterSet & pset )

explicit

Definition at line 22 of file PhotonOfflineClient.cc.

References muonRecoAnalyzer_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, photonAnalyzer_cfi::etBin, cosmicPhotonAnalyzer_cfi::etMax, photonAnalyzer_cfi::etMin, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), BeamMonitor_cff::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, and muonDTDigis_cfi::pset.

                                                                     {
   //dbe_ = 0;
   //dbe_ = edm::Service<DQMStore>().operator->();
   // dbe_->setVerbose(0);
   parameters_ = pset;
 
   analyzerName_ = pset.getParameter<string>("analyzerName");
   cutStep_ = pset.getParameter<double>("cutStep");
   numberOfSteps_ = pset.getParameter<int>("numberOfSteps");
 
   etMin = pset.getParameter<double>("etMin");
   etMax = pset.getParameter<double>("etMax");
   etBin = pset.getParameter<int>("etBin");
   etaMin = pset.getParameter<double>("etaMin");
   etaMax = pset.getParameter<double>("etaMax");
   etaBin = pset.getParameter<int>("etaBin");
   phiMin = pset.getParameter<double>("phiMin");
   phiMax = pset.getParameter<double>("phiMax");
   phiBin = pset.getParameter<int>("phiBin");
 
   standAlone_ = pset.getParameter<bool>("standAlone");
   batch_ = pset.getParameter<bool>("batch");
   excludeBkgHistos_ = pset.getParameter<bool>("excludeBkgHistos");
 
   outputFileName_ = pset.getParameter<string>("OutputFileName");
   inputFileName_ = pset.getUntrackedParameter<string>("InputFileName");
 
   histo_index_photons_ = 0;
   histo_index_conversions_ = 0;
   histo_index_efficiency_ = 0;
   histo_index_invMass_ = 0;
 
   types_.push_back("All");
   types_.push_back("GoodCandidate");
   if (!excludeBkgHistos_)
     types_.push_back("Background");
 
   parts_.push_back("AllEcal");
   parts_.push_back("Barrel");
   parts_.push_back("Endcaps");
 }

PhotonOfflineClient::~PhotonOfflineClient ( )

override

Definition at line 64 of file PhotonOfflineClient.cc.

64 {}

Member Function Documentation

void PhotonOfflineClient::book2DHistoVector	(	DQMStore::IBooker &	iBooker,
		std::vector< std::vector< MonitorElement * > > &	vecOfHist,
		std::string	histoType,
		std::string	histoName,
		std::string	title,
		int	xbin,
		double	xmin,
		double	xmax,
		int	ybin = `1`,
		double	ymin = `1`,
		double	ymax = `2`
	)

Definition at line 418 of file PhotonOfflineClient.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, HltBtagPostValidation_cff::histoName, CalibrationSummaryClient_cfi::kind, dqm::dqmstoreimpl::DQMStore::IBooker::pwd(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), overlapproblemtsosanalyzer_cfi::title, parallelization::uint, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

                                                          {
   int histo_index = 0;
 
   vector<MonitorElement*> temp1DVector;
 
   //determining which folder we're in
   bool conversionPlot = false;
   if (iBooker.pwd().find("Conversions") != string::npos)
     conversionPlot = true;
 
   if (conversionPlot) {
     histo_index_conversions_++;
     histo_index = histo_index_conversions_;
   } else {
     histo_index_photons_++;
     histo_index = histo_index_photons_;
   }
 
   stringstream histo_number_stream;
   histo_number_stream << "h_";
   if (histo_index < 10)
     histo_number_stream << "0";
   histo_number_stream << histo_index << "_";
 
   for (int cut = 0; cut != numberOfSteps_; ++cut) {  //looping over Et cut values
 
     for (uint type = 0; type != types_.size(); ++type) {  //looping over isolation type
 
       currentFolder_.str("");
       currentFolder_ << "Egamma/" + analyzerName_ + "/" << types_[type] << "Photons/Et above " << (cut + 1) * cutStep_
                      << " GeV";
       if (conversionPlot)
         currentFolder_ << "/Conversions";
 
       iBooker.setCurrentFolder(currentFolder_.str());
 
       string kind;
       if (conversionPlot)
         kind = " Conversions: ";
       else
         kind = " Photons: ";
 
       if (histoType == "1D")
         temp1DVector.push_back(
             iBooker.book1D(histo_number_stream.str() + histoName, types_[type] + kind + title, xbin, xmin, xmax));
       else if (histoType == "2D")
         temp1DVector.push_back(iBooker.book2D(
             histo_number_stream.str() + histoName, types_[type] + kind + title, xbin, xmin, xmax, ybin, ymin, ymax));
       else if (histoType == "Profile")
         temp1DVector.push_back(iBooker.bookProfile(
             histo_number_stream.str() + histoName, types_[type] + kind + title, xbin, xmin, xmax, ybin, ymin, ymax, ""));
       //else cout << "bad histoType\n";
     }
 
     temp2DVector.push_back(temp1DVector);
     temp1DVector.clear();
   }
 }

void PhotonOfflineClient::book3DHistoVector	(	DQMStore::IBooker &	iBooker,
		std::vector< std::vector< std::vector< MonitorElement * > > > &	vecOfHist,
		std::string	histoType,
		std::string	histoName,
		std::string	title,
		int	xbin,
		double	xmin,
		double	xmax,
		int	ybin = `1`,
		double	ymin = `1`,
		double	ymax = `2`
	)

Definition at line 487 of file PhotonOfflineClient.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, CalibrationSummaryClient_cfi::kind, dqm::dqmstoreimpl::DQMStore::IBooker::pwd(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), overlapproblemtsosanalyzer_cfi::title, parallelization::uint, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

                                                          {
   int histo_index = 0;
 
   vector<MonitorElement*> temp1DVector;
   vector<vector<MonitorElement*> > temp2DVector;
 
   //determining which folder we're in
   bool conversionPlot = false;
   if (iBooker.pwd().find("Conversions") != string::npos)
     conversionPlot = true;
 
   if (conversionPlot) {
     histo_index_conversions_++;
     histo_index = histo_index_conversions_;
   } else {
     histo_index_photons_++;
     histo_index = histo_index_photons_;
   }
 
   stringstream histo_number_stream;
   histo_number_stream << "h_";
   if (histo_index < 10)
     histo_number_stream << "0";
   histo_number_stream << histo_index << "_";
 
   for (int cut = 0; cut != numberOfSteps_; ++cut) {  //looping over Et cut values
 
     for (uint type = 0; type != types_.size(); ++type) {  //looping over isolation type
 
       for (uint part = 0; part != parts_.size(); ++part) {  //looping over different parts of the ecal
 
         currentFolder_.str("");
         currentFolder_ << "Egamma/" + analyzerName_ + "/" << types_[type] << "Photons/Et above " << (cut + 1) * cutStep_
                        << " GeV";
         if (conversionPlot)
           currentFolder_ << "/Conversions";
 
         iBooker.setCurrentFolder(currentFolder_.str());
 
         string kind;
         if (conversionPlot)
           kind = " Conversions: ";
         else
           kind = " Photons: ";
 
         if (histoType == "1D")
           temp1DVector.push_back(iBooker.book1D(histo_number_stream.str() + histoName + parts_[part],
                                                 types_[type] + kind + parts_[part] + ": " + title,
                                                 xbin,
                                                 xmin,
                                                 xmax));
         else if (histoType == "2D")
           temp1DVector.push_back(iBooker.book2D(histo_number_stream.str() + histoName + parts_[part],
                                                 types_[type] + kind + parts_[part] + ": " + title,
                                                 xbin,
                                                 xmin,
                                                 xmax,
                                                 ybin,
                                                 ymin,
                                                 ymax));
         else if (histoType == "Profile")
           temp1DVector.push_back(iBooker.bookProfile(histo_number_stream.str() + histoName + parts_[part],
                                                      types_[type] + kind + parts_[part] + ": " + title,
                                                      xbin,
                                                      xmin,
                                                      xmax,
                                                      ybin,
                                                      ymin,
                                                      ymax,
                                                      ""));
         //else cout << "bad histoType\n";
       }
 
       temp2DVector.push_back(temp1DVector);
       temp1DVector.clear();
     }
 
     temp3DVector.push_back(temp2DVector);
     temp2DVector.clear();
   }
 }

PhotonOfflineClient::MonitorElement * PhotonOfflineClient::bookHisto	(	DQMStore::IBooker &	iBooker,
		std::string	histoName,
		std::string	title,
		int	bin,
		double	min,
		double	max
	)

Definition at line 396 of file PhotonOfflineClient.cc.

References newFWLiteAna::bin, dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), HltBtagPostValidation_cff::histoName, SiStripPI::max, min(), dqm::dqmstoreimpl::DQMStore::IBooker::pwd(), and overlapproblemtsosanalyzer_cfi::title.

                                                                                                {
   int histo_index = 0;
   stringstream histo_number_stream;
 
   //determining which folder we're in
   if (iBooker.pwd().find("InvMass") != string::npos) {
     histo_index_invMass_++;
     histo_index = histo_index_invMass_;
   }
   if (iBooker.pwd().find("Efficiencies") != string::npos) {
     histo_index_efficiency_++;
     histo_index = histo_index_efficiency_;
   }
   histo_number_stream << "h_";
   if (histo_index < 10)
     histo_number_stream << "0";
   histo_number_stream << histo_index;
 
   return iBooker.book1D(histo_number_stream.str() + "_" + histoName, title, bin, min, max);
 }

void PhotonOfflineClient::dividePlots	(	MonitorElement *	dividend,
		MonitorElement *	numerator,
		MonitorElement *	denominator
	)

private

Definition at line 357 of file PhotonOfflineClient.cc.

References runTheMatrix::err, dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getEntries(), dqm::impl::MonitorElement::getNbinsX(), dqmiolumiharvest::j, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), dqm::impl::MonitorElement::setEfficiencyFlag(), dqm::impl::MonitorElement::setEntries(), mathSSE::sqrt(), and relativeConstraints::value.

                                                                    {
   double value, err;
 
   dividend->setEfficiencyFlag();
   if (denominator->getEntries() == 0)
     return;
 
   for (int j = 1; j <= numerator->getNbinsX(); j++) {
     if (denominator->getBinContent(j) != 0) {
       value = ((double)numerator->getBinContent(j)) / ((double)denominator->getBinContent(j));
       err = sqrt(value * (1 - value) / ((double)denominator->getBinContent(j)));
       dividend->setBinContent(j, value);
       dividend->setBinError(j, err);
     } else {
       dividend->setBinContent(j, 0);
       dividend->setBinError(j, 0);
     }
     dividend->setEntries(numerator->getEntries());
   }
 }

void PhotonOfflineClient::dividePlots	(	MonitorElement *	dividend,
		MonitorElement *	numerator,
		double	denominator
	)

private

Definition at line 380 of file PhotonOfflineClient.cc.

References runTheMatrix::err, dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsX(), dqmiolumiharvest::j, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), dqm::impl::MonitorElement::setEfficiencyFlag(), mathSSE::sqrt(), and relativeConstraints::value.

                                                                                                              {
   double value, err;
 
   dividend->setEfficiencyFlag();
   for (int j = 1; j <= numerator->getNbinsX(); j++) {
     if (denominator != 0) {
       value = ((double)numerator->getBinContent(j)) / denominator;
       err = sqrt(value * (1 - value) / denominator);
       dividend->setBinContent(j, value);
       dividend->setBinError(j, err);
     } else {
       dividend->setBinContent(j, 0);
     }
   }
 }

void PhotonOfflineClient::dqmEndJob	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter
	)

overridevirtual

Implements DQMEDHarvester.

Definition at line 68 of file PhotonOfflineClient.cc.

                                                                                       {
   if (!standAlone_)
     runClient(iBooker, iGetter);
 }

PhotonOfflineClient::MonitorElement * PhotonOfflineClient::retrieveHisto	(	DQMStore::IGetter &	iGetter,
		std::string	dir,
		std::string	name
	)

Definition at line 579 of file PhotonOfflineClient.cc.

References dqm::dqmstoreimpl::DQMStore::IGetter::getContents(), and parallelization::uint.

                                                                                      {
   //cout << "dir = " << dir << endl;
   //cout << "name = " << name << endl;
   vector<MonitorElement*> histoVector;
   uint indexOfRelevantHistogram = 0;
   string fullMEName = "";
   histoVector = iGetter.getContents(dir);
   for (uint index = 0; index != histoVector.size(); index++) {
     string MEName = histoVector[index]->getName();
     if (MEName.find(name) != string::npos) {
       indexOfRelevantHistogram = index;
       break;
     }
   }
   return histoVector[indexOfRelevantHistogram];
 }

void PhotonOfflineClient::runClient	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter
	)

virtual

Definition at line 74 of file PhotonOfflineClient.cc.

References gather_cfg::cout, TkAlMuonSelectors_cfi::cut, HLTTauDQMOffline_cfi::denominator, dqm::dqmstoreimpl::DQMStore::IGetter::dirExists(), muonRecoAnalyzer_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, photonAnalyzer_cfi::etBin, cosmicPhotonAnalyzer_cfi::etMax, photonAnalyzer_cfi::etMin, dqm::dqmstoreimpl::DQMStore::IGetter::get(), dqm::impl::MonitorElement::getIntValue(), HLTTauDQMOffline_cfi::numerator, BeamMonitor_cff::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, dqm::dqmstoreimpl::DQMStore::IGetter::setCurrentFolder(), and parallelization::uint.

                                                                                       {
   //  if(!dbe_) return;
 
   //if(batch_)  dbe_->open(inputFileName_);
 
   //std::cout << " PostProcessing analyzer name " << analyzerName_ << std::endl;
   if (!iGetter.dirExists("Egamma/" + analyzerName_)) {
     std::cout << "Folder Egamma/" + analyzerName_ + " does not exist - Abort the efficiency calculation " << std::endl;
     return;
   }
 
   //find out how many histograms are in the various folders
   histo_index_photons_ = iGetter.get("Egamma/" + analyzerName_ + "/numberOfHistogramsInPhotonsFolder")->getIntValue();
   histo_index_conversions_ =
       iGetter.get("Egamma/" + analyzerName_ + "/numberOfHistogramsInConversionsFolder")->getIntValue();
   histo_index_efficiency_ =
       iGetter.get("Egamma/" + analyzerName_ + "/numberOfHistogramsInEfficiencyFolder")->getIntValue();
   histo_index_invMass_ = iGetter.get("Egamma/" + analyzerName_ + "/numberOfHistogramsInInvMassFolder")->getIntValue();
 
   iGetter.setCurrentFolder("Egamma/" + analyzerName_ + "/");
 
   string AllPath = "Egamma/" + analyzerName_ + "/AllPhotons/";
   string IsoPath = "Egamma/" + analyzerName_ + "/GoodCandidatePhotons/";
   string NonisoPath = "Egamma/" + analyzerName_ + "/BackgroundPhotons/";
   string EffPath = "Egamma/" + analyzerName_ + "/Efficiencies/";
 
   //booking efficiency histograms
   iGetter.setCurrentFolder(EffPath);
 
   p_efficiencyVsEtaLoose_ = bookHisto(iBooker,
                                       "EfficiencyVsEtaLoose",
                                       "Fraction of Photons passing Loose Isolation vs #eta;#eta",
                                       etaBin,
                                       etaMin,
                                       etaMax);
   p_efficiencyVsEtLoose_ = bookHisto(iBooker,
                                      "EfficiencyVsEtLoose",
                                      "Fraction of Photons passing Loose Isolation vs E_{T};E_{T} (GeV)",
                                      etBin,
                                      etMin,
                                      etMax);
   p_efficiencyVsEtaTight_ = bookHisto(iBooker,
                                       "EfficiencyVsEtaTight",
                                       "Fraction of Photons passing Tight Isolation vs #eta;#eta",
                                       etaBin,
                                       etaMin,
                                       etaMax);
   p_efficiencyVsEtTight_ = bookHisto(iBooker,
                                      "EfficiencyVsEtTight",
                                      "Fraction of Photons passing Tight Isolation vs E_{T};E_{T} (GeV)",
                                      etBin,
                                      etMin,
                                      etMax);
 
   p_efficiencyVsEtaHLT_ =
       bookHisto(iBooker, "EfficiencyVsEtaHLT", "Fraction of Photons firing HLT vs #eta;#eta", etaBin, etaMin, etaMax);
   p_efficiencyVsEtHLT_ = bookHisto(
       iBooker, "EfficiencyVsEtHLT", "Fraction of Photons firing HLT vs E_{T};E_{T} (GeV)", etBin, etMin, etMax);
 
   p_convFractionVsEtaLoose_ =
       bookHisto(iBooker,
                 "ConvFractionVsEtaLoose",
                 "Fraction of Loosely Isolated Photons which are matched to two tracks vs #eta;#eta",
                 etaBin,
                 etaMin,
                 etaMax);
   p_convFractionVsEtLoose_ =
       bookHisto(iBooker,
                 "ConvFractionVsEtLoose",
                 "Fraction of Loosely Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",
                 etBin,
                 etMin,
                 etMax);
   p_convFractionVsEtaTight_ =
       bookHisto(iBooker,
                 "ConvFractionVsEtaTight",
                 "Fraction of Tightly Isolated Photons which are matched to two tracks vs #eta;#eta",
                 etaBin,
                 etaMin,
                 etaMax);
   p_convFractionVsEtTight_ =
       bookHisto(iBooker,
                 "ConvFractionVsEtTight",
                 "Fraction of Tightly Isolated Photons which are matched to two tracks vs E_{T};E_{T} (GeV)",
                 etBin,
                 etMin,
                 etMax);
 
   p_vertexReconstructionEfficiencyVsEta_ =
       bookHisto(iBooker,
                 "VertexReconstructionEfficiencyVsEta",
                 "Fraction of Converted Photons which have a valid vertex vs #eta;#eta",
                 etaBin,
                 etaMin,
                 etaMax);
 
   //booking conversion fraction histograms
   iGetter.setCurrentFolder(AllPath + "Et above 20 GeV/Conversions");
   book2DHistoVector(iBooker,
                     p_convFractionVsEt_,
                     "1D",
                     "convFractionVsEt",
                     "Fraction of Converted Photons vs E_{T};E_{T} (GeV)",
                     etBin,
                     etMin,
                     etMax);
   book3DHistoVector(iBooker,
                     p_convFractionVsPhi_,
                     "1D",
                     "convFractionVsPhi",
                     "Fraction of Converted Photons vs #phi;#phi",
                     phiBin,
                     phiMin,
                     phiMax);
   book2DHistoVector(iBooker,
                     p_convFractionVsEta_,
                     "1D",
                     "convFractionVsEta",
                     "Fraction of Converted Photons vs #eta;#eta",
                     etaBin,
                     etaMin,
                     etaMax);
 
   //booking bad channel fraction histograms
   iGetter.setCurrentFolder(AllPath + "Et above 20 GeV/");
   book2DHistoVector(iBooker,
                     p_badChannelsFractionVsPhi_,
                     "1D",
                     "badChannelsFractionVsPhi",
                     "Fraction of Photons which have at least one bad channel vs #phi;#phi",
                     phiBin,
                     phiMin,
                     phiMax);
   book2DHistoVector(iBooker,
                     p_badChannelsFractionVsEta_,
                     "1D",
                     "badChannelsFractionVsEta",
                     "Fraction of Photons which have at least one bad channel vs #eta;#eta",
                     etaBin,
                     etaMin,
                     etaMax);
   book2DHistoVector(iBooker,
                     p_badChannelsFractionVsEt_,
                     "1D",
                     "badChannelsFractionVsEt",
                     "Fraction of Photons which have at least one bad channel vs E_{T};E_{T} (GeV)",
                     etBin,
                     etMin,
                     etMax);
 
   //making efficiency plots
   MonitorElement* dividend;
   MonitorElement* numerator;
   MonitorElement* denominator;
 
   currentFolder_.str("");
   currentFolder_ << AllPath << "Et above 20 GeV/";
 
   //HLT efficiency plots
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtaHLT");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtaPostHLT");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtaPreHLT");
   dividePlots(dividend, numerator, denominator);
 
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtHLT");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtPostHLT");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtPreHLT");
   dividePlots(dividend, numerator, denominator);
 
   //efficiencies vs Eta
   denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEta");
 
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtaLoose");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtaLoose");
   dividePlots(dividend, numerator, denominator);
 
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtaTight");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtaTight");
   dividePlots(dividend, numerator, denominator);
 
   //efficiencies vs Et
   denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEtAllEcal");
 
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtLoose");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtLoose");
   dividePlots(dividend, numerator, denominator);
 
   dividend = retrieveHisto(iGetter, EffPath, "EfficiencyVsEtTight");
   numerator = retrieveHisto(iGetter, EffPath, "phoEtTight");
   dividePlots(dividend, numerator, denominator);
 
   //conversion fractions vs Eta
   dividend = retrieveHisto(iGetter, EffPath, "ConvFractionVsEtaLoose");
   numerator = retrieveHisto(iGetter, EffPath, "convEtaLoose");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtaLoose");
   dividePlots(dividend, numerator, denominator);
 
   dividend = retrieveHisto(iGetter, EffPath, "ConvFractionVsEtaTight");
   numerator = retrieveHisto(iGetter, EffPath, "convEtaTight");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtaTight");
   dividePlots(dividend, numerator, denominator);
 
   //conversion fractions vs Et
   dividend = retrieveHisto(iGetter, EffPath, "ConvFractionVsEtLoose");
   numerator = retrieveHisto(iGetter, EffPath, "convEtLoose");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtLoose");
   dividePlots(dividend, numerator, denominator);
 
   dividend = retrieveHisto(iGetter, EffPath, "ConvFractionVsEtTight");
   numerator = retrieveHisto(iGetter, EffPath, "convEtTight");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtTight");
   dividePlots(dividend, numerator, denominator);
 
   //conversion vertex recontruction efficiency
   dividend = retrieveHisto(iGetter, EffPath, "VertexReconstructionEfficiencyVsEta");
   numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvEta");
   denominator = retrieveHisto(iGetter, EffPath, "phoEtaVertex");
   dividePlots(dividend, numerator, denominator);
 
   iGetter.setCurrentFolder(EffPath);
 
   for (uint type = 0; type != types_.size(); ++type) {
     for (int cut = 0; cut != numberOfSteps_; ++cut) {
       currentFolder_.str("");
       currentFolder_ << "Egamma/" + analyzerName_ + "/" << types_[type] << "Photons/Et above " << (cut + 1) * cutStep_
                      << " GeV/";
 
       //making bad channel histograms
 
       //vs Et
       dividend = retrieveHisto(iGetter, currentFolder_.str(), "badChannelsFractionVsEt");
       numerator = retrieveHisto(iGetter, currentFolder_.str(), "phoEtBadChannels");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEtAllEcal");
       dividePlots(dividend, numerator, denominator);
 
       //vs eta
       dividend = retrieveHisto(iGetter, currentFolder_.str(), "badChannelsFractionVsEta");
       numerator = retrieveHisto(iGetter, currentFolder_.str(), "phoEtaBadChannels");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEta");
       dividePlots(dividend, numerator, denominator);
 
       //vs phi
       dividend = retrieveHisto(iGetter, currentFolder_.str(), "badChannelsFractionVsPhi");
       numerator = retrieveHisto(iGetter, currentFolder_.str(), "phoPhiBadChannels");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoPhiAllEcal");
       dividePlots(dividend, numerator, denominator);
 
       //making conversion fraction histograms
 
       //vs Et
       dividend = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "convFractionVsEt");
       numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvEtAllEcal");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEtAllEcal");
       dividePlots(dividend, numerator, denominator);
 
       //vs eta
       dividend = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "convFractionVsEta");
       numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvEtaForEfficiency");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoEta");
       dividePlots(dividend, numerator, denominator);
 
       //vs phi
       dividend = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "convFractionVsPhiAllEcal");
       numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvPhiForEfficiencyAllEcal");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoPhiAllEcal");
       dividePlots(dividend, numerator, denominator);
       dividend = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "convFractionVsPhiBarrel");
       numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvPhiForEfficiencyBarrel");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoPhiBarrel");
       dividePlots(dividend, numerator, denominator);
       dividend = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "convFractionVsPhiEndcaps");
       numerator = retrieveHisto(iGetter, currentFolder_.str() + "Conversions/", "phoConvPhiForEfficiencyEndcaps");
       denominator = retrieveHisto(iGetter, currentFolder_.str(), "phoPhiEndcaps");
       dividePlots(dividend, numerator, denominator);
 
       iGetter.setCurrentFolder(currentFolder_.str() + "Conversions/");
     }
   }
 
   // if(standAlone_) dbe_->save(outputFileName_);
   //else if(batch_) dbe_->save(inputFileName_);
 }

Member Data Documentation

std::string PhotonOfflineClient::analyzerName_

private

Definition at line 105 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::batch_

private

Definition at line 151 of file PhotonOfflineClient.h.

std::stringstream PhotonOfflineClient::currentFolder_

private

Definition at line 156 of file PhotonOfflineClient.h.

double PhotonOfflineClient::cutStep_

private

Definition at line 135 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etaBin

private

Definition at line 145 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMax

private

Definition at line 144 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMin

private

Definition at line 143 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etBin

private

Definition at line 142 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMax

private

Definition at line 141 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMin

private

Definition at line 140 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::excludeBkgHistos_

private

Definition at line 138 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_conversions_

private

Definition at line 159 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_efficiency_

private

Definition at line 160 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_invMass_

private

Definition at line 161 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_photons_

private

Definition at line 158 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::inputFileName_

private

Definition at line 154 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::minimalSetOfHistos_

private

Definition at line 137 of file PhotonOfflineClient.h.

int PhotonOfflineClient::numberOfSteps_

private

Definition at line 136 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::outputFileName_

private

Definition at line 153 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsEt_

private

Definition at line 124 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsEta_

private

Definition at line 122 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_badChannelsFractionVsPhi_

private

Definition at line 123 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_convFractionVsEt_

private

Definition at line 120 of file PhotonOfflineClient.h.

std::vector<std::vector<MonitorElement*> > PhotonOfflineClient::p_convFractionVsEta_

private

Definition at line 118 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaLoose_

private

Definition at line 113 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaTight_

private

Definition at line 115 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtLoose_

private

Definition at line 114 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtTight_

private

Definition at line 116 of file PhotonOfflineClient.h.

std::vector<std::vector<std::vector<MonitorElement*> > > PhotonOfflineClient::p_convFractionVsPhi_

private

Definition at line 119 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaHLT_

private

Definition at line 110 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaLoose_

private

Definition at line 106 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaTight_

private

Definition at line 108 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtHLT_

private

Definition at line 111 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtLoose_

private

Definition at line 107 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtTight_

private

Definition at line 109 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_vertexReconstructionEfficiencyVsEta_

private

Definition at line 126 of file PhotonOfflineClient.h.

edm::ParameterSet PhotonOfflineClient::parameters_

private

Definition at line 133 of file PhotonOfflineClient.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Types.EDAlias::clone(), and Mixins._TypedParameterizable::copy().

std::vector<std::string> PhotonOfflineClient::parts_

private

Definition at line 164 of file PhotonOfflineClient.h.

int PhotonOfflineClient::phiBin

private

Definition at line 148 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMax

private

Definition at line 147 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMin

private

Definition at line 146 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::standAlone_

private

Definition at line 150 of file PhotonOfflineClient.h.

std::vector<std::string> PhotonOfflineClient::types_

private

Definition at line 163 of file PhotonOfflineClient.h.

int PhotonOfflineClient::verbosity_

private

Definition at line 131 of file PhotonOfflineClient.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation