#include <CalibCalorimetry/EcalCorrelatedNoiseAnalysisModules/src/EcnaAnalyzer.cc>

Inheritance diagram for EcnaAnalyzer:

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

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::string	digiProducer_

std::string	EBdigiCollection_

std::string	EEdigiCollection_

std::string	eventHeaderCollection_

std::string	eventHeaderProducer_

TString	fAnalysisName

Int_t	fANY_RUN

Int_t *	fBuildEventDistribBad

Int_t *	fBuildEventDistribGood

TString	fCfgAnalyzerParametersFileName

TString	fCfgAnalyzerParametersFilePath

Int_t	fChozenGainNumber

Int_t	fChozenRunTypeNumber

Int_t *	fCounterQuad

Int_t	fCurrentEventNumber

TString *	fDateFirst

TString *	fDateLast

Int_t	fDeeDS5Memo1

Int_t	fDeeDS5Memo2

Int_t *	fDeeFromFedTcc

Int_t	fDeeIndexBegin

Int_t	fDeeIndexStop

TString *	fDeeNumberString

TString	fDynBaseLineSub

Int_t *	fESFromFedTcc

Int_t	fEvtNumber

Int_t	fEvtNumberMemo

std::ifstream	fFcin_f

Int_t *	fFedDigiOK

Int_t	fFedId

Int_t *	fFedIdCounter

Int_t *	fFedNbOfTreatedEvents

Int_t *	fFedStatus

Int_t *	fFedStatusOrder

Int_t	fFedTcc

Int_t	fFirstReqEvent

Int_t	fgMaxCar

Int_t	fLASER_STD

Int_t	fLastReqEvent

Int_t	fMaxCounterQuad

Int_t	fMaxFedIdCounter

Int_t	fMaxFedUnitCounter

Int_t	fMaxMgpaGainCounter

Int_t	fMaxRunTypeCounter

Int_t	fMaxTreatedStexCounter

Int_t	fMemoCutOK

Int_t *	fMemoDateFirstEvent

Int_t *	fMgpaGainCounter

Int_t	fMgpaGainNumber

TEcnaRun **	fMyCnaEBSM

TEcnaRun **	fMyCnaEEDee

TEcnaParEcal *	fMyEBEcal

TEcnaNumbering *	fMyEBNumbering

TEcnaObject *	fMyEcnaEBObjectManager

TEcnaObject *	fMyEcnaEEObjectManager

TEcnaParEcal *	fMyEEEcal

TEcnaNumbering *	fMyEENumbering

Int_t	fNbOfSamples

Int_t	fNbOfSelectedEvents

Int_t *	fNbOfTreatedFedsInDee

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

Int_t *	fRunTypeCounter

Int_t	fRunTypeNumber

Int_t *	fSMFromFedTcc

Int_t	fSMIndexBegin

Int_t	fSMIndexStop

Int_t *	fStexDigiOK

Int_t	fStexIndexBegin

Int_t	fStexIndexStop

TString	fStexName

Int_t *	fStexNbOfTreatedEvents

Int_t	fStexNumber

Int_t *	fStexStatus

time_t *	fTimeFirst

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_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename 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)

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

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

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

Definition at line 132 of file EcnaAnalyzer.h.

Member Enumeration Documentation

anonymous enum

Enumerator
kChannels
kGains
kFirstGainId

Definition at line 136 of file EcnaAnalyzer.h.

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

EcnaAnalyzer::kGains

Definition: EcnaAnalyzer.h:136

EcnaAnalyzer::kFirstGainId

Definition: EcnaAnalyzer.h:136

EcnaAnalyzer::kChannels

Definition: EcnaAnalyzer.h:136

Constructor & Destructor Documentation

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

explicit

Definition at line 44 of file EcnaAnalyzer.cc.

