#include <EcalDigiDisplay.h>

Inheritance diagram for EcalDigiDisplay:

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

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	readEBDigis (edm::Handle< EBDigiCollection > digis, int Mode)

void	readEEDigis (edm::Handle< EEDigiCollection > digis, int Mode)

void	readPNDigis (edm::Handle< EcalPnDiodeDigiCollection > PNs, int Mode)

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 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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

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)

Protected Attributes
bool	cryDigi

const std::string	digiProducer_

const std::string	ebDigiCollection_

const edm::EDGetTokenT < EBDigiCollection >	ebDigiToken_

const EcalElectronicsMapping *	ecalElectronicsMap_

const edm::ESGetToken < EcalElectronicsMapping, EcalMappingRcd >	ecalMappingToken_

const std::string	eeDigiCollection_

const edm::EDGetTokenT < EEDigiCollection >	eeDigiToken_

EcalFedMap *	fedMap

bool	inputIsOk

std::vector< int >	listChannels

std::vector< int >	listPns

std::vector< int >	listTowers

int	mode

bool	pnDigi

const edm::EDGetTokenT < EcalPnDiodeDigiCollection >	pnDiodeDigiToken_

const edm::EDGetTokenT < EcalRawDataCollection >	rawDataToken_

std::vector< std::string >	requestedEbs_

std::vector< int >	requestedFeds_

bool	ttDigi

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

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

void	endJob () override

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

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 20 of file EcalDigiDisplay.h.

Constructor & Destructor Documentation

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

Definition at line 25 of file EcalDigiDisplay.cc.

