#include <EcnaAnalyzer.h>

Inheritance diagram for EcnaAnalyzer:

Public Types
enum	{ kChannels = 1700, kGains = 3, kFirstGainId = 1 }

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
Bool_t	AnalysisOutcome (const TString &)

void	analyze (const edm::Event &, const edm::EventSetup &) override

void	CheckMsg (const Int_t &, const Int_t &)

void	CheckMsg (const Int_t &)

	EcnaAnalyzer (const edm::ParameterSet &)

Int_t	gainvalue (const Int_t &)

TString	runtype (const Int_t &)

	~EcnaAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	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)

Private Attributes
std::string	digiProducer_

std::string	EBdigiCollection_

edm::EDGetTokenT< EBDigiCollection >	EBdigiToken_

std::string	EEdigiCollection_

edm::EDGetTokenT< EEDigiCollection >	EEdigiToken_

std::string	eventHeaderCollection_

std::string	eventHeaderProducer_

edm::EDGetTokenT< EcalRawDataCollection >	eventHeaderToken_

TString	fAnalysisName

Int_t	fANY_RUN

std::vector< Int_t >	fBuildEventDistribBad

std::vector< Int_t >	fBuildEventDistribGood

TString	fCfgAnalyzerParametersFileName

TString	fCfgAnalyzerParametersFilePath

Int_t	fChozenGainNumber

Int_t	fChozenRunTypeNumber

Int_t	fCurrentEventNumber

std::vector< TString >	fDateFirst

std::vector< TString >	fDateLast

Int_t	fDeeDS5Memo1

Int_t	fDeeDS5Memo2

Int_t	fDeeIndexBegin

Int_t	fDeeIndexStop

std::vector< std::string >	fDeeNumberString

TString	fDynBaseLineSub

std::vector< Int_t >	fESFromFedTcc

Int_t	fEvtNumber

Int_t	fEvtNumberMemo

std::ifstream	fFcin_f

std::vector< Int_t >	fFedDigiOK

Int_t	fFedId

std::array< Int_t, fMaxFedIdCounter >	fFedIdCounter

std::vector< Int_t >	fFedNbOfTreatedEvents

std::vector< Int_t >	fFedStatus

std::vector< Int_t >	fFedStatusOrder

Int_t	fFedTcc

Int_t	fFirstReqEvent

Int_t	fLASER_STD

Int_t	fLastReqEvent

Int_t	fMaxFedUnitCounter = 0

Int_t	fMaxTreatedStexCounter = 0

Int_t	fMemoCutOK

std::vector< Int_t >	fMemoDateFirstEvent

std::array< Int_t, fMaxMgpaGainCounter >	fMgpaGainCounter

Int_t	fMgpaGainNumber

std::vector< std::unique_ptr< TEcnaRun > >	fMyCnaEBSM

std::vector< std::unique_ptr< TEcnaRun > >	fMyCnaEEDee

TEcnaParEcal	fMyEBEcal

TEcnaNumbering	fMyEBNumbering

TEcnaObject	fMyEcnaEBObjectManager

TEcnaObject	fMyEcnaEEObjectManager

TEcnaParEcal	fMyEEEcal

TEcnaNumbering	fMyEENumbering

Int_t	fNbOfSamples

Int_t	fNbOfSelectedEvents

std::vector< Int_t >	fNbOfTreatedFedsInDee

std::vector< Int_t >	fNbOfTreatedFedsInStex

Int_t	fNbOfTreatedStexs

Bool_t	fOutcomeError

Int_t	fPEDESTAL_GAP

Int_t	fPEDESTAL_STD

Int_t	fPEDSIM

Int_t	fPHYSICS_GLOBAL

Int_t	fRecNumber

Int_t	fReqNbOfEvts

Int_t	fRunNumber

std::array< Int_t, fMaxRunTypeCounter >	fRunTypeCounter

Int_t	fRunTypeNumber

std::vector< Int_t >	fSMFromFedTcc

Int_t	fSMIndexBegin

Int_t	fSMIndexStop

std::vector< Int_t >	fStexDigiOK

Int_t	fStexIndexBegin

Int_t	fStexIndexStop

TString	fStexName

std::vector< Int_t >	fStexNbOfTreatedEvents

Int_t	fStexNumber

std::vector< Int_t >	fStexStatus

std::vector< time_t >	fTimeFirst

std::vector< time_t >	fTimeLast

Int_t	fTreatedFedOrder

TString	fTTBELL

Int_t	iEvent_

Int_t	nChannels_

TString	sAnalysisName_

TString	sFirstReqEvent_

TString	sLastReqEvent_

TString	sNbOfSamples_

TString	sReqNbOfEvts_

TString	sStexName_

TString	sStexNumber_

unsigned int	verbosity_

Static Private Attributes
static constexpr Int_t	fgMaxCar = 512

static constexpr Int_t	fMaxFedIdCounter = 54

static constexpr Int_t	fMaxMgpaGainCounter = 4

static constexpr Int_t	fMaxRunTypeCounter = 26

Additional Inherited Members
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

EcnaAnalyzer.h Update: 11/0'/2011 Authors: B.Fabbro (berna.nosp@m.rd.f.nosp@m.abbro.nosp@m.@cea.nosp@m..fr) DSM/IRFU/SPP CEA-Saclay Copyright: Those valid for CEA sofware

ECNA web page: http://cms-fabbro.web.cern.ch/cms-fabbro/

cna_new/Correlated_Noise_Analysis/ECNA_cna_1.htm

--------------------------------— Analysis name codes

  TString  AnalysisName: code for the analysis. According to this code,
                         the analyzer selects the event type
                         (PEDESTAL_STD, PEDESTAL_GAP, LASER_STD, etc...)
                         and some other event characteristics
                         (example: the gain in pedestal runs:
                          AnalysisName = "Ped1" or "Ped6" or "Ped12")
                         The string AnalysisName is automatically
                         included in the name of the results files

              AnalysisName  RunType         Gain    DBLS (Dynamic
              BaseLine Substraction)
              ..........................................

              AdcAny        any run type       0    No

              AdcPed1       fPEDESTAL_STD      3    No
              AdcPed6       fPEDESTAL_STD      2    No
              AdcPed12      fPEDESTAL_STD      1    No

              AdcPeg12      fPEDESTAL_GAP      1    No

              AdcLaser      fLASER_STD         0    No
              AdcPes12      fPEDSIM            0    No

              AdcPhys       fPHYSICS_GLOBAL    0    No


              AdcSPed1      fPEDESTAL_STD      3    Yes
              AdcSPed6      fPEDESTAL_STD      2    Yes
              AdcSPed12     fPEDESTAL_STD      1    Yes

              AdcSPeg12     fPEDESTAL_GAP      1    Yes

              AdcSLaser     fLASER_STD         0    Yes
              AdcSPes12     fPEDSIM            0    Yes

CalibCalorimetry/EcalCorrelatedNoiseAnalysisModules/src/EcnaAnalyzer.cc

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

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

Definition at line 133 of file EcnaAnalyzer.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

Enumerator
kChannels
kGains
kFirstGainId

Definition at line 135 of file EcnaAnalyzer.h.

135 { kChannels = 1700, kGains = 3, kFirstGainId = 1 };

EcnaAnalyzer::kChannels

Definition: EcnaAnalyzer.h:135

EcnaAnalyzer::kGains

Definition: EcnaAnalyzer.h:135

EcnaAnalyzer::kFirstGainId

Definition: EcnaAnalyzer.h:135

Constructor & Destructor Documentation

◆ EcnaAnalyzer()

EcnaAnalyzer::EcnaAnalyzer ( const edm::ParameterSet & pSet )

explicit

Definition at line 45 of file EcnaAnalyzer.cc.

