#include <EcalPedHists.h>

Inheritance diagram for EcalPedHists:

Public Member Functions
	EcalPedHists (const edm::ParameterSet &ps)

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

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

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

void	endJob (void) override

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

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)

Private Member Functions
void	initHists (int FED)

std::string	intToString (int num)

void	readEBdigis (edm::Handle< EBDigiCollection > digis)

void	readEEdigis (edm::Handle< EEDigiCollection > digis)

Private Attributes
bool	allFEDsSelected_

const edm::InputTag	barrelDigiCollection_

const edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_

const EcalElectronicsMapping *	ecalElectronicsMap_

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	ecalMappingToken_

const edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

const edm::InputTag	endcapDigiCollection_

EcalFedMap *	fedMap_

std::map< int, stringHistMap >	FEDsAndHistMaps_

std::string	fileName_

const edm::InputTag	headerProducer_

bool	histsFilled_

bool	inputIsOk_

std::vector< int >	listChannels_

std::vector< std::string >	listEBs_

std::vector< int >	listFEDs_

std::vector< int >	listSamples_

const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_

TFile *	root_file_

int	runNum_

std::set< int >	theRealFedSet_

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)

Detailed Description

Definition at line 35 of file EcalPedHists.h.

Constructor & Destructor Documentation

◆ EcalPedHists()

EcalPedHists::EcalPedHists ( const edm::ParameterSet & ps )

Module which outputs a root file of ADC counts (all three gains) of user-selected channels (defaults to channel 1) for user-selected samples (defaults to samples 1,2,3) for user-selected supermodules.

Author: S. Cooper

Definition at line 13 of file EcalPedHists.cc.

References allFEDsSelected_, gather_cfg::cout, fedMap_, EcalFedMap::getFedFromSlice(), edm::ParameterSet::getUntrackedParameter(), histsFilled_, mps_fire::i, inputIsOk_, listChannels_, listEBs_, listFEDs_, listSamples_, and theRealFedSet_.

     : runNum_(-1),
       fileName_(ps.getUntrackedParameter<std::string>("fileName", std::string("ecalPedDigiDump"))),
       barrelDigiCollection_(ps.getParameter<edm::InputTag>("EBdigiCollection")),
       endcapDigiCollection_(ps.getParameter<edm::InputTag>("EEdigiCollection")),
       headerProducer_(ps.getParameter<edm::InputTag>("headerProducer")),
       rawDataToken_(consumes<EcalRawDataCollection>(headerProducer_)),
       ebDigiToken_(consumes<EBDigiCollection>(barrelDigiCollection_)),
       eeDigiToken_(consumes<EEDigiCollection>(endcapDigiCollection_)),
       ecalMappingToken_(esConsumes<edm::Transition::BeginRun>()) {
   using namespace std;
 
   fedMap_ = new EcalFedMap();
   histsFilled_ = false;
   //for(int i=601; i<655; ++i)
   //{
   //  listDefaults.push_back(i);
   //}
   listFEDs_ = ps.getUntrackedParameter<vector<int> >("listFEDs");
   listEBs_ = ps.getUntrackedParameter<vector<string> >("listEBs");
 
   if (listFEDs_.empty()) {
     allFEDsSelected_ = false;
     //if "actual" EB id given, then convert to FEDid and put in listFEDs_
     if (!listEBs_.empty()) {
       listFEDs_.clear();
       for (vector<string>::const_iterator itr = listEBs_.begin(); itr != listEBs_.end(); ++itr) {
         listFEDs_.push_back(fedMap_->getFedFromSlice(*itr));
       }
     }
   } else if (listFEDs_[0] == -1) {
     // Apply no selection if -1 is passed in FED list
     allFEDsSelected_ = true;
     //debug
     //cout << "no selection on FEDs!" << endl;
     //inputIsOk_=false;
     //return;
     //listFEDs_ = listDefaults;
   } else {
     //in this case, listFEDs should be populated
     allFEDsSelected_ = false;
   }
 
   if (!allFEDsSelected_) {
     // Verify FED numbers are valid
     for (vector<int>::const_iterator intIter = listFEDs_.begin(); intIter != listFEDs_.end(); intIter++) {
       if (((*intIter) < 601) || (654 < (*intIter))) {
         cout << "[EcalPedHists] FED value: " << (*intIter) << " found in listFEDs. "
              << " Valid range is 601-654. Returning." << endl;
         inputIsOk_ = false;
         return;
       } else
         theRealFedSet_.insert(*intIter);
     }
   }
 
   vector<int> listDefaults = vector<int>();
   listDefaults.clear();
   for (int i = 1; i < 1701; ++i) {
     listDefaults.push_back(i);
   }
   listChannels_ = ps.getUntrackedParameter<vector<int> >("listChannels", listDefaults);
   listDefaults.clear();
   // Get samples to plot (default to 1,2,3)
   listDefaults.push_back(0);
   listDefaults.push_back(1);
   listDefaults.push_back(2);
   listSamples_ = ps.getUntrackedParameter<vector<int> >("listSamples", listDefaults);
 
   inputIsOk_ = true;
   vector<int>::iterator intIter;
 
   // Verify crystal numbers are valid
   for (intIter = listChannels_.begin(); intIter != listChannels_.end(); ++intIter) {
     //TODO: Fix crystal index checking?
     //if ( ((*intIter) < 1)||(1700 < (*intIter)) )
     //{
     //  cout << "[EcalPedHists] ic value: " << (*intIter) << " found in listChannels. "
     //        << " Valid range is 1-1700. Returning." << endl;
     //  inputIsOk_ = false;
     //  return;
     //}
   }
   // Verify sample numbers are valid
   for (intIter = listSamples_.begin(); intIter != listSamples_.end(); intIter++) {
     if (((*intIter) < 1) || (10 < (*intIter))) {
       cout << "[EcalPedHists] sample number: " << (*intIter) << " found in listSamples. "
            << " Valid range is 1-10. Returning." << endl;
       inputIsOk_ = false;
       return;
     }
   }
 }

