#include <EcalURecHitHists.h>

Inheritance diagram for EcalURecHitHists:

Public Member Functions
	EcalURecHitHists (const edm::ParameterSet &)

	~EcalURecHitHists () 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 Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

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

void	endJob () override

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

void	initHists (int)

std::string	intToString (int num)

Private Attributes
TH1F *	allFedsHist_

TH1F *	allFedsTimingHist_

const edm::InputTag	ebUncalibratedRecHitCollection_

const edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	ebUncalibRecHitsToken_

const EcalElectronicsMapping *	ecalElectronicsMap_

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	ecalMappingToken_

const edm::InputTag	eeUncalibratedRecHitCollection_

const edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	eeUncalibRecHitsToken_

EcalFedMap *	fedMap_

std::map< int, TH1F * >	FEDsAndHists_

std::map< int, TH1F * >	FEDsAndTimingHists_

TFile *	file

std::string	fileName_

const double	histRangeMax_

const double	histRangeMin_

std::vector< int >	maskedChannels_

std::vector< std::string >	maskedEBs_

std::vector< int >	maskedFEDs_

int	runNum_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 54 of file EcalURecHitHists.h.

Constructor & Destructor Documentation

◆ EcalURecHitHists()

EcalURecHitHists::EcalURecHitHists ( const edm::ParameterSet & iConfig )

explicit

Definition at line 36 of file EcalURecHitHists.cc.

References allFedsHist_, allFedsTimingHist_, edm::BeginRun, fedMap_, EcalFedMap::getFedFromSlice(), edm::ParameterSet::getUntrackedParameter(), histRangeMax_, histRangeMin_, createfilelist::int, maskedChannels_, maskedEBs_, and maskedFEDs_.

     : ebUncalibratedRecHitCollection_(iConfig.getParameter<edm::InputTag>("EBUncalibratedRecHitCollection")),
       eeUncalibratedRecHitCollection_(iConfig.getParameter<edm::InputTag>("EEUncalibratedRecHitCollection")),
       ebUncalibRecHitsToken_(consumes<EcalUncalibratedRecHitCollection>(ebUncalibratedRecHitCollection_)),
       eeUncalibRecHitsToken_(consumes<EcalUncalibratedRecHitCollection>(eeUncalibratedRecHitCollection_)),
       ecalMappingToken_(esConsumes<edm::Transition::BeginRun>()),
       runNum_(-1),
       histRangeMax_(iConfig.getUntrackedParameter<double>("histogramMaxRange", 200.0)),
       histRangeMin_(iConfig.getUntrackedParameter<double>("histogramMinRange", -10.0)),
       fileName_(iConfig.getUntrackedParameter<std::string>("fileName", std::string("ecalURechHitHists"))) {
   vector<int> listDefaults;
   listDefaults.push_back(-1);
 
   maskedChannels_ = iConfig.getUntrackedParameter<vector<int> >("maskedChannels", listDefaults);
   maskedFEDs_ = iConfig.getUntrackedParameter<vector<int> >("maskedFEDs", listDefaults);
 
   vector<string> defaultMaskedEBs;
   defaultMaskedEBs.push_back("none");
   maskedEBs_ = iConfig.getUntrackedParameter<vector<string> >("maskedEBs", defaultMaskedEBs);
 
   fedMap_ = new EcalFedMap();
   string title1 = "Uncalib Rec Hits (ADC counts)";
   string name1 = "URecHitsAllFEDs";
   int numBins = (int)round(histRangeMax_ - histRangeMin_) + 1;
   allFedsHist_ = new TH1F(name1.c_str(), title1.c_str(), numBins, histRangeMin_, histRangeMax_);
   title1 = "Jitter for all FEDs";
   name1 = "JitterAllFEDs";
   allFedsTimingHist_ = new TH1F(name1.c_str(), title1.c_str(), 14, -7, 7);
 
   // load up the maskedFED list with the proper FEDids
   if (maskedFEDs_[0] == -1) {
     //if "actual" EB id given, then convert to FEDid and put in listFEDs_
     if (maskedEBs_[0] != "none") {
       maskedFEDs_.clear();
       for (vector<string>::const_iterator ebItr = maskedEBs_.begin(); ebItr != maskedEBs_.end(); ++ebItr) {
         maskedFEDs_.push_back(fedMap_->getFedFromSlice(*ebItr));
       }
     }
   }
 }

◆ ~EcalURecHitHists()

EcalURecHitHists::~EcalURecHitHists ( )

override

Definition at line 77 of file EcalURecHitHists.cc.

77 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 84 of file EcalURecHitHists.cc.