References AnalysisOutcome(), gather_cfg::cout, digiProducer_, EBdigiCollection_, EEdigiCollection_, eventHeaderCollection_, eventHeaderProducer_, 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, fLASER_STD, fLastReqEvent, fMaxFedIdCounter, fMaxFedUnitCounter, fMaxMgpaGainCounter, fMaxRunTypeCounter, fMaxTreatedStexCounter, fMemoCutOK, fMemoDateFirstEvent, fMgpaGainCounter, fMgpaGainNumber, fMyCnaEBSM, fMyCnaEEDee, fMyEBEcal, fMyEBNumbering, fMyEcnaEBObjectManager, fMyEcnaEEObjectManager, fMyEEEcal, fMyEENumbering, 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(), TEcnaParPaths::GetPathForResultsRootFiles(), mps_fire::i, TEcnaParEcal::MaxDeeInEE(), TEcnaParEcal::MaxDSInEE(), TEcnaParEcal::MaxSMInEB(), TEcnaParPaths::ResultsRootFilePath(), sAnalysisName_, sFirstReqEvent_, sLastReqEvent_, sNbOfSamples_, sReqNbOfEvts_, sStexName_, sStexNumber_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                       : 
   verbosity_(pSet.getUntrackedParameter("verbosity", 1U)),
   nChannels_(0), iEvent_(0)
 {
   //now do what ever initialization is needed
 
   using namespace edm;
   using namespace std;
 
   fMyEcnaEBObjectManager = new TEcnaObject();
   fMyEcnaEEObjectManager = new TEcnaObject();
 
   TEcnaParPaths* myPathEB = new TEcnaParPaths(fMyEcnaEBObjectManager);
   TEcnaParPaths* myPathEE = new TEcnaParPaths(fMyEcnaEEObjectManager);
 
   std::cout << "*EcnaAnalyzer-constructor> Check path for resultsq Root files." << std::endl;
 
   if( myPathEB->GetPathForResultsRootFiles() == kFALSE )
     {
       std::cout << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found." << std::endl;
       kill(getpid(),SIGUSR2);
     }
   else
     {
       std::cout << "*EcnaAnalyzer-constructor> Path for result files found = " << myPathEB->ResultsRootFilePath() << std::endl;
     }
 
   if( myPathEE->GetPathForResultsRootFiles() == kFALSE )
     {
       std::cout << "*EcnaAnalyzer-constructor> *** ERROR *** Path for result files not found." << std::endl;
       kill(getpid(),SIGUSR2);
     }
   else
     {
       std::cout << "*EcnaAnalyzer-constructor> Path for result files found = " << myPathEE->ResultsRootFilePath() << std::endl;
     }
 
 
   std::cout << "*EcnaAnalyzer-constructor> Parameter initialization." << std::endl;
 
   fgMaxCar = (Int_t)512;
   fTTBELL = '\007';
   fOutcomeError = kFALSE;
 
   fMyEBEcal      = new TEcnaParEcal(fMyEcnaEBObjectManager, "EB");
   fMyEBNumbering = new TEcnaNumbering(fMyEcnaEBObjectManager, "EB");
 
   fMyEEEcal      = new TEcnaParEcal(fMyEcnaEEObjectManager, "EE");
   fMyEENumbering = new TEcnaNumbering(fMyEcnaEEObjectManager, "EE");
 
   //========================================================================================== 
   //.................................. Get parameter values from python file
   eventHeaderProducer_   = pSet.getParameter<std::string>("eventHeaderProducer");
   digiProducer_          = pSet.getParameter<std::string>("digiProducer");
 
   eventHeaderCollection_ = pSet.getParameter<std::string>("eventHeaderCollection");
   EBdigiCollection_      = pSet.getParameter<std::string>("EBdigiCollection");
   EEdigiCollection_      = pSet.getParameter<std::string>("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;
   //===========================================================================================
 
   fRunTypeCounter = nullptr;
   fMaxRunTypeCounter = 26;
   fRunTypeCounter = new Int_t[fMaxRunTypeCounter];
   for(Int_t i=0; i<fMaxRunTypeCounter; i++){fRunTypeCounter[i] = 0;}
 
   fMgpaGainCounter = nullptr;
   fMaxMgpaGainCounter = 4;  // Because chozen gain = 0,1,2,3
   fMgpaGainCounter = new Int_t[fMaxMgpaGainCounter];
   for(Int_t i=0; i<fMaxMgpaGainCounter; i++){fMgpaGainCounter[i] = 0;}
 
   fFedIdCounter = nullptr;
   fMaxFedIdCounter = 54;
   fFedIdCounter = new Int_t[fMaxFedIdCounter];
   for(Int_t i=0; i<fMaxFedIdCounter; i++){fFedIdCounter[i] = 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 = new Int_t[fMaxFedUnitCounter];
   for(Int_t i=0; i<fMaxFedUnitCounter; i++){fFedDigiOK[i] = 0;}
 
   fFedNbOfTreatedEvents = new Int_t[fMaxFedUnitCounter];
   for(Int_t i=0; i<fMaxFedUnitCounter; i++){fFedNbOfTreatedEvents[i] = 0;}
 
   fFedStatus = new Int_t[fMaxFedUnitCounter];
   for(Int_t i=0; i<fMaxFedUnitCounter; i++){fFedStatus[i] = 0;}
 
   fFedStatusOrder = new Int_t[fMaxFedUnitCounter];
   for(Int_t i=0; i<fMaxFedUnitCounter; i++){fFedStatusOrder[i] = 0;}
 
   fDeeNumberString = new TString[fMaxFedUnitCounter];
   for(Int_t i=0; i<fMaxFedUnitCounter; i++){fDeeNumberString[i] = "SM";}
 
   if( fStexName == "Dee" )
     { 
       fDeeNumberString[ 0] = "Sector1 Dee4";
       fDeeNumberString[ 1] = "Sector2 Dee4";
       fDeeNumberString[ 2] = "Sector3 Dee4";
       fDeeNumberString[ 3] = "Sector4 Dee4";
       fDeeNumberString[ 4] = "Sector5 Dee4-Dee3";
       fDeeNumberString[ 5] = "Sector6 Dee3";
       fDeeNumberString[ 6] = "Sector7 Dee3";
       fDeeNumberString[ 7] = "Sector8 Dee3";
       fDeeNumberString[ 8] = "Sector9 Dee3";
       fDeeNumberString[ 9] = "Sector1 Dee1";
       fDeeNumberString[10] = "Sector2 Dee1";
       fDeeNumberString[11] = "Sector3 Dee1";
       fDeeNumberString[12] = "Sector4 Dee1";
       fDeeNumberString[13] = "Sector5 Dee1-Dee2";
       fDeeNumberString[14] = "Sector6 Dee2";
       fDeeNumberString[15] = "Sector7 Dee2";
       fDeeNumberString[16] = "Sector8 Dee2";
       fDeeNumberString[17] = "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)
   
   Int_t MaxSMAndDS = fMyEBEcal->MaxSMInEB()+fMyEEEcal->MaxDSInEE();
 
   fSMFromFedTcc = new Int_t[MaxSMAndDS];
   fESFromFedTcc = new Int_t[MaxSMAndDS];
   for(Int_t nFedTcc= 1; nFedTcc<=MaxSMAndDS; nFedTcc++){fESFromFedTcc[nFedTcc-1] = -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 = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fStexNbOfTreatedEvents[i] = 0;}
 
   fTimeFirst = new time_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fTimeFirst[i] = 0;}
   fTimeLast = new time_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fTimeLast[i] = 0;}
 
   fMemoDateFirstEvent = new Int_t[fMaxTreatedStexCounter];;
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fMemoDateFirstEvent[i] = 0;}
 
   Int_t MaxCar = fgMaxCar; 
   fDateFirst = new TString[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++)
     {
       fDateFirst[i].Resize(MaxCar);
       fDateFirst[i] = "*1st event date not found*";
     }
 
   MaxCar = fgMaxCar; 
   fDateLast = new TString[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++)
     {
       fDateLast[i].Resize(MaxCar);
       fDateLast[i] = "*last event date not found*";
     }
 
   fStexStatus = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fStexStatus[i] = 0;}
 
   fStexDigiOK = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fStexDigiOK[i] = 0;}
 
   fNbOfTreatedFedsInDee = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fNbOfTreatedFedsInDee[i] = 0;}
 
   fNbOfTreatedFedsInStex = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fNbOfTreatedFedsInStex[i] = 0;}
 
 
   //.......................... counters of events for GetSampleAdcValues
   fBuildEventDistribBad = nullptr;
   fBuildEventDistribBad = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fBuildEventDistribBad[i] = 0;}
 
   fBuildEventDistribGood = nullptr;
   fBuildEventDistribGood = new Int_t[fMaxTreatedStexCounter];
   for(Int_t i=0; i<fMaxTreatedStexCounter; i++){fBuildEventDistribGood[i] = 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;
 
   fMyCnaEBSM  = nullptr;
   fMyCnaEEDee = nullptr;
 
   fRunTypeNumber  = -1;
   fMgpaGainNumber = -1;
 
   fFedId  = -1;
   fFedTcc = -1;
 
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fAnalysisName        = " << fAnalysisName << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fNbOfSamples         = " << fNbOfSamples << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fFirstReqEvent       = " << fFirstReqEvent << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fLastReqEvent        = " << fLastReqEvent << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fReqNbOfEvts         = " << fReqNbOfEvts << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fStexName            = " << fStexName << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fStexNumber          = " << fStexNumber << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fChozenRunTypeNumber = " << fChozenRunTypeNumber << std::endl;
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> fChozenGainNumber    = " << fChozenGainNumber  << std::endl << std::endl;
 
   std::cout << "*EcnaAnalyzer::EcnaAnalyzer-constructor> Init done. " << std::endl;
 }