◆ ~EcalPedHists()

EcalPedHists::~EcalPedHists ( )

override

Definition at line 107 of file EcalPedHists.cc.

107 {}

Member Function Documentation

◆ analyze()

void EcalPedHists::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 174 of file EcalPedHists.cc.

References allFEDsSelected_, barrelDigiCollection_, mtdUncalibratedRecHits_cfi::barrelDigis, edm::SortedCollection< T, SORT >::begin(), MillePedeFileConverter_cfg::e, ebDigiToken_, eeDigiToken_, edm::SortedCollection< T, SORT >::end(), endcapDigiCollection_, mtdUncalibratedRecHits_cfi::endcapDigis, FEDsAndHistMaps_, headerProducer_, initHists(), inputIsOk_, edm::HandleBase::isValid(), rawDataToken_, readEBdigis(), readEEdigis(), runNum_, and theRealFedSet_.

                                                                     {
   using namespace std;
   using namespace edm;
 
   if (!inputIsOk_)
     return;
 
   // loop over the headers, this is to detect missing FEDs if all are selected
   if (allFEDsSelected_) {
     edm::Handle<EcalRawDataCollection> DCCHeaders;
     e.getByToken(rawDataToken_, DCCHeaders);
     if (!DCCHeaders.isValid()) {
       edm::LogError("EcalPedHists") << "Error! can't get the product " << headerProducer_;
       return;
     }
 
     for (EcalRawDataCollection::const_iterator headerItr = DCCHeaders->begin(); headerItr != DCCHeaders->end();
          ++headerItr) {
       int FEDid = 600 + headerItr->id();
       theRealFedSet_.insert(FEDid);
     }
   }
 
   // loop over fed list and make sure that there are histo maps
   for (set<int>::const_iterator fedItr = theRealFedSet_.begin(); fedItr != theRealFedSet_.end(); ++fedItr) {
     if (FEDsAndHistMaps_.find(*fedItr) == FEDsAndHistMaps_.end())
       initHists(*fedItr);
   }
 
   //debug
   //cout << "analyze...input is ok? " << inputIsOk_ << endl;
 
   bool barrelDigisFound = true;
   bool endcapDigisFound = true;
   // get the barrel digis
   // (one digi for each crystal)
   // TODO; SIC: fix this behavior
   Handle<EBDigiCollection> barrelDigis;
   e.getByToken(ebDigiToken_, barrelDigis);
   if (!barrelDigis.isValid()) {
     edm::LogError("EcalPedOffset") << "Error! can't get the product " << barrelDigiCollection_;
     barrelDigisFound = false;
   }
   // get the endcap digis
   // (one digi for each crystal)
   // TODO; SIC: fix this behavior
   Handle<EEDigiCollection> endcapDigis;
   e.getByToken(eeDigiToken_, endcapDigis);
   if (!endcapDigis.isValid()) {
     edm::LogError("EcalPedOffset") << "Error! can't get the product " << endcapDigiCollection_;
     endcapDigisFound = false;
   }
 
   if (barrelDigisFound)
     readEBdigis(barrelDigis);
   if (endcapDigisFound)
     readEEdigis(endcapDigis);
   if (!barrelDigisFound && !endcapDigisFound)
     edm::LogError("EcalPedOffset") << "No digis found in the event!";
 
   if (runNum_ == -1)
     runNum_ = e.id().run();
 }