References allFedsHist_, allFedsTimingHist_, edm::SortedCollection< T, SORT >::begin(), EcalElectronicsId::dccId(), ebUncalibRecHitsToken_, ecalElectronicsMap_, eeUncalibRecHitsToken_, edm::SortedCollection< T, SORT >::end(), FEDsAndHists_, FEDsAndTimingHists_, spr::find(), EcalElectronicsMapping::getElectronicsId(), ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), EBDetId::ic(), EEDetId::ic(), hit::id, iEvent, initHists(), LogDebug, maskedChannels_, maskedFEDs_, mps_fire::result, runNum_, and edm::SortedCollection< T, SORT >::size().

                                                                                   {
   int ievt = iEvent.id().event();
   Handle<EcalUncalibratedRecHitCollection> EBhits;
   Handle<EcalUncalibratedRecHitCollection> EEhits;
 
   iEvent.getByToken(ebUncalibRecHitsToken_, EBhits);
   LogDebug("EcalURecHitHists") << "event " << ievt << " hits collection size " << EBhits->size();
 
   iEvent.getByToken(eeUncalibRecHitsToken_, EEhits);
   LogDebug("EcalURecHitHists") << "event " << ievt << " hits collection size " << EEhits->size();
 
   for (EcalUncalibratedRecHitCollection::const_iterator hitItr = EBhits->begin(); hitItr != EBhits->end(); ++hitItr) {
     EcalUncalibratedRecHit hit = (*hitItr);
     EBDetId ebDet = hit.id();
     int ic = ebDet.ic();
     int hashedIndex = ebDet.hashedIndex();
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(ebDet);
     int FEDid = 600 + elecId.dccId();
     float ampli = hit.amplitude();
 
     vector<int>::iterator result;
     result = find(maskedFEDs_.begin(), maskedFEDs_.end(), FEDid);
     if (result != maskedFEDs_.end()) {
       LogWarning("EcalURecHitHists") << "skipping uncalRecHit for FED " << FEDid << " ; amplitude " << ampli;
       continue;
     }
 
     result = find(maskedChannels_.begin(), maskedChannels_.end(), hashedIndex);
     if (result != maskedChannels_.end()) {
       LogWarning("EcalURecHitHists") << "skipping uncalRecHit for channel: " << ic << " with amplitude " << ampli;
       continue;
     }
 
     // fill the proper hist
     TH1F* uRecHist = FEDsAndHists_[FEDid];
     TH1F* timingHist = FEDsAndTimingHists_[FEDid];
     if (uRecHist == nullptr) {
       initHists(FEDid);
       uRecHist = FEDsAndHists_[FEDid];
       timingHist = FEDsAndTimingHists_[FEDid];
     }
 
     uRecHist->Fill(ampli);
     allFedsHist_->Fill(ampli);
     timingHist->Fill(hit.jitter());
     allFedsTimingHist_->Fill(hit.jitter());
   }
 
   // Again for the endcap
   for (EcalUncalibratedRecHitCollection::const_iterator hitItr = EEhits->begin(); hitItr != EEhits->end(); ++hitItr) {
     EcalUncalibratedRecHit hit = (*hitItr);
     EEDetId eeDet = hit.id();
     int ic = eeDet.ic();
     int hashedIndex = eeDet.hashedIndex();
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(eeDet);
     int FEDid = 600 + elecId.dccId();
     float ampli = hit.amplitude();
 
     vector<int>::iterator result;
     result = find(maskedFEDs_.begin(), maskedFEDs_.end(), FEDid);
     if (result != maskedFEDs_.end()) {
       LogWarning("EcalURecHitHists") << "skipping uncalRecHit for FED " << FEDid << " ; amplitude " << ampli;
       continue;
     }
 
     result = find(maskedChannels_.begin(), maskedChannels_.end(), hashedIndex);
     if (result != maskedChannels_.end()) {
       LogWarning("EcalURecHitHists") << "skipping uncalRecHit for channel: " << ic << " with amplitude " << ampli;
       continue;
     }
 
     // fill the proper hist
     TH1F* uRecHist = FEDsAndHists_[FEDid];
     TH1F* timingHist = FEDsAndTimingHists_[FEDid];
     if (uRecHist == nullptr) {
       initHists(FEDid);
       uRecHist = FEDsAndHists_[FEDid];
       timingHist = FEDsAndTimingHists_[FEDid];
     }
 
     uRecHist->Fill(ampli);
     allFedsHist_->Fill(ampli);
     timingHist->Fill(hit.jitter());
     allFedsTimingHist_->Fill(hit.jitter());
   }
 
   if (runNum_ == -1) {
     runNum_ = iEvent.id().run();
   }
 }

