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

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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 11 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 57 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 245 of file GlobalHitsProdHistStripper.cc.

245 { return; }

◆ beginJob()

void GlobalHitsProdHistStripper::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 63 of file GlobalHitsProdHistStripper.cc.

63 { return; }

◆ beginRun()

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

override

Definition at line 72 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 65 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 120 of file GlobalHitsProdHistStripper.cc.

References dqm::implementation::IBooker::book1D(), dbe, edm::Run::getManyByType(), 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);
 
   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;
 }