#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 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 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 HltBtagPostValidation_cff::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_