◆ beginRun()

void EcalURecHitHists::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	c
	)

overrideprivate

Definition at line 199 of file EcalURecHitHists.cc.

References DummyCfis::c, ecalElectronicsMap_, and ecalMappingToken_.

                                                                      {
   ecalElectronicsMap_ = &c.getData(ecalMappingToken_);
 }

◆ endJob()

void EcalURecHitHists::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 206 of file EcalURecHitHists.cc.

References allFedsHist_, allFedsTimingHist_, DMR_cfg::cerr, ewkTauDQM_cfi::channels, DeadROC_duringRun::dir, fedMap_, FEDsAndHists_, FEDsAndTimingHists_, fileName_, EcalFedMap::getSliceFromFed(), compareTotals::hist, intToString(), maskedChannels_, runNum_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
   using namespace std;
   fileName_ += "-" + intToString(runNum_) + ".graph.root";
 
   TFile root_file_(fileName_.c_str(), "RECREATE");
 
   for (map<int, TH1F*>::const_iterator itr = FEDsAndHists_.begin(); itr != FEDsAndHists_.end(); ++itr) {
     string dir = fedMap_->getSliceFromFed(itr->first);
     TDirectory* FEDdir = gDirectory->mkdir(dir.c_str());
     FEDdir->cd();
 
     TH1F* hist = itr->second;
     if (hist != nullptr)
       hist->Write();
     else {
       cerr << "EcalPedHists: Error: This shouldn't happen!" << endl;
     }
     // Write out timing hist
     hist = FEDsAndTimingHists_[itr->first];
     if (hist != nullptr)
       hist->Write();
     else {
       cerr << "EcalPedHists: Error: This shouldn't happen!" << endl;
     }
     root_file_.cd();
   }
   allFedsHist_->Write();
   allFedsTimingHist_->Write();
   root_file_.Close();
 
   std::string channels;
   for (std::vector<int>::const_iterator itr = maskedChannels_.begin(); itr != maskedChannels_.end(); ++itr) {
     channels += intToString(*itr);
     channels += ",";
   }
 
   LogWarning("EcalMipGraphs") << "Masked channels are: " << channels << " and that is all!";
 }

◆ endRun()

void EcalURecHitHists::endRun	(	edm::Run const &	,
		edm::EventSetup const &	c
	)

overrideprivate

Definition at line 203 of file EcalURecHitHists.cc.

203 {}

◆ initHists()

void EcalURecHitHists::initHists ( int FED )

private

Definition at line 176 of file EcalURecHitHists.cc.

References fedMap_, FEDsAndHists_, FEDsAndTimingHists_, EcalFedMap::getSliceFromFed(), compareTotals::hist, histRangeMax_, histRangeMin_, createfilelist::int, and intToString().

Referenced by analyze().

                                         {
   using namespace std;
 
   string FEDid = intToString(FED);
   string title1 = "Uncalib Rec Hits (ADC counts) for ";
   title1.append(fedMap_->getSliceFromFed(FED));
   string name1 = "URecHitsFED";
   name1.append(intToString(FED));
   int numBins = (int)round(histRangeMax_ - histRangeMin_) + 1;
   TH1F* hist = new TH1F(name1.c_str(), title1.c_str(), numBins, histRangeMin_, histRangeMax_);
   FEDsAndHists_[FED] = hist;
   FEDsAndHists_[FED]->SetDirectory(nullptr);
 
   title1 = "Jitter for ";
   title1.append(fedMap_->getSliceFromFed(FED));
   name1 = "JitterFED";
   name1.append(intToString(FED));
   TH1F* timingHist = new TH1F(name1.c_str(), title1.c_str(), 14, -7, 7);
   FEDsAndTimingHists_[FED] = timingHist;
   FEDsAndTimingHists_[FED]->SetDirectory(nullptr);
 }

◆ intToString()

std::string EcalURecHitHists::intToString ( int num )

private

Definition at line 245 of file EcalURecHitHists.cc.

References EgammaValidation_cff::num.

Referenced by endJob(), and initHists().

                                                {
   using namespace std;
   ostringstream myStream;
   myStream << num << flush;
   return (myStream.str());  //returns the string form of the stringstream object
 }

Member Data Documentation

◆ allFedsHist_

TH1F* EcalURecHitHists::allFedsHist_

private

Definition at line 86 of file EcalURecHitHists.h.

Referenced by analyze(), EcalURecHitHists(), and endJob().

◆ allFedsTimingHist_

TH1F* EcalURecHitHists::allFedsTimingHist_

private

Definition at line 87 of file EcalURecHitHists.h.

Referenced by analyze(), EcalURecHitHists(), and endJob().

