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

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 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 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 64 of file PhotonOfflineClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file PhotonOfflineClient.cc.

References conversionPostprocessing_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, conversionPostprocessing_cfi::etBin, conversionPostprocessing_cfi::etMax, btagDijet_cfi::etMin, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), conversionPostprocessing_cfi::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 69 of file PhotonOfflineClient.cc.

70 {}

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 359 of file PhotonOfflineClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, callgraph::kind, DQMStore::IBooker::pwd(), DQMStore::IBooker::setCurrentFolder(), fftjetcommon_cfi::title, parallelization::uint(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::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 420 of file PhotonOfflineClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), TkAlMuonSelectors_cfi::cut, callgraph::kind, DQMStore::IBooker::pwd(), DQMStore::IBooker::setCurrentFolder(), fftjetcommon_cfi::title, parallelization::uint(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::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();
   }
 
 
 }

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

Definition at line 333 of file PhotonOfflineClient.cc.

References stringResolutionProvider_cfi::bin, DQMStore::IBooker::book1D(), SiStripPI::max, min(), DQMStore::IBooker::pwd(), and fftjetcommon_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 293 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setEfficiencyFlag(), 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 315 of file PhotonOfflineClient.cc.

References MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::setBinContent(), MonitorElement::setBinError(), 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 74 of file PhotonOfflineClient.cc.

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

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

Definition at line 486 of file PhotonOfflineClient.cc.

References 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 80 of file PhotonOfflineClient.cc.

References gather_cfg::cout, TkAlMuonSelectors_cfi::cut, pfDeepCMVADiscriminatorsJetTags_cfi::denominator, DQMStore::IGetter::dirExists(), conversionPostprocessing_cfi::etaBin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, conversionPostprocessing_cfi::etBin, conversionPostprocessing_cfi::etMax, btagDijet_cfi::etMin, DQMStore::IGetter::get(), MonitorElement::getIntValue(), pfDeepCMVADiscriminatorsJetTags_cfi::numerator, conversionPostprocessing_cfi::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, DQMStore::IGetter::removeElement(), 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_+"/");
   iGetter.removeElement("numberOfHistogramsInPhotonsFolder");
   iGetter.removeElement("numberOfHistogramsInConversionsFolder");
   iGetter.removeElement("numberOfHistogramsInEfficiencyFolder");
   iGetter.removeElement("numberOfHistogramsInInvMassFolder");
 
   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);
 
   iGetter.removeElement("phoEtaPreHLT");
   iGetter.removeElement("phoEtPreHLT");
   iGetter.removeElement("phoEtaPostHLT");
   iGetter.removeElement("phoEtPostHLT");
   iGetter.removeElement("phoEtaLoose");
   iGetter.removeElement("phoEtaTight");
   iGetter.removeElement("phoEtLoose");
   iGetter.removeElement("phoEtTight");
   iGetter.removeElement("phoEtaVertex");
   iGetter.removeElement("convEtaLoose");
   iGetter.removeElement("convEtaTight");
   iGetter.removeElement("convEtLoose");
   iGetter.removeElement("convEtTight");
 
   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/");
       iGetter.removeElement("phoConvEtaForEfficiency");
       iGetter.removeElement("phoConvPhiForEfficiencyAllEcal");
       iGetter.removeElement("phoConvPhiForEfficiencyBarrel");
       iGetter.removeElement("phoConvPhiForEfficiencyEndcaps");
     }
   }  
   
   // if(standAlone_) dbe_->save(outputFileName_);
   //else if(batch_) dbe_->save(inputFileName_);
 
 }

Member Data Documentation

std::string PhotonOfflineClient::analyzerName_

private

Definition at line 102 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::batch_

private

Definition at line 148 of file PhotonOfflineClient.h.

std::stringstream PhotonOfflineClient::currentFolder_

private

Definition at line 153 of file PhotonOfflineClient.h.

double PhotonOfflineClient::cutStep_

private

Definition at line 132 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etaBin

private

Definition at line 142 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMax

private

Definition at line 141 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etaMin

private

Definition at line 140 of file PhotonOfflineClient.h.

int PhotonOfflineClient::etBin

private

Definition at line 139 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMax

private

Definition at line 138 of file PhotonOfflineClient.h.

double PhotonOfflineClient::etMin

private

Definition at line 137 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::excludeBkgHistos_

private

Definition at line 135 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_conversions_

private

Definition at line 156 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_efficiency_

private

Definition at line 157 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_invMass_

private

Definition at line 158 of file PhotonOfflineClient.h.

int PhotonOfflineClient::histo_index_photons_

private

Definition at line 155 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::inputFileName_

private

Definition at line 151 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::minimalSetOfHistos_

private

Definition at line 134 of file PhotonOfflineClient.h.

int PhotonOfflineClient::numberOfSteps_

private

Definition at line 133 of file PhotonOfflineClient.h.

std::string PhotonOfflineClient::outputFileName_

private

Definition at line 150 of file PhotonOfflineClient.h.

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

private

Definition at line 121 of file PhotonOfflineClient.h.

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

private

Definition at line 119 of file PhotonOfflineClient.h.

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

private

Definition at line 120 of file PhotonOfflineClient.h.

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

private

Definition at line 117 of file PhotonOfflineClient.h.

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

private

Definition at line 115 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaLoose_

private

Definition at line 110 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtaTight_

private

Definition at line 112 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtLoose_

private

Definition at line 111 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_convFractionVsEtTight_

private

Definition at line 113 of file PhotonOfflineClient.h.

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

private

Definition at line 116 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaHLT_

private

Definition at line 107 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaLoose_

private

Definition at line 103 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtaTight_

private

Definition at line 105 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtHLT_

private

Definition at line 108 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtLoose_

private

Definition at line 104 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_efficiencyVsEtTight_

private

Definition at line 106 of file PhotonOfflineClient.h.

MonitorElement* PhotonOfflineClient::p_vertexReconstructionEfficiencyVsEta_

private

Definition at line 123 of file PhotonOfflineClient.h.

edm::ParameterSet PhotonOfflineClient::parameters_

private

Definition at line 130 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 161 of file PhotonOfflineClient.h.

int PhotonOfflineClient::phiBin

private

Definition at line 145 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMax

private

Definition at line 144 of file PhotonOfflineClient.h.

double PhotonOfflineClient::phiMin

private

Definition at line 143 of file PhotonOfflineClient.h.

bool PhotonOfflineClient::standAlone_

private

Definition at line 147 of file PhotonOfflineClient.h.

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

private

Definition at line 160 of file PhotonOfflineClient.h.

int PhotonOfflineClient::verbosity_

private

Definition at line 128 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