#include <EcalDigiDisplay.h>

Inheritance diagram for EcalDigiDisplay:

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

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

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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)

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

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...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Protected Attributes
bool	cryDigi

std::string	digiProducer_

std::string	ebDigiCollection_

const EcalElectronicsMapping *	ecalElectronicsMap_

std::string	eeDigiCollection_

EcalFedMap *	fedMap

bool	inputIsOk

std::vector< int >	listChannels

std::vector< int >	listPns

std::vector< int >	listTowers

int	mode

bool	pnDigi

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

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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Detailed Description

Definition at line 20 of file EcalDigiDisplay.h.

Constructor & Destructor Documentation

◆ EcalDigiDisplay()

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

Definition at line 33 of file EcalDigiDisplay.cc.

                                                           {
   //=========================================================================
  
   ebDigiCollection_ = ps.getParameter<std::string>("ebDigiCollection");
   eeDigiCollection_ = ps.getParameter<std::string>("eeDigiCollection");
   digiProducer_ = ps.getParameter<std::string>("digiProducer");
  
   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;
       }
     }
   }
 }

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

◆ ~EcalDigiDisplay()

EcalDigiDisplay::~EcalDigiDisplay ( )

override

Definition at line 168 of file EcalDigiDisplay.cc.

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

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 184 of file EcalDigiDisplay.cc.

                                                                        {
   //========================================================================
  
   if (!inputIsOk)
     return;
  
   //Get DCC headers
   edm::Handle<EcalRawDataCollection> dccHeader;
   try {
     e.getByLabel(digiProducer_, dccHeader);
   } catch (cms::Exception& ex) {
     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;
   try {
     e.getByLabel(digiProducer_, ebDigiCollection_, eb_digis);
     if (!eb_digis->empty())
       ebDigisFound = true;
   } catch (cms::Exception& ex) {
     edm::LogError("EcalDigiUnpackerModule") << "EB Digis were not found!";
   }
  
   //
   edm::Handle<EEDigiCollection> ee_digis;
   try {
     e.getByLabel(digiProducer_, eeDigiCollection_, ee_digis);
     if (!ee_digis->empty())
       eeDigisFound = true;
   } catch (cms::Exception& ex) {
     edm::LogError("EcalDigiUnpackerModule") << "EE Digis were not found!";
   }
  
   // retrieving crystal PN diodes from Event
   edm::Handle<EcalPnDiodeDigiCollection> pn_digis;
   try {
     e.getByLabel(digiProducer_, pn_digis);
     if (!pn_digis->empty())
       pnDigisFound = true;
   } catch (cms::Exception& ex) {
     edm::LogError("EcalDigiUnpackerModule") << "PNs were not found!";
   }
  
   //=============================
   //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);
   }
 }

References cryDigi, digiProducer_, MillePedeFileConverter_cfg::e, ebDigiCollection_, eeDigiCollection_, edm::SortedCollection< T, SORT >::empty(), edm::DataFrameContainer::empty(), inputIsOk, mode, pnDigi, readEBDigis(), readEEDigis(), readPNDigis(), and ttDigi.

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 174 of file EcalDigiDisplay.cc.

                                                                     {
   //========================================================================
   edm::LogInfo("EcalDigiDisplay") << "entering beginRun! ";
  
   edm::ESHandle<EcalElectronicsMapping> elecHandle;
   c.get<EcalMappingRcd>().get(elecHandle);
   ecalElectronicsMap_ = elecHandle.product();
 }

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

◆ endJob()

void EcalDigiDisplay::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 401 of file EcalDigiDisplay.cc.

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

Referenced by o2olib.O2ORunMgr::executeJob().

◆ readEBDigis()

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

protected

Definition at line 255 of file EcalDigiDisplay.cc.

                                                                              {
   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;
   }
 }

References edm::DataFrameContainer::begin(), gather_cfg::cout, EcalElectronicsId::dccId(), hgcalPerformanceValidation::df, ecalElectronicsMap_, edm::DataFrameContainer::end(), spr::find(), EcalElectronicsMapping::getElectronicsId(), mps_fire::i, LEDCalibrationChannels::ieta, inputIsOk, LEDCalibrationChannels::iphi, listChannels, listTowers, requestedFeds_, edm::DataFrameContainer::size(), hgcalTowerProducer_cfi::tower, and groupFilesInBlocks::tt.

Referenced by analyze().

◆ readEEDigis()

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

protected

Definition at line 308 of file EcalDigiDisplay.cc.

                                                                              {
   //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;
   }
 }

References edm::DataFrameContainer::begin(), gather_cfg::cout, EcalElectronicsId::dccId(), hgcalPerformanceValidation::df, ecalElectronicsMap_, edm::DataFrameContainer::end(), spr::find(), EcalElectronicsMapping::getElectronicsId(), mps_fire::i, listChannels, listTowers, requestedFeds_, and edm::DataFrameContainer::size().

Referenced by analyze().

◆ readPNDigis()

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

protected

Definition at line 353 of file EcalDigiDisplay.cc.

                                                                                     {
   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";
     }
   }
 }