TEcnaRun.h

Go to the documentation of this file.

#ifndef   CL_TEcnaRun_H
#define   CL_TEcnaRun_H

// system include files
#include <memory>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <string>
#include <vector>
#include <ctime>
#include <cmath>
#include <cstdio>
#include "Riostream.h"

// ROOT include files
#include "TObject.h"
#include "TSystem.h"
#include "TString.h"
#include "TVectorD.h"

// user include files
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaObject.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParEcal.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaHeader.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParCout.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParPaths.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaRead.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaWrite.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaRootFile.h"


class TEcnaRun: public TObject {
  
 private:

  //............ attributes

  Int_t fgMaxCar;      // Max nb of caracters for char* 

  Int_t fCnaCommand,  fCnaError;

  Int_t fCnew;          // flags for dynamical allocation
  Int_t fCdelete;       

  TString fTTBELL;

  Int_t* fMiscDiag;                          // Counters for miscellaneous diagnostics
  Int_t  fNbOfMiscDiagCounters;
  Int_t  fMaxMsgIndexForMiscDiag;

  TEcnaObject    *fObjectManager; // for ECNA object management
  TEcnaHeader    *fFileHeader;    // header for result type file
  TEcnaParEcal   *fEcal;          // for access to the Ecal current subdetector parameters
  TEcnaNumbering *fEcalNumbering; // for access to the Ecal channel, Stin and Stex numbering
  TEcnaParCout   *fCnaParCout;    // for comment/error messages
  TEcnaParPaths  *fCnaParPaths;   // for file access
  TEcnaWrite     *fCnaWrite;      // for access to the results files

  //  TEcnaRootFile *gCnaRootFile;

  TString fFlagSubDet;
  TString fStexName, fStinName;

  Bool_t  fOpenRootFile;   // flag open ROOT file (open = kTRUE, close = kFALSE)
  Int_t   fReadyToReadData;

  TString  fRootFileName;
  TString  fRootFileNameShort;
  TString  fNewRootFileName;
  TString  fNewRootFileNameShort;

  Int_t   fSpecialStexStinNotIndexed;  // management of event distribution building
  Int_t   fStinIndexBuilt;
  Int_t   fBuildEvtNotSkipped;

  Int_t   fNbSampForFic;
  Int_t   fNbSampForCalc;

  Int_t   fNumberOfEvents;

  Int_t   fMemoReadNumberOfEventsforSamples;

  Double_t*** fT3d_AdcValues;   // 3D array[channel][sample][event] ADC values
  Double_t**  fT3d2_AdcValues; 
  Double_t*   fT3d1_AdcValues;
  Int_t*      fTagAdcEvt;

  Int_t**     fT2d_NbOfEvts; // 2D array[channel][sample] max nb of evts read for a given (channel,sample) 
  Int_t*      fT1d_NbOfEvts;
  Int_t*      fTagNbOfEvts;

  Int_t*      fT1d_StexStinFromIndex; // 1D array[Stin] Stin Number as a function of the index Stin
  Int_t*      fTagStinNumbers;

  Double_t**  fT2d_ev;        // 2D array[channel][sample] for expectation values
  Double_t*   fT1d_ev;
  Int_t*      fTagMSp;

  Double_t**  fT2d_sig;       // 2D array[channel][sample] for sigmass
  Double_t*   fT1d_sig;
  Int_t*      fTagSSp;
 
  Double_t*** fT3d_cov_ss;    // 3D array[channel][sample][sample] for (sample,sample) covariances
  Double_t**  fT3d2_cov_ss;
  Double_t*   fT3d1_cov_ss;
  Int_t*      fTagCovCss;

  Double_t*** fT3d_cor_ss;    // 3D array[channel][sample][sample] for (sample,sample) correlations
  Double_t**  fT3d2_cor_ss;
  Double_t*   fT3d1_cor_ss;
  Int_t*      fTagCorCss;

//...........................................................................
  Double_t*   fT1d_ev_ev;     // 1D array[channel] for expectation values of the expectation values of the samples
  Int_t*      fTagPed;        // (PEDESTAL)
  Double_t*   fT1d_av_mped;   // 1D array[Stin] for expectation values of the Pesdestals of the Stins
  Int_t*      fTagAvPed;      // (AVERAGED PEDESTAL)

  Double_t*   fT1d_evsamp_of_sigevt; // 1D array[channel] for expectation values of the sigmas of the samples
  Int_t*      fTagTno;               // (TOTAL NOISE)
  Double_t*   fT1d_av_totn;          // 1D array[Stin] for expectation values of the total noise
  Int_t*      fTagAvTno;             //(AVERAGED TOTAL NOISE)

