#include <GlobalHitsProdHistStripper.h>

Inheritance diagram for GlobalHitsProdHistStripper:

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

void	beginJob (void) override

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

void	endJob () override

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

	GlobalHitsProdHistStripper (const edm::ParameterSet &)

	~GlobalHitsProdHistStripper () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
unsigned int	count

DQMStore *	dbe

bool	doOutput

std::string	fName

int	frequency

bool	getAllProvenances

std::vector< MonitorElement * >	me

MonitorElement *	meCaloEcal [2]

MonitorElement *	meCaloEcalE [2]

MonitorElement *	meCaloEcalEta

MonitorElement *	meCaloEcalPhi

MonitorElement *	meCaloEcalToF [2]

MonitorElement *	meCaloHcal [2]

MonitorElement *	meCaloHcalE [2]

MonitorElement *	meCaloHcalEta

MonitorElement *	meCaloHcalPhi

MonitorElement *	meCaloHcalToF [2]

MonitorElement *	meCaloPreSh [2]

MonitorElement *	meCaloPreShE [2]

MonitorElement *	meCaloPreShEta

MonitorElement *	meCaloPreShPhi

MonitorElement *	meCaloPreShToF [2]

MonitorElement *	meGeantTrkE

MonitorElement *	meGeantTrkPt

MonitorElement *	meGeantVtxX [2]

MonitorElement *	meGeantVtxY [2]

MonitorElement *	meGeantVtxZ [2]

MonitorElement *	meMCG4Trk [2]

MonitorElement *	meMCG4Vtx [2]

MonitorElement *	meMCRGP [2]

MonitorElement *	meMuon [2]

MonitorElement *	meMuonCscToF [2]

MonitorElement *	meMuonCscZ

MonitorElement *	meMuonDtR

MonitorElement *	meMuonDtToF [2]

MonitorElement *	meMuonEta

MonitorElement *	meMuonPhi

MonitorElement *	meMuonRpcBR

MonitorElement *	meMuonRpcBToF [2]

MonitorElement *	meMuonRpcFToF [2]

MonitorElement *	meMuonRpcFZ

MonitorElement *	meTrackerPx [2]

MonitorElement *	meTrackerPxBR

MonitorElement *	meTrackerPxBToF

MonitorElement *	meTrackerPxEta

MonitorElement *	meTrackerPxFToF

MonitorElement *	meTrackerPxFZ

MonitorElement *	meTrackerPxPhi

MonitorElement *	meTrackerSi [2]

MonitorElement *	meTrackerSiBR

MonitorElement *	meTrackerSiBToF

MonitorElement *	meTrackerSiEta

MonitorElement *	meTrackerSiFToF

MonitorElement *	meTrackerSiFZ

MonitorElement *	meTrackerSiPhi

std::map< std::string, MonitorElement * >	monitorElements

std::string	outputfile

bool	printProvenanceInfo

int	verbosity

int	vtxunit

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Class to fill dqm monitor elements from existing EDM file

Author: M. Strang SUNY-Buffalo

Definition at line 41 of file GlobalHitsProdHistStripper.h.

Member Typedef Documentation

◆ DQMStore

typedef dqm::legacy::DQMStore GlobalHitsProdHistStripper::DQMStore

Definition at line 43 of file GlobalHitsProdHistStripper.h.

◆ MonitorElement

typedef dqm::legacy::MonitorElement GlobalHitsProdHistStripper::MonitorElement

Definition at line 44 of file GlobalHitsProdHistStripper.h.

Constructor & Destructor Documentation

◆ GlobalHitsProdHistStripper()

GlobalHitsProdHistStripper::GlobalHitsProdHistStripper ( const edm::ParameterSet & iPSet )

explicit

Definition at line 12 of file GlobalHitsProdHistStripper.cc.