EcnaAnalyzer::~EcnaAnalyzer ( )

override

Definition at line 417 of file EcnaAnalyzer.cc.

References CheckMsg(), gather_cfg::cout, TauDecayModes::dec, fBuildEventDistribBad, fBuildEventDistribGood, fDateFirst, fDateLast, fDeeIndexBegin, fDeeIndexStop, fFedId, fFedIdCounter, alignBH_cfg::fixed, fMaxFedIdCounter, fMaxMgpaGainCounter, fMaxRunTypeCounter, fMgpaGainCounter, fMgpaGainNumber, fMyCnaEBSM, fMyCnaEEDee, fMyEBEcal, fMyEBNumbering, fMyEEEcal, fMyEENumbering, fOutcomeError, fRunTypeCounter, fRunTypeNumber, fSMIndexBegin, fSMIndexStop, fStexIndexBegin, fStexIndexStop, fStexName, fStexNumber, fTimeFirst, fTimeLast, fTTBELL, gainvalue(), TEcnaRun::GetReadyToCompute(), mps_fire::i, n0, runtype(), TEcnaRun::SampleValues(), TEcnaRun::StartStopDate(), and TEcnaRun::StartStopTime().

 {
   // do anything here that needs to be done at destruction time
   // (e.g. close files, deallocate resources etc.)
 
   using namespace std;
   //..................................... format numerical values
   std::cout << std::setiosflags(std::ios::showpoint | std::ios::uppercase);
   std::cout << std::setprecision(3) << std::setw(6);
   cout.setf(std::ios::dec, std::ios::basefield);
   cout.setf(std::ios::fixed, std::ios::floatfield);
   cout.setf(std::ios::left, std::ios::adjustfield);
   cout.setf(std::ios::right, std::ios::adjustfield);
 
   std::cout << "EcnaAnalyzer::~EcnaAnalyzer()> destructor is going to be executed." << std::endl;
 
   if( fOutcomeError == kTRUE )return;
 
   //-------------------------------------------------------------------------------
 
   //....................................................... EB (SM)
   if( fMyCnaEBSM == nullptr && fStexName == "SM" )
     {
       std::cout << std::endl << "!EcnaAnalyzer-destructor> **** ERROR **** fMyCnaEBSM = " << fMyCnaEBSM
         << ". !===> ECNA HAS NOT BEEN INITIALIZED." << std::endl
         << "  Last event run type = " << runtype(fRunTypeNumber)
         << ", fRunTypeNumber = " << fRunTypeNumber
         << ", last event Mgpa gain = " << gainvalue(fMgpaGainNumber)
         << ", fMgpaGainNumber = " << fMgpaGainNumber
         << ", last event fFedId(+601) = " << fFedId+601 << std::endl << std::endl;
     }
   else
     {
       for(Int_t iSM = fSMIndexBegin; iSM < fSMIndexStop; iSM++)
     {
       if( fMyCnaEBSM[iSM] != 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 )
         {
           std::cout << "!EcnaAnalyzer-destructor> PROBLEM with write ROOT file for SM" << iSM+1
                 << fTTBELL << std::endl;
         }
         }
       else
         {
           std::cout << "*EcnaAnalyzer-destructor> Calculations and writing on file already done for SM "
             << iSM+1 << std::endl;
         }
     }
       delete fMyCnaEBSM;
     }
   //....................................................... EE (Dee)
 
   if( fMyCnaEEDee == nullptr && fStexName == "Dee" )
     {
       std::cout << std::endl << "!EcnaAnalyzer-destructor> **** ERROR **** fMyCnaEEDee = " << fMyCnaEEDee
         << ". !===> ECNA HAS NOT BEEN INITIALIZED." << std::endl
         << "  Last event run type = " << runtype(fRunTypeNumber)
         << ", fRunTypeNumber = " << fRunTypeNumber
         << ", last event Mgpa gain = " << gainvalue(fMgpaGainNumber)
         << ", fMgpaGainNumber = " << fMgpaGainNumber
         << ", last event fFedId(+601) = " << fFedId+601 << std::endl << std::endl;
     }
   else
     {
       for(Int_t iDee = fDeeIndexBegin; iDee < fDeeIndexStop; iDee++)
     {
       if( fMyCnaEEDee[iDee] != 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 )
         {
           std::cout << "!EcnaAnalyzer-destructor> PROBLEM with write ROOT file for Dee" << iDee+1
                 << fTTBELL << std::endl;
         }
         }
       else
         {
           std::cout << "*EcnaAnalyzer-destructor> Calculations and writing on file already done for Dee "
             << iDee+1 << std::endl;
         }
     }
       delete fMyCnaEEDee;
     }
   std::cout <<endl;
 
   //-----------------------------------------------------------------------------------
 
   std::cout << "*EcnaAnalyzer-destructor> Status of events returned by GetSampleAdcValues(): "
         << std::endl;
 
   for(Int_t i0Stex=fStexIndexBegin; i0Stex<fStexIndexStop; i0Stex++ )
     {
       std::cout << fStexName << i0Stex+1 << "> Status OK: " << fBuildEventDistribGood[i0Stex]
         << " / ERROR(S): " << fBuildEventDistribBad[i0Stex];
       if( fBuildEventDistribBad[i0Stex] > 0 )
     {std::cout << " <=== SHOULD BE EQUAL TO ZERO ! " << fTTBELL;}
       std::cout << std::endl;
     }
 
   std::cout << std::endl<< "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - " <<  std::endl;
   
   std::cout << "*EcnaAnalyzer-destructor> Run types seen in event headers before selection:" << std::endl;
 
   for(Int_t i=0; i<fMaxRunTypeCounter; i++)
     {
       std::cout << " => " << std::setw(10) << fRunTypeCounter[i]
         << " event header(s) with run type " << runtype(i) << std::endl; 
     }
 
   std::cout << std::endl<< "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - " <<  std::endl;
 
   std::cout << "*EcnaAnalyzer-destructor> Mgpa gains seen in event headers before selection:" << std::endl;
 
   for(Int_t i=0; i<fMaxMgpaGainCounter; i++)
     {
       std::cout << " => " << std::setw(10) << fMgpaGainCounter[i]
         << " event header(s) with gain " << gainvalue(i) << std::endl; 
     }
 
   std::cout << std::endl<< "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - " <<  std::endl;
 
   std::cout << "*EcnaAnalyzer-destructor> Numbers of selected events for each FED:" << std::endl;
 
   for(Int_t i=0; i<fMaxFedIdCounter; i++)
     {
       std::cout << " => FedId " << i+601 << ": "
             << std::setw(10) << fFedIdCounter[i] << " events" << std::endl;
     }
 
   std::cout << std::endl<< "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - " <<  std::endl;
 
   if( fStexNumber == 0 )
     {
       // std::cout << "*EcnaAnalyzer-destructor> fDateFirst = " << fDateFirst[0] << std::endl
       //           << "                          fDateLast  = " << fDateLast[fMaxTreatedStexCounter-1] << std::endl << std::endl;
     }
   if( fStexNumber > 0 )
     {
       std::cout << "*EcnaAnalyzer-destructor> fDateFirst = " << fDateFirst[fStexNumber-1] << std::endl
         << "                          fDateLast  = " << fDateLast[fStexNumber-1] << std::endl << std::endl;
     }
 
   std::cout << std::endl<< "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - " <<  std::endl;
 
   Int_t n0 =0; CheckMsg(n0);
 
   delete fMyEBNumbering;
   delete fMyEENumbering;
 
   delete fMyEBEcal;
   delete fMyEEEcal;
 
   std::cout << "*EcnaAnalyzer-destructor> End of execution." << std::endl;
 }