  Double_t*   fT1d_ev_cor_ss;   // 1D array[channel] for expectation values of the cor(s,s)
  Int_t*      fTagMeanCorss;    // (MEAN COR(S,S))
  Double_t*   fT1d_av_ev_corss; // 1D array[Stin] for expectation values of the mean cor(s,s)
  Int_t*      fTagAvMeanCorss;  // (AVERAGED MEAN COR(S,S))

  Double_t*   fT1d_sigevt_of_evsamp;  // 1D array[channel] for sigmas of the expectation values of the samples
  Int_t*      fTagLfn;                // (LOW FREQUENCY NOISE)
  Double_t*   fT1d_av_lofn;           // 1D array[Stin]  the expectation values of the low frequency noise
  Int_t*      fTagAvLfn;              // (AVERAGED LOW FREQUENCY NOISE)

  Double_t*   fT1d_evevt_of_sigsamp; // 1D array[channel] for sigmas of the sigmas of the samples
  Int_t*      fTagHfn;               // (HIGH FREQUENCY NOISE)
  Double_t*   fT1d_av_hifn;          // 1D array[channel] for expectation values of the high frequency noise
  Int_t*      fTagAvHfn;             // (AVERAGED HIGH FREQUENCY NOISE)

  Double_t*   fT1d_sig_cor_ss;   // 1D array[channel] for sigmas of the cor(s,s)
  Int_t*      fTagSigCorss;      // (SIGMA OF COR(S,S))
  Double_t*   fT1d_av_sig_corss; // 1D array[channel] for expectation values of sigmas  the  of the cor(s,s)
  Int_t*      fTagAvSigCorss;    // (AVERAGED SIGMA OF COR(S,S))

  //...........................................................................
  Double_t**  fT2d_lf_cov; // 2D array[channel][channel] for (channel,channel) low frequency covariances
  Double_t*   fT2d1_lf_cov;
  Int_t*      fTagLfCov;

  Double_t**  fT2d_lf_cor; // 2D array[channel][channel] for (channel,channel) low frequency correlations
  Double_t*   fT2d1_lf_cor;
  Int_t*      fTagLfCor;

  //...........................................................................
  Double_t**  fT2d_hf_cov; // 2D array[channel][channel] for (channel,channel) low frequency covariances
  Double_t*   fT2d1_hf_cov;
  Int_t*      fTagHfCov;

  Double_t**  fT2d_hf_cor; // 2D array[channel][channel] for (channel,channel) low frequency correlations
  Double_t*   fT2d1_hf_cor;
  Int_t*      fTagHfCor;

  //------------------------- 2 tableaux (ci,cj)
  Double_t**  fT2d_lfcc_mostins;  // 2D array[Stin][Stin] for (Stin,Stin) mean cov(c,c)
  Double_t*   fT2d1_lfcc_mostins; // (relevant ones) averaged over samples
  Int_t*      fTagLFccMoStins;

  Double_t**  fT2d_hfcc_mostins;  // 2D array[Stin][Stin] for (Stin,Stin) mean cor(c,c)
  Double_t*   fT2d1_hfcc_mostins; // (relevant ones) averaged over samples
  Int_t*      fTagHFccMoStins;

  //------------------------------------------------------------------------------------

  Int_t**     fT2dCrysNumbersTable;
  Int_t*      fT1dCrysNumbersTable;

  std::ofstream    fFcout_f;

  Int_t       fFlagPrint;
  Int_t       fCodePrintComments, fCodePrintWarnings, fCodePrintAllComments, fCodePrintNoComment;


 public: 

  //................. constructors
  
  TEcnaRun();                            //  constructor without argument
  TEcnaRun(TEcnaObject*, const TString&);               //  constructors with argument (FOR USER'S DECLARATION)
  TEcnaRun(TEcnaObject*, const TString&, const Int_t&);

  //TEcnaRun(const TString&, const Int_t&, const TEcnaParPaths*, const TEcnaParCout*);
  //TEcnaRun(const TString&);               //  constructors with argument (FOR USER'S DECLARATION)
  //TEcnaRun(const TString&, const Int_t&);

  TEcnaRun(const TEcnaRun&);   //  copy constructor

  //.................... C++ methods

  //TEcnaRun&  operator=(const TEcnaRun&);  //  overloading of the operator=

  //................. destructor
  
  ~TEcnaRun() override;
  
  //...................................................... methods that will (should) be private

  void Init();

  void SetEcalSubDetector(const TString&);

  Bool_t GetPathForResults();

  Bool_t OpenRootFile(const Text_t *, const TString&);
  Bool_t CloseRootFile(const Text_t *);

  //======================================= methods for the user =========================================
  void GetReadyToReadData(const TString&, const Int_t&, const Int_t&, const Int_t&, const Int_t&, const Int_t&);
  void GetReadyToReadData(const TString&, const Int_t&, const Int_t&, const Int_t&, const Int_t&, const Int_t&, const Int_t&);
  