◆ beginRun()

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

overrideprotected

Definition at line 109 of file EcalPedHists.cc.

References DummyCfis::c, ecalElectronicsMap_, and ecalMappingToken_.

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

◆ endJob()

void EcalPedHists::endJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 115 of file EcalPedHists.cc.

References DMR_cfg::cerr, DeadROC_duringRun::dir, fedMap_, FEDsAndHistMaps_, fileName_, EcalFedMap::getSliceFromFed(), compareTotals::hist, histsFilled_, inputIsOk_, intToString(), listChannels_, hfnoseParametersInitialization_cfi::name2, root_file_, runNum_, and theRealFedSet_.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
   using namespace std;
   if (inputIsOk_) {
     //debug
     //cout << "endJob:creating root file!" << endl;
 
     fileName_ += "-" + intToString(runNum_) + ".graph.root";
 
     TFile root_file_(fileName_.c_str(), "RECREATE");
     //Loop over FEDs first
     for (set<int>::const_iterator FEDitr = theRealFedSet_.begin(); FEDitr != theRealFedSet_.end(); ++FEDitr) {
       if (!histsFilled_)
         break;
       string dir = fedMap_->getSliceFromFed(*FEDitr);
       TDirectory* FEDdir = gDirectory->mkdir(dir.c_str());
       FEDdir->cd();
       //root_file_.mkdir(dir.c_str());
       //root_file_.cd(dir.c_str());
       map<string, TH1F*> mapHistos = FEDsAndHistMaps_[*FEDitr];
 
       //Loop over channels; write histos and directory structure
       for (vector<int>::const_iterator itr = listChannels_.begin(); itr != listChannels_.end(); itr++) {
         //debug
         //cout << "loop over channels" << endl;
 
         TH1F* hist = nullptr;
         string chnl = intToString(*itr);
         string name1 = "Cry";
         name1.append(chnl + "Gain1");
         string name2 = "Cry";
         name2.append(chnl + "Gain6");
         string name3 = "Cry";
         name3.append(chnl + "Gain12");
         hist = mapHistos[name1];
         // This is a sanity check only
         if (hist != nullptr) {
           string cryDirName = "Cry_" + chnl;
           TDirectory* cryDir = FEDdir->mkdir(cryDirName.c_str());
           cryDir->cd();
           hist->SetDirectory(cryDir);
           hist->Write();
           hist = mapHistos[name2];
           hist->SetDirectory(cryDir);
           hist->Write();
           hist = mapHistos[name3];
           hist->SetDirectory(cryDir);
           hist->Write();
           //root_file_.cd(dir.c_str());
           root_file_.cd();
         } else {
           cerr << "EcalPedHists: Error: This shouldn't happen!" << endl;
         }
       }
       root_file_.cd();
     }
     root_file_.Close();
   }
 }

◆ endRun()

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

overrideprotected

Definition at line 113 of file EcalPedHists.cc.

113 {}

◆ initHists()

void EcalPedHists::initHists ( int FED )

private

Definition at line 239 of file EcalPedHists.cc.