References AnalysisOutcome(), digiProducer_, EBdigiCollection_, EBdigiToken_, EEdigiCollection_, EEdigiToken_, eventHeaderCollection_, eventHeaderProducer_, eventHeaderToken_, Exception, fAnalysisName, fANY_RUN, fBuildEventDistribBad, fBuildEventDistribGood, fChozenGainNumber, fChozenRunTypeNumber, fCurrentEventNumber, fDateFirst, fDateLast, fDeeDS5Memo1, fDeeDS5Memo2, fDeeIndexBegin, fDeeIndexStop, fDeeNumberString, fDynBaseLineSub, fESFromFedTcc, fEvtNumber, fEvtNumberMemo, fFedDigiOK, fFedId, fFedIdCounter, fFedNbOfTreatedEvents, fFedStatus, fFedStatusOrder, fFedTcc, fFirstReqEvent, fgMaxCar, ntuplemaker::fill, fLASER_STD, fLastReqEvent, fMaxFedUnitCounter, fMaxTreatedStexCounter, fMemoCutOK, fMemoDateFirstEvent, fMgpaGainCounter, fMgpaGainNumber, fMyEBEcal, fMyEcnaEBObjectManager, fMyEcnaEEObjectManager, fMyEEEcal, fNbOfSamples, fNbOfSelectedEvents, fNbOfTreatedFedsInDee, fNbOfTreatedFedsInStex, fNbOfTreatedStexs, fOutcomeError, fPEDESTAL_GAP, fPEDESTAL_STD, fPEDSIM, fPHYSICS_GLOBAL, fRecNumber, fReqNbOfEvts, fRunNumber, fRunTypeCounter, fRunTypeNumber, fSMFromFedTcc, fSMIndexBegin, fSMIndexStop, fStexDigiOK, fStexIndexBegin, fStexIndexStop, fStexName, fStexNbOfTreatedEvents, fStexNumber, fStexStatus, fTimeFirst, fTimeLast, fTreatedFedOrder, fTTBELL, edm::ParameterSet::getParameter(), mps_fire::i, ProducerED_cfi::InputTag, TEcnaParEcal::MaxDeeInEE(), TEcnaParEcal::MaxDSInEE(), TEcnaParEcal::MaxSMInEB(), sAnalysisName_, sFirstReqEvent_, sLastReqEvent_, sNbOfSamples_, sReqNbOfEvts_, sStexName_, sStexNumber_, AlCaHLTBitMon_QueryRunRegistry::string, and trackerHitRTTI::vector.

     : verbosity_(pSet.getUntrackedParameter("verbosity", 1U)),
       nChannels_(0),
       iEvent_(0),
       fMyEBNumbering(&fMyEcnaEBObjectManager, "EB"),
       fMyEBEcal(&fMyEcnaEBObjectManager, "EB"),
       fMyEENumbering(&fMyEcnaEEObjectManager, "EE"),
       fMyEEEcal(&fMyEcnaEEObjectManager, "EE") {
   // now do what ever initialization is needed
 
   std::unique_ptr<TEcnaParPaths> myPathEB = std::make_unique<TEcnaParPaths>(&fMyEcnaEBObjectManager);
   std::unique_ptr<TEcnaParPaths> myPathEE = std::make_unique<TEcnaParPaths>(&fMyEcnaEEObjectManager);
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-constructor> Check path for resultsq Root files.";
 
   if (myPathEB->GetPathForResultsRootFiles() == kFALSE) {
     edm::LogError("ecnaAnal") << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found.";
     throw cms::Exception("Calibration") << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found.";
   } else {
     edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-constructor> Path for result files found = "
                                  << myPathEB->ResultsRootFilePath();
   }
 
   if (myPathEE->GetPathForResultsRootFiles() == kFALSE) {
     edm::LogError("ecnaAnal") << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found.";
     throw cms::Exception("Calibration") << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found.";
   } else {
     edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-constructor> Path for result files found = "
                                  << myPathEE->ResultsRootFilePath();
   }
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-constructor> Parameter initialization.";
 
   fTTBELL = '\007';
   fOutcomeError = kFALSE;
 
   //==========================================================================================
   //.................................. Get parameter values from python file
   eventHeaderProducer_ = pSet.getParameter<std::string>("eventHeaderProducer");
   digiProducer_ = pSet.getParameter<std::string>("digiProducer");
 
   eventHeaderCollection_ = pSet.getParameter<std::string>("eventHeaderCollection");
   eventHeaderToken_ = consumes<EcalRawDataCollection>(edm::InputTag(digiProducer_, eventHeaderCollection_));
 
   EBdigiCollection_ = pSet.getParameter<std::string>("EBdigiCollection");
   EEdigiCollection_ = pSet.getParameter<std::string>("EEdigiCollection");
   EBdigiToken_ = consumes<EBDigiCollection>(edm::InputTag(digiProducer_, EBdigiCollection_));
   EEdigiToken_ = consumes<EEDigiCollection>(edm::InputTag(digiProducer_, EEdigiCollection_));
 
   sAnalysisName_ = pSet.getParameter<std::string>("sAnalysisName");
   sNbOfSamples_ = pSet.getParameter<std::string>("sNbOfSamples");
   sFirstReqEvent_ = pSet.getParameter<std::string>("sFirstReqEvent");
   sLastReqEvent_ = pSet.getParameter<std::string>("sLastReqEvent");
   sReqNbOfEvts_ = pSet.getParameter<std::string>("sReqNbOfEvts");
   sStexName_ = pSet.getParameter<std::string>("sStexName");
   sStexNumber_ = pSet.getParameter<std::string>("sStexNumber");
 
   fAnalysisName = sAnalysisName_.Data();
   fNbOfSamples = atoi(sNbOfSamples_.Data());
   fFirstReqEvent = atoi(sFirstReqEvent_.Data());
   fLastReqEvent = atoi(sLastReqEvent_.Data());
   fReqNbOfEvts = atoi(sReqNbOfEvts_.Data());
   fStexName = sStexName_.Data();
   fStexNumber = atoi(sStexNumber_.Data());
 
   //------------------------------- ERRORS in requested evts numbers
   if (fFirstReqEvent < 1) {
     fOutcomeError = AnalysisOutcome("ERR_FNEG");
   }
 
   if ((fLastReqEvent >= fFirstReqEvent) && (fReqNbOfEvts > fLastReqEvent - fFirstReqEvent + 1)) {
     fOutcomeError = AnalysisOutcome("ERR_LREQ");
   }
 
   if (fOutcomeError == kTRUE)
     return;
   //===========================================================================================
 
   std::fill(fRunTypeCounter.begin(), fRunTypeCounter.end(), 0);
   std::fill(fMgpaGainCounter.begin(), fMgpaGainCounter.end(), 0);
   std::fill(fFedIdCounter.begin(), fFedIdCounter.end(), 0);
 
   fEvtNumber = 0;
   fEvtNumberMemo = -1;
   fRecNumber = 0;
 
   fDeeDS5Memo1 = 0;
   fDeeDS5Memo2 = 0;
 
   fCurrentEventNumber = 0;
   fNbOfSelectedEvents = 0;
 
   fMemoCutOK = 0;
   fTreatedFedOrder = 0;
   fNbOfTreatedStexs = 0;
 
   //-------------- Fed
   if (fStexName == "SM") {
     fMaxFedUnitCounter = fMyEBEcal.MaxSMInEB();
   }  // EB: FED Unit = SM
   if (fStexName == "Dee") {
     fMaxFedUnitCounter = fMyEEEcal.MaxDSInEE();
   }  // EE: FED Unit = Data Sector
 
   fFedDigiOK = std::vector<Int_t>(fMaxFedUnitCounter, 0);
   fFedNbOfTreatedEvents = std::vector<Int_t>(fMaxFedUnitCounter, 0);
   fFedStatus = std::vector<Int_t>(fMaxFedUnitCounter, 0);
   fFedStatusOrder = std::vector<Int_t>(fMaxFedUnitCounter, 0);
 
   if (fStexName != "Dee") {
     fDeeNumberString = std::vector<std::string>(fMaxFedUnitCounter, "SM");
   } else {
     fDeeNumberString = {"Sector1 Dee4",
                         "Sector2 Dee4",
                         "Sector3 Dee4",
                         "Sector4 Dee4",
                         "Sector5 Dee4-Dee3",
                         "Sector6 Dee3",
                         "Sector7 Dee3",
                         "Sector8 Dee3",
                         "Sector9 Dee3",
                         "Sector1 Dee1",
                         "Sector2 Dee1",
                         "Sector3 Dee1",
                         "Sector4 Dee1",
                         "Sector5 Dee1-Dee2",
                         "Sector6 Dee2",
                         "Sector7 Dee2",
                         "Sector8 Dee2",
                         "Sector9 Dee2"};
   }
   //............................... arrays fSMFromFedDcc and fESFromFedTcc
   //
   //  FED-TCC:   1   2   3   4   5   6   7   8   9
   //      Dee:   3   3   3   4   4   4   4  4-3  3
   //       DS:   7   8   9   1   2   3   4   5   6
   //       ES:   7   8   9   1   2   3   4   5   6  (ES = DS)
   //
   //  FED-TCC:  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25
   //  26  27
   //       SM:  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34
   //       35  36 SM:  -1  -2  -3  -4  -5  -6  -7  -8  -9 -10 -11 -12 -13 -14
   //       -15 -16 -17 -18
   //
   //  FED-TCC:  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43
   //  44  45
   //       SM:   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16
   //       17  18
   //
   //  FED-TCC:  46  47  48  49  50  51  52  53  54
   //      Dee:   2   2   2   1   1   1   1  1-2  2
   //       DS:   7   8   9   1   2   3   4   5   6
   //       ES:  16  17  18  10  11  12  13  14  15  (ES = DS + 9)
 
   const Int_t MaxSMAndDS = fMyEBEcal.MaxSMInEB() + fMyEEEcal.MaxDSInEE();
 
   fSMFromFedTcc = std::vector<Int_t>(MaxSMAndDS, -1);
   fESFromFedTcc = std::vector<Int_t>(MaxSMAndDS, -1);
 
   for (Int_t nFedTcc = 1; nFedTcc <= 3; nFedTcc++) {
     fESFromFedTcc[nFedTcc - 1] = nFedTcc + 6;
   }  // Dee3, ES 7,8,9
   for (Int_t nFedTcc = 4; nFedTcc <= 9; nFedTcc++) {
     fESFromFedTcc[nFedTcc - 1] = nFedTcc - 3;
   }  // Dee4, ES 1,2,3,4,5; Dee3, DS 5,6
 
   for (Int_t nFedTcc = 10; nFedTcc <= 27; nFedTcc++) {
     fSMFromFedTcc[nFedTcc - 1] = nFedTcc + 9;
   }  // EB-  SM 19 to 36
   for (Int_t nFedTcc = 28; nFedTcc <= 45; nFedTcc++) {
     fSMFromFedTcc[nFedTcc - 1] = nFedTcc - 27;
   }  // EB+  SM  1 to 18
 
   for (Int_t nFedTcc = 46; nFedTcc <= 48; nFedTcc++) {
     fESFromFedTcc[nFedTcc - 1] = nFedTcc - 30;
   }  // Dee2, ES 16,17,18
   for (Int_t nFedTcc = 49; nFedTcc <= 54; nFedTcc++) {
     fESFromFedTcc[nFedTcc - 1] = nFedTcc - 39;
   }  // Dee1, ES 10,11,12,13,14; Dee2, ES 14,15
 
   //............................... Nb of treated events for "AdcPeg12" and
   //"AdcSPeg12" analysis
   //-------------- Stex
   if (fStexName == "SM") {
     fMaxTreatedStexCounter = fMyEBEcal.MaxSMInEB();
   }  // EB: Stex = SM
   if (fStexName == "Dee") {
     fMaxTreatedStexCounter = fMyEEEcal.MaxDeeInEE();
   }  // EE: Stex = Dee
 
   fStexNbOfTreatedEvents = std::vector<Int_t>(fMaxTreatedStexCounter, 0);
 
   fTimeFirst = std::vector<time_t>(fMaxTreatedStexCounter, 0);
   fTimeLast = std::vector<time_t>(fMaxTreatedStexCounter, 0);
 
   fMemoDateFirstEvent = std::vector<Int_t>(fMaxTreatedStexCounter, 0);
 
   Int_t MaxCar = fgMaxCar;
   fDateFirst = std::vector<TString>(fMaxTreatedStexCounter);
   for (Int_t i = 0; i < fMaxTreatedStexCounter; i++) {
     fDateFirst[i].Resize(MaxCar);
     fDateFirst[i] = "*1st event date not found*";
   }
 
   MaxCar = fgMaxCar;
   fDateLast = std::vector<TString>(fMaxTreatedStexCounter);
   for (Int_t i = 0; i < fMaxTreatedStexCounter; i++) {
     fDateLast[i].Resize(MaxCar);
     fDateLast[i] = "*last event date not found*";
   }
 
   fStexStatus = std::vector(fMaxTreatedStexCounter, 0);
   fStexDigiOK = std::vector(fMaxTreatedStexCounter, 0);
 
   fNbOfTreatedFedsInDee = std::vector(fMaxTreatedStexCounter, 0);
   fNbOfTreatedFedsInStex = std::vector(fMaxTreatedStexCounter, 0);
 
   //.......................... counters of events for GetSampleAdcValues
   fBuildEventDistribBad = std::vector(fMaxTreatedStexCounter, 0);
 
   fBuildEventDistribGood = std::vector(fMaxTreatedStexCounter, 0);
 
   //----------------------------------- Analysis name codes
   //------------------------------------------
   //
   //                  AnalysisName  RunType         Gain    DBLS (Dynamic
   //                  BaseLine Substraction)
   //
   //                  AdcAny        any run type       0    no
   //
   //                  AdcPed1       fPEDESTAL_STD      3    No
   //                  AdcPed6       fPEDESTAL_STD      2    No
   //                  AdcPed12      fPEDESTAL_STD      1    No
   //
   //                  AdcPeg12      fPEDESTAL_GAP      1    No
   //
   //                  AdcLaser      fLASER_STD         0    No
   //                  AdcPes12      fPEDSIM            0    No
   //
   //                  AdcPhys       fPHYSICS_GLOBAL    0    No
   //
   //
   //                  AdcSPed1      fPEDESTAL_STD      3    Yes
   //                  AdcSPed6      fPEDESTAL_STD      2    Yes
   //                  AdcSPed12     fPEDESTAL_STD      1    Yes
   //
   //                  AdcSPeg12     fPEDESTAL_GAP      1    Yes
   //
   //                  AdcSLaser     fLASER_STD         0    Yes
   //                  AdcSPes12     fPEDSIM            0    Yes
   //
   //--------------------------------------------------------------------------------------------------
 
   //................ Run type list
 
   fLASER_STD = 4;
   fPEDESTAL_STD = 9;
   fPHYSICS_GLOBAL = 13;
   fPEDESTAL_GAP = 18;
   fPEDSIM = 24;
 
   fANY_RUN = 25;
 
   //................ Chozen run type from analysis name
   fChozenRunTypeNumber = fANY_RUN;  // default
   if (fAnalysisName == "AdcAny") {
     fChozenRunTypeNumber = fANY_RUN;
   }
   if (fAnalysisName == "AdcPed1" || fAnalysisName == "AdcPed6" || fAnalysisName == "AdcPed12" ||
       fAnalysisName == "AdcSPed1" || fAnalysisName == "AdcSPed6" || fAnalysisName == "AdcSPed12") {
     fChozenRunTypeNumber = fPEDESTAL_STD;
   }
   if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12") {
     fChozenRunTypeNumber = fPEDESTAL_GAP;
   }
   if (fAnalysisName == "AdcLaser" || fAnalysisName == "AdcSLaser") {
     fChozenRunTypeNumber = fLASER_STD;
   }
   if (fAnalysisName == "AdcPhys") {
     fChozenRunTypeNumber = fPHYSICS_GLOBAL;
   }
   if (fAnalysisName == "AdcPes12 " || fAnalysisName == "AdcSPes12 ") {
     fChozenRunTypeNumber = fPEDSIM;
   }
 
   //................ Gains from analysis name
   fChozenGainNumber = 0;  // default => event always accepted if fChozenGainNumber = 0 ( see
                           // USER's Analysis cut in ::analyze(...) )
   if (fAnalysisName == "AdcAny") {
     fChozenGainNumber = 0;
   }
   if (fAnalysisName == "AdcPed1" || fAnalysisName == "AdcSPed1") {
     fChozenGainNumber = 3;
   }
   if (fAnalysisName == "AdcPed6" || fAnalysisName == "AdcSPed6") {
     fChozenGainNumber = 2;
   }
   if (fAnalysisName == "AdcPed12" || fAnalysisName == "AdcSPed12") {
     fChozenGainNumber = 1;
   }
   if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12") {
     fChozenGainNumber = 0;
   }
   if (fAnalysisName == "AdcLaser" || fAnalysisName == "AdcSLaser") {
     fChozenGainNumber = 0;
   }
   if (fAnalysisName == "AdcPes12 " || fAnalysisName == "AdcSPes12 ") {
     fChozenGainNumber = 0;
   }
   if (fAnalysisName == "AdcPhys") {
     fChozenGainNumber = 0;
   }
 
   //............... Flag for Dynamic BaseLine Substraction from analysis name
   fDynBaseLineSub = "no";  // default
   if (fAnalysisName == "AdcAny" || fAnalysisName == "AdcPed1" || fAnalysisName == "AdcPed6" ||
       fAnalysisName == "AdcPed12" || fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcLaser" ||
       fAnalysisName == "AdcPhys" || fAnalysisName == "AdcPes12 ") {
     fDynBaseLineSub = "no";
   }
   if (fAnalysisName == "AdcSPed1" || fAnalysisName == "AdcSPed6" || fAnalysisName == "AdcSPed12" ||
       fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcSLaser" || fAnalysisName == "AdcSPes12 ") {
     fDynBaseLineSub = "yes";
   }
 
   //....................... Index range for ECNA init and for loop on
   // GetSampleAdcValues calls
   if (fStexNumber == 0) {
     if (fStexName == "SM") {
       fSMIndexBegin = 0;
       fSMIndexStop = fMyEBEcal.MaxSMInEB();
       fStexIndexBegin = fSMIndexBegin;
       fStexIndexStop = fSMIndexStop;
       fDeeIndexBegin = 0;
       fDeeIndexStop = 0;
     }
     if (fStexName == "Dee") {
       fSMIndexBegin = 0;
       fSMIndexStop = 0;
       fDeeIndexBegin = 0;
       fDeeIndexStop = fMyEEEcal.MaxDeeInEE();
       fStexIndexBegin = fDeeIndexBegin;
       fStexIndexStop = fDeeIndexStop;
     }
   } else {
     if (fStexName == "SM") {
       fSMIndexBegin = fStexNumber - 1;
       fSMIndexStop = fStexNumber;
       fStexIndexBegin = fSMIndexBegin;
       fStexIndexStop = fSMIndexStop;
       fDeeIndexBegin = 0;
       fDeeIndexStop = 0;
     }
     if (fStexName == "Dee") {
       fSMIndexBegin = 0;
       fSMIndexStop = 0;
       fDeeIndexBegin = fStexNumber - 1;
       fDeeIndexStop = fStexNumber;
       fStexIndexBegin = fDeeIndexBegin;
       fStexIndexStop = fDeeIndexStop;
     }
   }
 
   //......... DATA DEPENDENT PARAMETERS
   fRunNumber = 0;
 
   fRunTypeNumber = -1;
   fMgpaGainNumber = -1;
 
   fFedId = -1;
   fFedTcc = -1;
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fAnalysisName        = " << fAnalysisName;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fNbOfSamples         = " << fNbOfSamples;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fFirstReqEvent       = " << fFirstReqEvent;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fLastReqEvent        = " << fLastReqEvent;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fReqNbOfEvts         = " << fReqNbOfEvts;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fStexName            = " << fStexName;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fStexNumber          = " << fStexNumber;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fChozenRunTypeNumber = "
                                << fChozenRunTypeNumber;
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fChozenGainNumber    = "
                                << fChozenGainNumber << std::endl;
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer::EcnaAnalyzer-constructor> Init done. ";
 }