References cryDigi, fedMap, EcalFedMap::getFedFromSlice(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, inputIsOk, listChannels, listPns, listTowers, mode, pnDigi, requestedEbs_, requestedFeds_, and ttDigi.

     : ebDigiCollection_(ps.getParameter<std::string>("ebDigiCollection")),
       eeDigiCollection_(ps.getParameter<std::string>("eeDigiCollection")),
       digiProducer_(ps.getParameter<std::string>("digiProducer")),
       rawDataToken_(consumes<EcalRawDataCollection>(edm::InputTag(digiProducer_))),
       ebDigiToken_(consumes<EBDigiCollection>(edm::InputTag(digiProducer_, ebDigiCollection_))),
       eeDigiToken_(consumes<EEDigiCollection>(edm::InputTag(digiProducer_, eeDigiCollection_))),
       pnDiodeDigiToken_(consumes<EcalPnDiodeDigiCollection>(edm::InputTag(digiProducer_))),
       ecalMappingToken_(esConsumes<edm::Transition::BeginRun>()) {
   //=========================================================================
   requestedFeds_ = ps.getUntrackedParameter<std::vector<int> >("requestedFeds");
   requestedEbs_ = ps.getUntrackedParameter<std::vector<std::string> >("requestedEbs");
 
   cryDigi = ps.getUntrackedParameter<bool>("cryDigi");
   ttDigi = ps.getUntrackedParameter<bool>("ttDigi");
   pnDigi = ps.getUntrackedParameter<bool>("pnDigi");
 
   mode = ps.getUntrackedParameter<int>("mode");
   listChannels = ps.getUntrackedParameter<std::vector<int> >("listChannels");
   listTowers = ps.getUntrackedParameter<std::vector<int> >("listTowers");
   listPns = ps.getUntrackedParameter<std::vector<int> >("listPns");
 
   std::vector<int> listDefaults;
   listDefaults.push_back(-1);
   requestedFeds_ = ps.getUntrackedParameter<std::vector<int> >("requestedFeds", listDefaults);
   bool fedIsGiven = false;
 
   std::vector<std::string> ebDefaults;
   ebDefaults.push_back("none");
   requestedEbs_ = ps.getUntrackedParameter<std::vector<std::string> >("requestedEbs", ebDefaults);
   // FEDs and EBs
   if (requestedFeds_[0] != -1) {
     edm::LogInfo("EcalDigiDisplay") << "FED id is given! Goining to beginRun! ";
     fedIsGiven = true;
   } else {
     if (requestedEbs_[0] != "none") {
       //EB id is given and convert to FED id
       requestedFeds_.clear();
       fedMap = new EcalFedMap();
       for (std::vector<std::string>::const_iterator ebItr = requestedEbs_.begin(); ebItr != requestedEbs_.end();
            ++ebItr) {
         requestedFeds_.push_back(fedMap->getFedFromSlice(*ebItr));
       }
       delete fedMap;
     } else {
       //Select all FEDs in the Event
       for (int i = 601; i < 655; ++i) {
         requestedFeds_.push_back(i);
       }
     }
   }
 
   //Channel list
   listChannels = ps.getUntrackedParameter<std::vector<int> >("listChannels", listDefaults);
   //Tower list
   listTowers = ps.getUntrackedParameter<std::vector<int> >("listTowers", listDefaults);
 
   //Consistancy checks:
   std::vector<int>::iterator fedIter;
   std::vector<int>::iterator intIter;
   inputIsOk = true;
 
   if (fedIsGiven) {
     for (fedIter = requestedFeds_.begin(); fedIter != requestedFeds_.end(); ++fedIter) {
       if ((*fedIter) < 601 || (*fedIter) > 655) {
         edm::LogError("EcalDigiDisplay") << " FED value: " << (*fedIter) << " found in requetsedFeds. "
                                          << " Valid range is 601-654. Returning.";
         inputIsOk = false;
         return;
       }  //Loop over requetsed FEDS
     }
   }
   bool barrelSM = false;
   //Loop over and Check if Barrel SM is picked up
   for (fedIter = requestedFeds_.begin(); fedIter != requestedFeds_.end(); ++fedIter) {
     if ((*fedIter) > 609 && (*fedIter) < 646 && inputIsOk)  // if EB SM is being picked up
       barrelSM = true;
   }
 
   if (barrelSM) {
     if (cryDigi) {
       // Check with channels in Barrel
       for (intIter = listChannels.begin(); intIter != listChannels.end(); intIter++) {
         if (((*intIter) < 1) || (1700 < (*intIter))) {
           edm::LogError("EcalDigiDisplay") << " ic value: " << (*intIter) << " found in listChannels. "
                                            << " Valid range is 1-1700. Returning.";
           inputIsOk = false;
           return;
         }
       }
     }
     //Check with Towers in Barrel
     if (ttDigi) {
       for (intIter = listTowers.begin(); intIter != listTowers.end(); intIter++) {
         if (((*intIter) < 1) || (70 < (*intIter))) {
           edm::LogError("EcalDigiDisplay") << " TT value: " << (*intIter) << " found in listTowers. "
                                            << " Valid range for EB SM is 1-70. Returning.";
           inputIsOk = false;
           return;
         }
       }
     }
   } else  //if EE DCC is being picked up
       if (ttDigi) {
     //Check with Towers in Endcap
     for (intIter = listTowers.begin(); intIter != listTowers.end(); intIter++) {
       if ((*intIter) > 34) {
         edm::LogError("EcalDigiDisplay") << " TT value: " << (*intIter) << " found in listTowers. "
                                          << " Valid range for EE DCC is 1-34. Returning.";
         inputIsOk = false;
         return;
       }
     }
   }
 
   //PNs
   listPns = ps.getUntrackedParameter<std::vector<int> >("listPns", listDefaults);
   /*
   if ( listPns[0] != -1 ) pnDigi = true;
   else {
     listPns.clear();
     for ( int i=1; i < 11; ++i ) {
       listPns.push_back(i);
     }
   }
   */
   if (pnDigi) {
     for (intIter = listPns.begin(); intIter != listPns.end(); intIter++) {
       if (((*intIter) < 1) || (10 < (*intIter))) {
         edm::LogError("EcalDigiDisplay") << " Pn number : " << (*intIter) << " found in listPns. "
                                          << " Valid range is 1-10. Returning.";
         inputIsOk = false;
         return;
       }
     }
   }
 }

EcalDigiDisplay::~EcalDigiDisplay ( )

override

Definition at line 163 of file EcalDigiDisplay.cc.

                                   {
   //=========================================================================
   //delete *;
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 181 of file EcalDigiDisplay.cc.

References cryDigi, ebDigiToken_, eeDigiToken_, edm::Event::getByToken(), inputIsOk, edm::HandleBase::isValid(), mode, pnDigi, pnDiodeDigiToken_, rawDataToken_, readEBDigis(), readEEDigis(), readPNDigis(), and ttDigi.

                                                                        {
   //========================================================================
 
   if (!inputIsOk)
     return;
 
   //Get DCC headers
   edm::Handle<EcalRawDataCollection> dccHeader;
   e.getByToken(rawDataToken_, dccHeader);
   if (!dccHeader.isValid()) {
     edm::LogError("EcalDigiUnpackerModule") << "Can't get DCC Headers!";
   }
 
   //
   bool ebDigisFound = false;
   bool eeDigisFound = false;
   bool pnDigisFound = false;
   // retrieving crystal data from Event
   edm::Handle<EBDigiCollection> eb_digis;
   e.getByToken(ebDigiToken_, eb_digis);
   if (!eb_digis.isValid()) {
     edm::LogError("EcalDigiUnpackerModule") << "EB Digis were not found!";
   }
   if (!eb_digis->empty())
     ebDigisFound = true;
 
   //
   edm::Handle<EEDigiCollection> ee_digis;
   e.getByToken(eeDigiToken_, ee_digis);
   if (!ee_digis.isValid()) {
     edm::LogError("EcalDigiUnpackerModule") << "EE Digis were not found!";
   }
   if (!ee_digis->empty())
     eeDigisFound = true;
 
   // retrieving crystal PN diodes from Event
   edm::Handle<EcalPnDiodeDigiCollection> pn_digis;
   e.getByToken(pnDiodeDigiToken_, pn_digis);
   if (!pn_digis.isValid()) {
     edm::LogError("EcalDigiUnpackerModule") << "PNs were not found!";
   }
   if (!pn_digis->empty())
     pnDigisFound = true;
 
   //=============================
   //Call for funcitons
   //=============================
   if (cryDigi || ttDigi) {
     if (ebDigisFound)
       readEBDigis(eb_digis, mode);
     if (eeDigisFound)
       readEEDigis(ee_digis, mode);
     if (!(ebDigisFound || eeDigisFound)) {
       edm::LogWarning("EcalDigiUnpackerModule") << "No Digis were found! Returning..";
       return;
     }
   }
   if (pnDigi) {
     if (pnDigisFound)
       readPNDigis(pn_digis, mode);
   }
 }

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

overrideprivate

Definition at line 169 of file EcalDigiDisplay.cc.

References ecalElectronicsMap_, ecalMappingToken_, and edm::EventSetup::getData().

                                                                     {
   //========================================================================
   edm::LogInfo("EcalDigiDisplay") << "entering beginRun! ";
 
   ecalElectronicsMap_ = &c.getData(ecalMappingToken_);
 }

void EcalDigiDisplay::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 394 of file EcalDigiDisplay.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

                              {
   //==================================================
   edm::LogInfo("EcalDigiDisplay") << "DONE!.... ";
 }

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

overrideprivate

Definition at line 177 of file EcalDigiDisplay.cc.

177 {}

void EcalDigiDisplay::readEBDigis	(	edm::Handle< EBDigiCollection >	digis,
		int	Mode
	)

protected

Definition at line 248 of file EcalDigiDisplay.cc.

References gather_cfg::cout, EcalElectronicsId::dccId(), isotrackApplyRegressor::df, ecalElectronicsMap_, spr::find(), EcalElectronicsMapping::getElectronicsId(), mps_fire::i, inputIsOk, listChannels, listTowers, requestedFeds_, EcalDataFrame::sample(), hgcalTowerProducer_cfi::tower, and groupFilesInBlocks::tt.

Referenced by analyze().

                                                                              {
   for (EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
     EBDetId detId = EBDetId((*digiItr).id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
 
     int FEDid = elecId.dccId() + 600;
     std::vector<int>::iterator fedIter = find(requestedFeds_.begin(), requestedFeds_.end(), FEDid);
     if (fedIter == requestedFeds_.end())
       continue;
 
     int ic = EBDetId((*digiItr).id()).ic();
     int tt = EBDetId((*digiItr).id()).tower().iTT();
 
     //Check if Mode is set 1 or 2
     if (Mode == 1) {
       edm::LogInfo("EcalDigiDisplay") << "\n\n^^^^^^^^^^^^^^^^^^ [EcalDigiDisplay]  digi cry collection size "
                                       << digis->size();
       edm::LogInfo("EcalDigiDisplay") << "                       [EcalDigiDisplay]  dumping first " << listChannels[0]
                                       << " crystals\n";
       //It will break if all required digis are dumpped
       if (ic > listChannels[0])
         continue;
     } else if (Mode == 2) {
       std::vector<int>::iterator icIterCh;
       std::vector<int>::iterator icIterTt;
       icIterCh = find(listChannels.begin(), listChannels.end(), ic);
       icIterTt = find(listTowers.begin(), listTowers.end(), tt);
       if (icIterCh == listChannels.end() && icIterTt == listTowers.end())
         continue;
       edm::LogInfo("EcalDigiDisplay") << "\n\n^^^^^^^^^^^^^^^^^^ [EcalDigiDisplay]  digi cry collection size "
                                       << digis->size();
     } else {
       edm::LogInfo("EcalDigiDisplay") << "[EcalDigiDisplay] parameter mode set to: " << Mode
                                       << ". Only mode 1 and 2 are allowed. Returning...";
       inputIsOk = false;
       return;
     }
     std::cout << "FEDID: " << FEDid << std::endl;
     std::cout << "Tower: " << EBDetId((*digiItr).id()).tower().iTT() << " ic-cry: " << EBDetId((*digiItr).id()).ic()
               << " i-phi: " << EBDetId((*digiItr).id()).iphi() << " j-eta: " << EBDetId((*digiItr).id()).ieta()
               << std::endl;
     //Get Samples
     for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
       EBDataFrame df(*digiItr);
       if (!(i % 3))
         std::cout << "\n\t";
       std::cout << "sId: " << (i + 1) << " " << df.sample(i) << "\t";
     }
     std::cout << " " << std::endl;
   }
 }

void EcalDigiDisplay::readEEDigis	(	edm::Handle< EEDigiCollection >	digis,
		int	Mode
	)

protected

Definition at line 301 of file EcalDigiDisplay.cc.

References gather_cfg::cout, EcalElectronicsId::dccId(), isotrackApplyRegressor::df, ecalElectronicsMap_, spr::find(), EcalElectronicsMapping::getElectronicsId(), mps_fire::i, listChannels, listTowers, requestedFeds_, and EcalDataFrame::sample().

Referenced by analyze().

                                                                              {
   //For Endcap so far works only  Mode 2
   if (Mode != 2) {
     std::cout << "For Endcap mode needs to be set to 2" << std::endl;
     return;
   }
 
   for (EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
     //Make sure that digis are form requested place
     EEDetId detId = EEDetId((*digiItr).id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
 
     int FEDid = elecId.dccId() + 600;
     std::vector<int>::iterator fedIter = find(requestedFeds_.begin(), requestedFeds_.end(), FEDid);
     if (fedIter == requestedFeds_.end())
       continue;
 
     edm::LogInfo("EcalDigiDisplay") << "\n\n^^^^^^^^^^^^^^^^^^ [EcalDigiDisplay]  digi cry collection size "
                                     << digis->size();
 
     int crystalId = 10000 * FEDid + 100 * elecId.towerId() + 5 * (elecId.stripId() - 1) + elecId.xtalId();
     int chId = elecId.towerId();  // this is a channel in Endcap DCC, sometimes also called as Super Crystal
 
     std::vector<int>::iterator icIterCh;
     std::vector<int>::iterator icIterTt;
     icIterCh = find(listChannels.begin(), listChannels.end(), crystalId);
     icIterTt = find(listTowers.begin(), listTowers.end(), chId);
     if (icIterCh == listChannels.end() && icIterTt == listTowers.end())
       continue;
 
     std::cout << "FEDID: " << FEDid << std::endl;
     std::cout << "Tower: " << elecId.towerId() << "crystalId: " << crystalId
               << " i-x: " << EEDetId((*digiItr).id()).ix() << " j-y: " << EEDetId((*digiItr).id()).iy() << std::endl;
 
     //Get samples
     for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
       EEDataFrame df(*digiItr);
       if (!(i % 3))
         std::cout << "\n\t";
       std::cout << "sId: " << (i + 1) << " " << df.sample(i) << "\t";
     }
     std::cout << " " << std::endl;
   }
 }

void EcalDigiDisplay::readPNDigis	(	edm::Handle< EcalPnDiodeDigiCollection >	PNs,
		int	Mode
	)

protected

Definition at line 346 of file EcalDigiDisplay.cc.

References gather_cfg::cout, spr::find(), EcalPnDiodeDetId::iDCCId(), inputIsOk, listPns, and requestedFeds_.

Referenced by analyze().

                                                                                     {
   int pnDigiCounter = 0;
 
   //Loop over PN digis
   for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNs->begin(); pnItr != PNs->end(); ++pnItr) {
     EcalPnDiodeDetId pnDetId = EcalPnDiodeDetId((*pnItr).id());
     //Make sure that we look at the requested place
     int FEDid = pnDetId.iDCCId() + 600;
     std::vector<int>::iterator fedIter = find(requestedFeds_.begin(), requestedFeds_.end(), FEDid);
     if (fedIter == requestedFeds_.end())
       continue;
     int pnNum = (*pnItr).id().iPnId();
 
     if (Mode == 1) {
       edm::LogInfo("EcalDigiDisplay") << "\n\n^^^^^^^^^^^^^^^^^^ EcalDigiDisplay  digi PN collection.  Size: "
                                       << PNs->size();
       edm::LogInfo("EcalDigiDisplay") << "                       [EcalDigiDisplay]  dumping first " << listPns[0]
                                       << " PNs ";
 
       if ((pnDigiCounter++) >= listPns[0])
         break;
     } else if (Mode == 2) {
       edm::LogInfo("EcalDigiDisplay") << "\n\n^^^^^^^^^^^^^^^^^^ EcalDigiDisplay  digi PN collection.  Size: "
                                       << PNs->size();
 
       // Check that we look at PN from the given list
       std::vector<int>::iterator pnIter;
       pnIter = find(listPns.begin(), listPns.end(), pnNum);
       if (pnIter == listPns.end())
         continue;
     } else {
       edm::LogError("EcalDigiDisplay") << "[EcalDigiDisplay] parameter mode set to: " << Mode
                                        << ". Only mode 1 and 2 are allowed. Returning...";
       inputIsOk = false;
       return;
     }
 
     std::cout << "DCCID: " << pnDetId.iDCCId() << std::endl;
     std::cout << "\nPN num: " << (*pnItr).id().iPnId();
     for (int samId = 0; samId < (*pnItr).size(); samId++) {
       if (!(samId % 3))
         std::cout << "\n\t";
       std::cout << "sId: " << (samId + 1) << " " << (*pnItr).sample(samId) << "\t";
     }
   }
 }