References FEDsAndHistMaps_, intToString(), listChannels_, and hfnoseParametersInitialization_cfi::name2.

Referenced by analyze().

                                     {
   using namespace std;
   //using namespace edm;
 
   std::map<string, TH1F*> histMap;
   //debug
   //cout << "Initializing map for FED:" << *FEDitr << endl;
   for (vector<int>::const_iterator intIter = listChannels_.begin(); intIter != listChannels_.end(); ++intIter) {
     //Put 3 histos (1 per gain) for the channel into the map
     string FEDid = intToString(FED);
     string chnl = intToString(*intIter);
     string title1 = "Gain1 ADC Counts for channel ";
     title1.append(chnl);
     string name1 = "Cry";
     name1.append(chnl + "Gain1");
     string title2 = "Gain6 ADC Counts for channel ";
     title2.append(chnl);
     string name2 = "Cry";
     name2.append(chnl + "Gain6");
     string title3 = "Gain12 ADC Counts for channel ";
     title3.append(chnl);
     string name3 = "Cry";
     name3.append(chnl + "Gain12");
     histMap.insert(make_pair(name1, new TH1F(name1.c_str(), title1.c_str(), 75, 175.0, 250.0)));
     histMap[name1]->SetDirectory(nullptr);
     histMap.insert(make_pair(name2, new TH1F(name2.c_str(), title2.c_str(), 75, 175.0, 250.0)));
     histMap[name2]->SetDirectory(nullptr);
     histMap.insert(make_pair(name3, new TH1F(name3.c_str(), title3.c_str(), 75, 175.0, 250.0)));
     histMap[name3]->SetDirectory(nullptr);
   }
   FEDsAndHistMaps_.insert(make_pair(FED, histMap));
 }

◆ intToString()

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

private

Definition at line 377 of file EcalPedHists.cc.

References EgammaValidation_cff::num.

Referenced by endJob(), initHists(), readEBdigis(), and readEEdigis().

                                            {
   using namespace std;
   //
   // outputs the number into the string stream and then flushes
   // the buffer (makes sure the output is put into the stream)
   //
   ostringstream myStream;  //creates an ostringstream object
   myStream << num << flush;
   return (myStream.str());  //returns the string form of the stringstream object
 }

◆ readEBdigis()

void EcalPedHists::readEBdigis ( edm::Handle< EBDigiCollection > digis )

private

Definition at line 272 of file EcalPedHists.cc.

References gpuClustering::adc, edm::DataFrameContainer::begin(), DMR_cfg::cerr, EcalElectronicsId::dccId(), hcalRecHitTable_cff::detId, ecalElectronicsMap_, edm::DataFrameContainer::end(), FEDsAndHistMaps_, spr::find(), ecalLiteDTU::gainId(), EcalElectronicsMapping::getElectronicsId(), compareTotals::hist, histsFilled_, intToString(), listChannels_, listSamples_, hfnoseParametersInitialization_cfi::name2, ecalGpuTask_cfi::sample, and theRealFedSet_.

Referenced by analyze().

                                                                 {
   using namespace std;
   using namespace edm;
   //debug
   //cout << "readEBdigis" << endl;
 
   // Loop over digis
   for (EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
     EBDetId detId = EBDetId(digiItr->id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
     int FEDid = 600 + elecId.dccId();
     int crystalId = detId.ic();
 
     //debug
     //cout << "FEDid:" << FEDid << " cryId:" << crystalId << endl;
     //cout << "FEDid:" << FEDid << endl;
     //Select desired supermodules only
     set<int>::const_iterator fedIter = find(theRealFedSet_.begin(), theRealFedSet_.end(), FEDid);
     if (fedIter == theRealFedSet_.end())
       continue;
 
     // Select desired channels only
     vector<int>::iterator icIter;
     icIter = find(listChannels_.begin(), listChannels_.end(), crystalId);
     if (icIter == listChannels_.end())
       continue;
 
     // Get the adc counts from the selected samples and fill the corresponding histogram
     // Must subtract 1 from user-given sample list (e.g., user's sample 1 -> sample 0)
     for (vector<int>::iterator itr = listSamples_.begin(); itr != listSamples_.end(); itr++) {
       histsFilled_ = true;
       map<string, TH1F*> mapHistos = FEDsAndHistMaps_[FEDid];
       string chnl = intToString(crystalId);
       string name1 = "Cry";
       name1.append(chnl + "Gain1");
       string name2 = "Cry";
       name2.append(chnl + "Gain6");
       string name3 = "Cry";
       name3.append(chnl + "Gain12");
       TH1F* hist = nullptr;
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 3)
         hist = mapHistos[name1];
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 2)
         hist = mapHistos[name2];
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 1)
         hist = mapHistos[name3];
       if (hist != nullptr)
         hist->Fill(((EBDataFrame)(*digiItr)).sample(*itr - 1).adc());
       else
         cerr << "EcalPedHistDumper: Error: This shouldn't happen!" << endl;
     }
   }
 }