Member Function Documentation

Bool_t EcnaAnalyzer::AnalysisOutcome ( const TString & s_opt )

Definition at line 1311 of file EcnaAnalyzer.cc.

References CheckMsg(), gather_cfg::cout, 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) )
       )
     {
       std::cout << std::endl
             << "**************************** ANALYSIS REPORT > OK **************************************" << std::endl
             << "*EcnaAnalyzer::AnalysisOutcome(...)> The maximum requested number of events and the maximum" << std::endl
             << "                                     number of treated " << fStexName
             << "'s have been reached." << std::endl
             << "                                     Analysis successfully ended from EcnaAnalyzer " << std::endl
             << "                                     by SIGNAL: kill(getpid(),SIGUSR2)." << std::endl
             << "                                     Number of selected events   = " << fNbOfSelectedEvents << std::endl
             << "                                     Last requested event number = " << fLastReqEvent << std::endl
             << "                                     Current event number        = " << fCurrentEventNumber << std::endl;
 
       Int_t n0 = 0; CheckMsg(n0);
 
       std::cout << "****************************************************************************************" << std::endl
             << std::endl;
 
       result = kTRUE;
       kill(getpid(),SIGUSR2);
     }
 
       if( fLastReqEvent >= fFirstReqEvent && fCurrentEventNumber > fLastReqEvent &&
       ! ( (fStexNumber > 0 && fNbOfTreatedStexs == 1) ||
           (fStexNumber == 0 && fNbOfTreatedStexs == MaxNbOfStex) ) )
     {
       std::cout << std::endl
             << "**************************** ANALYSIS REPORT >>> *** WARNING *** WARNING *** WARNING ***" << std::endl
             << "*EcnaAnalyzer::AnalysisOutcome(...)> Last event reached before completion of analysis." << std::endl
             << "                                     Analysis ended from EcnaAnalyzer " << std::endl
             << "                                     by SIGNAL: kill(getpid(),SIGUSR2)." << std::endl
             << "                                     Number of selected events   = " << fNbOfSelectedEvents << std::endl
             << "                                     Last requested event number = " << fLastReqEvent << std::endl
             << "                                     Current event number        = " << fCurrentEventNumber << std::endl;
 
       Int_t n0 = 0; CheckMsg(n0);
 
       std::cout << "****************************************************************************************" << std::endl
             << std::endl;
       
       result = kTRUE;
       kill(getpid(),SIGUSR2);
     }
     }
   else
     {
       if( s_opt == "ERR_FNEG" )
     {
       std::cout << std::endl
             << "**************************** ANALYSIS REPORT >>> **** ERROR **** ERROR **** ERROR ******" << std::endl
             << "*EcnaAnalyzer::AnalysisOutcome(...)> First event number = " << fFirstReqEvent
             << ". Should be strictly potitive." << std::endl
             << "                             Analysis ended from EcnaAnalyzer " << std::endl
             << "                             by SIGNAL: kill(getpid(),SIGUSR2)." << std::endl;
 
       std::cout << "****************************************************************************************" << std::endl
             << std::endl;
 
       result = kTRUE;
       kill(getpid(),SIGUSR2);
     }
       if( s_opt == "ERR_LREQ" )
     {
       std::cout << std::endl
             << "**************************** ANALYSIS REPORT >>> **** ERROR **** ERROR **** ERROR ******" << std::endl
             << "*EcnaAnalyzer::analyze(...)> Requested number of events = " << fReqNbOfEvts << "." << std::endl
             << "                             Too large compared to the event range: "
             << fFirstReqEvent << " - " << fLastReqEvent << std::endl
             << "                             Analysis ended from EcnaAnalyzer " << std::endl
             << "                             by SIGNAL: kill(getpid(),SIGUSR2)." << std::endl;
 
       std::cout << "****************************************************************************************" << std::endl
             << std::endl;
 
       result = kTRUE;
       kill(getpid(),SIGUSR2);
     }
     }
   return result;
 } // end of EcnaAnalyzer::AnalysisOutcome(const Int_t& n_option)

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