References dbe, doOutput, fName, frequency, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), Utilities::operator, outputfile, printProvenanceInfo, AlCaHLTBitMon_QueryRunRegistry::string, verbosity, and vtxunit.

     : fName(""),
       verbosity(0),
       frequency(0),
       vtxunit(0),
       getAllProvenances(false),
       printProvenanceInfo(false),
       outputfile(""),
       count(0) {
   std::string MsgLoggerCat = "GlobalHitsProdHistStripper_GlobalHitsProdHistStripper";
 
   // get information from parameter set
   fName = iPSet.getUntrackedParameter<std::string>("Name");
   verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
   frequency = iPSet.getUntrackedParameter<int>("Frequency");
   vtxunit = iPSet.getUntrackedParameter<int>("VtxUnit");
   outputfile = iPSet.getParameter<std::string>("OutputFile");
   doOutput = iPSet.getParameter<bool>("DoOutput");
   edm::ParameterSet m_Prov = iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
   getAllProvenances = m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
   printProvenanceInfo = m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
 
   // use value of first digit to determine default output level (inclusive)
   // 0 is none, 1 is basic, 2 is fill output, 3 is gather output
   verbosity %= 10;
 
   // get dqm info
   dbe = nullptr;
   dbe = edm::Service<DQMStore>().operator->();
 
   // print out Parameter Set information being used
   if (verbosity >= 0) {
     edm::LogInfo(MsgLoggerCat) << "\n===============================\n"
                                << "Initialized as EDAnalyzer with parameter values:\n"
                                << "    Name           = " << fName << "\n"
                                << "    Verbosity      = " << verbosity << "\n"
                                << "    Frequency      = " << frequency << "\n"
                                << "    VtxUnit        = " << vtxunit << "\n"
                                << "    OutputFile     = " << outputfile << "\n"
                                << "    DoOutput      = " << doOutput << "\n"
                                << "    GetProv        = " << getAllProvenances << "\n"
                                << "    PrintProv      = " << printProvenanceInfo << "\n"
                                << "===============================\n";
   }
 }

◆ ~GlobalHitsProdHistStripper()

GlobalHitsProdHistStripper::~GlobalHitsProdHistStripper ( )

override

Definition at line 58 of file GlobalHitsProdHistStripper.cc.

References dbe, doOutput, outputfile, and dqm::implementation::DQMStore::save().

                                                         {
   if (doOutput)
     if (!outputfile.empty() && dbe)
       dbe->save(outputfile);
 }

Member Function Documentation

◆ analyze()

void GlobalHitsProdHistStripper::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 251 of file GlobalHitsProdHistStripper.cc.

251 { return; }

◆ beginJob()

void GlobalHitsProdHistStripper::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 64 of file GlobalHitsProdHistStripper.cc.

64 { return; }

◆ beginRun()

void GlobalHitsProdHistStripper::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 73 of file GlobalHitsProdHistStripper.cc.