◆ ~EcnaAnalyzer()

EcnaAnalyzer::~EcnaAnalyzer ( )

override

Definition at line 433 of file EcnaAnalyzer.cc.

References CheckMsg(), fBuildEventDistribBad, fBuildEventDistribGood, fDateFirst, fDateLast, fDeeIndexBegin, fDeeIndexStop, fFedId, fFedIdCounter, fMaxFedIdCounter, fMaxMgpaGainCounter, fMaxRunTypeCounter, fMgpaGainCounter, fMgpaGainNumber, fMyCnaEBSM, fMyCnaEEDee, fOutcomeError, fRunTypeCounter, fRunTypeNumber, fSMIndexBegin, fSMIndexStop, fStexIndexBegin, fStexIndexStop, fStexName, fStexNumber, fTimeFirst, fTimeLast, fTTBELL, gainvalue(), mps_fire::i, n0, and runtype().

                             {
   // do anything here that needs to be done at destruction time
   // (e.g. close files, deallocate resources etc.)
 
   //..................................... format numerical values
   edm::LogVerbatim("ecnaAnal") << "EcnaAnalyzer::~EcnaAnalyzer()> destructor is going to be executed." << std::endl;
 
   if (fOutcomeError == kTRUE)
     return;
 
   //-------------------------------------------------------------------------------
 
   //....................................................... EB (SM)
   if (fMyCnaEBSM.empty() && fStexName == "SM") {
     edm::LogVerbatim("ecnaAnal") << "\n!EcnaAnalyzer-destructor> **** ERROR **** fMyCnaEBSM is empty"
                                  << ". !===> ECNA HAS NOT BEEN INITIALIZED."
                                  << "\n  Last event run type = " << runtype(fRunTypeNumber)
                                  << ", fRunTypeNumber = " << fRunTypeNumber
                                  << ", last event Mgpa gain = " << gainvalue(fMgpaGainNumber)
                                  << ", fMgpaGainNumber = " << fMgpaGainNumber
                                  << ", last event fFedId(+601) = " << fFedId + 601 << std::endl;
   } else {
     for (Int_t iSM = fSMIndexBegin; iSM < fSMIndexStop; iSM++) {
       if (fMyCnaEBSM[iSM].get() != nullptr) {
         //........................................ register dates 1 and 2
         fMyCnaEBSM[iSM]->StartStopDate(fDateFirst[iSM], fDateLast[iSM]);
         fMyCnaEBSM[iSM]->StartStopTime(fTimeFirst[iSM], fTimeLast[iSM]);
 
         //........................................ Init .root file
         fMyCnaEBSM[iSM]->GetReadyToCompute();
         fMyCnaEBSM[iSM]->SampleValues();
 
         //........................................ write the sample values in
         //.root file
         if (fMyCnaEBSM[iSM]->WriteRootFile() == kFALSE) {
           edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer-destructor> PROBLEM with write ROOT file for SM" << iSM + 1
                                        << fTTBELL;
         }
       } else {
         edm::LogVerbatim("ecnaAnal")
             << "*EcnaAnalyzer-destructor> Calculations and writing on file already done for SM " << iSM + 1;
       }
     }
   }
   //....................................................... EE (Dee)
 
   if (fMyCnaEEDee.empty() && fStexName == "Dee") {
     edm::LogVerbatim("ecnaAnal") << "\n!EcnaAnalyzer-destructor> **** ERROR **** fMyCnaEEDee is empty"
                                  << ". !===> ECNA HAS NOT BEEN INITIALIZED."
                                  << "\n  Last event run type = " << runtype(fRunTypeNumber)
                                  << ", fRunTypeNumber = " << fRunTypeNumber
                                  << ", last event Mgpa gain = " << gainvalue(fMgpaGainNumber)
                                  << ", fMgpaGainNumber = " << fMgpaGainNumber
                                  << ", last event fFedId(+601) = " << fFedId + 601 << std::endl;
   } else {
     for (Int_t iDee = fDeeIndexBegin; iDee < fDeeIndexStop; iDee++) {
       if (fMyCnaEEDee[iDee].get() != nullptr) {
         //........................................ register dates 1 and 2
         fMyCnaEEDee[iDee]->StartStopDate(fDateFirst[iDee], fDateLast[iDee]);
         fMyCnaEEDee[iDee]->StartStopTime(fTimeFirst[iDee], fTimeLast[iDee]);
 
         //........................................ Init .root file
         fMyCnaEEDee[iDee]->GetReadyToCompute();
         fMyCnaEEDee[iDee]->SampleValues();
 
         //........................................ write the sample values in
         //.root file
         if (fMyCnaEEDee[iDee]->WriteRootFile() == kFALSE) {
           edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer-destructor> PROBLEM with write ROOT file for Dee" << iDee + 1
                                        << " " << fTTBELL;
         }
       } else {
         edm::LogVerbatim("ecnaAnal")
             << "*EcnaAnalyzer-destructor> Calculations and writing on file already done for Dee " << iDee + 1;
       }
     }
   }
   edm::LogVerbatim("ecnaAnal") << std::endl;
 
   //-----------------------------------------------------------------------------------
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> Status of events returned by GetSampleAdcValues(): ";
 
   for (Int_t i0Stex = fStexIndexBegin; i0Stex < fStexIndexStop; i0Stex++) {
     edm::LogVerbatim("ecnaAnal") << fStexName << i0Stex + 1 << "> Status OK: " << fBuildEventDistribGood[i0Stex]
                                  << " / ERROR(S): " << fBuildEventDistribBad[i0Stex];
     if (fBuildEventDistribBad[i0Stex] > 0) {
       edm::LogVerbatim("ecnaAnal") << " <=== SHOULD BE EQUAL TO ZERO ! " << fTTBELL;
     }
     edm::LogVerbatim("ecnaAnal") << std::endl;
   }
 
   edm::LogVerbatim("ecnaAnal") << "\n- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ";
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> Run types seen in event headers before selection:";
 
   for (Int_t i = 0; i < fMaxRunTypeCounter; i++) {
     edm::LogVerbatim("ecnaAnal") << " => " << std::setw(10) << fRunTypeCounter[i] << " event header(s) with run type "
                                  << runtype(i);
   }
 
   edm::LogVerbatim("ecnaAnal") << "\n- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ";
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> Mgpa gains seen in event headers before selection:";
 
   for (Int_t i = 0; i < fMaxMgpaGainCounter; i++) {
     edm::LogVerbatim("ecnaAnal") << " => " << std::setw(10) << fMgpaGainCounter[i] << " event header(s) with gain "
                                  << gainvalue(i);
   }
 
   edm::LogVerbatim("ecnaAnal") << "\n- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ";
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> Numbers of selected events for each FED:";
 
   for (Int_t i = 0; i < fMaxFedIdCounter; i++) {
     edm::LogVerbatim("ecnaAnal") << " => FedId " << i + 601 << ": " << std::setw(10) << fFedIdCounter[i] << " events";
   }
 
   edm::LogVerbatim("ecnaAnal") << "\n- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ";
 
   if (fStexNumber > 0) {
     edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> fDateFirst = " << fDateFirst[fStexNumber - 1]
                                  << "\n                          fDateLast  = " << fDateLast[fStexNumber - 1]
                                  << std::endl;
   }
 
   edm::LogVerbatim("ecnaAnal") << "\n- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ";
 
   Int_t n0 = 0;
   CheckMsg(n0);
 
   edm::LogVerbatim("ecnaAnal") << "*EcnaAnalyzer-destructor> End of execution.";
 }

Member Function Documentation

◆ AnalysisOutcome()

Bool_t EcnaAnalyzer::AnalysisOutcome ( const TString & s_opt )

Definition at line 1332 of file EcnaAnalyzer.cc.

References CheckMsg(), fCurrentEventNumber, fFirstReqEvent, fLastReqEvent, fMyEBEcal, fMyEEEcal, fNbOfSelectedEvents, fNbOfTreatedStexs, fReqNbOfEvts, fStexName, fStexNumber, TEcnaParEcal::MaxDeeInEE(), TEcnaParEcal::MaxSMInEB(), n0, and mps_fire::result.