◆ ebUncalibratedRecHitCollection_

const edm::InputTag EcalURecHitHists::ebUncalibratedRecHitCollection_

private

Definition at line 69 of file EcalURecHitHists.h.

◆ ebUncalibRecHitsToken_

const edm::EDGetTokenT<EcalUncalibratedRecHitCollection> EcalURecHitHists::ebUncalibRecHitsToken_

private

Definition at line 72 of file EcalURecHitHists.h.

Referenced by analyze().

◆ ecalElectronicsMap_

const EcalElectronicsMapping* EcalURecHitHists::ecalElectronicsMap_

private

Definition at line 91 of file EcalURecHitHists.h.

Referenced by analyze(), and beginRun().

◆ ecalMappingToken_

const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalURecHitHists::ecalMappingToken_

private

Definition at line 74 of file EcalURecHitHists.h.

Referenced by beginRun().

◆ eeUncalibratedRecHitCollection_

const edm::InputTag EcalURecHitHists::eeUncalibratedRecHitCollection_

private

Definition at line 70 of file EcalURecHitHists.h.

◆ eeUncalibRecHitsToken_

const edm::EDGetTokenT<EcalUncalibratedRecHitCollection> EcalURecHitHists::eeUncalibRecHitsToken_

private

Definition at line 73 of file EcalURecHitHists.h.

Referenced by analyze().

◆ fedMap_

EcalFedMap* EcalURecHitHists::fedMap_

private

Definition at line 90 of file EcalURecHitHists.h.

Referenced by EcalURecHitHists(), endJob(), and initHists().

◆ FEDsAndHists_

std::map<int, TH1F*> EcalURecHitHists::FEDsAndHists_

private

Definition at line 83 of file EcalURecHitHists.h.

Referenced by analyze(), endJob(), and initHists().

◆ FEDsAndTimingHists_

std::map<int, TH1F*> EcalURecHitHists::FEDsAndTimingHists_

private

Definition at line 84 of file EcalURecHitHists.h.

Referenced by analyze(), endJob(), and initHists().

◆ file

TFile* EcalURecHitHists::file

private

Definition at line 89 of file EcalURecHitHists.h.

Referenced by python.rootplot.utilities.RootFile::cd(), ztee.GZipLog::finish(), ztee.GZipLog::flush_block(), python.rootplot.utilities.RootFile::get(), tfile.TFileService::stop(), core.VertexHistograms.VertexHistograms::write(), core.TreeAnalyzerNumpy.TreeAnalyzerNumpy::write(), and ztee.GZipLog::write_block().

◆ fileName_

std::string EcalURecHitHists::fileName_

private

Definition at line 78 of file EcalURecHitHists.h.

Referenced by endJob().

◆ histRangeMax_

const double EcalURecHitHists::histRangeMax_

private

Definition at line 77 of file EcalURecHitHists.h.

Referenced by EcalURecHitHists(), and initHists().

◆ histRangeMin_

const double EcalURecHitHists::histRangeMin_

private

Definition at line 77 of file EcalURecHitHists.h.

Referenced by EcalURecHitHists(), and initHists().

◆ maskedChannels_

std::vector<int> EcalURecHitHists::maskedChannels_

private

Definition at line 80 of file EcalURecHitHists.h.

Referenced by analyze(), EcalURecHitHists(), and endJob().

◆ maskedEBs_

std::vector<std::string> EcalURecHitHists::maskedEBs_

private

Definition at line 82 of file EcalURecHitHists.h.

Referenced by EcalURecHitHists().

◆ maskedFEDs_

std::vector<int> EcalURecHitHists::maskedFEDs_

private

Definition at line 81 of file EcalURecHitHists.h.

Referenced by analyze(), and EcalURecHitHists().

◆ runNum_

int EcalURecHitHists::runNum_

private

Definition at line 76 of file EcalURecHitHists.h.

Referenced by analyze(), and endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalURecHitHists()

◆ ~EcalURecHitHists()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ initHists()

◆ intToString()

Member Data Documentation

◆ allFedsHist_

◆ allFedsTimingHist_

◆ ebUncalibratedRecHitCollection_

◆ ebUncalibRecHitsToken_

◆ ecalElectronicsMap_

◆ ecalMappingToken_

◆ eeUncalibratedRecHitCollection_

◆ eeUncalibRecHitsToken_

◆ fedMap_

◆ FEDsAndHists_

◆ FEDsAndTimingHists_

◆ file

◆ fileName_

◆ histRangeMax_

◆ histRangeMin_

◆ maskedChannels_

◆ maskedEBs_

◆ maskedFEDs_

◆ runNum_