  Bool_t GetSampleAdcValues(const Int_t&, const Int_t&, const Int_t&, const Int_t&, const Double_t&);
  Bool_t ReadSampleAdcValues();
  Bool_t ReadSampleAdcValues(const Int_t&);

  void StartStopDate(const TString&, const TString&);
  void StartStopTime(time_t, time_t);

  //................... Calculation methods ( associated to GetReadyToReadData(...) )
  //-------------------- Standard Calculations
  void StandardCalculations();       // see list in the method itself (.cc file)

  void SampleMeans();                // Calculation of the expectation values over the events
                                     // for each sample and for each channel
  void SampleSigmas();               // Calculation of the variances over the events
                                     // for each sample and for each channel
  //...........................................
  void CovariancesBetweenSamples();  // Calculation of the (s,s) covariances over the events
                                     // for each channel
  void CorrelationsBetweenSamples(); // Calculation of the (s,s) correlations over the events
                                     // for each channel
  //..........................................................
  void Pedestals();                        
  void TotalNoise();
  void LowFrequencyNoise();
  void HighFrequencyNoise();
  void MeanCorrelationsBetweenSamples();
  void SigmaOfCorrelationsBetweenSamples();

  //..........................................................
  void AveragePedestals();
  void AverageTotalNoise();
  void AverageLowFrequencyNoise();
  void AverageHighFrequencyNoise();
  void AverageMeanCorrelationsBetweenSamples();
  void AverageSigmaOfCorrelationsBetweenSamples();

  //---------- Calculations involving cov and cor between channels
  //
  //   Recommended calling sequences: expert1, expert1 + expert2,  expert2
  //
  //   NOT recommended: expert2 + expert1  (lost of time and place)
  //
  //-------------------- Expert 1 Calculations
  void Expert1Calculations();       // see list in the method itself (.cc file)

  void LowFrequencyCovariancesBetweenChannels();
  void HighFrequencyCovariancesBetweenChannels();

  void LowFrequencyCorrelationsBetweenChannels();
  void HighFrequencyCorrelationsBetweenChannels();

  //-------------------- Expert 2 Calculations
  void Expert2Calculations();       // see list in the method itself (.cc file)

  void LowFrequencyMeanCorrelationsBetweenTowers();
  void HighFrequencyMeanCorrelationsBetweenTowers();

  void LowFrequencyMeanCorrelationsBetweenSCs();
  void HighFrequencyMeanCorrelationsBetweenSCs();

  //===================================== "technical" methods ==========================================

  void GetReadyToCompute(); // Make result root file name and check events
  void SampleValues();      // 3D histo of the sample ADC value for each triple (channel,sample,event)

  //.................................... Technical calculation methods (Stin = Tower or SC)
  void LowFrequencyMeanCorrelationsBetweenStins();
  void HighFrequencyMeanCorrelationsBetweenStins();

  //...................................... ROOT file methods
  TString GetRootFileName();
  TString GetRootFileNameShort();
  TString GetNewRootFileName();
  TString GetNewRootFileNameShort();

  Bool_t WriteRootFile();
  Bool_t WriteNewRootFile(const TString&);
  Bool_t WriteRootFile(const Text_t *, Int_t&);

  void TRootStinNumbers();
  void TRootNbOfEvts(const Int_t&);

  void TRootAdcEvt(const Int_t&, const Int_t&);

  void TRootMSp(const Int_t&);
  void TRootSSp(const Int_t&);

  void TRootCovCss(const Int_t&, const Int_t&);
  void TRootCorCss(const Int_t&, const Int_t&);

  void TRootLfCov();
  void TRootLfCor();

  void TRootHfCov();
  void TRootHfCor();

  void TRootLFccMoStins();
  void TRootHFccMoStins();

  void TRootPed();
  void TRootTno();
  void TRootMeanCorss();

  void TRootLfn();
  void TRootHfn();
  void TRootSigCorss();

  void TRootAvPed();
  void TRootAvEvCorss();
  void TRootAvSigCorss();
  void TRootAvTno();
  void TRootAvLfn();
  void TRootAvHfn();

  //................................ Flags Print Comments/Debug
  void PrintNoComment();   // (default) Set flags to forbid the printing of all the comments
                           // except ERRORS
  void PrintWarnings();    // Set flags to authorize printing of some warnings
  void PrintComments();    // Set flags to authorize printing of infos and some comments
                           // concerning initialisations
  void PrintAllComments(); // Set flags to authorize printing of all the comments

ClassDefOverride(TEcnaRun,1) // Calculation of correlated noises from data
};  

#endif    //  CL_TEcnaRun_H