Referenced by analyze(), and EcnaAnalyzer().

                                                          {
   //---- STOP if end of analysis
 
   Bool_t result = kFALSE;
 
   if (s_opt == "EVT") {
     Int_t MaxNbOfStex = 0;
     if (fStexName == "SM") {
       MaxNbOfStex = fMyEBEcal.MaxSMInEB();
     }
     if (fStexName == "Dee") {
       MaxNbOfStex = fMyEEEcal.MaxDeeInEE();
     }
 
     if (((fStexNumber > 0 && fNbOfTreatedStexs == 1) || (fStexNumber == 0 && fNbOfTreatedStexs == MaxNbOfStex)) &&
         ((fLastReqEvent < fFirstReqEvent) ||
          (fLastReqEvent >= fFirstReqEvent && fCurrentEventNumber <= fLastReqEvent))) {
       edm::LogVerbatim("ecnaAnal")
           << "\n**************************** ANALYSIS REPORT > OK **************************************"
           << "\n*EcnaAnalyzer::AnalysisOutcome(...)> The maximum requested number of events and the maximum"
           << "\n                                     number of treated " << fStexName << "'s have been reached."
           << "\n                                     Analysis successfully ended from EcnaAnalyzer "
           << "\n                                     Number of selected events   = " << fNbOfSelectedEvents
           << "\n                                     Last requested event number = " << fLastReqEvent
           << "\n                                     Current event number        = " << fCurrentEventNumber;
 
       Int_t n0 = 0;
       CheckMsg(n0);
 
       edm::LogVerbatim("ecnaAnal")
           << "****************************************************************************************\n";
 
       result = kTRUE;
       return result;
     }
 
     if (fLastReqEvent >= fFirstReqEvent && fCurrentEventNumber > fLastReqEvent &&
         !((fStexNumber > 0 && fNbOfTreatedStexs == 1) || (fStexNumber == 0 && fNbOfTreatedStexs == MaxNbOfStex))) {
       edm::LogVerbatim("ecnaAnal") << "\n**************************** ANALYSIS REPORT >>> *** "
                                       "WARNING *** WARNING *** WARNING ***"
                                    << "\n*EcnaAnalyzer::AnalysisOutcome(...)> Last event reached "
                                       "before completion of analysis."
                                    << "\n                                     Analysis ended from EcnaAnalyzer "
                                    << "\n                                     Number of selected events   = "
                                    << fNbOfSelectedEvents
                                    << "\n                                     Last requested event number = "
                                    << fLastReqEvent
                                    << "\n                                     Current event number        = "
                                    << fCurrentEventNumber;
 
       Int_t n0 = 0;
       CheckMsg(n0);
 
       edm::LogVerbatim("ecnaAnal")
           << "****************************************************************************************" << std::endl;
 
       result = kTRUE;
       return result;
     }
   } else {
     if (s_opt == "ERR_FNEG") {
       edm::LogVerbatim("ecnaAnal")
           << "\n**************************** ANALYSIS REPORT >>> **** ERROR **** ERROR **** ERROR ******"
           << "\n*EcnaAnalyzer::AnalysisOutcome(...)> First event number = " << fFirstReqEvent
           << ". Should be strictly potitive."
           << "\n                             Analysis ended from EcnaAnalyzer ";
 
       edm::LogVerbatim("ecnaAnal")
           << "****************************************************************************************" << std::endl;
 
       result = kTRUE;
       return result;
     }
     if (s_opt == "ERR_LREQ") {
       edm::LogVerbatim("ecnaAnal")
           << "\n**************************** ANALYSIS REPORT >>> **** ERROR **** ERROR **** ERROR ******"
           << "\n*EcnaAnalyzer::analyze(...)> Requested number of events = " << fReqNbOfEvts << "."
           << "\n                             Too large compared to the event range: " << fFirstReqEvent << " - "
           << fLastReqEvent << "\n                             Analysis ended from EcnaAnalyzer ";
 
       edm::LogVerbatim("ecnaAnal")
           << "****************************************************************************************" << std::endl;
 
       result = kTRUE;
       return result;
     }
   }
   return result;
 }  // end of EcnaAnalyzer::AnalysisOutcome(const Int_t& n_option)

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 573 of file EcnaAnalyzer.cc.

                                                                               {
   //..................................... format numerical values
   std::cout << std::setiosflags(std::ios::showpoint | std::ios::uppercase);
   std::cout << std::setprecision(3) << std::setw(6);
   std::cout.setf(std::ios::dec, std::ios::basefield);
   std::cout.setf(std::ios::fixed, std::ios::floatfield);
   std::cout.setf(std::ios::left, std::ios::adjustfield);
   std::cout.setf(std::ios::right, std::ios::adjustfield);
 
   fRecNumber++;
 
   Int_t iFreq = (fLastReqEvent - fFirstReqEvent + 1) / 5;
   if (iFreq <= 0) {
     iFreq = 10000;
   }
 
   Int_t MaxSMAndDS = fMyEBEcal.MaxSMInEB() + fMyEEEcal.MaxDSInEE();
 
   //********************************************* EVENT TREATMENT
   //********************************
   const edm::Handle<EcalRawDataCollection> &pEventHeader = iEvent.getHandle(eventHeaderToken_);
   const EcalRawDataCollection *myEventHeader = nullptr;
   if (pEventHeader.isValid()) {
     myEventHeader = pEventHeader.product();
   } else {
     edm::LogError("ecnaAnal") << "Error! can't get the product " << eventHeaderCollection_.c_str();
   }
   //........... Decode myEventHeader infos
   for (EcalRawDataCollection::const_iterator headerItr = myEventHeader->begin(); headerItr != myEventHeader->end();
        ++headerItr) {
     //===> fRunNumber, fRunTypeNumber, fMgpaGainNumber, fFedId, fEvtNumber
     //     will be used in AnalysisOutcome(...) below
     fRunNumber = (Int_t)headerItr->getRunNumber();
     if (fRunNumber <= 0) {
       fRunNumber = (Int_t)iEvent.id().run();
     }
     fRunTypeNumber = (Int_t)headerItr->getRunType();
     fMgpaGainNumber = (Int_t)headerItr->getMgpaGain();
     fFedId = (Int_t)headerItr->fedId() - 601;  // 1st Fed = 601, FedId = Fed number - 1
     fEvtNumber = (Int_t)headerItr->getLV1();
     if (fEvtNumber <= 0) {
       fEvtNumber = (Int_t)iEvent.id().event();
     }
 
     if (fEvtNumber != fEvtNumberMemo) {
       fEvtNumberMemo = fEvtNumber;
 
       //============================================
       //  cmsRun INTERRUPTION if analysis complete
       //  or if fCurrentEventNumber >= LastReqEvent
       //============================================
       if (AnalysisOutcome("EVT") == kTRUE) {
         return;
       }
 
       // no interruption => event has to be analyzed
 
       fCurrentEventNumber++;
 
       if (fRecNumber == 1 || fRecNumber == 50 || fRecNumber == 100 || fRecNumber == 500 || fRecNumber == 1000 ||
           fRecNumber % iFreq == 0) {
         Int_t n1 = 1;
         CheckMsg(n1);
       }
 
       if (fCurrentEventNumber < fFirstReqEvent)
         return;  // skip events before fFirstReqEvent
     }
 
     //.................. Increment Run type and MgpaGain counters
     if (fRunTypeNumber >= 0 && fRunTypeNumber < fMaxRunTypeCounter) {
       fRunTypeCounter[fRunTypeNumber]++;
     }
     if (fMgpaGainNumber >= 0 && fMgpaGainNumber < fMaxMgpaGainCounter) {
       fMgpaGainCounter[fMgpaGainNumber]++;
     }
 
     //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% User's analysis cut
     //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
     if (!(fRunNumber > 0 && (fRunTypeNumber == fChozenRunTypeNumber || fChozenRunTypeNumber == fANY_RUN) &&
           (fMgpaGainNumber == fChozenGainNumber || fChozenGainNumber == 0)))
       return;
 
     //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
     if (fMemoCutOK == 0) {
       fMemoCutOK = 1;
     }
 
     //---- Accelerating selection with "FED-TCC" number [ from
     // headerItr->getDccInTCCCommand() ]
     //     Arrays fSMFromFedTcc[] and fESFromFedTcc[] are initialised in Init()
 
     if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
         fAnalysisName == "AdcAny") {
       fFedTcc = (Int_t)headerItr->getDccInTCCCommand();
 
       if (fFedTcc >= 1 && fFedTcc <= MaxSMAndDS) {
         if (fStexName == "SM") {
           if (fFedTcc < 10 || fFedTcc > 45)
             return;
 
           if (fSMFromFedTcc[fFedTcc - 1] >= 1 && fSMFromFedTcc[fFedTcc - 1] <= fMyEBEcal.MaxSMInEB() &&
               fStexNbOfTreatedEvents[fSMFromFedTcc[fFedTcc - 1] - 1] >= fReqNbOfEvts)
             return;
         }
 
         if (fStexName == "Dee") {
           if (fFedTcc >= 10 && fFedTcc <= 45)
             return;
 
           if (fESFromFedTcc[fFedTcc - 1] >= 1 && fESFromFedTcc[fFedTcc - 1] <= fMyEEEcal.MaxDSInEE() &&
               fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1] >= fReqNbOfEvts)
             return;
         }
       }  // end of if( fFedTcc >= 1 && fFedTcc <= MaxSMAndDS )
     }    // end of if( fAnalysisName == "AdcPeg12"  || fAnalysisName == "AdcSPeg12"
     // ...)
 
     //.................. Increment FedId counters
     if (fFedId >= 0 && fFedId < fMaxFedIdCounter) {
       fFedIdCounter[fFedId]++;
     }
 
   }  // end of for(EcalRawDataCollection::const_iterator
   // headerItr=myEventHeader->begin(); headerItr !=
   // myEventHeader->end();++headerItr)
 
   if (fMemoCutOK == 0)
     return;  // return if no event passed the user's analysis cut
 
   //========================== SELECTED EVENTS ================================
   fNbOfSelectedEvents++;
   if (fNbOfSelectedEvents == 1) {
     Int_t n2 = 2;
     CheckMsg(n2);
   }
 
   //============================ Ecna init for the pointers array
   //=================================
   //.................................................................. EB (SM)
   if (fMyCnaEBSM.empty() && fStexName == "SM") {
     fMyCnaEBSM.resize(fMyEBEcal.MaxSMInEB());
   }
   //.................................................................. EE (Dee)
   if (fMyCnaEEDee.empty() && fStexName == "Dee") {
     fMyCnaEEDee.resize(fMyEEEcal.MaxDeeInEE());
   }
 
   //============================ EVENT TREATMENT ==============================
   Int_t MaxNbOfStex = 0;
   if (fStexName == "SM") {
     MaxNbOfStex = fMyEBEcal.MaxSMInEB();
   }
   if (fStexName == "Dee") {
     MaxNbOfStex = fMyEEEcal.MaxDeeInEE();
   }
 
   if ((fStexNumber > 0 && fNbOfTreatedStexs == 0) || (fStexNumber == 0 && fNbOfTreatedStexs < MaxNbOfStex)) {
     //================================================================= Record
     // type EB (SM)
     if (fStexName == "SM" && fSMIndexBegin < fSMIndexStop) {
       //......................................... Get digisEB
       const edm::Handle<EBDigiCollection> &pdigisEB = iEvent.getHandle(EBdigiToken_);
       const EBDigiCollection *digisEB = nullptr;
       if (pdigisEB.isValid()) {
         digisEB = pdigisEB.product();
       } else {
         edm::LogError("ecnaAnal") << "Error! can't get the product " << EBdigiCollection_.c_str();
       }
 
       // Initialize vectors if not already done
       if (int(digisEB->size()) > nChannels_) {
         nChannels_ = digisEB->size();
       }
 
       // Int_t print_count = 0;
       if (Int_t(digisEB->end() - digisEB->begin()) >= 0 &&
           Int_t(digisEB->end() - digisEB->begin()) <= Int_t(digisEB->size())) {
         //..........................................EB
         //===============================================================================
         //
         //                    Loop over Ecal barrel digisEB (Xtals)
         //
         //===============================================================================
 
         for (EBDigiCollection::const_iterator digiItr = digisEB->begin(); digiItr != digisEB->end(); ++digiItr) {
           EBDetId id_crystal(digiItr->id());
           // Int_t HashedIndex = id_crystal.hashedIndex();
 
           Int_t i0SM = id_crystal.ism() - 1;  //   <============== GET the SM number - 1 here
 
           if (i0SM >= 0 && i0SM < fMaxTreatedStexCounter) {
             if (fMyCnaEBSM[i0SM].get() == nullptr && fStexStatus[i0SM] != 2) {
               //=============================== Init Ecna EB
               //===============================
               fMyCnaEBSM[i0SM] = std::make_unique<TEcnaRun>(&fMyEcnaEBObjectManager, "EB", fNbOfSamples);
               fMyCnaEBSM[i0SM]->GetReadyToReadData(
                   fAnalysisName, fRunNumber, fFirstReqEvent, fLastReqEvent, fReqNbOfEvts, i0SM + 1, fRunTypeNumber);
 
               edm::LogVerbatim("ecnaAnal")
                   << "*EcnaAnalyzer::analyze(...)> ********* INIT ECNA EB ********* "
                   << "\n                                   fAnalysisName = " << fAnalysisName
                   << "\n                                      fRunNumber = " << fRunNumber
                   << "\n                                  fFirstReqEvent = " << fFirstReqEvent
                   << "\n                                   fLastReqEvent = " << fLastReqEvent
                   << "\n                                    fReqNbOfEvts = " << fReqNbOfEvts
                   << "\n                                              SM = " << i0SM + 1
                   << "\n                                        run type = " << runtype(fRunTypeNumber);
               //============================================================================
             }
 
             if (fStexStatus[i0SM] < 2)  // nothing to do if status=2 reached
             {
               fStexDigiOK[i0SM]++;
               if (fStexDigiOK[i0SM] == 1) {
                 fStexNbOfTreatedEvents[i0SM]++;
               }
 
               if (fStexNbOfTreatedEvents[i0SM] >= 1 && fStexNbOfTreatedEvents[i0SM] <= fReqNbOfEvts) {
                 //......................................... date of first event
                 //(in real time)
                 edm::Timestamp Time = iEvent.time();
                 edm::TimeValue_t t_current_ev_time = (cond::Time_t)Time.value();
                 time_t i_current_ev_time = (time_t)(t_current_ev_time >> 32);
                 const time_t *p_current_ev_time = &i_current_ev_time;
                 char *astime = ctime(p_current_ev_time);
 
                 if (fStexDigiOK[i0SM] == 1 && fStexNbOfTreatedEvents[i0SM] == 1 &&
                     (fStexNumber == 0 || i0SM + 1 == fStexNumber)) {
                   fTimeFirst[i0SM] = i_current_ev_time;
                   fDateFirst[i0SM] = astime;
                   fTimeLast[i0SM] = i_current_ev_time;
                   fDateLast[i0SM] = astime;
                   edm::LogVerbatim("ecnaAnal") << "*----> beginning of analysis for " << fStexName << i0SM + 1
                                                << ". First analyzed event date : " << astime;
                 }
 
                 if (i_current_ev_time < fTimeFirst[i0SM]) {
                   fTimeFirst[i0SM] = i_current_ev_time;
                   fDateFirst[i0SM] = astime;
                 }
                 if (i_current_ev_time > fTimeLast[i0SM]) {
                   fTimeLast[i0SM] = i_current_ev_time;
                   fDateLast[i0SM] = astime;
                 }
 
                 //=============================================> CUT on i0SM
                 // value
                 if ((fStexNumber > 0 && i0SM == fStexNumber - 1) || (fStexNumber == 0)) {
                   Int_t iEta = id_crystal.ietaSM();  // ietaSM() : range = [1,85]
                   Int_t iPhi = id_crystal.iphiSM();  // iphiSM() : range = [1,20]
 
                   Int_t n1SMCrys = (iEta - 1) * (fMyEBEcal.MaxTowPhiInSM() * fMyEBEcal.MaxCrysPhiInTow()) +
                                    iPhi;                                              // range = [1,1700]
                   Int_t n1SMTow = fMyEBNumbering.Get1SMTowFrom1SMCrys(n1SMCrys);      // range = [1,68]
                   Int_t i0TowEcha = fMyEBNumbering.Get0TowEchaFrom1SMCrys(n1SMCrys);  // range = [0,24]
 
                   Int_t NbOfSamplesFromDigis = digiItr->size();
 
                   EBDataFrame df(*digiItr);
 
                   if (NbOfSamplesFromDigis > 0 && NbOfSamplesFromDigis <= fMyEBEcal.MaxSampADC()) {
                     Double_t adcDBLS = (Double_t)0;
                     // Three 1st samples mean value for Dynamic Base Line
                     // Substraction (DBLS)
                     if (fDynBaseLineSub == "yes") {
                       for (Int_t i0Sample = 0; i0Sample < 3; i0Sample++) {
                         adcDBLS += (Double_t)(df.sample(i0Sample).adc());
                       }
                       adcDBLS /= (Double_t)3;
                     }
                     // Loop over the samples
                     for (Int_t i0Sample = 0; i0Sample < fNbOfSamples; i0Sample++) {
                       Double_t adc = (Double_t)(df.sample(i0Sample).adc()) - adcDBLS;
                       //................................................. Calls
                       // to GetSampleAdcValues
                       if (fMyCnaEBSM[i0SM]->GetSampleAdcValues(
                               fStexNbOfTreatedEvents[i0SM], n1SMTow, i0TowEcha, i0Sample, adc) == kTRUE) {
                         fBuildEventDistribGood[i0SM]++;
                       } else {
                         fBuildEventDistribBad[i0SM]++;
                       }
                     }
                   } else {
                     edm::LogVerbatim("ecnaAnal")
                         << "EcnaAnalyzer::analyze(...)> NbOfSamplesFromDigis out of bounds = " << NbOfSamplesFromDigis;
                   }
                 }  // end of if( (fStexNumber > 0 && i0SM == fStexNumber-1) ||
                    // (fStexNumber == 0) )
               }    // end of if( fStexNbOfTreatedEvents[i0SM] >= 1 &&
                    // fStexNbOfTreatedEvents[i0SM] <= fReqNbOfEvts )
             }      // end of if( fStexStatus[i0SM] < 2 )
           }        // end of if( i0SM >= 0 && i0SM<fMaxTreatedStexCounter  )
         }          // end of for (EBDigiCollection::const_iterator digiItr =
                    // digisEB->begin();
                    //             digiItr != digisEB->end(); ++digiItr)
 
         for (Int_t i0SM = 0; i0SM < fMaxTreatedStexCounter; i0SM++) {
           fStexDigiOK[i0SM] = 0;  // reset fStexDigiOK[i0SM] after loop on digis
         }
 
       }  // end of if( Int_t(digisEB->end()-digisEB->begin()) >= 0 &&
          // Int_t(digisEB->end()-digisEB->begin()) <=  Int_t(digisEB->size()) )
     }    // end of if( fStexName == "SM" && fSMIndexBegin < fSMIndexStop )
 
     //=============================================================== Record
     // type EE (Dee)
     if (fStexName == "Dee" && fDeeIndexBegin < fDeeIndexStop) {
       //......................................... Get digisEE
       const edm::Handle<EEDigiCollection> &pdigisEE = iEvent.getHandle(EEdigiToken_);
       const EEDigiCollection *digisEE = nullptr;
       if (pdigisEE.isValid()) {
         digisEE = pdigisEE.product();
       } else {
         edm::LogError("ecnaAnal") << "Error! can't get the product " << EEdigiCollection_.c_str();
       }
 
       // Initialize vectors if not already done
       if (int(digisEE->size()) > nChannels_) {
         nChannels_ = digisEE->size();
       }
 
       // Int_t print_count = 0;
       if (Int_t(digisEE->end() - digisEE->begin()) >= 0 &&
           Int_t(digisEE->end() - digisEE->begin()) <= Int_t(digisEE->size())) {
         //======================================================================================
         //
         //                           Loop over Ecal endcap digisEE (Xtals)
         //
         //======================================================================================
 
         for (EEDigiCollection::const_iterator digiItr = digisEE->begin(); digiItr != digisEE->end(); ++digiItr) {
           EEDetId id_crystal(digiItr->id());
 
           Int_t iX_data = id_crystal.ix();        // iX_data : range = [1,100]
           Int_t iY_data = id_crystal.iy();        // iY_data : range = [1,100]
           Int_t i_quad = id_crystal.iquadrant();  // iquadrant() : range = [1,4]
           Int_t i_sgnZ = id_crystal.zside();      //     zside() : values = -1,+1
 
           Int_t iX = iX_data;
           Int_t iY = iY_data;  // iY : range = [1,100]
 
           //.......... See
           // CMSSW/DataFormats/EcalDetId/src/EEDetId.cc::ixQuadrantOne()  [ in
           // which ix() = iX_data ]
           if (i_quad == 1 || i_quad == 4) {
             iX = iX_data - 50;
           }  // iX_data : range = [51,100], iX : range = [1,50]
           if (i_quad == 3 || i_quad == 2) {
             iX = 51 - iX_data;
           }  // iX_data : range = [50,1],   iX : range = [1,50]
 
           Int_t n1DeeCrys =
               (iX - 1) * (fMyEEEcal.MaxSCIYInDee() * fMyEEEcal.MaxCrysIYInSC()) + iY;  // n1DeeCrys: range = [1,5000]
 
           Int_t n1DeeNumber = 0;
           if (i_quad == 1 && i_sgnZ == 1) {
             n1DeeNumber = 2;
           }
           if (i_quad == 1 && i_sgnZ == -1) {
             n1DeeNumber = 3;
           }
           if (i_quad == 2 && i_sgnZ == 1) {
             n1DeeNumber = 1;
           }
           if (i_quad == 2 && i_sgnZ == -1) {
             n1DeeNumber = 4;
           }
           if (i_quad == 3 && i_sgnZ == 1) {
             n1DeeNumber = 1;
           }
           if (i_quad == 3 && i_sgnZ == -1) {
             n1DeeNumber = 4;
           }
           if (i_quad == 4 && i_sgnZ == 1) {
             n1DeeNumber = 2;
           }
           if (i_quad == 4 && i_sgnZ == -1) {
             n1DeeNumber = 3;
           }
 
           Int_t i0Dee = n1DeeNumber - 1;  //   <============== GET the Dee number - 1 here
 
           if (i0Dee >= 0 && i0Dee < fMaxTreatedStexCounter) {
             if (fMyCnaEEDee[i0Dee].get() == nullptr && fStexStatus[i0Dee] != 2) {
               //=============================== Init Ecna EE
               //===============================
               fMyCnaEEDee[i0Dee] = std::make_unique<TEcnaRun>(&fMyEcnaEEObjectManager, "EE", fNbOfSamples);
               fMyCnaEEDee[i0Dee]->GetReadyToReadData(
                   fAnalysisName, fRunNumber, fFirstReqEvent, fLastReqEvent, fReqNbOfEvts, i0Dee + 1, fRunTypeNumber);
 
               edm::LogVerbatim("ecnaAnal")
                   << "*EcnaAnalyzer::analyze(...)> ********* INIT ECNA EE ********* "
                   << "\n                                   fAnalysisName = " << fAnalysisName
                   << "\n                                      fRunNumber = " << fRunNumber
                   << "\n                                  fFirstReqEvent = " << fFirstReqEvent
                   << "\n                                   fLastReqEvent = " << fLastReqEvent
                   << "\n                                    fReqNbOfEvts = " << fReqNbOfEvts
                   << "\n                                             Dee = " << i0Dee + 1
                   << "\n                                        run type = " << runtype(fRunTypeNumber);
               //============================================================================
             }
 
             if (fStexStatus[i0Dee] < 2)  // nothing to do if status=2 reached
             {
               Bool_t cOKForTreatment = kFALSE;
 
               if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
                   fAnalysisName == "AdcAny") {
                 if (fFedTcc >= 1 && fFedTcc <= MaxSMAndDS) {
                   fFedDigiOK[fESFromFedTcc[fFedTcc - 1] - 1]++;
 
                   if (!(fESFromFedTcc[fFedTcc - 1] == 5 || fESFromFedTcc[fFedTcc - 1] == 14)) {
                     if (fFedDigiOK[fESFromFedTcc[fFedTcc - 1] - 1] == 1) {
                       fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1]++;
                     }
                     if (fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1] >= 1 &&
                         fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1] <= fReqNbOfEvts) {
                       fStexNbOfTreatedEvents[i0Dee] = fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1];
                       cOKForTreatment = kTRUE;
                     }
                   }
                   if (fESFromFedTcc[fFedTcc - 1] == 5 || fESFromFedTcc[fFedTcc - 1] == 14) {
                     if (fFedDigiOK[fESFromFedTcc[fFedTcc - 1] - 1] == 1) {
                       fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1]++;
                       fDeeDS5Memo1 = n1DeeNumber;
                       fStexNbOfTreatedEvents[i0Dee] = fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1];
                     } else {
                       if (fDeeDS5Memo2 == 0) {
                         if (n1DeeNumber != fDeeDS5Memo1) {
                           // change of Dee in Data sector 5
                           fDeeDS5Memo2 = n1DeeNumber;
                           fStexNbOfTreatedEvents[i0Dee] = fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1];
                         }
                       }
                     }
                     if (fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1] >= 1 &&
                         fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc - 1] - 1] <= fReqNbOfEvts) {
                       cOKForTreatment = kTRUE;
                     }
                   }
                 }  // end of if( fFedTcc >= 1 && fFedTcc <= MaxSMAndDS )
               }    // end of if( fAnalysisName == "AdcPeg12"  || fAnalysisName ==
                    // "AdcSPeg12" .... )
               else {
                 fStexDigiOK[i0Dee]++;
                 if (fStexDigiOK[i0Dee] == 1) {
                   fStexNbOfTreatedEvents[i0Dee]++;
                 }
                 if (fStexNbOfTreatedEvents[i0Dee] >= 1 && fStexNbOfTreatedEvents[i0Dee] <= fReqNbOfEvts) {
                   cOKForTreatment = kTRUE;
                 }
               }
 
               if (cOKForTreatment == kTRUE) {
                 //......................................... date of first event
                 //(in real time)
                 edm::Timestamp Time = iEvent.time();
                 edm::TimeValue_t t_current_ev_time = (cond::Time_t)Time.value();
                 time_t i_current_ev_time = (time_t)(t_current_ev_time >> 32);
                 const time_t *p_current_ev_time = &i_current_ev_time;
                 char *astime = ctime(p_current_ev_time);
 
                 if ((!(fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
                        fAnalysisName == "AdcAny") &&
                      fStexDigiOK[i0Dee] == 1 && fStexNbOfTreatedEvents[i0Dee] == 1) ||
                     ((fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
                       fAnalysisName == "AdcAny") &&
                      fFedDigiOK[fESFromFedTcc[fFedTcc - 1] - 1] == 1 && fStexNbOfTreatedEvents[i0Dee] == 1 &&
                      fMemoDateFirstEvent[i0Dee] == 0)) {
                   fTimeFirst[i0Dee] = i_current_ev_time;
                   fDateFirst[i0Dee] = astime;
                   fTimeLast[i0Dee] = i_current_ev_time;
                   fDateLast[i0Dee] = astime;
                   edm::LogVerbatim("ecnaAnal")
                       << "----- beginning of analysis for " << fStexName << i0Dee + 1 << "-------"
                       << "\n First event date  = " << astime << "\n Nb of selected evts = " << fNbOfSelectedEvents
                       << "\n---------------------------------------------------------------";
                   fMemoDateFirstEvent[i0Dee]++;
                 }
 
                 if (i_current_ev_time < fTimeFirst[i0Dee]) {
                   fTimeFirst[i0Dee] = i_current_ev_time;
                   fDateFirst[i0Dee] = astime;
                 }
                 if (i_current_ev_time > fTimeLast[i0Dee]) {
                   fTimeLast[i0Dee] = i_current_ev_time;
                   fDateLast[i0Dee] = astime;
                 }
 
                 //=============================================> cut on i0Dee
                 // value
                 if ((fStexNumber > 0 && i0Dee == fStexNumber - 1) || (fStexNumber == 0)) {
                   TString sDir = fMyEENumbering.GetDeeDirViewedFromIP(n1DeeNumber);
                   Int_t n1DeeSCEcna = fMyEENumbering.Get1DeeSCEcnaFrom1DeeCrys(n1DeeCrys, sDir);
                   Int_t i0SCEcha = fMyEENumbering.Get1SCEchaFrom1DeeCrys(n1DeeCrys, sDir) - 1;
 
                   Int_t NbOfSamplesFromDigis = digiItr->size();
 
                   EEDataFrame df(*digiItr);
 
                   if (NbOfSamplesFromDigis > 0 && NbOfSamplesFromDigis <= fMyEEEcal.MaxSampADC()) {
                     Double_t adcDBLS = (Double_t)0;
                     // Three 1st samples mean value for Dynamic Base Line
                     // Substraction (DBLS)
                     if (fDynBaseLineSub == "yes") {
                       for (Int_t i0Sample = 0; i0Sample < 3; i0Sample++) {
                         adcDBLS += (Double_t)(df.sample(i0Sample).adc());
                       }
                       adcDBLS /= (Double_t)3;
                     }
                     // Loop over the samples
                     for (Int_t i0Sample = 0; i0Sample < fNbOfSamples; i0Sample++) {
                       Double_t adc = (Double_t)(df.sample(i0Sample).adc()) - adcDBLS;
                       //................................................. Calls
                       // to GetSampleAdcValues
                       if (fMyCnaEEDee[i0Dee]->GetSampleAdcValues(
                               fStexNbOfTreatedEvents[i0Dee], n1DeeSCEcna, i0SCEcha, i0Sample, adc) == kTRUE) {
                         fBuildEventDistribGood[i0Dee]++;
                       } else {
                         fBuildEventDistribBad[i0Dee]++;
                       }
                     }
                   } else {
                     edm::LogVerbatim("ecnaAnal")
                         << "EcnaAnalyzer::analyze(...)> NbOfSamplesFromDigis out of bounds = " << NbOfSamplesFromDigis;
                   }
                 }  // end of if( (fStexNumber > 0 && i0Dee == fStexNumber-1) ||
                    // (fStexNumber == 0) )
               }    // end of if( fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc-1]-1]
                    // >= 1 &&
               // fFedNbOfTreatedEvents[fESFromFedTcc[fFedTcc-1]-1] <=
               // fReqNbOfEvts )
             }  // end of if( fStexStatus[i0Dee] < 2 )
           }    // end of if( i0Dee >= 0 && i0Dee<fMaxTreatedStexCounter )
         }      // end of for (EBDigiCollection::const_iterator digiItr =
                // digisEB->begin();
                //             digiItr != digisEB->end(); ++digiItr)
 
         // reset fStexDigiOK[i0Dee] or fFedDigiOK[i0Dee] to zero after loop on
         // digis
         if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
             fAnalysisName == "AdcAny") {
           for (Int_t i0FedES = 0; i0FedES < fMaxFedUnitCounter; i0FedES++) {
             fFedDigiOK[i0FedES] = 0;
           }
 
           // reset fDeeDS5Memo1 and fDeeDS5Memo2 (for Data sector 5 versus Dee
           // number  management)
           fDeeDS5Memo1 = 0;
           fDeeDS5Memo2 = 0;
         } else {
           for (Int_t i0Dee = 0; i0Dee < fMaxTreatedStexCounter; i0Dee++) {
             fStexDigiOK[i0Dee] = 0;
           }
         }
 
       }  // end of if( Int_t(digisEB->end()-digisEB->begin()) >= 0 &&
          // Int_t(digisEB->end()-digisEB->begin()) <=  Int_t(digisEB->size()) )
 
     }  // end of if( fStexName == "Dee" && fDeeIndexBegin < fDeeIndexStop )
   }    // end of if( (fStexNumber > 0 && fNbOfTreatedStexs == 0) || (fStexNumber ==
   // 0 && fNbOfTreatedStexs < MaxNbOfStex) )
 
   //=============================================================================================
   //
   //                    Number of treated events. Setting Stex and Fed status.
   //
   //=============================================================================================
 
   // (take into account the "Accelerating selection with FED number" section -
   // see above -)
   if (fStexName == "SM" || (fStexName == "Dee" &&
                             !(fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" ||
                               fAnalysisName == "AdcPhys" || fAnalysisName == "AdcAny")))  // one FED = one SM = one Stex
   {
     for (Int_t i0Stex = fStexIndexBegin; i0Stex < fStexIndexStop; i0Stex++) {
       if (fStexStatus[i0Stex] != 2)  // do not change fStexStatus[i0Stex] if already set to 2
                                      // even if fStexNbOfTreatedEvents[i0Stex] == fReqNbOfEvts
       {
         if (fStexNbOfTreatedEvents[i0Stex] == fReqNbOfEvts) {
           fStexStatus[i0Stex] = 1;
         }
         if (fStexNbOfTreatedEvents[i0Stex] > fReqNbOfEvts) {
           fStexStatus[i0Stex] = 2;
         }
       }
     }
   }
 
   // one FED = one Data Sector (DS or ES)
   if (fStexName == "Dee" && (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" ||
                              fAnalysisName == "AdcPhys" || fAnalysisName == "AdcAny")) {
     for (Int_t i0FedES = 0; i0FedES < fMaxFedUnitCounter; i0FedES++) {
       if (fFedStatus[i0FedES] != 2)  // do not change fFedStatus[i0FedES] if already set to 2
                                      // even if fFedNbOfTreatedEvents[i0FedES] == fReqNbOfEvts
       {
         if (fFedNbOfTreatedEvents[i0FedES] == fReqNbOfEvts) {
           fFedStatus[i0FedES] = 1;
           fTreatedFedOrder++;
           fFedStatusOrder[i0FedES] = fTreatedFedOrder;
         }
         if (fFedNbOfTreatedEvents[i0FedES] > fReqNbOfEvts) {
           fFedStatus[i0FedES] = 2;
         }
       }
     }
 
     Int_t j0Fed = 4;
     //..................................................... Dee 4 (DS 1,2,3,4 ;
     // ES 1,2,3,4)
     for (Int_t i0FedES = 0; i0FedES <= 3; i0FedES++) {
       if (fFedStatus[i0FedES] == 1) {
         fNbOfTreatedFedsInDee[3]++;
         fFedStatus[i0FedES] = 2;
       }
     }
 
     //..................................................... Dee 3, Dee 4 (DS 5 ;
     // ES 5)
     j0Fed = 4;
     if (fFedStatus[j0Fed] == 1) {
       fNbOfTreatedFedsInDee[3]++;
       fNbOfTreatedFedsInDee[2]++;
       fFedStatus[j0Fed] = 2;
     }
 
     //.................................................... Dee 3 (DS 6,7,8,9 ;
     // ES 6,7,8,9)
     for (Int_t i0FedES = 5; i0FedES <= 8; i0FedES++) {
       if (fFedStatus[i0FedES] == 1) {
         fNbOfTreatedFedsInDee[2]++;
         fFedStatus[i0FedES] = 2;
       }
     }
 
     //..................................................... Dee 1 (DS 1,2,3,4 ;
     // ES 10,11,12,13)
     for (Int_t i0FedES = 9; i0FedES <= 12; i0FedES++) {
       if (fFedStatus[i0FedES] == 1) {
         fNbOfTreatedFedsInDee[0]++;
         fFedStatus[i0FedES] = 2;
       }
     }
 
     //..................................................... Dee 1, Dee 2 (DS 5 ;
     // ES 5)
     j0Fed = 13;
     if (fFedStatus[j0Fed] == 1) {
       fNbOfTreatedFedsInDee[0]++;
       fNbOfTreatedFedsInDee[1]++;
       fFedStatus[j0Fed] = 2;
     }
 
     //..................................................... Dee 2 (DS 6,7,8,9 ;
     // ES 15,16,17,18)
     for (Int_t i0FedES = 14; i0FedES <= 17; i0FedES++) {
       if (fFedStatus[i0FedES] == 1) {
         fNbOfTreatedFedsInDee[1]++;
         fFedStatus[i0FedES] = 2;
       }
     }
 
     //-----------------------------------------------------
     for (Int_t i0Dee = 0; i0Dee < 4; i0Dee++) {
       if (fNbOfTreatedFedsInStex[i0Dee] >= 0 && fNbOfTreatedFedsInStex[i0Dee] < 5) {
         fNbOfTreatedFedsInStex[i0Dee] = fNbOfTreatedFedsInDee[i0Dee];
       }
       if (fNbOfTreatedFedsInDee[i0Dee] == 5) {
         fStexStatus[i0Dee] = 1;
         fNbOfTreatedFedsInDee[i0Dee] = 0;
       }
     }
 
   }  // end of if( fStexName == "Dee" &&
   // ( fAnalysisName == "AdcPeg12"  || fAnalysisName == "AdcSPeg12" ... ) )
 
   //----------------------------------------------------------------------------------------------
   for (Int_t i0Stex = fStexIndexBegin; i0Stex < fStexIndexStop; i0Stex++) {
     if (fStexStatus[i0Stex] == 1) {
       fNbOfTreatedStexs++;  // increase nb of treated Stex's only if
                             // fStexStatus[i0Stex] == 1
       //....................................................... date of last
       // event edm::Timestamp Time = iEvent.time(); edm::TimeValue_t
       // t_current_ev_time = (cond::Time_t)Time.value(); time_t
       // i_current_ev_time = (time_t)(t_current_ev_time>>32); const time_t*
       // p_current_ev_time = &i_current_ev_time; char*          astime =
       // ctime(p_current_ev_time); fTimeLast[i0Stex] = i_current_ev_time;
       // fDateLast[i0Stex] = astime;
 
       // if( i_current_ev_time > fTimeLast[i0Stex] )
       // {fTimeLast[i0Stex] = i_current_ev_time; fDateLast[i0Stex] = astime;}
 
       edm::LogVerbatim("ecnaAnal") << "---------- End of analysis for " << fStexName << i0Stex + 1 << " -----------";
       Int_t n3 = 3;
       CheckMsg(n3, i0Stex);
       edm::LogVerbatim("ecnaAnal") << " Number of selected events = " << fNbOfSelectedEvents;
       edm::LogVerbatim("ecnaAnal") << std::endl
                                    << fNbOfTreatedStexs << " " << fStexName << "'s with " << fReqNbOfEvts
                                    << " events analyzed."
                                    << "\n---------------------------------------------------------";
 
       //================================= WRITE RESULTS FILE
       if (fStexName == "SM") {
         if (fMyCnaEBSM[i0Stex].get() != nullptr) {
           //........................................ register dates 1 and 2
           fMyCnaEBSM[i0Stex]->StartStopDate(fDateFirst[i0Stex], fDateLast[i0Stex]);
           fMyCnaEBSM[i0Stex]->StartStopTime(fTimeFirst[i0Stex], fTimeLast[i0Stex]);
 
           //........................................ Init .root file
           fMyCnaEBSM[i0Stex]->GetReadyToCompute();
           fMyCnaEBSM[i0Stex]->SampleValues();
 
           //........................................ write the sample values in
           //.root file
           if (fMyCnaEBSM[i0Stex]->WriteRootFile() == kFALSE) {
             edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer::analyze> PROBLEM with write ROOT file for SM" << i0Stex + 1
                                          << fTTBELL;
           }
         }
         // set pointer to zero in order to avoid recalculation and rewriting at
         // the destructor level
         fMyCnaEBSM[i0Stex].reset();
         edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer::analyze> Set memory free: delete done for SM " << i0Stex + 1;
       }
 
       if (fStexName == "Dee") {
         if (fMyCnaEEDee[i0Stex].get() != nullptr) {
           //........................................ register dates 1 and 2
           fMyCnaEEDee[i0Stex]->StartStopDate(fDateFirst[i0Stex], fDateLast[i0Stex]);
           fMyCnaEEDee[i0Stex]->StartStopTime(fTimeFirst[i0Stex], fTimeLast[i0Stex]);
 
           //........................................ Init .root file
           fMyCnaEEDee[i0Stex]->GetReadyToCompute();
           fMyCnaEEDee[i0Stex]->SampleValues();
 
           //........................................ write the sample values in
           //.root file
           if (fMyCnaEEDee[i0Stex]->WriteRootFile() == kFALSE) {
             edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer::analyze> PROBLEM with write ROOT file for Dee" << i0Stex + 1
                                          << " " << fTTBELL;
           }
         }
         // set pointer to zero in order to avoid recalculation and rewriting at
         // the destructor level
         fMyCnaEEDee[i0Stex].reset();
         edm::LogVerbatim("ecnaAnal") << "!EcnaAnalyzer::analyze> Set memory free: delete done for Dee " << i0Stex + 1;
       }
 
       fStexStatus[i0Stex] = 2;  // set fStexStatus[i0Stex] to 2 definitively
       edm::LogVerbatim("ecnaAnal") << "*---------------------------------------------------------------------------- ";
 
     }  // end of if( fStexStatus[i0Stex] == 1 )
   }    // end of for(Int_t i0Stex=fStexIndexBegin; i0Stex<fStexIndexStop; i0Stex++)
 }