override

Definition at line 597 of file EcnaAnalyzer.cc.

References ecalMGPA::adc(), EcalMGPASample::adc(), AnalysisOutcome(), edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), MessageLogger_cfi::cerr, CheckMsg(), gather_cfg::cout, TauDecayModes::dec, digiProducer_, EBdigiCollection_, EEdigiCollection_, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), edm::EventID::event(), eventHeaderCollection_, eventHeaderProducer_, cppFunctionSkipper::exception, 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(), edm::Event::getByLabel(), TEcnaNumbering::GetDeeDirViewedFromIP(), TEcnaRun::GetReadyToCompute(), TEcnaRun::GetReadyToReadData(), edm::EventBase::id(), EBDetId::ism(), EEDetId::ix(), TEcnaParEcal::MaxCrysIYInSC(), TEcnaParEcal::MaxCrysPhiInTow(), TEcnaParEcal::MaxDeeInEE(), TEcnaParEcal::MaxDSInEE(), TEcnaParEcal::MaxSCIYInDee(), TEcnaParEcal::MaxSMInEB(), TEcnaParEcal::MaxTowPhiInSM(), nChannels_, edm::Handle< T >::product(), edm::EventID::run(), runtype(), EcalDataFrame::sample(), TEcnaRun::SampleValues(), edm::DataFrameContainer::size(), TEcnaRun::StartStopDate(), TEcnaRun::StartStopTime(), edm::EventBase::time(), and edm::Timestamp::value().

 {
   using namespace std;
   //..................................... format numerical values
   std::cout << std::setiosflags(std::ios::showpoint | std::ios::uppercase);
   std::cout << std::setprecision(3) << std::setw(6);
   cout.setf(std::ios::dec, std::ios::basefield);
   cout.setf(std::ios::fixed, std::ios::floatfield);
   cout.setf(std::ios::left, std::ios::adjustfield);
   cout.setf(std::ios::right, std::ios::adjustfield);
 
   using namespace edm;
 
   fRecNumber++; 
 
   Int_t iFreq = (fLastReqEvent - fFirstReqEvent + 1)/5;
   if( iFreq <= 0 ){iFreq = 10000;}
 
   Int_t MaxSMAndDS = fMyEBEcal->MaxSMInEB()+fMyEEEcal->MaxDSInEE();
 
   //********************************************* EVENT TREATMENT ********************************
   Handle<EcalRawDataCollection> pEventHeader;
   const EcalRawDataCollection* myEventHeader = nullptr;
   try{
     iEvent.getByLabel(eventHeaderProducer_, eventHeaderCollection_, pEventHeader);
     myEventHeader = pEventHeader.product();
   }catch (std::exception& ex ){
     std::cerr << "Error! can't get the product " << eventHeaderCollection_.c_str() << std::endl;
   }
   //........... 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 == nullptr && fStexName == "SM" )
     {
       fMyCnaEBSM = new TEcnaRun*[fMyEBEcal->MaxSMInEB()];
       for(Int_t i0SM = 0; i0SM < fMyEBEcal->MaxSMInEB(); i0SM++)
     {fMyCnaEBSM[i0SM] = nullptr;}
     }
   //.................................................................. EE (Dee)
   if( fMyCnaEEDee == nullptr && fStexName == "Dee" )
     {
       fMyCnaEEDee = new TEcnaRun*[fMyEEEcal->MaxDeeInEE()];
       for(Int_t iDee = 0; iDee < fMyEEEcal->MaxDeeInEE(); iDee++)
     {fMyCnaEEDee[iDee] = nullptr;}
     }
   
   //============================ 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
       Handle<EBDigiCollection> pdigisEB;
       const EBDigiCollection* digisEB = nullptr;
       try{
         iEvent.getByLabel(digiProducer_, EBdigiCollection_, pdigisEB); 
         digisEB = pdigisEB.product();
       }catch (std::exception& ex ){
         std::cerr << "Error! can't get the product " << EBdigiCollection_.c_str() << std::endl;
       }
       
       // 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] == nullptr && fStexStatus[i0SM] != 2 )
             {
               //=============================== Init Ecna EB ===============================
               fMyCnaEBSM[i0SM] = new TEcnaRun(fMyEcnaEBObjectManager, "EB", fNbOfSamples);
               fMyCnaEBSM[i0SM]->GetReadyToReadData(fAnalysisName,  fRunNumber,
                                    fFirstReqEvent, fLastReqEvent, fReqNbOfEvts,
                                    i0SM+1,         fRunTypeNumber);
               
               std::cout << "*EcnaAnalyzer::analyze(...)> ********* INIT ECNA EB ********* " << std::endl
                     << "                                   fAnalysisName = " << fAnalysisName << std::endl
                     << "                                      fRunNumber = " << fRunNumber << std::endl
                     << "                                  fFirstReqEvent = " << fFirstReqEvent << std::endl
                     << "                                   fLastReqEvent = " << fLastReqEvent << std::endl
                     << "                                    fReqNbOfEvts = " << fReqNbOfEvts << std::endl
                     << "                                              SM = " << i0SM+1 << std::endl
                     << "                                        run type = " << runtype(fRunTypeNumber)
                     << std::endl;
               //============================================================================
             }
               
               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;
                   std::cout << "*----> beginning of analysis for " << fStexName << i0SM+1
                         << ". First analyzed event date : " << astime << std::endl;
                   //      << " t_current_ev_time = " << t_current_ev_time  << std::endl
                   //      << " i_current_ev_time = " << i_current_ev_time  << std::endl
                   //      << " p_current_ev_time = " << p_current_ev_time  << std::endl
                 }
 
                   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
                     {
                       std::cout << "EcnaAnalyzer::analyze(...)> NbOfSamplesFromDigis out of bounds = "
                         << NbOfSamplesFromDigis << std::endl;
                     }
                 } // 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
       Handle<EEDigiCollection> pdigisEE;
       const EEDigiCollection* digisEE = nullptr;
       try{
         iEvent.getByLabel(digiProducer_, EEdigiCollection_, pdigisEE); 
         digisEE = pdigisEE.product();
       }catch (std::exception& ex ){
         std::cerr << "Error! can't get the product " << EEdigiCollection_.c_str() << std::endl;
       }
       
       // 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] == nullptr  && fStexStatus[i0Dee] != 2 )
             {
               //=============================== Init Ecna EE ===============================
               fMyCnaEEDee[i0Dee] = new TEcnaRun(fMyEcnaEEObjectManager, "EE", fNbOfSamples);
               fMyCnaEEDee[i0Dee]->GetReadyToReadData(fAnalysisName,  fRunNumber,
                                  fFirstReqEvent, fLastReqEvent, fReqNbOfEvts,
                                  i0Dee+1,         fRunTypeNumber);
               
               std::cout << "*EcnaAnalyzer::analyze(...)> ********* INIT ECNA EE ********* " << std::endl
                     << "                                   fAnalysisName = " << fAnalysisName << std::endl
                     << "                                      fRunNumber = " << fRunNumber << std::endl
                     << "                                  fFirstReqEvent = " << fFirstReqEvent << std::endl
                     << "                                   fLastReqEvent = " << fLastReqEvent << std::endl
                     << "                                    fReqNbOfEvts = " << fReqNbOfEvts << std::endl
                     << "                                             Dee = " << i0Dee+1 << std::endl
                     << "                                        run type = " << runtype(fRunTypeNumber)
                     << std::endl;
               //============================================================================
             }
 
               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;
                   std::cout << "----- beginning of analysis for " << fStexName << i0Dee+1 << "-------"  << std::endl
                     //<< " t_current_ev_time = " << t_current_ev_time  << std::endl
                     //<< " i_current_ev_time = " << i_current_ev_time  << std::endl
                     //<< " p_current_ev_time = " << p_current_ev_time  << std::endl
                         << " First event date  = " << astime << std::endl
                         << " Nb of selected evts = " << fNbOfSelectedEvents << std::endl
                         << "---------------------------------------------------------------"  << std::endl;
                   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
                     {
                       std::cout << "EcnaAnalyzer::analyze(...)> NbOfSamplesFromDigis out of bounds = "
                         << NbOfSamplesFromDigis << std::endl;
                     }
                 } // 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;}
       
       std::cout << "---------- End of analysis for " << fStexName << i0Stex+1 << " -----------" << std::endl;
       Int_t n3 = 3; CheckMsg(n3, i0Stex);
       // std::cout     << " t_current_ev_time = " << t_current_ev_time  << std::endl
       //<< " i_current_ev_time = " << i_current_ev_time  << std::endl
       //<< " p_current_ev_time = " << p_current_ev_time  << std::endl
       // std::cout      << " Last analyzed event date  = " << astime << std::endl;
       std::cout         << " Number of selected events = " << fNbOfSelectedEvents << std::endl;
       std::cout << std::endl << fNbOfTreatedStexs << " " << fStexName
             << "'s with " << fReqNbOfEvts << " events analyzed." << std::endl
             << "---------------------------------------------------------"  << std::endl;
       
       //================================= WRITE RESULTS FILE
       if( fStexName == "SM" )
         {
           if( fMyCnaEBSM[i0Stex] != 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 )
             {
               std::cout << "!EcnaAnalyzer::analyze> PROBLEM with write ROOT file for SM" << i0Stex+1
                 << fTTBELL << std::endl;
             }
         }
           // set pointer to zero in order to avoid recalculation and rewriting at the destructor level
           delete fMyCnaEBSM[i0Stex];  fMyCnaEBSM[i0Stex] = nullptr;
           std::cout << "!EcnaAnalyzer::analyze> Set memory free: delete done for SM " << i0Stex+1 << std::endl;
         }
 
       if( fStexName == "Dee" )
         {
           if( fMyCnaEEDee[i0Stex] != 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 )
             {
               std::cout << "!EcnaAnalyzer::analyze> PROBLEM with write ROOT file for Dee" << i0Stex+1
                 << fTTBELL << std::endl;
             }
         }
           // set pointer to zero in order to avoid recalculation and rewriting at the destructor level
           delete fMyCnaEEDee[i0Stex]; fMyCnaEEDee[i0Stex] = nullptr;
           std::cout << "!EcnaAnalyzer::analyze> Set memory free: delete done for Dee " << i0Stex+1 << std::endl;
         }
 
       fStexStatus[i0Stex] = 2;        // set fStexStatus[i0Stex] to 2 definitively
       std::cout << "*---------------------------------------------------------------------------- " << std::endl;
 
     } // end of if( fStexStatus[i0Stex] == 1 )
     } // end of for(Int_t i0Stex=fStexIndexBegin; i0Stex<fStexIndexStop; i0Stex++)
 }

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