◆ readEEdigis()

void EcalPedHists::readEEdigis ( edm::Handle< EEDigiCollection > digis )

private

Definition at line 326 of file EcalPedHists.cc.

References gpuClustering::adc, edm::DataFrameContainer::begin(), DMR_cfg::cerr, EcalElectronicsId::dccId(), hcalRecHitTable_cff::detId, ecalElectronicsMap_, edm::DataFrameContainer::end(), FEDsAndHistMaps_, spr::find(), ecalLiteDTU::gainId(), EcalElectronicsMapping::getElectronicsId(), compareTotals::hist, histsFilled_, intToString(), listChannels_, listSamples_, hfnoseParametersInitialization_cfi::name2, ecalGpuTask_cfi::sample, and theRealFedSet_.

Referenced by analyze().

                                                                 {
   using namespace std;
   using namespace edm;
   //debug
   //cout << "readEEdigis" << endl;
 
   // Loop over digis
   for (EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
     EEDetId detId = EEDetId(digiItr->id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
     int FEDid = 600 + elecId.dccId();
     int crystalId = 10000 * FEDid + 100 * elecId.towerId() + 5 * (elecId.stripId() - 1) + elecId.xtalId();
 
     //Select desired FEDs only
     set<int>::const_iterator fedIter = find(theRealFedSet_.begin(), theRealFedSet_.end(), FEDid);
     if (fedIter == theRealFedSet_.end())
       continue;
 
     // Select desired channels only
     vector<int>::iterator icIter;
     icIter = find(listChannels_.begin(), listChannels_.end(), crystalId);
     if (icIter == listChannels_.end())
       continue;
 
     // Get the adc counts from the selected samples and fill the corresponding histogram
     // Must subtract 1 from user-given sample list (e.g., user's sample 1 -> sample 0)
     for (vector<int>::iterator itr = listSamples_.begin(); itr != listSamples_.end(); itr++) {
       histsFilled_ = true;
       map<string, TH1F*> mapHistos = FEDsAndHistMaps_[FEDid];
       string chnl = intToString(crystalId);
       string name1 = "Cry";
       name1.append(chnl + "Gain1");
       string name2 = "Cry";
       name2.append(chnl + "Gain6");
       string name3 = "Cry";
       name3.append(chnl + "Gain12");
       TH1F* hist = nullptr;
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 3)
         hist = mapHistos[name1];
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 2)
         hist = mapHistos[name2];
       if (((EBDataFrame)(*digiItr)).sample(*itr - 1).gainId() == 1)
         hist = mapHistos[name3];
       if (hist != nullptr)
         hist->Fill(((EBDataFrame)(*digiItr)).sample(*itr - 1).adc());
       else
         cerr << "EcalPedHistDumper: Error: This shouldn't happen!" << endl;
     }
   }
 }

Member Data Documentation

◆ allFEDsSelected_

bool EcalPedHists::allFEDsSelected_

private

Definition at line 54 of file EcalPedHists.h.

Referenced by analyze(), and EcalPedHists().

◆ barrelDigiCollection_

const edm::InputTag EcalPedHists::barrelDigiCollection_

private

Definition at line 57 of file EcalPedHists.h.

Referenced by analyze().

◆ ebDigiToken_

const edm::EDGetTokenT<EBDigiCollection> EcalPedHists::ebDigiToken_

private

Definition at line 70 of file EcalPedHists.h.