◆ CheckMsg() [1/2]

void EcnaAnalyzer::CheckMsg	(	const Int_t &	MsgNum,
		const Int_t &	i0Stex
	)

Definition at line 1427 of file EcnaAnalyzer.cc.

References fAnalysisName, fChozenGainNumber, fChozenRunTypeNumber, fCurrentEventNumber, fDeeNumberString, fEvtNumber, fFedId, fFedNbOfTreatedEvents, fFedStatus, fFedStatusOrder, fFedTcc, fMaxFedUnitCounter, fMgpaGainNumber, fNbOfSelectedEvents, fNbOfTreatedFedsInStex, fNbOfTreatedStexs, fRecNumber, fReqNbOfEvts, fRunNumber, fRunTypeNumber, fStexIndexBegin, fStexIndexStop, fStexName, fStexNbOfTreatedEvents, fStexStatus, gainvalue(), and runtype().

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), and ~EcnaAnalyzer().

                                                                     {
   //------ Cross-check messages
 
   if (MsgNum == 1) {
     edm::LogVerbatim("ecnaAnal") << "---------------- CROSS-CHECK A ------------------ "
                                  << "\n**************** CURRENT EVENT ****************** ";
   }
   if (MsgNum == 2) {
     edm::LogVerbatim("ecnaAnal") << "---------------- CROSS-CHECK B ------------------ "
                                  << "\n**** FIRST EVENT PASSING USER'S ANALYSIS CUT **** ";
   }
   if (MsgNum == 3) {
     edm::LogVerbatim("ecnaAnal") << "---------------- CROSS-CHECK C ------------------ "
                                  << "\n*** CURRENT VALUES BEFORE RESULT FILE WRITING *** ";
   }
   if (MsgNum == 3 || MsgNum == 4) {
     edm::LogVerbatim("ecnaAnal") << "          fRecNumber = " << fRecNumber
                                  << "\n          fEvtNumber = " << fEvtNumber;
   }
 
   edm::LogVerbatim("ecnaAnal") << " fCurrentEventNumber = " << fCurrentEventNumber
                                << "\n fNbOfSelectedEvents = " << fNbOfSelectedEvents
                                << "\n          fRunNumber = " << fRunNumber
                                << "\n     Chozen run type = " << runtype(fChozenRunTypeNumber)
                                << "\n            Run type = " << runtype(fRunTypeNumber)
                                << "\n             fFedTcc = " << fFedTcc << "\n        fFedId(+601) = " << fFedId + 601
                                << "\n           fStexName = " << fStexName
                                << "\n         Chozen gain = " << gainvalue(fChozenGainNumber)
                                << "\n           Mgpa Gain = " << gainvalue(fMgpaGainNumber) << std::endl;
 
   if (fAnalysisName == "AdcPeg12" || fAnalysisName == "AdcSPeg12" || fAnalysisName == "AdcPhys" ||
       fAnalysisName == "AdcAny") {
     if (fStexName == "SM") {
       for (Int_t j0Stex = fStexIndexBegin; j0Stex < fStexIndexStop; j0Stex++) {
         Int_t nStexNbOfTreatedEvents = fStexNbOfTreatedEvents[j0Stex];
         if (fStexStatus[j0Stex] == 1) {
           nStexNbOfTreatedEvents = fStexNbOfTreatedEvents[j0Stex];
         }
         if (fStexStatus[j0Stex] == 2) {
           nStexNbOfTreatedEvents = fStexNbOfTreatedEvents[j0Stex];
         }
 
         edm::LogVerbatim("ecnaAnal") << fStexName << std::setw(3) << j0Stex + 1 << ": " << std::setw(5)
                                      << nStexNbOfTreatedEvents << " events. " << fStexName
                                      << " status: " << fStexStatus[j0Stex];
         if (j0Stex == i0Stex) {
           edm::LogVerbatim("ecnaAnal") << " (going to write file for this " << fStexName << ").";
         }
       }
     }
 
     if (fStexName == "Dee") {
       for (Int_t i0FedES = 0; i0FedES < fMaxFedUnitCounter; i0FedES++) {
         Int_t nFedNbOfTreatedEvents = fFedNbOfTreatedEvents[i0FedES];
         if (fFedStatus[i0FedES] == 1) {
           nFedNbOfTreatedEvents = fFedNbOfTreatedEvents[i0FedES];
         }
         if (fFedStatus[i0FedES] == 2) {
           nFedNbOfTreatedEvents = fFedNbOfTreatedEvents[i0FedES];
         }
 
         edm::LogVerbatim("ecnaAnal") << "Fed (ES) " << std::setw(3) << i0FedES + 1 << ": " << std::setw(5)
                                      << nFedNbOfTreatedEvents << " events."
                                      << " Fed status: " << fFedStatus[i0FedES] << ", order: " << std::setw(3)
                                      << fFedStatusOrder[i0FedES] << " (" << fDeeNumberString[i0FedES] << ")";
       }
 
       for (Int_t j0Stex = fStexIndexBegin; j0Stex < fStexIndexStop; j0Stex++) {
         edm::LogVerbatim("ecnaAnal") << fStexName << std::setw(3) << j0Stex + 1 << ": " << std::setw(5)
                                      << fNbOfTreatedFedsInStex[j0Stex] << " analyzed Fed(s). " << fStexName
                                      << " status: " << fStexStatus[j0Stex];
         if (j0Stex == i0Stex) {
           edm::LogVerbatim("ecnaAnal") << " (going to write file for this " << fStexName << ").";
         }
       }
     }
 
     edm::LogVerbatim("ecnaAnal") << "Number of " << fStexName << "'s with " << fReqNbOfEvts
                                  << " events analyzed: " << fNbOfTreatedStexs;
   }
 
   if (MsgNum == 1 || MsgNum == 2) {
     edm::LogVerbatim("ecnaAnal") << "*---------------------------------------------------------------------------- ";
   }
   if (MsgNum == 3) {
     edm::LogVerbatim("ecnaAnal") << "*............................................................................ ";
   }
 
 }  // end of EcnaAnalyzer::CheckMsg(const Int_t& MsgNum, const Int_t& i0Stex)

