#include <GlobalHitsProdHistStripper.h>

Inheritance diagram for GlobalHitsProdHistStripper:

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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

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

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)

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
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Class to fill dqm monitor elements from existing EDM file

Author: M. Strang SUNY-Buffalo

Definition at line 42 of file GlobalHitsProdHistStripper.h.

Constructor & Destructor Documentation

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, DQMStore::setVerbose(), DQMStore::showDirStructure(), 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->();
   if (dbe) {
     if (verbosity > 0 ) {
       dbe->setVerbose(1);
     } else {
       dbe->setVerbose(0);
     }
   }
 
   if (dbe) {
     if (verbosity > 0 ) dbe->showDirStructure();
   }
 
   // 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 ( )

override

Definition at line 74 of file GlobalHitsProdHistStripper.cc.

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

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

Member Function Documentation

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

override

Definition at line 275 of file GlobalHitsProdHistStripper.cc.

 {
   return;
 }

void GlobalHitsProdHistStripper::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 80 of file GlobalHitsProdHistStripper.cc.

 {
   return;
 }

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

override

Definition at line 94 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;
 }

void GlobalHitsProdHistStripper::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 85 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;
 }

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

override

Definition at line 148 of file GlobalHitsProdHistStripper.cc.

References DQMStore::book1D(), dbe, edm::Run::getManyByType(), mps_fire::i, edm::HandleBase::isValid(), me, edm::HandleBase::provenance(), 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;
 }