Referenced by analyze().

◆ ecalElectronicsMap_

const EcalElectronicsMapping* EcalPedHists::ecalElectronicsMap_

private

Definition at line 74 of file EcalPedHists.h.

Referenced by beginRun(), readEBdigis(), and readEEdigis().

◆ ecalMappingToken_

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

private

Definition at line 72 of file EcalPedHists.h.

Referenced by beginRun().

◆ eeDigiToken_

const edm::EDGetTokenT<EEDigiCollection> EcalPedHists::eeDigiToken_

private

Definition at line 71 of file EcalPedHists.h.

Referenced by analyze().

◆ endcapDigiCollection_

const edm::InputTag EcalPedHists::endcapDigiCollection_

private

Definition at line 58 of file EcalPedHists.h.

Referenced by analyze().

◆ fedMap_

EcalFedMap* EcalPedHists::fedMap_

private

Definition at line 66 of file EcalPedHists.h.

Referenced by EcalPedHists(), and endJob().

◆ FEDsAndHistMaps_

std::map<int, stringHistMap> EcalPedHists::FEDsAndHistMaps_

private

Definition at line 64 of file EcalPedHists.h.

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

◆ fileName_

std::string EcalPedHists::fileName_

private

Definition at line 56 of file EcalPedHists.h.

Referenced by endJob().

◆ headerProducer_

const edm::InputTag EcalPedHists::headerProducer_

private

Definition at line 59 of file EcalPedHists.h.

Referenced by analyze().

◆ histsFilled_

bool EcalPedHists::histsFilled_

private

Definition at line 55 of file EcalPedHists.h.

Referenced by EcalPedHists(), endJob(), readEBdigis(), and readEEdigis().

◆ inputIsOk_

bool EcalPedHists::inputIsOk_

private

Definition at line 53 of file EcalPedHists.h.

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

◆ listChannels_

std::vector<int> EcalPedHists::listChannels_

private

Definition at line 60 of file EcalPedHists.h.

Referenced by EcalPedHists(), endJob(), initHists(), readEBdigis(), and readEEdigis().

◆ listEBs_

std::vector<std::string> EcalPedHists::listEBs_

private

Definition at line 63 of file EcalPedHists.h.

Referenced by EcalPedHists().

◆ listFEDs_

std::vector<int> EcalPedHists::listFEDs_

private

Definition at line 62 of file EcalPedHists.h.

Referenced by EcalPedHists().

◆ listSamples_

std::vector<int> EcalPedHists::listSamples_

private

Definition at line 61 of file EcalPedHists.h.

Referenced by EcalPedHists(), readEBdigis(), and readEEdigis().

◆ rawDataToken_

const edm::EDGetTokenT<EcalRawDataCollection> EcalPedHists::rawDataToken_

private

Definition at line 69 of file EcalPedHists.h.

Referenced by analyze().

◆ root_file_

TFile* EcalPedHists::root_file_

private

Definition at line 67 of file EcalPedHists.h.

Referenced by endJob().

◆ runNum_

int EcalPedHists::runNum_

private

Definition at line 52 of file EcalPedHists.h.

Referenced by analyze(), and endJob().

◆ theRealFedSet_

std::set<int> EcalPedHists::theRealFedSet_

private

Definition at line 65 of file EcalPedHists.h.

Referenced by analyze(), EcalPedHists(), endJob(), readEBdigis(), and readEEdigis().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalPedHists()

◆ ~EcalPedHists()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ initHists()

◆ intToString()

◆ readEBdigis()

◆ readEEdigis()

Member Data Documentation

◆ allFEDsSelected_

◆ barrelDigiCollection_

◆ ebDigiToken_

◆ ecalElectronicsMap_

◆ ecalMappingToken_

◆ eeDigiToken_

◆ endcapDigiCollection_

◆ fedMap_

◆ FEDsAndHistMaps_

◆ fileName_

◆ headerProducer_

◆ histsFilled_

◆ inputIsOk_

◆ listChannels_

◆ listEBs_

◆ listFEDs_

◆ listSamples_

◆ rawDataToken_

◆ root_file_

◆ runNum_

◆ theRealFedSet_