◆ CheckMsg() [2/2]

void EcnaAnalyzer::CheckMsg ( const Int_t & MsgNum )

Definition at line 1422 of file EcnaAnalyzer.cc.

References CheckMsg().

                                                {
   Int_t nm1 = -1;
   CheckMsg(MsgNum, nm1);
 }

◆ gainvalue()

Int_t EcnaAnalyzer::gainvalue ( const Int_t & numgain )

Definition at line 1607 of file EcnaAnalyzer.cc.

References relativeConstraints::value.

Referenced by CheckMsg(), and ~EcnaAnalyzer().

                                                   {
   Int_t value = 0;
 
   if (numgain == 1) {
     value = 12;
   }
   if (numgain == 2) {
     value = 6;
   }
   if (numgain == 3) {
     value = 1;
   }
 
   return value;
 }

◆ runtype()

TString EcnaAnalyzer::runtype ( const Int_t & numtype )

Definition at line 1517 of file EcnaAnalyzer.cc.

Referenced by analyze(), CheckMsg(), and ~EcnaAnalyzer().

                                                   {
   TString cType = "?";
 
   if (numtype == 0) {
     cType = "COSMICS";
   }
   if (numtype == 1) {
     cType = "BEAMH4";
   }
   if (numtype == 2) {
     cType = "BEAMH2";
   }
   if (numtype == 3) {
     cType = "MTCC";
   }
   if (numtype == 4) {
     cType = "LASER_STD";
   }
   if (numtype == 5) {
     cType = "LASER_POWER_SCAN";
   }
   if (numtype == 6) {
     cType = "LASER_DELAY_SCAN";
   }
   if (numtype == 7) {
     cType = "TESTPULSE_SCAN_MEM";
   }
   if (numtype == 8) {
     cType = "TESTPULSE_MGPA";
   }
   if (numtype == 9) {
     cType = "PEDESTAL_STD";
   }
   if (numtype == 10) {
     cType = "PEDESTAL_OFFSET_SCAN";
   }
   if (numtype == 11) {
     cType = "PEDESTAL_25NS_SCAN";
   }
   if (numtype == 12) {
     cType = "LED_STD";
   }
 
   if (numtype == 13) {
     cType = "PHYSICS_GLOBAL";
   }
   if (numtype == 14) {
     cType = "COSMICS_GLOBAL";
   }
   if (numtype == 15) {
     cType = "HALO_GLOBAL";
   }
 
   if (numtype == 16) {
     cType = "LASER_GAP";
   }
   if (numtype == 17) {
     cType = "TESTPULSE_GAP";
   }
   if (numtype == 18) {
     cType = "PEDESTAL_GAP";
   }
   if (numtype == 19) {
     cType = "LED_GAP";
   }
 
   if (numtype == 20) {
     cType = "PHYSICS_LOCAL";
   }
   if (numtype == 21) {
     cType = "COSMICS_LOCAL";
   }
   if (numtype == 22) {
     cType = "HALO_LOCAL";
   }
   if (numtype == 23) {
     cType = "CALIB_LOCAL";
   }
 
   //.......................................... non-CMS types
   if (numtype == 24) {
     cType = "PEDSIM";
   }  // SIMULATION
   if (numtype == 25) {
     cType = "ANY_RUN";
   }  // ANY RUN (ALL TYPES ACCEPTED)
 
   return cType;
 }

