#include <EcalURecHitHists.h>

Inheritance diagram for EcalURecHitHists:

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

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

ESProxyIndex const *	esGetTokenIndices (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::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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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	initHists (int)

std::string	intToString (int num)

Private Attributes
TH1F *	allFedsHist_

TH1F *	allFedsTimingHist_

edm::InputTag	EBUncalibratedRecHitCollection_

const EcalElectronicsMapping *	ecalElectronicsMap_

edm::InputTag	EEUncalibratedRecHitCollection_

EcalFedMap *	fedMap_

std::map< int, TH1F * >	FEDsAndHists_

std::map< int, TH1F * >	FEDsAndTimingHists_

TFile *	file

std::string	fileName_

double	histRangeMax_

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

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

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

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

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

Definition at line 54 of file EcalURecHitHists.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file EcalURecHitHists.cc.

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

     : EBUncalibratedRecHitCollection_(iConfig.getParameter<edm::InputTag>("EBUncalibratedRecHitCollection")),
       EEUncalibratedRecHitCollection_(iConfig.getParameter<edm::InputTag>("EEUncalibratedRecHitCollection")),
       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 ( )

override

Definition at line 74 of file EcalURecHitHists.cc.

74 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 81 of file EcalURecHitHists.cc.

References allFedsHist_, allFedsTimingHist_, EcalUncalibratedRecHit::amplitude(), edm::SortedCollection< T, SORT >::begin(), EcalElectronicsId::dccId(), EBUncalibratedRecHitCollection_, ecalElectronicsMap_, EEUncalibratedRecHitCollection_, edm::SortedCollection< T, SORT >::end(), edm::EventID::event(), FEDsAndHists_, FEDsAndTimingHists_, spr::find(), edm::Event::getByLabel(), EcalElectronicsMapping::getElectronicsId(), ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), EBDetId::ic(), EEDetId::ic(), EcalUncalibratedRecHit::id(), edm::EventBase::id(), initHists(), EcalUncalibratedRecHit::jitter(), LogDebug, maskedChannels_, maskedFEDs_, mps_fire::result, edm::EventID::run(), runNum_, and edm::SortedCollection< T, SORT >::size().

                                                                                   {
   int ievt = iEvent.id().event();
   Handle<EcalUncalibratedRecHitCollection> EBhits;
   Handle<EcalUncalibratedRecHitCollection> EEhits;
 
   iEvent.getByLabel(EBUncalibratedRecHitCollection_, EBhits);
   LogDebug("EcalURecHitHists") << "event " << ievt << " hits collection size " << EBhits->size();
 
   iEvent.getByLabel(EEUncalibratedRecHitCollection_, 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();
   }
 }

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 196 of file EcalURecHitHists.cc.

References ecalElectronicsMap_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                                      {
   edm::ESHandle<EcalElectronicsMapping> elecHandle;
   c.get<EcalMappingRcd>().get(elecHandle);
   ecalElectronicsMap_ = elecHandle.product();
 }

void EcalURecHitHists::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 203 of file EcalURecHitHists.cc.

References allFedsHist_, allFedsTimingHist_, MessageLogger_cfi::cerr, dir, fedMap_, FEDsAndHists_, FEDsAndTimingHists_, fileName_, EcalFedMap::getSliceFromFed(), create_public_lumi_plots::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!";
 }

void EcalURecHitHists::initHists ( int FED )

private

Definition at line 173 of file EcalURecHitHists.cc.

References fedMap_, FEDsAndHists_, FEDsAndTimingHists_, EcalFedMap::getSliceFromFed(), create_public_lumi_plots::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);
 }

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

private

Definition at line 242 of file EcalURecHitHists.cc.

Referenced by endJob(), and initHists().

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