References count, frequency, getAllProvenances, edm::Run::getAllStableProvenance(), mps_fire::i, printProvenanceInfo, edm::RunBase::run(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

                                                                                          {
   std::string MsgLoggerCat = "GlobalHitsProdHistStripper_beginRun";
   // keep track of number of runs processed
   ++count;
 
   int nrun = iRun.run();
 
   if (verbosity > 0) {
     edm::LogInfo(MsgLoggerCat) << "Processing run " << nrun << " (" << count << " runs total)";
   } else if (verbosity == 0) {
     if (nrun % frequency == 0 || count == 1) {
       edm::LogInfo(MsgLoggerCat) << "Processing run " << nrun << " (" << count << " runs total)";
     }
   }
 
   if (getAllProvenances) {
     std::vector<const edm::StableProvenance *> AllProv;
     iRun.getAllStableProvenance(AllProv);
 
     if (verbosity >= 0)
       edm::LogInfo(MsgLoggerCat) << "Number of Provenances = " << AllProv.size();
 
     if (printProvenanceInfo && (verbosity >= 0)) {
       TString eventout("\nProvenance info:\n");
 
       for (unsigned int i = 0; i < AllProv.size(); ++i) {
         eventout += "\n       ******************************";
         eventout += "\n       Module       : ";
         eventout += AllProv[i]->moduleLabel();
         eventout += "\n       ProductID    : ";
         eventout += AllProv[i]->productID().id();
         eventout += "\n       ClassName    : ";
         eventout += AllProv[i]->className();
         eventout += "\n       InstanceName : ";
         eventout += AllProv[i]->productInstanceName();
         eventout += "\n       BranchName   : ";
         eventout += AllProv[i]->branchName();
       }
       eventout += "\n       ******************************\n";
       edm::LogInfo(MsgLoggerCat) << eventout << "\n";
       printProvenanceInfo = false;
     }
     getAllProvenances = false;
   }
 
   return;
 }

◆ endJob()

void GlobalHitsProdHistStripper::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 66 of file GlobalHitsProdHistStripper.cc.

References count, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

                                         {
   std::string MsgLoggerCat = "GlobalHitsProdHistStripper_endJob";
   if (verbosity >= 0)
     edm::LogInfo(MsgLoggerCat) << "Terminating having processed " << count << " runs.";
   return;
 }

◆ endRun()

void GlobalHitsProdHistStripper::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 121 of file GlobalHitsProdHistStripper.cc.

References dqm::implementation::IBooker::book1D(), dbe, Exception, mps_fire::i, edm::HandleBase::isValid(), me, HerwigMaxPtPartonFilter_cfi::moduleLabel, edm::HandleBase::provenance(), dqm::implementation::DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, parallelization::uint, and x.

                                                                                        {
   std::string MsgLoggerCat = "GlobalHitsProdHistStripper_endRun";
 
   edm::Handle<TH1F> histogram1D;
   std::vector<edm::Handle<TH1F>> allhistogram1D;
 
   //iRun.getManyByType(allhistogram1D);
   throw cms::Exception("UnsupportedFunction") << "GlobalHitsProdHistStripper::endRun: "
                                               << "getManyByType has not been supported by the Framework since 2015. "
                                               << "This module has been broken since then. Maybe it should be deleted. "
                                               << "Another possibility is to upgrade to use GetterOfProducts instead.";
 
   me.resize(allhistogram1D.size());
 
   for (uint i = 0; i < allhistogram1D.size(); ++i) {
     histogram1D = allhistogram1D[i];
     if (!histogram1D.isValid()) {
       edm::LogWarning(MsgLoggerCat) << "Invalid histogram extracted from event.";
       continue;
     }
 
     me[i] = nullptr;
 
     /*
     std::cout << "Extracting histogram: " << std::endl
               << "       Module       : "
               << (histogram1D.provenance()->branchDescription()).moduleLabel()
               << std::endl
               << "       ProductID    : "
               <<
     (histogram1D.provenance()->branchDescription()).productID().id()
               << std::endl
               << "       ClassName    : "
               << (histogram1D.provenance()->branchDescription()).className()
               << std::endl
               << "       InstanceName : "
               <<
     (histogram1D.provenance()->branchDescription()).productInstanceName()
               << std::endl
               << "       BranchName   : "
               << (histogram1D.provenance()->branchDescription()).branchName()
               << std::endl;
     */
 
     if ((histogram1D.provenance()->branchDescription()).moduleLabel() != "globalhitsprodhist")
       continue;
 
     std::string histname = histogram1D->GetName();
 
     std::string subhist1 = histname.substr(1, 5);
     std::string subhist2 = histname.substr(1, 4);
 
     if (dbe) {
       if (subhist1 == "CaloE" || subhist1 == "CaloP") {
         dbe->setCurrentFolder("GlobalHitsV/ECal");
       } else if (subhist1 == "CaloH") {
         dbe->setCurrentFolder("GlobalHitsV/HCal");
       } else if (subhist1 == "Geant" || subhist2 == "MCG4" || subhist1 == "MCRGP") {
         dbe->setCurrentFolder("GlobalHitsV/MCGeant");
       } else if (subhist2 == "Muon") {
         dbe->setCurrentFolder("GlobalHitsV/Muon");
       } else if (subhist1 == "Track") {
         dbe->setCurrentFolder("GlobalHitsV/Tracker");
       }
 
       me[i] = dbe->book1D(histname,
                           histogram1D->GetTitle(),
                           histogram1D->GetXaxis()->GetNbins(),
                           histogram1D->GetXaxis()->GetXmin(),
                           histogram1D->GetXaxis()->GetXmax());
       me[i]->setAxisTitle(histogram1D->GetXaxis()->GetTitle(), 1);
       me[i]->setAxisTitle(histogram1D->GetYaxis()->GetTitle(), 2);
     }
 
     std::string mename = me[i]->getName();
 
     // std::cout << "Extracting histogram " << histname
     //        << " into MonitorElement " << mename
     //        << std::endl;
 
     for (Int_t x = 1; x <= histogram1D->GetXaxis()->GetNbins(); ++x) {
       Double_t binx = histogram1D->GetBinCenter(x);
       Double_t value = histogram1D->GetBinContent(x);
       me[i]->Fill(binx, value);
     }
   }
   return;
 }