Definition at line 1411 of file EcnaAnalyzer.cc.

References gather_cfg::cout, 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 )
     {std::cout << "---------------- CROSS-CHECK A ------------------ " << std::endl
            << "**************** CURRENT EVENT ****************** " << std::endl;}
   if( MsgNum == 2 )
     {std::cout << "---------------- CROSS-CHECK B ------------------ " << std::endl
            << "**** FIRST EVENT PASSING USER'S ANALYSIS CUT **** " << std::endl;}
   if( MsgNum == 3 )
     {std::cout << "---------------- CROSS-CHECK C ------------------ " << std::endl
            << "*** CURRENT VALUES BEFORE RESULT FILE WRITING *** " << std::endl;}
   if( MsgNum == 3 || MsgNum == 4 )
     {std::cout << "          fRecNumber = " << fRecNumber << std::endl
            << "          fEvtNumber = " << fEvtNumber << std::endl;}
   
   std::cout << " fCurrentEventNumber = " << fCurrentEventNumber << std::endl
         << " fNbOfSelectedEvents = " << fNbOfSelectedEvents << std::endl
         << "          fRunNumber = " << fRunNumber << std::endl
         << "     Chozen run type = " << runtype(fChozenRunTypeNumber) << std::endl
         << "            Run type = " << runtype(fRunTypeNumber) << std::endl
         << "             fFedTcc = " << fFedTcc << std::endl
         << "        fFedId(+601) = " << fFedId+601 << std::endl
         << "           fStexName = " << fStexName << std::endl
         << "         Chozen gain = " << gainvalue(fChozenGainNumber) << std::endl 
         << "           Mgpa Gain = " << gainvalue(fMgpaGainNumber) << std::endl << 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];}
           
           std::cout << fStexName << std::setw(3) << j0Stex+1 << ": "
             << std::setw(5) << nStexNbOfTreatedEvents << " events. "
             << fStexName << " status: " << fStexStatus[j0Stex];
           if( j0Stex == i0Stex ){std::cout << " (going to write file for this " << fStexName << ").";}
           std::cout << std::endl; 
         }
     }
 
       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];}
           
           std::cout << "Fed (ES) " << std::setw(3) << i0FedES+1 << ": "
             << std::setw(5) << nFedNbOfTreatedEvents << " events."
             << " Fed status: " << fFedStatus[i0FedES]
             << ", order: " << std::setw(3) << fFedStatusOrder[i0FedES]
             << " (" << fDeeNumberString[i0FedES] << ")" << std::endl;
         }
       
       for(Int_t j0Stex=fStexIndexBegin; j0Stex<fStexIndexStop; j0Stex++)
         {
           std::cout << fStexName << std::setw(3) << j0Stex+1 << ": "
             << std::setw(5) << fNbOfTreatedFedsInStex[j0Stex] << " analyzed Fed(s). "
             << fStexName << " status: " << fStexStatus[j0Stex];
           if( j0Stex == i0Stex ){std::cout << " (going to write file for this " << fStexName << ").";}
           std::cout << std::endl; 
         }
     }
 
       std::cout << "Number of " << fStexName << "'s with "
         << fReqNbOfEvts << " events analyzed: " << fNbOfTreatedStexs << std::endl;
     }
   
   if( MsgNum == 1 || MsgNum == 2 )
     {std::cout << "*---------------------------------------------------------------------------- " << std::endl;}
   if( MsgNum == 3 )
     {std::cout << "*............................................................................ " << std::endl;}
 
 } // end of EcnaAnalyzer::CheckMsg(const Int_t& MsgNum, const Int_t& i0Stex)