Member Data Documentation

◆ digiProducer_

std::string EcnaAnalyzer::digiProducer_

private

Definition at line 156 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ EBdigiCollection_

std::string EcnaAnalyzer::EBdigiCollection_

private

Definition at line 158 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ EBdigiToken_

edm::EDGetTokenT<EBDigiCollection> EcnaAnalyzer::EBdigiToken_

private

Definition at line 161 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ EEdigiCollection_

std::string EcnaAnalyzer::EEdigiCollection_

private

Definition at line 159 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ EEdigiToken_

edm::EDGetTokenT<EEDigiCollection> EcnaAnalyzer::EEdigiToken_

private

Definition at line 162 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ eventHeaderCollection_

std::string EcnaAnalyzer::eventHeaderCollection_

private

Definition at line 157 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ eventHeaderProducer_

std::string EcnaAnalyzer::eventHeaderProducer_

private

Definition at line 155 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ eventHeaderToken_

edm::EDGetTokenT<EcalRawDataCollection> EcnaAnalyzer::eventHeaderToken_

private

Definition at line 160 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fAnalysisName

TString EcnaAnalyzer::fAnalysisName

private

Definition at line 188 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fANY_RUN

Int_t EcnaAnalyzer::fANY_RUN

private

Definition at line 236 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fBuildEventDistribBad

std::vector<Int_t> EcnaAnalyzer::fBuildEventDistribBad

private

Definition at line 180 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fBuildEventDistribGood

std::vector<Int_t> EcnaAnalyzer::fBuildEventDistribGood

private

Definition at line 181 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fCfgAnalyzerParametersFileName

TString EcnaAnalyzer::fCfgAnalyzerParametersFileName

private

Definition at line 185 of file EcnaAnalyzer.h.

◆ fCfgAnalyzerParametersFilePath

TString EcnaAnalyzer::fCfgAnalyzerParametersFilePath

private

Definition at line 183 of file EcnaAnalyzer.h.

◆ fChozenGainNumber

Int_t EcnaAnalyzer::fChozenGainNumber

private

Definition at line 189 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fChozenRunTypeNumber

Int_t EcnaAnalyzer::fChozenRunTypeNumber

private

Definition at line 190 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fCurrentEventNumber

Int_t EcnaAnalyzer::fCurrentEventNumber

private

Definition at line 177 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fDateFirst

std::vector<TString> EcnaAnalyzer::fDateFirst

private

Definition at line 245 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fDateLast

std::vector<TString> EcnaAnalyzer::fDateLast

private

Definition at line 246 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fDeeDS5Memo1

Int_t EcnaAnalyzer::fDeeDS5Memo1

private

Definition at line 220 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fDeeDS5Memo2

Int_t EcnaAnalyzer::fDeeDS5Memo2

private

Definition at line 221 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fDeeIndexBegin

Int_t EcnaAnalyzer::fDeeIndexBegin

