EcnaAnalyzer Class Reference

#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 final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
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.

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

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.

References gpuClustering::adc, AnalysisOutcome(), edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), CheckMsg(), gather_cfg::cout, TauDecayModes::dec, hgcalPerformanceValidation::df, EBdigiCollection_, EBdigiToken_, EEdigiCollection_, EEdigiToken_, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), eventHeaderCollection_, eventHeaderToken_, fAnalysisName, fANY_RUN, fBuildEventDistribBad, fBuildEventDistribGood, fChozenGainNumber, fChozenRunTypeNumber, fCurrentEventNumber, fDateFirst, fDateLast, fDeeDS5Memo1, fDeeDS5Memo2, fDeeIndexBegin, fDeeIndexStop, fDynBaseLineSub, fESFromFedTcc, fEvtNumber, fEvtNumberMemo, fFedDigiOK, fFedId, fFedIdCounter, fFedNbOfTreatedEvents, fFedStatus, fFedStatusOrder, fFedTcc, fFirstReqEvent, alignBH_cfg::fixed, fLastReqEvent, fMaxFedIdCounter, fMaxFedUnitCounter, fMaxMgpaGainCounter, fMaxRunTypeCounter, fMaxTreatedStexCounter, fMemoCutOK, fMemoDateFirstEvent, fMgpaGainCounter, fMgpaGainNumber, fMyCnaEBSM, fMyCnaEEDee, fMyEBEcal, fMyEBNumbering, fMyEcnaEBObjectManager, fMyEcnaEEObjectManager, fMyEEEcal, fMyEENumbering, fNbOfSamples, fNbOfSelectedEvents, fNbOfTreatedFedsInDee, fNbOfTreatedFedsInStex, fNbOfTreatedStexs, fRecNumber, fReqNbOfEvts, fRunNumber, fRunTypeCounter, fRunTypeNumber, fSMFromFedTcc, fSMIndexBegin, fSMIndexStop, fStexDigiOK, fStexIndexBegin, fStexIndexStop, fStexName, fStexNbOfTreatedEvents, fStexNumber, fStexStatus, fTimeFirst, fTimeLast, fTreatedFedOrder, fTTBELL, TEcnaNumbering::Get0TowEchaFrom1SMCrys(), TEcnaNumbering::Get1DeeSCEcnaFrom1DeeCrys(), TEcnaNumbering::Get1SCEchaFrom1DeeCrys(), TEcnaNumbering::Get1SMTowFrom1SMCrys(), TEcnaNumbering::GetDeeDirViewedFromIP(), l1tPhase2CaloJetEmulator_cfi::iEta, iEvent, EBDetId::ism(), edm::HandleBase::isValid(), EEDetId::ix(), TEcnaParEcal::MaxCrysIYInSC(), TEcnaParEcal::MaxCrysPhiInTow(), TEcnaParEcal::MaxDeeInEE(), TEcnaParEcal::MaxDSInEE(), TEcnaParEcal::MaxSampADC(), TEcnaParEcal::MaxSCIYInDee(), TEcnaParEcal::MaxSMInEB(), TEcnaParEcal::MaxTowPhiInSM(), nChannels_, edm::Handle< T >::product(), runtype(), b2gDoubleLeptonHLTEventValidation_cfi::sDir, edm::DataFrameContainer::size(), and RecoSummaryTask_cfi::Time.

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