void EcnaAnalyzer::CheckMsg ( const Int_t & MsgNum )

Definition at line 1408 of file EcnaAnalyzer.cc.

References CheckMsg().

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

EcnaAnalyzer::CheckMsg

void CheckMsg(const Int_t &, const Int_t &)

Definition: EcnaAnalyzer.cc:1411

Int_t EcnaAnalyzer::gainvalue ( const Int_t & numgain )

Definition at line 1533 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;
 }

TString EcnaAnalyzer::runtype ( const Int_t & numtype )

Definition at line 1494 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

std::string EcnaAnalyzer::digiProducer_

private

Definition at line 158 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

std::string EcnaAnalyzer::EBdigiCollection_

private

Definition at line 160 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

std::string EcnaAnalyzer::EEdigiCollection_

private

Definition at line 161 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

std::string EcnaAnalyzer::eventHeaderCollection_

private

Definition at line 159 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

std::string EcnaAnalyzer::eventHeaderProducer_

private

Definition at line 157 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

TString EcnaAnalyzer::fAnalysisName

private

Definition at line 186 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fANY_RUN

private

Definition at line 236 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fBuildEventDistribBad

private

Definition at line 179 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fBuildEventDistribGood

private

Definition at line 180 of file EcnaAnalyzer.h.

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

TString EcnaAnalyzer::fCfgAnalyzerParametersFileName

private

Definition at line 183 of file EcnaAnalyzer.h.

TString EcnaAnalyzer::fCfgAnalyzerParametersFilePath

private

Definition at line 182 of file EcnaAnalyzer.h.

Int_t EcnaAnalyzer::fChozenGainNumber

private

Definition at line 187 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fChozenRunTypeNumber

private

Definition at line 188 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fCounterQuad

private

Definition at line 275 of file EcnaAnalyzer.h.

Int_t EcnaAnalyzer::fCurrentEventNumber

private

Definition at line 176 of file EcnaAnalyzer.h.

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

TString* EcnaAnalyzer::fDateFirst

private

Definition at line 245 of file EcnaAnalyzer.h.

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

TString* EcnaAnalyzer::fDateLast

private

Definition at line 246 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fDeeDS5Memo1

private

Definition at line 219 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fDeeDS5Memo2

private

Definition at line 220 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fDeeFromFedTcc

private

Definition at line 212 of file EcnaAnalyzer.h.

Int_t EcnaAnalyzer::fDeeIndexBegin

private

Definition at line 204 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fDeeIndexStop

private

Definition at line 205 of file EcnaAnalyzer.h.

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

TString* EcnaAnalyzer::fDeeNumberString