private

Definition at line 206 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fDeeIndexStop

Int_t EcnaAnalyzer::fDeeIndexStop

private

Definition at line 207 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fDeeNumberString

std::vector<std::string> EcnaAnalyzer::fDeeNumberString

private

Definition at line 217 of file EcnaAnalyzer.h.

Referenced by CheckMsg(), and EcnaAnalyzer().

◆ fDynBaseLineSub

TString EcnaAnalyzer::fDynBaseLineSub

private

Definition at line 191 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fESFromFedTcc

std::vector<Int_t> EcnaAnalyzer::fESFromFedTcc

private

Definition at line 213 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fEvtNumber

Int_t EcnaAnalyzer::fEvtNumber

private

Definition at line 174 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fEvtNumberMemo

Int_t EcnaAnalyzer::fEvtNumberMemo

private

Definition at line 175 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fFcin_f

std::ifstream EcnaAnalyzer::fFcin_f

private

Definition at line 186 of file EcnaAnalyzer.h.

◆ fFedDigiOK

std::vector<Int_t> EcnaAnalyzer::fFedDigiOK

private

Definition at line 228 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fFedId

Int_t EcnaAnalyzer::fFedId

private

Definition at line 216 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fFedIdCounter

std::array<Int_t, fMaxFedIdCounter> EcnaAnalyzer::fFedIdCounter

private

Definition at line 272 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fFedNbOfTreatedEvents

std::vector<Int_t> EcnaAnalyzer::fFedNbOfTreatedEvents

private

Definition at line 229 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fFedStatus

std::vector<Int_t> EcnaAnalyzer::fFedStatus

private

Definition at line 227 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fFedStatusOrder

std::vector<Int_t> EcnaAnalyzer::fFedStatusOrder

private

Definition at line 215 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fFedTcc

Int_t EcnaAnalyzer::fFedTcc

private

Definition at line 211 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fFirstReqEvent

Int_t EcnaAnalyzer::fFirstReqEvent

private

Definition at line 196 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), and EcnaAnalyzer().

◆ fgMaxCar

constexpr Int_t EcnaAnalyzer::fgMaxCar = 512

staticprivate

Definition at line 148 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fLASER_STD

Int_t EcnaAnalyzer::fLASER_STD

private

Definition at line 239 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fLastReqEvent

Int_t EcnaAnalyzer::fLastReqEvent

private

Definition at line 197 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), and EcnaAnalyzer().

◆ fMaxFedIdCounter

constexpr Int_t EcnaAnalyzer::fMaxFedIdCounter = 54

staticprivate

Definition at line 271 of file EcnaAnalyzer.h.

Referenced by analyze(), and ~EcnaAnalyzer().

◆ fMaxFedUnitCounter

Int_t EcnaAnalyzer::fMaxFedUnitCounter = 0

private

Definition at line 226 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fMaxMgpaGainCounter

constexpr Int_t EcnaAnalyzer::fMaxMgpaGainCounter = 4

staticprivate

Definition at line 268 of file EcnaAnalyzer.h.

Referenced by analyze(), and ~EcnaAnalyzer().

◆ fMaxRunTypeCounter

constexpr Int_t EcnaAnalyzer::fMaxRunTypeCounter = 26

staticprivate

Definition at line 265 of file EcnaAnalyzer.h.

Referenced by analyze(), and ~EcnaAnalyzer().

◆ fMaxTreatedStexCounter

Int_t EcnaAnalyzer::fMaxTreatedStexCounter = 0

private

Definition at line 219 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fMemoCutOK

Int_t EcnaAnalyzer::fMemoCutOK

private

Definition at line 231 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fMemoDateFirstEvent

std::vector<Int_t> EcnaAnalyzer::fMemoDateFirstEvent

private

Definition at line 248 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fMgpaGainCounter

std::array<Int_t, fMaxMgpaGainCounter> EcnaAnalyzer::fMgpaGainCounter

private

Definition at line 269 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fMgpaGainNumber

Int_t EcnaAnalyzer::fMgpaGainNumber

private

Definition at line 202 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fMyCnaEBSM

std::vector<std::unique_ptr<TEcnaRun> > EcnaAnalyzer::fMyCnaEBSM

private

Definition at line 253 of file EcnaAnalyzer.h.

Referenced by analyze(), and ~EcnaAnalyzer().

◆ fMyCnaEEDee

std::vector<std::unique_ptr<TEcnaRun> > EcnaAnalyzer::fMyCnaEEDee

private

Definition at line 254 of file EcnaAnalyzer.h.

Referenced by analyze(), and ~EcnaAnalyzer().

◆ fMyEBEcal

TEcnaParEcal EcnaAnalyzer::fMyEBEcal

private

Definition at line 257 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), and EcnaAnalyzer().

◆ fMyEBNumbering

TEcnaNumbering EcnaAnalyzer::fMyEBNumbering

private

Definition at line 256 of file EcnaAnalyzer.h.

Referenced by analyze().

◆ fMyEcnaEBObjectManager

TEcnaObject EcnaAnalyzer::fMyEcnaEBObjectManager

private

Definition at line 250 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fMyEcnaEEObjectManager

TEcnaObject EcnaAnalyzer::fMyEcnaEEObjectManager

private

Definition at line 251 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fMyEEEcal

TEcnaParEcal EcnaAnalyzer::fMyEEEcal

private

Definition at line 260 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), and EcnaAnalyzer().

◆ fMyEENumbering

TEcnaNumbering EcnaAnalyzer::fMyEENumbering

private

Definition at line 259 of file EcnaAnalyzer.h.

Referenced by analyze().

◆ fNbOfSamples

Int_t EcnaAnalyzer::fNbOfSamples

private

Definition at line 193 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fNbOfSelectedEvents

Int_t EcnaAnalyzer::fNbOfSelectedEvents

private

Definition at line 178 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fNbOfTreatedFedsInDee

std::vector<Int_t> EcnaAnalyzer::fNbOfTreatedFedsInDee

private

Definition at line 233 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fNbOfTreatedFedsInStex

std::vector<Int_t> EcnaAnalyzer::fNbOfTreatedFedsInStex

private

Definition at line 234 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fNbOfTreatedStexs

Int_t EcnaAnalyzer::fNbOfTreatedStexs

private

Definition at line 232 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fOutcomeError

Bool_t EcnaAnalyzer::fOutcomeError

private

Definition at line 172 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fPEDESTAL_GAP

Int_t EcnaAnalyzer::fPEDESTAL_GAP

private

Definition at line 238 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fPEDESTAL_STD

Int_t EcnaAnalyzer::fPEDESTAL_STD

private

Definition at line 237 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fPEDSIM

Int_t EcnaAnalyzer::fPEDSIM

private

Definition at line 241 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fPHYSICS_GLOBAL

Int_t EcnaAnalyzer::fPHYSICS_GLOBAL

private

Definition at line 240 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ fRecNumber

Int_t EcnaAnalyzer::fRecNumber

private

Definition at line 176 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fReqNbOfEvts

Int_t EcnaAnalyzer::fReqNbOfEvts

private

Definition at line 201 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fRunNumber

Int_t EcnaAnalyzer::fRunNumber

private

Definition at line 194 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fRunTypeCounter

std::array<Int_t, fMaxRunTypeCounter> EcnaAnalyzer::fRunTypeCounter

private

Definition at line 266 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fRunTypeNumber

Int_t EcnaAnalyzer::fRunTypeNumber

private

Definition at line 195 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fSMFromFedTcc

std::vector<Int_t> EcnaAnalyzer::fSMFromFedTcc

private

Definition at line 212 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fSMIndexBegin

Int_t EcnaAnalyzer::fSMIndexBegin

private

Definition at line 204 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fSMIndexStop

Int_t EcnaAnalyzer::fSMIndexStop

private

Definition at line 205 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fStexDigiOK

std::vector<Int_t> EcnaAnalyzer::fStexDigiOK

private

Definition at line 222 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fStexIndexBegin

Int_t EcnaAnalyzer::fStexIndexBegin

private

Definition at line 208 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fStexIndexStop

Int_t EcnaAnalyzer::fStexIndexStop

private

Definition at line 209 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fStexName

TString EcnaAnalyzer::fStexName

private

Definition at line 198 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), CheckMsg(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fStexNbOfTreatedEvents

std::vector<Int_t> EcnaAnalyzer::fStexNbOfTreatedEvents

private

Definition at line 223 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fStexNumber

Int_t EcnaAnalyzer::fStexNumber

private

Definition at line 199 of file EcnaAnalyzer.h.

Referenced by AnalysisOutcome(), analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fStexStatus

std::vector<Int_t> EcnaAnalyzer::fStexStatus

private

Definition at line 224 of file EcnaAnalyzer.h.

Referenced by analyze(), CheckMsg(), and EcnaAnalyzer().

◆ fTimeFirst

std::vector<time_t> EcnaAnalyzer::fTimeFirst

private

Definition at line 243 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fTimeLast

std::vector<time_t> EcnaAnalyzer::fTimeLast

private

Definition at line 244 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ fTreatedFedOrder

Int_t EcnaAnalyzer::fTreatedFedOrder

private

Definition at line 214 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

◆ fTTBELL

TString EcnaAnalyzer::fTTBELL

private

Definition at line 149 of file EcnaAnalyzer.h.

Referenced by analyze(), EcnaAnalyzer(), and ~EcnaAnalyzer().

◆ iEvent_

Int_t EcnaAnalyzer::iEvent_

private

Definition at line 154 of file EcnaAnalyzer.h.

◆ nChannels_

Int_t EcnaAnalyzer::nChannels_

private

Definition at line 153 of file EcnaAnalyzer.h.

Referenced by analyze().

◆ sAnalysisName_

TString EcnaAnalyzer::sAnalysisName_

private

Definition at line 164 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sFirstReqEvent_

TString EcnaAnalyzer::sFirstReqEvent_

private

Definition at line 166 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sLastReqEvent_

TString EcnaAnalyzer::sLastReqEvent_

private

Definition at line 167 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sNbOfSamples_

TString EcnaAnalyzer::sNbOfSamples_

private

Definition at line 165 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sReqNbOfEvts_

TString EcnaAnalyzer::sReqNbOfEvts_

private

Definition at line 168 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sStexName_

TString EcnaAnalyzer::sStexName_

private

Definition at line 169 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ sStexNumber_

TString EcnaAnalyzer::sStexNumber_

private

Definition at line 170 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

◆ verbosity_

unsigned int EcnaAnalyzer::verbosity_

private

Definition at line 152 of file EcnaAnalyzer.h.

Public Types

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

cna_new/Correlated_Noise_Analysis/ECNA_cna_1.htm

--------------------------------— Analysis name codes

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ EcnaAnalyzer()

◆ ~EcnaAnalyzer()

Member Function Documentation

◆ AnalysisOutcome()

◆ analyze()

◆ CheckMsg() [1/2]

◆ CheckMsg() [2/2]

◆ gainvalue()

◆ runtype()

Member Data Documentation

◆ digiProducer_

◆ EBdigiCollection_

◆ EBdigiToken_

◆ EEdigiCollection_

◆ EEdigiToken_

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ eventHeaderToken_

◆ fAnalysisName

◆ fANY_RUN

◆ fBuildEventDistribBad

◆ fBuildEventDistribGood

◆ fCfgAnalyzerParametersFileName

◆ fCfgAnalyzerParametersFilePath

◆ fChozenGainNumber

◆ fChozenRunTypeNumber

◆ fCurrentEventNumber

◆ fDateFirst

◆ fDateLast

◆ fDeeDS5Memo1

◆ fDeeDS5Memo2

◆ fDeeIndexBegin

◆ fDeeIndexStop

◆ fDeeNumberString

◆ fDynBaseLineSub

◆ fESFromFedTcc

◆ fEvtNumber

◆ fEvtNumberMemo

◆ fFcin_f

◆ fFedDigiOK

◆ fFedId

◆ fFedIdCounter

◆ fFedNbOfTreatedEvents

◆ fFedStatus

◆ fFedStatusOrder

◆ fFedTcc

◆ fFirstReqEvent

◆ fgMaxCar

◆ fLASER_STD

◆ fLastReqEvent

◆ fMaxFedIdCounter

◆ fMaxFedUnitCounter

◆ fMaxMgpaGainCounter

◆ fMaxRunTypeCounter

◆ fMaxTreatedStexCounter

◆ fMemoCutOK

◆ fMemoDateFirstEvent

◆ fMgpaGainCounter

◆ fMgpaGainNumber

◆ fMyCnaEBSM

◆ fMyCnaEEDee

◆ fMyEBEcal

◆ fMyEBNumbering

◆ fMyEcnaEBObjectManager

◆ fMyEcnaEEObjectManager

◆ fMyEEEcal

◆ fMyEENumbering

◆ fNbOfSamples

◆ fNbOfSelectedEvents