private

Definition at line 216 of file EcnaAnalyzer.h.

Referenced by CheckMsg(), and EcnaAnalyzer().

TString EcnaAnalyzer::fDynBaseLineSub

private

Definition at line 189 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fESFromFedTcc

private

Definition at line 211 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fEvtNumber

private

Definition at line 173 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fEvtNumberMemo

private

Definition at line 174 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

std::ifstream EcnaAnalyzer::fFcin_f

private

Definition at line 184 of file EcnaAnalyzer.h.

Int_t* EcnaAnalyzer::fFedDigiOK

private

Definition at line 227 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fFedId

private

Definition at line 215 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fFedIdCounter

private

Definition at line 272 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fFedNbOfTreatedEvents

private

Definition at line 228 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fFedStatus

private

Definition at line 226 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fFedStatusOrder

private

Definition at line 214 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fFedTcc

private

Definition at line 209 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fFirstReqEvent

private

Definition at line 194 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fgMaxCar

private

Definition at line 150 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fLASER_STD

private

Definition at line 239 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fLastReqEvent

private

Definition at line 195 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMaxCounterQuad

private

Definition at line 274 of file EcnaAnalyzer.h.

Int_t EcnaAnalyzer::fMaxFedIdCounter

private

Definition at line 271 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMaxFedUnitCounter

private

Definition at line 225 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMaxMgpaGainCounter

private

Definition at line 268 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMaxRunTypeCounter

private

Definition at line 265 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMaxTreatedStexCounter

private

Definition at line 218 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fMemoCutOK

private

Definition at line 230 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fMemoDateFirstEvent

private

Definition at line 248 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fMgpaGainCounter

private

Definition at line 269 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fMgpaGainNumber

private

Definition at line 200 of file EcnaAnalyzer.h.

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

TEcnaRun** EcnaAnalyzer::fMyCnaEBSM

private

Definition at line 253 of file EcnaAnalyzer.h.

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

TEcnaRun** EcnaAnalyzer::fMyCnaEEDee

private

Definition at line 254 of file EcnaAnalyzer.h.

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

TEcnaParEcal* EcnaAnalyzer::fMyEBEcal

private

Definition at line 257 of file EcnaAnalyzer.h.

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

TEcnaNumbering* EcnaAnalyzer::fMyEBNumbering

private

Definition at line 256 of file EcnaAnalyzer.h.

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

TEcnaObject* EcnaAnalyzer::fMyEcnaEBObjectManager

private

Definition at line 250 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

TEcnaObject* EcnaAnalyzer::fMyEcnaEEObjectManager

private

Definition at line 251 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

TEcnaParEcal* EcnaAnalyzer::fMyEEEcal

private

Definition at line 260 of file EcnaAnalyzer.h.

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

TEcnaNumbering* EcnaAnalyzer::fMyEENumbering

private

Definition at line 259 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fNbOfSamples

private

Definition at line 191 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fNbOfSelectedEvents

private

Definition at line 177 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fNbOfTreatedFedsInDee

private

Definition at line 232 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t* EcnaAnalyzer::fNbOfTreatedFedsInStex

private

Definition at line 233 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fNbOfTreatedStexs

private

Definition at line 231 of file EcnaAnalyzer.h.

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

Bool_t EcnaAnalyzer::fOutcomeError

private

Definition at line 171 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fPEDESTAL_GAP

private

Definition at line 238 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fPEDESTAL_STD

private

Definition at line 237 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fPEDSIM

private

Definition at line 241 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fPHYSICS_GLOBAL

private

Definition at line 240 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

Int_t EcnaAnalyzer::fRecNumber

private

Definition at line 175 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fReqNbOfEvts

private

Definition at line 199 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fRunNumber

private

Definition at line 192 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fRunTypeCounter

private

Definition at line 266 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fRunTypeNumber

private

Definition at line 193 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fSMFromFedTcc

private

Definition at line 210 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fSMIndexBegin

private

Definition at line 202 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fSMIndexStop

private

Definition at line 203 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fStexDigiOK

private

Definition at line 221 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

Int_t EcnaAnalyzer::fStexIndexBegin

private

Definition at line 206 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fStexIndexStop

private

Definition at line 207 of file EcnaAnalyzer.h.

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

TString EcnaAnalyzer::fStexName

private

Definition at line 196 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fStexNbOfTreatedEvents

private

Definition at line 222 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fStexNumber

private

Definition at line 197 of file EcnaAnalyzer.h.

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

Int_t* EcnaAnalyzer::fStexStatus

private

Definition at line 223 of file EcnaAnalyzer.h.

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

time_t* EcnaAnalyzer::fTimeFirst

private

Definition at line 243 of file EcnaAnalyzer.h.

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

time_t* EcnaAnalyzer::fTimeLast

private

Definition at line 244 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::fTreatedFedOrder

private

Definition at line 213 of file EcnaAnalyzer.h.

Referenced by analyze(), and EcnaAnalyzer().

TString EcnaAnalyzer::fTTBELL

private

Definition at line 151 of file EcnaAnalyzer.h.

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

Int_t EcnaAnalyzer::iEvent_

private

Definition at line 156 of file EcnaAnalyzer.h.

Int_t EcnaAnalyzer::nChannels_

private

Definition at line 155 of file EcnaAnalyzer.h.

Referenced by analyze().

TString EcnaAnalyzer::sAnalysisName_

private

Definition at line 163 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sFirstReqEvent_

private

Definition at line 165 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sLastReqEvent_

private

Definition at line 166 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sNbOfSamples_

private

Definition at line 164 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sReqNbOfEvts_

private

Definition at line 167 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sStexName_

private

Definition at line 168 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

TString EcnaAnalyzer::sStexNumber_

private

Definition at line 169 of file EcnaAnalyzer.h.

Referenced by EcnaAnalyzer().

unsigned int EcnaAnalyzer::verbosity_

private

Definition at line 154 of file EcnaAnalyzer.h.

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

cna_new/Correlated_Noise_Analysis/ECNA_cna_1.htm

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation