TEcnaRun.cc

Go to the documentation of this file.

//----------Author's Names: B.Fabbro, FX Gentit DSM/IRFU/SPP CEA-Saclay
//----------Copyright: Those valid for CEA sofware
//----------Modified: 30/01/2014

#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaRun.h"

//--------------------------------------
//  TEcnaRun.cc
//  Class creation: 03 Dec 2002
//  Documentation: see TEcnaRun.h
//--------------------------------------

R__EXTERN TEcnaRootFile *gCnaRootFile;

ClassImp(TEcnaRun)
//___________________________________________________________________________
//

TEcnaRun::TEcnaRun()
{
//Constructor without argument: nothing special done

 // std::cout << "[Info Management] CLASS: TEcnaRun.           CREATE OBJECT: this = " << this << std::endl;
}

TEcnaRun::TEcnaRun(TEcnaObject* pObjectManager, const TString& SubDet)
{
//Constructor with argument: call to Init() and declare fEcal according to SubDet value ("EB" or "EE")

 // std::cout << "[Info Management] CLASS: TEcnaRun.           CREATE OBJECT: this = " << this << std::endl;

  Init();
  fObjectManager = (TEcnaObject*)pObjectManager;
  Long_t i_this = (Long_t)this;
  pObjectManager->RegisterPointer("TEcnaRun", i_this);

  //............................ fCnaParCout
  fCnaParCout = nullptr;
  Long_t iCnaParCout = pObjectManager->GetPointerValue("TEcnaParCout");
  if( iCnaParCout == 0 )
    {fCnaParCout = new TEcnaParCout(pObjectManager); /*fCnew++*/}
  else
    {fCnaParCout = (TEcnaParCout*)iCnaParCout;}

  //............................ fCnaParPaths
  fCnaParPaths = nullptr;
  Long_t iCnaParPaths = pObjectManager->GetPointerValue("TEcnaParPaths");
  if( iCnaParPaths == 0 )
    {fCnaParPaths = new TEcnaParPaths(pObjectManager); /*fCnew++*/}
  else
    {fCnaParPaths = (TEcnaParPaths*)iCnaParPaths;}

  //fCfgResultsRootFilePath    = fCnaParPaths->ResultsRootFilePath();
  //fCfgHistoryRunListFilePath = fCnaParPaths->HistoryRunListFilePath();

  //ffFileHeader = 0;
  //fconst Text_t *h_name  = "CnaHeader";   //==> voir cette question avec FXG
  //fconst Text_t *h_title = "CnaHeader";   //==> voir cette question avec FXG
  //ffFileHeader = new TEcnaHeader(h_name, h_title);     //fCnew++;

  //............................ fFileHeader
  const Text_t *h_name  = "CnaHeader";  //==> voir cette question avec FXG
  const Text_t *h_title = "CnaHeader";  //==> voir cette question avec FXG

  fFileHeader = nullptr;
  //Int_t iFileHeader = pObjectManager->GetPointerValue("TEcnaHeader");
  Long_t iFileHeader = 0;  // one TEcnaHeader object for each file since they can be open simultaneously 
  if( iFileHeader == 0 )
    {fFileHeader = new TEcnaHeader(pObjectManager, h_name, h_title); /*fCnew++*/}
  else
    {fFileHeader = (TEcnaHeader*)iFileHeader;}

  SetEcalSubDetector(SubDet.Data());
  fNbSampForFic = fEcal->MaxSampADC(); // DEFAULT Number of samples for ROOT file
}

TEcnaRun::TEcnaRun(TEcnaObject* pObjectManager, const TString& SubDet, const Int_t& NbOfSamples)
{
  //fCnaParPaths = 0; fCnaParPaths = new TEcnaParPaths();  //fCnew++;
  //fCnaParCout = 0;  fCnaParCout  = new TEcnaParCout();   //fCnew++;

  Init();
  fObjectManager = (TEcnaObject*)pObjectManager;
  Long_t i_this = (Long_t)this;
  pObjectManager->RegisterPointer("TEcnaRun", i_this);

  //............................ fCnaParCout
  fCnaParCout = nullptr;
  Long_t iCnaParCout = pObjectManager->GetPointerValue("TEcnaParCout");
  if( iCnaParCout == 0 )
    {fCnaParCout = new TEcnaParCout(pObjectManager); /*fCnew++*/}
  else
    {fCnaParCout = (TEcnaParCout*)iCnaParCout;}

  //............................ fCnaParPaths
  fCnaParPaths = nullptr;
  Long_t iCnaParPaths = pObjectManager->GetPointerValue("TEcnaParPaths");
  if( iCnaParPaths == 0 )
    {fCnaParPaths = new TEcnaParPaths(pObjectManager); /*fCnew++*/}
  else
    {fCnaParPaths = (TEcnaParPaths*)iCnaParPaths;}

  //fCfgResultsRootFilePath    = fCnaParPaths->ResultsRootFilePath();
  //fCfgHistoryRunListFilePath = fCnaParPaths->HistoryRunListFilePath();

  //............................ fFileHeader
  const Text_t *h_name  = "CnaHeader";  //==> voir cette question avec FXG
  const Text_t *h_title = "CnaHeader";  //==> voir cette question avec FXG

  fFileHeader = nullptr;
  //Long_t iFileHeader = pObjectManager->GetPointerValue("TEcnaHeader");
  Long_t iFileHeader = 0;  // one TEcnaHeader object for each file since they can be open simultaneously 
  if( iFileHeader == 0 )
    {fFileHeader = new TEcnaHeader(pObjectManager, h_name, h_title); /*fCnew++*/}
  else
    {fFileHeader = (TEcnaHeader*)iFileHeader;}

  SetEcalSubDetector(SubDet.Data());
  if( NbOfSamples>0 && NbOfSamples<=fEcal->MaxSampADC() )
    {
      fNbSampForFic = NbOfSamples;
    }
  else
    {
      std::cout << "TEcnaRun/CONSTRUCTOR> Number of required samples = " << NbOfSamples
       << ": OUT OF RANGE. Set to the default value (= " << fEcal->MaxSampADC() << ")."
       << fTTBELL << std::endl;
      fNbSampForFic = fEcal->MaxSampADC(); // DEFAULT Number of samples for file reading
    }

}

//.... return true or false according to the existence of the path. The path itself is in an attribute of fCnaParPaths.
Bool_t TEcnaRun::GetPathForResults(){return fCnaParPaths->GetPathForResultsRootFiles();}

void TEcnaRun::Init()
{
//Initialisation

  fCnew       = 0;
  fCdelete    = 0;
  fCnaCommand = 0;
  fCnaError   = 0;

  fTTBELL = '\007';

  //........................... TString file names init
  fgMaxCar  = (Int_t)512;

  //................ MiscDiag counters .................
  fMaxMsgIndexForMiscDiag = (Int_t)10;
  fNbOfMiscDiagCounters   = (Int_t)50;
  fMiscDiag = nullptr;

  fNumberOfEvents = 0;
  //............................. init pointers  ( Init() )
  fT3d_AdcValues  = nullptr;
  fT3d2_AdcValues = nullptr;
  fT3d1_AdcValues = nullptr;

  fT1d_StexStinFromIndex = nullptr;

  fT2d_NbOfEvts = nullptr;
  fT1d_NbOfEvts = nullptr; 

  fT2d_ev        = nullptr;
  fT1d_ev        = nullptr;
  fT2d_sig       = nullptr;
  fT1d_sig       = nullptr;

  fT3d_cov_ss    = nullptr;
  fT3d2_cov_ss   = nullptr;
  fT3d1_cov_ss   = nullptr;

  fT3d_cor_ss    = nullptr;
  fT3d2_cor_ss   = nullptr;
  fT3d1_cor_ss   = nullptr;

  fT2d_lf_cov  = nullptr;
  fT2d1_lf_cov = nullptr;

  fT2d_lf_cor  = nullptr;
  fT2d1_lf_cor = nullptr;

  fT2d_hf_cov  = nullptr;
  fT2d1_hf_cov = nullptr;

  fT2d_hf_cor  = nullptr;
  fT2d1_hf_cor = nullptr;

  fT2d_lfcc_mostins  = nullptr;
  fT2d1_lfcc_mostins = nullptr;

  fT2d_hfcc_mostins  = nullptr;
  fT2d1_hfcc_mostins = nullptr;

  fT1d_ev_ev            = nullptr;
  fT1d_evsamp_of_sigevt = nullptr;
  fT1d_ev_cor_ss        = nullptr;
  fT1d_av_mped          = nullptr;
  fT1d_av_totn          = nullptr;
  fT1d_av_lofn          = nullptr;
  fT1d_av_hifn          = nullptr;
  fT1d_av_ev_corss      = nullptr;
  fT1d_av_sig_corss     = nullptr;

  fT1d_sigevt_of_evsamp = nullptr;
  fT1d_evevt_of_sigsamp = nullptr;
  fT1d_sig_cor_ss       = nullptr;

  fT2dCrysNumbersTable = nullptr;
  fT1dCrysNumbersTable = nullptr;

  //................................ tags   ( Init() )
  fTagStinNumbers = nullptr;

  fTagNbOfEvts = nullptr;

  fTagAdcEvt = nullptr;

  fTagMSp    = nullptr;
  fTagSSp    = nullptr;

  fTagCovCss = nullptr;
  fTagCorCss = nullptr;

  fTagHfCov  = nullptr;
  fTagHfCor  = nullptr;
  fTagLfCov  = nullptr;
  fTagLfCor  = nullptr;

  fTagLFccMoStins = nullptr;
  fTagHFccMoStins = nullptr;

  fTagPed = nullptr;
  fTagTno = nullptr;
  fTagMeanCorss = nullptr;

  fTagLfn = nullptr;
  fTagHfn = nullptr;
  fTagSigCorss = nullptr;

  fTagAvPed = nullptr;
  fTagAvTno = nullptr;
  fTagAvLfn = nullptr;
  fTagAvHfn = nullptr;

  fTagAvMeanCorss = nullptr;
  fTagAvSigCorss  = nullptr;

  //................................................... Code Print  ( Init() )
  fCodePrintNoComment   = fCnaParCout->GetCodePrint("NoComment");
  fCodePrintWarnings    = fCnaParCout->GetCodePrint("Warnings ");      // => default value
  fCodePrintComments    = fCnaParCout->GetCodePrint("Comments");
  fCodePrintAllComments = fCnaParCout->GetCodePrint("AllComments");

  fFlagPrint = fCodePrintWarnings;

  //...................................................
  gCnaRootFile = nullptr;
  fOpenRootFile = kFALSE;
  fReadyToReadData = 0;

  //.............................................. Miscellaneous
  fSpecialStexStinNotIndexed = -1;

  fStinIndexBuilt     = 0;
  fBuildEvtNotSkipped = 0;

  fMemoReadNumberOfEventsforSamples = 0;

}// end of Init()

//========================================================================
void  TEcnaRun::SetEcalSubDetector(const TString& SubDet)
{
 // Set Subdetector (EB or EE)

  Int_t MaxCar = fgMaxCar;
  fFlagSubDet.Resize(MaxCar);
  fFlagSubDet = SubDet.Data();


  fEcal          = nullptr; fEcal          = new TEcnaParEcal(fFlagSubDet.Data());           //fCnew++;
  fEcalNumbering = nullptr; fEcalNumbering = new TEcnaNumbering(fFlagSubDet.Data(), fEcal);  //fCnew++;
  fCnaWrite      = nullptr;

  fCnaWrite =
    new TEcnaWrite(fFlagSubDet.Data(), fCnaParPaths, fCnaParCout, fEcal, fEcalNumbering); //fCnew++;

  if( fFlagSubDet == "EB" ){fStexName = "SM "; fStinName = "tower";}
  if( fFlagSubDet == "EE" ){fStexName = "Dee"; fStinName = " SC  ";}
}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//
//                    copy constructor
//
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

TEcnaRun::TEcnaRun(const TEcnaRun& dcop)
{
  std::cout << "*TEcnaRun::TEcnaRun(const TEcnaRun& dcop)> "
       << " Now is the time to write a copy constructor" 
       << std::endl;
  
  //{ Int_t cintoto;  cin >> cintoto; }
  
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//
//                    overloading of the operator=
//
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

//TEcnaRun& TEcnaRun::operator=(const TEcnaRun& dcop)
//{
//Overloading of the operator=
//
//  fCopy(dcop);
//  return *this;
//}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//
//                            destructor
//
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

TEcnaRun::~TEcnaRun()
{
  //Destructor
  
  if(fFlagPrint == fCodePrintAllComments)
    {
      std::cout << "*TEcnaRun::~TEcnaRun()> Entering destructor." << std::endl;
    }

  if(fFlagPrint != fCodePrintNoComment || fFlagPrint == fCodePrintWarnings )
    {
      if( fBuildEvtNotSkipped > 0 )
    {
      std::cout << "************************************************************************************* "
           << std::endl;
      std::cout << "*TEcnaRun::~TEcnaRun()> Nb of calls to GetSampleAdcValues by cmsRun: "
           << fBuildEvtNotSkipped << std::endl;
      std::cout << "************************************************************************************* "
           << std::endl;
    }
    }

  if(fFlagPrint == fCodePrintAllComments)
    {
      Int_t misc_czero = 0;
      for(Int_t i = 0; i < fNbOfMiscDiagCounters; i++)
    {
      if( fMiscDiag[i] != 0 )
        {
          std::cout << "                          fMiscDiag Counter "
           << std::setw(3) << i << " = " << std::setw(9) << fMiscDiag[i]
           << " (INFO: alloc on non zero freed zone) " << std::endl;
        }
      else
        {
          misc_czero++;
        }
    }
      std::cout << "                          Nb of fMiscDiag counters at zero: "
       << misc_czero << " (total nb of counters: "
       << fNbOfMiscDiagCounters << ")" << std::endl;
    }
  
  if (fMiscDiag                != nullptr){delete [] fMiscDiag;                 fCdelete++;}
  
  //if (fFileHeader              != 0){delete fFileHeader;                  fCdelete++;}
  //if (fEcal                    != 0){delete fEcal;                        fCdelete++;}
  //if (fEcalNumbering           != 0){delete fEcalNumbering;               fCdelete++;}
  //if (fCnaParCout              != 0){delete fCnaParCout;                  fCdelete++;}
  //if (fCnaParPaths             != 0){delete fCnaParPaths;                 fCdelete++;}
  //if (fCnaWrite                != 0){delete fCnaWrite;                    fCdelete++;}
  
  if (fT1d_StexStinFromIndex   != nullptr){delete [] fT1d_StexStinFromIndex;    fCdelete++;}

  if (fT2d_NbOfEvts        != nullptr){delete [] fT2d_NbOfEvts;         fCdelete++;}
  if (fT1d_NbOfEvts        != nullptr){delete [] fT1d_NbOfEvts;         fCdelete++;}

  if (fT3d_AdcValues        != nullptr){delete [] fT3d_AdcValues;         fCdelete++;}
  if (fT3d2_AdcValues       != nullptr){delete [] fT3d2_AdcValues;        fCdelete++;}
  if (fT3d1_AdcValues       != nullptr){delete [] fT3d1_AdcValues;        fCdelete++;}

  if (fT2d_ev              != nullptr){delete [] fT2d_ev;               fCdelete++;}
  if (fT1d_ev              != nullptr){delete [] fT1d_ev;               fCdelete++;}

  if (fT2d_sig             != nullptr){delete [] fT2d_sig;              fCdelete++;}
  if (fT1d_sig             != nullptr){delete [] fT1d_sig;              fCdelete++;}

  if (fT3d_cov_ss          != nullptr){delete [] fT3d_cov_ss;           fCdelete++;}
  if (fT3d2_cov_ss         != nullptr){delete [] fT3d2_cov_ss;          fCdelete++;}
  if (fT3d1_cov_ss         != nullptr){delete [] fT3d1_cov_ss;          fCdelete++;}

  if (fT3d_cor_ss          != nullptr){delete [] fT3d_cor_ss;           fCdelete++;}
  if (fT3d2_cor_ss         != nullptr){delete [] fT3d2_cor_ss;          fCdelete++;}
  if (fT3d1_cor_ss         != nullptr){delete [] fT3d1_cor_ss;          fCdelete++;}

  if (fT2d_lf_cov          != nullptr){delete [] fT2d_lf_cov;           fCdelete++;}
  if (fT2d1_lf_cov         != nullptr){delete [] fT2d1_lf_cov;          fCdelete++;}

  if (fT2d_lf_cor          != nullptr){delete [] fT2d_lf_cor;           fCdelete++;}
  if (fT2d1_lf_cor         != nullptr){delete [] fT2d1_lf_cor;          fCdelete++;}

  if (fT2d_hf_cov          != nullptr){delete [] fT2d_hf_cov;           fCdelete++;}
  if (fT2d1_hf_cov         != nullptr){delete [] fT2d1_hf_cov;          fCdelete++;}

  if (fT2d_hf_cor          != nullptr){delete [] fT2d_hf_cor;           fCdelete++;}
  if (fT2d1_hf_cor         != nullptr){delete [] fT2d1_hf_cor;          fCdelete++;}

  if (fT2d_lfcc_mostins    != nullptr){delete [] fT2d_lfcc_mostins;     fCdelete++;}
  if (fT2d1_lfcc_mostins   != nullptr){delete [] fT2d1_lfcc_mostins ;   fCdelete++;}

  if (fT2d_hfcc_mostins    != nullptr){delete [] fT2d_hfcc_mostins ;    fCdelete++;}
  if (fT2d1_hfcc_mostins   != nullptr){delete [] fT2d1_hfcc_mostins;    fCdelete++;}

  if (fT1d_ev_ev             != nullptr){delete [] fT1d_ev_ev;              fCdelete++;}
  if (fT1d_evsamp_of_sigevt  != nullptr){delete [] fT1d_evsamp_of_sigevt;   fCdelete++;}
  if (fT1d_ev_cor_ss         != nullptr){delete [] fT1d_ev_cor_ss;          fCdelete++;}
  if (fT1d_av_mped           != nullptr){delete [] fT1d_av_mped;            fCdelete++;}
  if (fT1d_av_totn           != nullptr){delete [] fT1d_av_totn;            fCdelete++;}
  if (fT1d_av_lofn           != nullptr){delete [] fT1d_av_lofn;            fCdelete++;}
  if (fT1d_av_hifn           != nullptr){delete [] fT1d_av_hifn;            fCdelete++;}
  if (fT1d_av_ev_corss       != nullptr){delete [] fT1d_av_ev_corss;        fCdelete++;}
  if (fT1d_av_sig_corss      != nullptr){delete [] fT1d_av_sig_corss;       fCdelete++;}

  if (fT1d_sigevt_of_evsamp    != nullptr){delete [] fT1d_sigevt_of_evsamp;     fCdelete++;}
  if (fT1d_evevt_of_sigsamp    != nullptr){delete [] fT1d_evevt_of_sigsamp;     fCdelete++;}
  if (fT1d_sig_cor_ss          != nullptr){delete [] fT1d_sig_cor_ss;           fCdelete++;}

  if (fT2dCrysNumbersTable     != nullptr){delete [] fT2dCrysNumbersTable;      fCdelete++;}
  if (fT1dCrysNumbersTable     != nullptr){delete [] fT1dCrysNumbersTable;      fCdelete++;}

  if (fTagStinNumbers          != nullptr){delete [] fTagStinNumbers;           fCdelete++;}
  if (fTagNbOfEvts             != nullptr){delete [] fTagNbOfEvts;              fCdelete++;}
  if (fTagAdcEvt               != nullptr){delete [] fTagAdcEvt;                fCdelete++;}
  if (fTagMSp                  != nullptr){delete [] fTagMSp;                   fCdelete++;}
  if (fTagSSp                  != nullptr){delete [] fTagSSp;                   fCdelete++;}

  if (fTagCovCss               != nullptr){delete [] fTagCovCss;                fCdelete++;}
  if (fTagCorCss               != nullptr){delete [] fTagCorCss;                fCdelete++;}

  if (fTagHfCov                != nullptr){delete [] fTagHfCov;                 fCdelete++;}
  if (fTagHfCor                != nullptr){delete [] fTagHfCor;                 fCdelete++;}
  if (fTagLfCov                != nullptr){delete [] fTagLfCov;                 fCdelete++;}
  if (fTagLfCor                != nullptr){delete [] fTagLfCor;                 fCdelete++;}

  if (fTagLFccMoStins          != nullptr){delete [] fTagLFccMoStins;           fCdelete++;}
  if (fTagHFccMoStins          != nullptr){delete [] fTagHFccMoStins;           fCdelete++;}

  if (fTagPed                 != nullptr){delete [] fTagPed;                  fCdelete++;}
  if (fTagTno                 != nullptr){delete [] fTagTno;                  fCdelete++;}
  if (fTagMeanCorss           != nullptr){delete [] fTagMeanCorss;            fCdelete++;}

  if (fTagLfn                != nullptr){delete [] fTagLfn;                 fCdelete++;}
  if (fTagHfn                != nullptr){delete [] fTagHfn;                 fCdelete++;}
  if (fTagSigCorss           != nullptr){delete [] fTagSigCorss;            fCdelete++;}

  if (fTagAvPed          != nullptr){delete [] fTagAvPed;           fCdelete++;}
  if (fTagAvTno          != nullptr){delete [] fTagAvTno;           fCdelete++;}
  if (fTagAvLfn          != nullptr){delete [] fTagAvLfn;           fCdelete++;}
  if (fTagAvHfn          != nullptr){delete [] fTagAvHfn;           fCdelete++;}
  if (fTagAvMeanCorss    != nullptr){delete [] fTagAvMeanCorss;     fCdelete++;}
  if (fTagAvSigCorss     != nullptr){delete [] fTagAvSigCorss;      fCdelete++;}

  if ( fCnew != fCdelete )
    {
      std::cout << "!TEcnaRun::~TEcnaRun()> WRONG MANAGEMENT OF MEMORY ALLOCATIONS: fCnew = "
       << fCnew << ", fCdelete = " << fCdelete << fTTBELL << std::endl;
    }
  else
    {
      // std::cout << "*TEcnaRun::~TEcnaRun()> Management of memory allocations: OK. fCnew = "
      //   << fCnew << ", fCdelete = " << fCdelete << std::endl;
    }

  if(fFlagPrint == fCodePrintAllComments)
    {
      std::cout << "*TEcnaRun::~TEcnaRun()> Exiting destructor (this = " << this << ")." << std::endl
       << "~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#~#"
       << std::endl;
    }

 // std::cout << "[Info Management] CLASS: TEcnaRun.           DESTROY OBJECT: this = " << this << std::endl;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//
//                             M  E  T  H  O  D  S
//
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

//============================================================================
//
//                           GetReadyToReadData(...)
//                  
//  Preparation of the result file name + tags allocations
//  + ADC event distribution array allocation + nb of events array allocation
//
//============================================================================

void TEcnaRun::GetReadyToReadData(const TString&     typ_ana, const Int_t& run_number,
                 const Int_t& nfirst,  const Int_t& nlast, const Int_t& nbevts,
                 const Int_t& Stex)
{
//Preparation of the data reading. Set part of the header. No RunType as argument.
//Use default value = 99999999 and call method with all the arguments (see below)

  Int_t RunType = 99999999;
  GetReadyToReadData(typ_ana, run_number, nfirst, nlast, nbevts, Stex, RunType);
}
//--------------------------------------------------------------------------------
void TEcnaRun::GetReadyToReadData( const TString&       typ_ana, const Int_t& run_number,
                   const Int_t&  nfirst,  const Int_t& nlast, const Int_t& nbevts,
                   const Int_t&  Stex,    const Int_t& run_type)
{
//Preparation of the data reading. Set part of the header
//
//   [nfirst, nlast] = [1,50] (50 events) or [151,300] (150 events) or etc...

  Int_t nrangeevts = nlast - nfirst + 1;  // number of events in range

  if( nrangeevts < nbevts )
    {
      if( nlast >= nfirst )
    {
      std::cout << "*TEcnaRun::GetReadyToReadData(...)> --- WARNING ---> number of events = " << nbevts
           << ", out of range (range = " << nfirst << "," << nlast << ")" << std::endl
           << "                                    The number of found events will be less " << std::endl
           << "                                    than the number of requested events." << std::endl;
    }
      if( nlast < nfirst )
    {
      std::cout << "*TEcnaRun::GetReadyToReadData(...)> --- INFO ---> last requested event number = " << nlast
           << ", less than first requested event number (= " << nfirst << ")" << std::endl
           << "                                    File will be read until EOF if the number of found events" << std::endl
           << "                                    remains less than the number of requested events." << std::endl;
    }

    }
  
  //............. allocation for counters
  fMiscDiag = new Int_t[fNbOfMiscDiagCounters];       fCnew++;
  for (Int_t iz=0; iz<fNbOfMiscDiagCounters; iz++){fMiscDiag[iz] = (Int_t)0;}

  //************** CHECK OF ARGUMENTS: nfirst_arg and nbevts_arg
  Int_t nentries = 99999999;   // => to be reintroduced as argument (like run_type) (?)
  if ( nfirst <= nentries )
    {
      //--------------------- test positivity of nfirst_arg        (GetReadyToReadData)
      if ( nfirst > 0 )
    {
      //-------- test compatibility between the last requested event number
      //         and the number of entries
      if( nlast <= nentries )
        { 
          const Text_t *h_name  = "CnaHeader";   //==> voir cette question avec FXG
          const Text_t *h_title = "CnaHeader";   //==> voir cette question avec FXG
          //fFileHeader->HeaderParameters(h_name,     h_title ,
          //                      typ_ana,    fNbSampForFic,
          //                  run_number, nfirst,  nlast,  nbevts,
          //                  Stex,       nentries);                 fCnew++;

          if ( fEcal->MaxStinEcnaInStex() > 0  &&  fEcal->MaxCrysInStin() > 0  &&  fNbSampForFic > 0 ) 
        {
          if( fFileHeader == nullptr ){fFileHeader = new TEcnaHeader(fObjectManager, h_name, h_title);} // fCnew++;

          fFileHeader->HeaderParameters(typ_ana,    fNbSampForFic,
                        run_number, nfirst, nlast, nbevts,
                        Stex,       run_type);

          // After this call to TEcnaHeader, we have:               (GetReadyToReadData)
          //     fFileHeader->fTypAna            = typ_ana
          //     fFileHeader->fNbOfSamples       = fNbSampForFic
          //     fFileHeader->fRunNumber         = run_number
          //     fFileHeader->fFirstReqEvtNumber = nfirst
          //     fFileHeader->fLastReqEvtNumber  = nlast
          //     fFileHeader->fReqNbOfEvts       = nbevts
          //     fFileHeader->fStex              = Stex number
          //     fFileHeader->fRunType           = run_type
          
          // fFileHeader->Print();
          
          // {Int_t cintoto; std::cout << "taper 0 pour continuer" << std::endl; cin >> cintoto;}
          
          //  fFileHeader->SetName("CnaHeader");              *======> voir FXG
          //  fFileHeader->SetTitle("CnaHeader");
          
          //......................................... allocation tags + init of them (GetReadyToReadData)
          
          fTagStinNumbers = new Int_t[1]; fCnew++; fTagStinNumbers[0] = (Int_t)0;
          fTagNbOfEvts    = new Int_t[1]; fCnew++; fTagNbOfEvts[0]    = (Int_t)0;
          
          fTagAdcEvt      = new Int_t[fEcal->MaxCrysEcnaInStex()]; fCnew++;
          for (Int_t iz=0; iz<fEcal->MaxCrysEcnaInStex(); iz++){fTagAdcEvt[iz] = (Int_t)0;}
          
          fTagMSp         = new Int_t[1]; fCnew++; fTagMSp[0] = (Int_t)0;
          fTagSSp         = new Int_t[1]; fCnew++; fTagSSp[0] = (Int_t)0;
          
          fTagCovCss      = new Int_t[fEcal->MaxCrysEcnaInStex()]; fCnew++;
          for (Int_t iz=0; iz<fEcal->MaxCrysEcnaInStex(); iz++){fTagCovCss[iz] = (Int_t)0;}
          
          fTagCorCss      = new Int_t[fEcal->MaxCrysEcnaInStex()]; fCnew++;
          for (Int_t iz=0; iz<fEcal->MaxCrysEcnaInStex(); iz++){fTagCorCss[iz] = (Int_t)0;}
          
          fTagLfCov       = new Int_t[1]; fCnew++; fTagLfCov[0] = (Int_t)0;
          fTagLfCor       = new Int_t[1]; fCnew++; fTagLfCor[0] = (Int_t)0;
          
          fTagHfCov       = new Int_t[1]; fCnew++; fTagHfCov[0] = (Int_t)0;
          fTagHfCor       = new Int_t[1]; fCnew++; fTagHfCor[0] = (Int_t)0;
          
          fTagLFccMoStins = new Int_t[1]; fCnew++; fTagLFccMoStins[0] = (Int_t)0;
          fTagHFccMoStins = new Int_t[1]; fCnew++; fTagHFccMoStins[0] = (Int_t)0;
          
          fTagPed         = new Int_t[1]; fCnew++; fTagPed[0]         = (Int_t)0;
          fTagTno         = new Int_t[1]; fCnew++; fTagTno[0]         = (Int_t)0;
          fTagMeanCorss   = new Int_t[1]; fCnew++; fTagMeanCorss[0]   = (Int_t)0; 
          
          fTagLfn         = new Int_t[1]; fCnew++; fTagLfn[0]         = (Int_t)0; 
          fTagHfn         = new Int_t[1]; fCnew++; fTagHfn[0]         = (Int_t)0; 
          fTagSigCorss    = new Int_t[1]; fCnew++; fTagSigCorss[0]    = (Int_t)0; 
          
          fTagAvPed       = new Int_t[1]; fCnew++; fTagAvPed[0]       = (Int_t)0; 
          fTagAvTno       = new Int_t[1]; fCnew++; fTagAvTno[0]       = (Int_t)0; 
          fTagAvLfn       = new Int_t[1]; fCnew++; fTagAvLfn[0]       = (Int_t)0; 
          fTagAvHfn       = new Int_t[1]; fCnew++; fTagAvHfn[0]       = (Int_t)0; 
          fTagAvMeanCorss = new Int_t[1]; fCnew++; fTagAvMeanCorss[0] = (Int_t)0; 
          fTagAvSigCorss  = new Int_t[1]; fCnew++; fTagAvSigCorss[0]  = (Int_t)0; 
          
          //====================================================================================
          //
          //   allocation for fT1d_StexStinFromIndex[] and init to fSpecialStexStinNotIndexed
          //
          //====================================================================================
          
          if(fT1d_StexStinFromIndex == nullptr)
            {
              fT1d_StexStinFromIndex = new Int_t[fEcal->MaxStinEcnaInStex()];          fCnew++;
            }
          for ( Int_t i0_Stin = 0; i0_Stin < fEcal->MaxStinEcnaInStex(); i0_Stin++ )
            {
              fT1d_StexStinFromIndex[i0_Stin] = fSpecialStexStinNotIndexed;
            }
          
          //-------------------------------------------------------------  (GetReadyToReadData)
      
          //====================================================================================
          //
          //   allocation of the 3D array fT3d_AdcValues[channel][sample][events] (ADC values)
          //
          //   This array is filled in the GetSampleAdcValues(...) method
          //
          //====================================================================================
          
          if(fT3d_AdcValues == nullptr)
            {
              //............ Allocation for the 3d array 
              std::cout << "*TEcnaRun::GetReadyToReadData(...)> Allocation of 3D array for ADC distributions."
               << " Nb of requested evts = " << fFileHeader->fReqNbOfEvts << std::endl
               << "                                    This number must not be too large"
               << " (no failure after this message means alloc OK)." << std::endl;

              fT3d_AdcValues = new Double_t**[fEcal->MaxCrysEcnaInStex()];        fCnew++;  

              fT3d2_AdcValues  =
            new  Double_t*[fEcal->MaxCrysEcnaInStex()*
                       fNbSampForFic];                          fCnew++;  

              fT3d1_AdcValues  =
            new   Double_t[fEcal->MaxCrysEcnaInStex()*
                       fNbSampForFic*
                       fFileHeader->fReqNbOfEvts];                       fCnew++;

              for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++){
            fT3d_AdcValues[i0StexEcha] = &fT3d2_AdcValues[0] + i0StexEcha*fNbSampForFic;
            for(Int_t j0Sample=0; j0Sample<fNbSampForFic; j0Sample++){
              fT3d2_AdcValues[fNbSampForFic*i0StexEcha + j0Sample] = &fT3d1_AdcValues[0]+
                fFileHeader->fReqNbOfEvts*(fNbSampForFic*i0StexEcha+j0Sample);}}
            }
          //................................. Init to zero
          for (Int_t iza=0; iza<fEcal->MaxCrysEcnaInStex(); iza++)
            {
              for (Int_t izb=0; izb<fNbSampForFic; izb++)
            {
              for (Int_t izc=0; izc<fFileHeader->fReqNbOfEvts; izc++)
                {
                  if( fT3d_AdcValues[iza][izb][izc] != (Double_t)0 )
                {
                  fMiscDiag[0]++;
                  fT3d_AdcValues[iza][izb][izc] = (Double_t)0;
                }
                }
            }
            }

          //--------------------------------------------------------- (GetReadyToReadData)    
          //====================================================================================
          //
          //   allocation of the 2D array fT2d_NbOfEvts[channel][sample] (Max nb of evts)
          //
          //====================================================================================
          
          if (fT2d_NbOfEvts == nullptr)
            {
              fT2d_NbOfEvts  = new Int_t*[fEcal->MaxCrysEcnaInStex()];           fCnew++;
              fT1d_NbOfEvts  = new  Int_t[fEcal->MaxCrysEcnaInStex()*
                              fNbSampForFic];            fCnew++;
              
              for(Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
            {
              fT2d_NbOfEvts[i0StexEcha] =
                &fT1d_NbOfEvts[0] + i0StexEcha*fNbSampForFic;
            }
              
              //................ Init the array to 0
              for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
            {
              for(Int_t i0Sample=0; i0Sample<fNbSampForFic; i0Sample++)
                {
                  fT2d_NbOfEvts[i0StexEcha][i0Sample] = 0;
                }
            }
            }
          else
            {
              std::cerr << "!TEcnaRun::GetReadyToReadData(...)> *** ERROR *** No allocation for fT2d_NbOfEvts!"
               << " Pointer already not NULL " << fTTBELL << std::endl;
              // {Int_t cintoto; std::cout << "Enter: 0 and RETURN to continue or: CTRL C to exit"
              //           << std::endl; std::cin >> cintoto;}
            }
        }
          else
        {
          std::cerr << std::endl
               << "!TEcnaRun::GetReadyToReadData(...)> "
               << " *** ERROR *** " << std::endl
               << " --------------------------------------------------"
               << std::endl
               << " NULL or NEGATIVE values for arguments" << std::endl
               << " with expected positive values:"        << std::endl
               << " Number of Stins in Stex = " << fEcal->MaxStinEcnaInStex()  << std::endl
               << " Number of crystals in Stin     = " << fEcal->MaxCrysInStin() << std::endl
               << " Number of samples by channel    = " << fNbSampForFic << std::endl
               << std::endl
               << std::endl
               << " hence, no memory allocation for array member has been performed." << std::endl;
          
          std::cout << "Enter: 0 and RETURN to continue or: CTRL C to exit";
          Int_t toto;
          std::cin >> toto;
        }
          //----------------------------------------------------------- (GetReadyToReadData)
          if(fFlagPrint == fCodePrintAllComments ){
        std::cout << std::endl;
        std::cout << "*TEcnaRun::GetReadyToReadData(...)>" << std::endl
             << "          The method has been called with the following argument values:"  << std::endl
             << "          Analysis name                = "
             << fFileHeader->fTypAna << std::endl
             << "          Run number                   = "
             << fFileHeader->fRunNumber << std::endl
             << "          Run type                     = "
             << fFileHeader->fRunType << std::endl
             << "          First requested event number = "
             << fFileHeader->fFirstReqEvtNumber << std::endl
             << "          Last requested event number  = "
             << fFileHeader->fLastReqEvtNumber << std::endl
             << "          " << fStexName.Data() << " number                  = "
             << fFileHeader->fStex << std::endl
             << "          Number of " << fStinName.Data()
             << " in " << fStexName.Data() << "       = "
             << fEcal->MaxStinEcnaInStex()  << std::endl
             << "          Number of crystals in " << fStinName.Data() << "  = "
             << fEcal->MaxCrysInStin()  << std::endl
             << "          Number of samples by channel = "
             << fNbSampForFic  << std::endl
             << std::endl;}
          
          fReadyToReadData = 1;                        // set flag
        }
      else
        {
          if (fFlagPrint != fCodePrintNoComment){
        std::cout << "!TEcnaRun::GetReadyToReadData(...) > WARNING/CORRECTION:" << std::endl
             << "! The fisrt requested event number is not positive (nfirst = " << nfirst << ") "
             << fTTBELL << std::endl;}
        }
    }
      else
    {
      if (fFlagPrint != fCodePrintNoComment){
        std::cout << std::endl << "!TEcnaRun::GetReadyToReadData(...)> WARNING/CORRECTION:" << std::endl
         << "! The number of requested events (nbevts = " << nbevts << ") is too large." << std::endl
         << "! Last event number = " << nlast << " > number of entries = " << nentries << ". "
         << fTTBELL << std::endl << std::endl;}
    }
    }
  else
    {
      std::cout << "!TEcnaRun::GetReadyToReadData(...) *** ERROR ***> "
       << " The first requested event number is greater than the number of entries."
       << fTTBELL << std::endl;
    }
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "*TEcnaRun::GetReadyToReadData(...)> Leaving the method. fReadyToReadData = "
     << fReadyToReadData << std::endl; }

} // end of GetReadyToReadData

//====================================================================================================
//
//     GetSampleAdcValues: method called by the CMSSW analyzer (cmsRun)
//
//     At each event, put the Sample ADC value in the 3D array: fT3d_AdcValues[i0StexEcha][i0Sample][i0EventIndex]
//
//         |============================================================|
//         |                                                            |
//         | (Stex,Stin) means: (SM,tower) for EB and: (Dee,SC) for EE  |
//         |                                                            |
//         |============================================================|
//
//    (Stin number <-> Stin index correspondance, ADC sample values)
//
//        THIS METHOD IS CALLED INSIDE THE LOOPS OVER:
//          ( EVENTS (tower or SC (CRYSTAL in tower or SC (SAMPLES)))) 
//
//  Arguments: Event    = event number.              Range = [ 1, fFileHeader->fReqNbOfEvts ]
//             n1StexStin = Stin number in Stex.     Range = [ 1, fEcal->MaxStinEcnaInStex() ]
//             i0StinEcha = channel number in Stin.  Range = [ 0, fEcal->MaxCrysInStin() - 1 ]     
//             sample   = ADC sample number.         Range = [ 0, fNbSampForFic - 1 ]
//             adcvalue = ADC sample value.
//
//====================================================================================================
Bool_t TEcnaRun::GetSampleAdcValues(const Int_t&    n1EventNumber, const Int_t& n1StexStin,
                    const Int_t&    i0StinEcha,    const Int_t& i0Sample,
                    const Double_t& adcvalue)
{
  //Building of the arrays fT1d_StexStinFromIndex[] and fT3d_AdcValues[][][]
  
  fBuildEvtNotSkipped++;      // event not skipped by cmsRun

  Bool_t ret_code = kFALSE;

  Int_t i0EventIndex   = n1EventNumber - 1;  // INDEX FOR Event number  
  Int_t i0StexStinEcna = n1StexStin - 1;     // INDEX FOR StexStin = Number_of_the_Stin_in_Stex - 1

  Int_t i_trouve = 0;
  //.................................................................. (GetSampleAdcValues)
  if(fReadyToReadData == 1)  
    {
      if( n1StexStin>= 1 && n1StexStin <= fEcal->MaxStinEcnaInStex() )
    {      
      if( i0StinEcha >= 0 && i0StinEcha < fEcal->MaxCrysInStin() )
        {
          if( i0Sample >= 0 && i0Sample < fEcal->MaxSampADC() )
        {
          //..... Put the StexStin number in 1D array fT1d_StexStinFromIndex[] = Stin index + 1
          if( fT1d_StexStinFromIndex != nullptr )  // table fT1d_StexStinFromIndex[index] already allocated
            {
              ret_code = kTRUE;

              // StexStin already indexed
              if( n1StexStin == fT1d_StexStinFromIndex[i0StexStinEcna] ){i_trouve = 1;}
          
              // StexStin index not found: set index for new StexStin
              if (i_trouve != 1 )
            {
              if( fT1d_StexStinFromIndex[i0StexStinEcna] == fSpecialStexStinNotIndexed )
                {
                  fT1d_StexStinFromIndex[i0StexStinEcna] = n1StexStin;
                  fFileHeader->fStinNumbersCalc = 1;
                  fTagStinNumbers[0] = 1;
                  fStinIndexBuilt++;                      //  number of found Stins
                  
                  if(fFlagPrint == fCodePrintAllComments)
                {
                  if( fStinIndexBuilt == 1 )
                    {
                      std::cout << std::endl << "*TEcnaRun::GetSampleAdcValues(...)> event " << n1EventNumber
                       << " : first event for " << fStexName.Data() << " " << fFileHeader->fStex
                       << "; " << fStinName.Data() << "s : ";
                    }
                  if( fFlagSubDet == "EB" )
                    {std::cout << fT1d_StexStinFromIndex[i0StexStinEcna] << ", ";}
                  if( fFlagSubDet == "EE" )
                    {std::cout << fEcalNumbering->
                       GetDeeSCConsFrom1DeeSCEcna(fFileHeader->fStex, fT1d_StexStinFromIndex[i0StexStinEcna])
                      << ", ";}
                }
                  //.................................................... (GetSampleAdcValues)
                  if(fFlagPrint == fCodePrintAllComments)
                {
                  std::cout << " (" << fStinIndexBuilt << " " << fStinName.Data()
                       << " found), channel " << i0StinEcha << ", i0Sample " << i0Sample << std::endl;
                }
                  ret_code = kTRUE;
                } // if ( fT1d_StexStinFromIndex[i0StexStinEcna] == fSpecialStexStinNotIndexed )
              else
                {
                  std::cout << "!TEcnaRun::GetSampleAdcValues(...)> *** ERROR ***> NOT ALLOWED if RESULT. " 
                   << " n1StexStin = " << n1StexStin << ", fT1d_StexStinFromIndex["
                   << i0StexStinEcna << "] = "
                   << fT1d_StexStinFromIndex[i0StexStinEcna]
                   << ", fStinIndexBuilt = " << fStinIndexBuilt
                   << fTTBELL << std::endl;
                  ret_code = kFALSE;
                }
            }  //  if (i_trouve != 1 )

            } //  if( fT1d_StexStinFromIndex != 0 ) 
          else
            {
              std::cout << "!TEcnaRun, GetSampleAdcValues *** ERROR ***> "
               << " fT1d_StexStinFromIndex = " << fT1d_StexStinFromIndex
               << " fT1d_StexStinFromIndex[] ALLOCATION NOT DONE" << fTTBELL << std::endl;
              ret_code = kFALSE;
            } //.................................................................. (GetSampleAdcValues)
        } // end of if( i0Sample >= 0 && i0Sample < fNbSampForFic )
          else
        {
          //.......Reading data => Message and error only if sample >= fEcal->MaxSampADC()
          //       (not fNbSampForFic, the later is used only for calculations)
          if( i0Sample >= fEcal->MaxSampADC() )
            {
              std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
               << " sample number = " << i0Sample << ". OUT OF BOUNDS"
               << " (max = " << fNbSampForFic << ")"
               << fTTBELL << std::endl;
              ret_code = kFALSE;
            }
          else
            {
              ret_code = kTRUE;
            }
        }// else of if( i0Sample >= 0 && i0Sample < fNbSampForFic )
        } // end of if( i0StinEcha >= 0 && i0StinEcha < fEcal->MaxCrysInStin() )
      else
        {
          std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
           << " channel number in " << fStinName.Data() << " = " << i0StinEcha << ". OUT OF BOUNDS"
           << " (max = " << fEcal->MaxCrysInStin() << ")" 
           << fTTBELL << std::endl;
          ret_code = kFALSE;
        } // else of if( i0StinEcha >= 0 && i0StinEcha < fEcal->MaxCrysInStin() )
    }
      else
    {
      std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
           << fStinName.Data() << " number in " << fStexName.Data() << " = " << n1StexStin << ". OUT OF BOUNDS"
           << " (max = " << fEcal->MaxStinEcnaInStex() << ")"
           << fTTBELL << std::endl;
      ret_code = kFALSE;
    }

      //.................................................................. (GetSampleAdcValues)
      //........ Filling of the 2D array of the event numbers in the data reading loop and 
      //         filling of the 3D array of the ADC sample values
      //
      //                           ONLY if ret_code == kTRUE

      if( ret_code == kTRUE )
    {
      //............ 1) Conversion (Stin,i0StinEcha) -> i0StexEcha  (same numbering for EB and EE)
      //==========================================================================================
      //   n1StexStin (Tower or SC):      1            2            3
      //   iStexStin                      0            1            2 
      //
      //   i0StinEcha:                 0......24    0......24    0......24
      //
      //   i0StexEcha         :        0......24   25......49   50......74  grouped by StexStin's
      //   i0StexEcha+1 (Xtal):        1......25   26......50   51......75       
      //
      //==========================================================================================

      Int_t i0StexEcha = i0StexStinEcna*fEcal->MaxCrysInStin() + i0StinEcha;
          
      //--------------------------------------------------------- (GetSampleAdcValues)
      if( i0StexEcha >= 0 && i0StexEcha < fEcal->MaxCrysEcnaInStex() )
        {
          //............ 2) Increase of the nb of evts for (StexEcha,sample) (events found in the data)
          (fT2d_NbOfEvts[i0StexEcha][i0Sample])++;     // value after first incrementation = 1
          fTagNbOfEvts[0] = 1;
          fFileHeader->fNbOfEvtsCalc = 1;
          
          //............ 3) Filling of the 3D array of the ADC values
          if ( i0EventIndex >= 0 && i0EventIndex < fFileHeader->fReqNbOfEvts )
        {  
          if( i0Sample >= 0 && i0Sample < fNbSampForFic )
            {
              fT3d_AdcValues[i0StexEcha][i0Sample][i0EventIndex] = adcvalue;
            }
          else
            {
              std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
               << " sample index = " << i0Sample << ". OUT OF BOUNDS"
               << " (max = " << fNbSampForFic << ")"
               << fTTBELL << std::endl;
            }
        }
          else
        {
          std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
               << " event number = " << n1EventNumber << ". OUT OF BOUNDS"
               << " (max = " << fFileHeader->fReqNbOfEvts << ")"
               << fTTBELL << std::endl;
          ret_code = kFALSE;
        }
        }
      else
        {
          std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> "
           << " CHANNEL NUMBER OUT OF BOUNDS" << std::endl
           << " i0StexEcha number = " << i0StexEcha
           << " , n1StexStin = " << n1StexStin
           << " , i0StinEcha = " << i0StinEcha
           << " , fEcal->MaxCrysEcnaInStex() = " << fEcal->MaxCrysEcnaInStex() 
           << fTTBELL << std::endl; 
          ret_code = kFALSE;
          // {Int_t cintoto; std::cout << "TAPER 0 POUR CONTINUER" << std::endl; cin >> cintoto;}
        }
    } // end of if( ret_code == kTRUE )
      else
    {
      std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> ret_code = kFALSE "
           << fTTBELL << std::endl;
    }
    } // end of if(fReadyToReadData == 1)
  else
    {
      std::cout << "!TEcnaRun::GetSampleAdcValues(...) *** ERROR ***> GetReadyToReadData(...) not called."
       << fTTBELL << std::endl;
      ret_code = kFALSE;
    }
  //.................................................................. (GetSampleAdcValues)
  if (ret_code == kFALSE)
    {
      std::cout << "!TEcnaRun::GetSampleAdcValues(...)> *** ERROR ***> ret_code = " << ret_code
       << " (FALSE). Event: " << n1EventNumber
       << ", " << fStexName.Data() << ": " << fFileHeader->fStex
       << ", " << fStinName.Data() << ": " << n1StexStin
       << ", channel: " << i0StinEcha
       << ", Sample: "   << i0Sample
       << ", ADC value: " << adcvalue << std::endl;
    } 
  return ret_code;
}
//------------- ( end of GetSampleAdcValues ) -----------------------
//====================================================================================================
//
//  ReadSampleAdcValues: Get the Sample ADC values from file by using TEcnaRead.
//
//====================================================================================================
Bool_t TEcnaRun::ReadSampleAdcValues(){return ReadSampleAdcValues(fEcal->MaxSampADC());}

Bool_t TEcnaRun::ReadSampleAdcValues(const Int_t& nb_samp_for_calc)
{
  // read the Sample ADC values from "ADC" result root files                     (ReadSampleAdcValues)

  // put the number of sample for calculations in attribute fNbSampForCalc
  // and call the method without arguments
  // We must have: nb_samp_for_calc <= fFileHeader->fNbOfSamples (= nb of samples in ROOT file)

  fNbSampForCalc = nb_samp_for_calc;

//   TEcnaRead* MyRootFile = new TEcnaRead(fFlagSubDet.Data(), fCnaParPaths, fCnaParCout,
//                     fFileHeader, fEcalNumbering, fCnaWrite);          //  fCnew++;
  
  TEcnaRead* MyRootFile = new TEcnaRead(fObjectManager, fFlagSubDet.Data());              //  fCnew++;

  MyRootFile->PrintNoComment();

  MyRootFile->FileParameters(fFileHeader->fTypAna, fFileHeader->fNbOfSamples, fFileHeader->fRunNumber,
                 fFileHeader->fFirstReqEvtNumber, fFileHeader->fLastReqEvtNumber,
                 fFileHeader->fReqNbOfEvts,       fFileHeader->fStex,
                 fCnaParPaths->ResultsRootFilePath().Data());

  Bool_t ok_read = MyRootFile->LookAtRootFile();

  fFileHeader->fStartTime = MyRootFile->GetStartTime();
  fFileHeader->fStopTime  = MyRootFile->GetStopTime();
  fFileHeader->fStartDate = MyRootFile->GetStartDate();
  fFileHeader->fStopDate  = MyRootFile->GetStopDate();
  
  if ( ok_read == kTRUE )
    {
      fRootFileName      = MyRootFile->GetRootFileName();
      fRootFileNameShort = MyRootFile->GetRootFileNameShort();
      std::cout << "*TEcnaRun::ReadSampleAdcValues> Reading sample ADC values from file: " << std::endl
       << "           " << fRootFileName << std::endl;
      
      size_t i_no_data = 0;

      //.......... Read the StinNumbers in the old file                     (ReadSampleAdcValues)
      TVectorD vec(fEcal->MaxStinEcnaInStex());
      for(Int_t i=0; i<fEcal->MaxStinEcnaInStex(); i++){vec(i)=(Double_t)0.;}
      vec = MyRootFile->ReadStinNumbers(fEcal->MaxStinEcnaInStex());
      if( MyRootFile->DataExist() == kTRUE )
    {
      fTagStinNumbers[0] = 1;
      fFileHeader->fStinNumbersCalc = 1;
      for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
        {fT1d_StexStinFromIndex[i0StexStinEcna] = (Int_t)vec(i0StexStinEcna);}
    }
      else
    {
      i_no_data++;
    }
      //.......... Read the Numbers of Events in the old file                      (ReadSampleAdcValues)
      TMatrixD partial_matrix(fEcal->MaxCrysInStin(), fFileHeader->fNbOfSamples);
      for(Int_t i=0; i<fEcal->MaxCrysInStin(); i++)
    {for(Int_t j=0; j<fFileHeader->fNbOfSamples; j++){partial_matrix(i,j)=(Double_t)0.;}}

      for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStin = MyRootFile->GetStexStinFromIndex(i0StexStinEcna);
      if(n1StexStin != -1)
        {
          partial_matrix =
        MyRootFile->ReadNumberOfEventsForSamples
        (n1StexStin, fEcal->MaxCrysInStin(), fFileHeader->fNbOfSamples);

          if( MyRootFile->DataExist() == kTRUE )
        {
          fTagNbOfEvts[0] = 1;
          fFileHeader->fNbOfEvtsCalc = 1;
          for(Int_t i0StinCrys=0; i0StinCrys<fEcal->MaxCrysInStin(); i0StinCrys++)
            {
              Int_t i0StexEcha = (n1StexStin-1)*fEcal->MaxCrysInStin() + i0StinCrys;
              for(Int_t i0Sample=0; i0Sample<fFileHeader->fNbOfSamples; i0Sample++)
            {fT2d_NbOfEvts[i0StexEcha][i0Sample] = (Int_t)partial_matrix(i0StinCrys,i0Sample);}
            }
        }
          else
        {
          i_no_data++;
        }
        }
    }

      //.......... Read the Sample ADC values in the old file                     (ReadSampleAdcValues)
      Double_t*** fT3d_read_AdcValues =
    MyRootFile->ReadSampleAdcValuesSameFile
    (fEcal->MaxCrysEcnaInStex(), fFileHeader->fNbOfSamples, fFileHeader->fReqNbOfEvts);
      
      if( MyRootFile->DataExist() == kTRUE )
    {             
      for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
        {
          for(Int_t i0Sample=0; i0Sample<fFileHeader->fNbOfSamples;i0Sample++)
        {      
          for(Int_t i_event=0; i_event<fFileHeader->fReqNbOfEvts; i_event++)
            {fT3d_AdcValues[i0StexEcha][i0Sample][i_event] = 
               fT3d_read_AdcValues[i0StexEcha][i0Sample][i_event];}
        }
        }
    }
      else
    {
      i_no_data++;
    }
      if(i_no_data)
    {
      std::cout << "!TEcnaRun::ReadSampleAdcValues(...)> *ERROR* =====> "
           << " Read failure. i_no_data = " << i_no_data << fTTBELL << std::endl;  
    }
    }
  else
    {
      std::cout << "!TEcnaRun::ReadSampleAdcValues(...)> *ERROR* =====> "
       << " ROOT file not found" << fTTBELL << std::endl;     
    }
  delete MyRootFile;        //  fCdelete++;
  return ok_read;
}
//------------- ( end of ReadSampleAdcValues ) -----------------------
//-------------------------------------------------------------------------
//
//    Get the ROOT file name (long and short)
//
//-------------------------------------------------------------------------
TString TEcnaRun::GetRootFileName(){return fRootFileName;}
TString TEcnaRun::GetRootFileNameShort(){return fRootFileNameShort;}

//###################################################################################################
//
// THE FOLLOWING METHODS ARE CALLED AFTER THE LOOPS OVER EVENTS, STINS, CRYSTALS AND SAMPLES
//
//###################################################################################################
//=========================================================================
//
//     Set start time, stop time, StartDate, StopDate
//
//=========================================================================
void TEcnaRun::StartStopTime(time_t t_startime, time_t t_stoptime)
{
// Put the start an stop time (if they exist) in fFileHeader class attributes.
  
  fFileHeader->fStartTime = t_startime;
  fFileHeader->fStopTime  = t_stoptime;
}

void TEcnaRun::StartStopDate(const TString& c_startdate, const TString& c_stopdate)
{
// Put the start an stop date (if they exist) in fFileHeader class attributes.

  fFileHeader->fStartDate = c_startdate;
  fFileHeader->fStopDate  = c_stopdate;
}

//=========================================================================
//
//                         GetReadyToCompute()   (technical)
//
//=========================================================================
void TEcnaRun::GetReadyToCompute()
{
//
  // MAKE THE RESULTS FILE NAME and
  // CHECK OF THE NUMBER OF FOUND EVENTS AND init fNumberOfEvents
  // (number used to compute the average values over the events)
  // The number of events fNumberOfEvents is extracted from the array fT2d_NbOfEvts[] 
  // which has been filled by the GetSampleAdcValues(...) method

  //..................... Making of the Root File name that will be written
  fCnaWrite->RegisterFileParameters(fFileHeader->fTypAna.Data(), fFileHeader->fNbOfSamples,
                    fFileHeader->fRunNumber,
                    fFileHeader->fFirstReqEvtNumber, fFileHeader->fLastReqEvtNumber,
                    fFileHeader->fReqNbOfEvts,       fFileHeader->fStex);
  
  fCnaWrite->fMakeResultsFileName();  // set fRootFileName, fRootFileNameShort

  //..................... Checking numbers of found events channel by channel
  if( fT2d_NbOfEvts != nullptr )
    {
      fNumberOfEvents = fCnaWrite->NumberOfEventsAnalysis(fT2d_NbOfEvts, fEcal->MaxCrysEcnaInStex(),
                              fNbSampForFic, fFileHeader->fReqNbOfEvts);
    }
  else
    {
      std::cout << "*TEcnaRun::GetReadyToCompute()> no data? fT2d_NbOfEvts = " << fT2d_NbOfEvts << std::endl;
    }
}  
//  end of GetReadyToCompute()

//-------------------------------------------------------------------
//
//                      SampleValues()      (technical)
//
//  Written in .root file corresponding to analysis name
//  beginning with: "Adc" (see EcnaAnalyzer.cc in package "Modules")
//
//-------------------------------------------------------------------
void TEcnaRun::SampleValues()
{
//3D histo of the sample ADC values for all the triples (StexEcha, sample, event)

  // The histo is already in fT3d_AdcValues[][][]
  // this method sets the "Tag", increment the "f...Calc" (and must be kept for that)
  // f...Calc > 0  => allow writing on file.
 
  if( fFileHeader->fAdcEvtCalc > 0 ){fFileHeader->fAdcEvtCalc = 0;}
  for( Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {fTagAdcEvt[i0StexEcha] = 1;        fFileHeader->fAdcEvtCalc++;}
}

//=========================================================================
//
//               C A L C U L A T I O N    M E T H O D S
//
//     fTag... = 1 => Calculation done. OK for writing on result file
//     ...Calc++   => Incrementation for result file size. 
//
//=========================================================================
void TEcnaRun::StandardCalculations()
{
  SampleMeans();
  SampleSigmas();
  CorrelationsBetweenSamples();
  
  Pedestals();                          // => mean over Xtal's
  TotalNoise();
  LowFrequencyNoise();
  HighFrequencyNoise();
  MeanCorrelationsBetweenSamples();
  SigmaOfCorrelationsBetweenSamples();
  
  AveragePedestals();                  // Average => mean over Stin's (Tower if EB, SC if EE)
  AverageTotalNoise();
  AverageLowFrequencyNoise();
  AverageHighFrequencyNoise();
  AverageMeanCorrelationsBetweenSamples();
  AverageSigmaOfCorrelationsBetweenSamples();
}

void TEcnaRun::Expert1Calculations()
{
  // long time, big file

  LowFrequencyCorrelationsBetweenChannels();
  HighFrequencyCorrelationsBetweenChannels();
}

void TEcnaRun::Expert2Calculations()
{
  // long time, no big file
  // expert 1 is called (if not called before) without writing in file.
  // Results are used only in memory to compute expert2 calculations

  LowFrequencyMeanCorrelationsBetweenStins();
  HighFrequencyMeanCorrelationsBetweenStins();
}
//====================================================================
//
//       E X P E C T A T I O N   V A L U E S  ,  V A R I A N C E S
// 
//====================================================================
//----------------------------------------------------------------
//  Calculation of the expectation values of the samples
//  for all the StexEchas
//
//  SMean(c,s)  = E_e[A(c,s,e*)]
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  E_e : average over the events
//----------------------------------------------------------------
void TEcnaRun::SampleMeans()
{
// Calculation of the expectation values over events
// for the samples 0 to fNbSampForCalc and for all the StexEchas

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::SampleMeans() " << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout  << "           Calculation: sample expectation values over the events"
      << " for each channel." << std::endl;}
  
  //................... Allocation fT2d_ev
  if ( fT2d_ev == nullptr ){
    Int_t n_samp = fNbSampForCalc;
    Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT2d_ev = new Double_t*[n_StexEcha];             fCnew++;  
    fT1d_ev = new  Double_t[n_StexEcha*n_samp];      fCnew++;   
    for(Int_t i = 0 ; i < n_StexEcha ; i++){
      fT2d_ev[i] = &fT1d_ev[0] + i*n_samp;}
  }
  //................... init fT2d_ev to zero
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
    {
      if( fT2d_ev[i0StexEcha][i0Sample] != (Double_t)0 )
        {fMiscDiag[1]++; fT2d_ev[i0StexEcha][i0Sample] = (Double_t)0;}
    } 
    }
  
  //................... Calculation
  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t i0Sample = 0 ; i0Sample < fNbSampForCalc ; i0Sample++)
    {
      for( Int_t i_event = 0; i_event < fNumberOfEvents; i_event++ )
        {
          fT2d_ev[i0StexEcha][i0Sample] += fT3d_AdcValues[i0StexEcha][i0Sample][i_event];
        }
      fT2d_ev[i0StexEcha][i0Sample] /= fNumberOfEvents;
    }
    }
  fTagMSp[0] = 1;        fFileHeader->fMSpCalc++;
}

//--------------------------------------------------------
//  Calculation of the sigmas of the samples
//  for all the StexEchas
//
//  SSigma(c,s) = sqrt{ Cov_e[A(c,s,e*),A(c,s,e*)] }
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  Cov_e : covariance over the events
//--------------------------------------------------------
void TEcnaRun::SampleSigmas() 
{
//Calculation of the sigmas of the samples for all the StexEchas
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::SampleSigmas()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation: sample ADC sigmas over the events "
     << " for each channel." << std::endl;}
  
  //... preliminary calculation of the expectation values if not done yet.
  //    The tag is set to 1 after call to the method. It is reset to 0
  //    because the expectation values must not be written in the result ROOT file
  //    (since the tag was equal to 0)
  if(fTagMSp[0] != 1){SampleMeans(); fTagMSp[0] = 0;}

  //................... Allocation fT2d_sig
  if( fT2d_sig == nullptr){
    Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    Int_t n_samp = fNbSampForCalc;
    fT2d_sig = new Double_t*[n_StexEcha];                fCnew++;        
    fT1d_sig = new  Double_t[n_StexEcha*n_samp];         fCnew++;  
    for(Int_t i0StexEcha = 0 ; i0StexEcha < n_StexEcha ; i0StexEcha++){
      fT2d_sig[i0StexEcha] = &fT1d_sig[0] + i0StexEcha*n_samp;}
  }
  // ................... init fT2d_sig to zero
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
    {
      if( fT2d_sig[i0StexEcha][i0Sample] != (Double_t)0 )
        {fMiscDiag[2]++; fT2d_sig[i0StexEcha][i0Sample] = (Double_t)0;}
    } 
    }
  
  //................... Calculation
  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t i0Sample = 0 ; i0Sample < fNbSampForCalc ; i0Sample++)
    {
      Double_t variance = (Double_t)0.;
      for( Int_t i_event = 0; i_event < fNumberOfEvents; i_event++ )
        {
          Double_t ecart = fT3d_AdcValues[i0StexEcha][i0Sample][i_event] - fT2d_ev[i0StexEcha][i0Sample];
          variance += ecart*ecart;
        }
      variance /= fNumberOfEvents;
      fT2d_sig[i0StexEcha][i0Sample] = sqrt(variance);
    }
    }
  fTagSSp[0] = 1;        fFileHeader->fSSpCalc++;
}

//====================================================================
//
//       C O V A R I A N C E S   &   C O R R E L A T I O N S
//
//                 B E T W E E N   S A M P L E S
//
//====================================================================
//-----------------------------------------------------------
//  Calculation of the covariances between samples
//  for all the StexEchas
//  Cov(c;s,s') = Cov_e[ A(c,s,e*) , A(c,s',e*) ]
//              = E_e[ ( A(c,s,e*) - E_e[A(c,s,e*)] )*
//                     ( A(c,s',e*) - E_e[A(c,s',e*)] ) ]
//  A(c,s,e)    : ADC value for channel c, sample s, event e
//  E_e , Cov_e : average, covariance over the events
//-----------------------------------------------------------
void TEcnaRun::CovariancesBetweenSamples()
{
  //Calculation of the covariances between samples for all the StexEchas
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::CovariancesBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation: covariances between samples"
         << " for each channel." << std::endl;}
  
  //................... Allocations cov_ss
  if( fT3d_cov_ss == nullptr ){
    const Int_t n_samp = fNbSampForCalc;
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT3d_cov_ss  = new Double_t**[n_StexEcha];                fCnew++;  
    fT3d2_cov_ss = new  Double_t*[n_StexEcha*n_samp];         fCnew++;  
    fT3d1_cov_ss = new   Double_t[n_StexEcha*n_samp*n_samp];  fCnew++;  
    for(Int_t i = 0 ; i < n_StexEcha ; i++){
      fT3d_cov_ss[i] = &fT3d2_cov_ss[0] + i*n_samp;
      for(Int_t j = 0 ; j < n_samp ; j++){
    fT3d2_cov_ss[n_samp*i+j] = &fT3d1_cov_ss[0]+n_samp*(n_samp*i+j);}}
  }
  
  //.................. Calculation (= init)
  //.................. computation of half of the matrix, diagonal included)

  //... preliminary calculation of the expectation values if not done yet.
  //    The tag is set to 1 after call to the method. It is reset to 0
  //    because the expectation values must not be written in the result ROOT file
  //    (since the tag was equal to 0)
  //    Results in array  fT2d_ev[j0StexEcha][i0Sample]
  if(fTagMSp[0] != 1){SampleMeans(); fTagMSp[0] = 0;}


  for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {
      for (Int_t i0Sample = 0 ; i0Sample < fNbSampForCalc ; i0Sample++)
    {
      for (Int_t j0Sample = 0 ; j0Sample <= i0Sample; j0Sample++)
        {
          fT3d_cov_ss[j0StexEcha][i0Sample][j0Sample] = (Double_t)0;
          for( Int_t i_event = 0; i_event < fNumberOfEvents; i_event++ )
        {
          fT3d_cov_ss[j0StexEcha][i0Sample][j0Sample] +=
            (fT3d_AdcValues[j0StexEcha][i0Sample][i_event] - fT2d_ev[j0StexEcha][i0Sample])
            *(fT3d_AdcValues[j0StexEcha][j0Sample][i_event] - fT2d_ev[j0StexEcha][j0Sample]);
        }
          fT3d_cov_ss[j0StexEcha][i0Sample][j0Sample] /= (Double_t)fNumberOfEvents;
          fT3d_cov_ss[j0StexEcha][j0Sample][i0Sample] = fT3d_cov_ss[j0StexEcha][i0Sample][j0Sample];
        }
    }
      fTagCovCss[j0StexEcha] = 1;     fFileHeader->fCovCssCalc++;
    }
}

//-----------------------------------------------------------
//
//  Calculation of the correlations between samples
//  for all the StexEchas
//  Cor(c;s,s') = Cov(c;s,s')/sqrt{ Cov(c;s,s)*Cov(c;s',s') }
//-----------------------------------------------------------
void TEcnaRun::CorrelationsBetweenSamples()
{
//Calculation of the correlations between samples for all the StexEchas

  //... preliminary calculation of the covariances if not done yet.
  //    Test only the first tag since the cov are computed globaly
  //    but set all the tags to 0 because we don't want to write
  //    the covariances in the result ROOT file    
  if ( fTagCovCss[0] != 1 ){CovariancesBetweenSamples();
  for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {fTagCovCss[j0StexEcha] = 0;}}
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::CorrelationsBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation: correlations between samples"
     << " for each channel." << std::endl;}

  //................... Allocations cor_ss
  if( fT3d_cor_ss == nullptr){
    const Int_t n_samp = fNbSampForCalc;
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT3d_cor_ss  = new Double_t**[n_StexEcha];                fCnew++;  
    fT3d2_cor_ss = new  Double_t*[n_StexEcha*n_samp];         fCnew++;  
    fT3d1_cor_ss = new   Double_t[n_StexEcha*n_samp*n_samp];  fCnew++;  
    for(Int_t i = 0 ; i < n_StexEcha ; i++){
      fT3d_cor_ss[i] = &fT3d2_cor_ss[0] + i*n_samp;
      for(Int_t j = 0 ; j < n_samp ; j++){
    fT3d2_cor_ss[n_samp*i+j] = &fT3d1_cor_ss[0]+n_samp*(n_samp*i+j);}}
  }
  
  //..................... calculation of the correlations (=init)
  //......................computation of half of the matrix, diagonal included (verif = 1)
    
 for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
   {
     for (Int_t i0Sample = 0 ; i0Sample < fNbSampForCalc ; i0Sample++)
       {
     for (Int_t j0Sample = 0 ; j0Sample <= i0Sample ; j0Sample++)
       {
         if( (fT3d_cov_ss[j0StexEcha][i0Sample][i0Sample] > 0) &&
         (fT3d_cov_ss[j0StexEcha][j0Sample][j0Sample] > 0) )
           {
         fT3d_cor_ss[j0StexEcha][i0Sample][j0Sample] = fT3d_cov_ss[j0StexEcha][i0Sample][j0Sample]/
           ( sqrt(fT3d_cov_ss[j0StexEcha][i0Sample][i0Sample])*sqrt(fT3d_cov_ss[j0StexEcha][j0Sample][j0Sample]) );
           }
         else
           {
         (fT3d_cor_ss)[j0StexEcha][i0Sample][j0Sample] = (Double_t)0; // prevoir compteur + fTTBELL
           }
         fT3d_cor_ss[j0StexEcha][j0Sample][i0Sample] = fT3d_cor_ss[j0StexEcha][i0Sample][j0Sample];
       }
       }
     fTagCorCss[j0StexEcha] = 1;          fFileHeader->fCorCssCalc++;
   }
}

//===========================================================================
//
//     M E A N   P E D E S T A L S ,   T O T A L    N O I S E ,
//     L O W     F R E Q U E N C Y    N O I S E , 
//     H I G H   F R E Q U E N C Y    N O I S E 
//     M E A N   O F  C O R ( S , S ),  S I G M A   O F  C O R ( S , S )
//
//===========================================================================
//-------------------------------------------------------------------------
//
//  Calculation of the Pedestals for each channel in Stex
//  tag: Ped
//  Pedestal(c ) = E_e[ E_s[A(c ,s*,e*)] ]
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  E_e : average over the events
//  E_s : average over the samples
//      
//-------------------------------------------------------------------------
void TEcnaRun::Pedestals()
{
// Calculation, for each channel, of the expectation values
// (over the samples 0 to fNbSampForCalc-1) of the ADC expectation values 
// (over the events)

  //... preliminary calculation of the expectation values if not done yet
  if ( fTagMSp[0] != 1 ){SampleMeans(); fTagMSp[0]=0;}

  //................... Allocation ev_ev + init to zero (mandatory)
  if( fT1d_ev_ev == nullptr ){fT1d_ev_ev = new Double_t[fEcal->MaxCrysEcnaInStex()]; fCnew++;}
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {if( fT1d_ev_ev[i0StexEcha] != (Double_t)0 )
      {fMiscDiag[11]++; fT1d_ev_ev[i0StexEcha] = (Double_t)0;}}
  
  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::Pedestals()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the channels, of the expectation values (over the samples 1 to "
     << fNbSampForCalc << ")" << std::endl
     << "          of the ADC expectation values (over the events)." << std::endl;}
  
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {     
      for (Int_t i0Sample = 0; i0Sample < fNbSampForCalc; i0Sample++)
    {
      fT1d_ev_ev[i0StexEcha] += fT2d_ev[i0StexEcha][i0Sample];
    }
      fT1d_ev_ev[i0StexEcha] /= fNbSampForCalc;
    }
  fTagPed[0] = 1;                      fFileHeader->fPedCalc++;
}
//------------------------ (end of Pedestals) ----------------------------

//------------------------------------------------------------------------------
//
//  Calculation of the TN (Total Noise)
//  tag: Tno
//    
//  TotalNoise(c)  = E_s[ sqrt{ E_e[ ( A(c ,s*,e*) - E_e[A(c ,s*,e*)] )^2 ] } ]
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  E_e : average over the events
//  E_s : average over the samples
//
//------------------------------------------------------------------------------
void TEcnaRun::TotalNoise()
{
// Calculation, for each channel, of the expectation values
// (over the samples 0 to fNbSampForCalc-1) of the sigmas 
// (over the events)

  //... preliminary calculation of the sigmas if not done yet
  if ( fTagSSp[0] != 1 ){SampleSigmas(); fTagSSp[0]=0;}

  //................... Allocation ev_ev + init to zero (mandatory)
  if( fT1d_evsamp_of_sigevt == nullptr ){fT1d_evsamp_of_sigevt = new Double_t[fEcal->MaxCrysEcnaInStex()]; fCnew++;}
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {if( fT1d_evsamp_of_sigevt[i0StexEcha] != (Double_t)0 )
      {fMiscDiag[12]++; fT1d_evsamp_of_sigevt[i0StexEcha] = (Double_t)0;}}
  
  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::TotalNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the channels, of the expectation values (over the samples 1 to "
     << fNbSampForCalc << ")" << std::endl
     << "          of the ADC expectation values (over the events)." << std::endl;}

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {     
      for (Int_t i0Sample = 0; i0Sample < fNbSampForCalc; i0Sample++)
    {
      if( fT2d_sig[i0StexEcha][i0Sample] < 0)
        {
          std::cout << "!TEcnaRun::TotalNoise() *** ERROR ***> Negative sigma!"
           << fTTBELL << std::endl;
        }
      else
        {
          fT1d_evsamp_of_sigevt[i0StexEcha] += fT2d_sig[i0StexEcha][i0Sample];
        }
    }
      fT1d_evsamp_of_sigevt[i0StexEcha] /= fNbSampForCalc;
    }
  fTagTno[0] = 1;                     fFileHeader->fTnoCalc++;
}
//------------------------ (end of TotalNoise) ----------------------------

//---------------------------------------------------------------------------------
//
//  Calculation of the LFN  (Low Frequency Noise)
//  tag: Lfn
//
//  LowFqNoise(c) = sqrt{ E_e[ ( E_s[A(c ,s*,e*)] - E_e[ E_s[A(c ,s*,e*)] ] )^2 ] }
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  E_e : average over the events
//  E_s : average over the samples
//
//---------------------------------------------------------------------------------
void TEcnaRun::LowFrequencyNoise()
{
// Calculation, for each channel, of the sigma (over the events)
// of the ADC expectation values (over the samples 0 to fNbSampForCalc-1)

  //................... Allocation fT1d_sigevt_of_evsamp + init to zero (mandatory)
  if( fT1d_sigevt_of_evsamp == nullptr ){
    fT1d_sigevt_of_evsamp = new Double_t[fEcal->MaxCrysEcnaInStex()];              fCnew++;  
  }
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( fT1d_sigevt_of_evsamp[i0StexEcha] != (Double_t)0 )
    {fMiscDiag[13]++; fT1d_sigevt_of_evsamp[i0StexEcha] = (Double_t)0;}
    }

  //................... Allocation mean_over_samples
  TVectorD  mean_over_samples(fNumberOfEvents);
  for(Int_t i=0; i<fNumberOfEvents; i++){mean_over_samples(i)=(Double_t)0.;}

  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::LowFrequencyNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for each channel, of the sigma (over the events)" << std::endl
     << "           of the ADC expectation values (over the samples 1 to "
     <<  fNbSampForCalc << ")." << std::endl;}

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      // Calculation of the mean over the events of the mean over the samples
      Double_t mean_over_events = (Double_t)0;
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
    {
      // Calculation, for each event, of the mean over the samples  
      mean_over_samples(n_event) = (Double_t)0.;
      for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
        {
          mean_over_samples(n_event) += fT3d_AdcValues[i0StexEcha][i0Sample][n_event];
        }
      mean_over_samples(n_event) /= (Double_t)fNbSampForCalc;

      mean_over_events += mean_over_samples(n_event);
    }
       mean_over_events /= (Double_t)fNumberOfEvents;
       
       // Calculation of the sigma over the events of the mean over the samples
       Double_t var = (Double_t)0;
       for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
     {
      Double_t ecart =  mean_over_samples(n_event) - mean_over_events;
      var += ecart*ecart;
     }
       var /= (Double_t)fNumberOfEvents;
       
       fT1d_sigevt_of_evsamp[i0StexEcha] = sqrt(var);
    }
  fTagLfn[0] = 1;                      fFileHeader->fLfnCalc++;
}
//------------------------ (end of LowFrequencyNoise) ----------------------------

//---------------------------------------------------------------------------------
//
//  Calculation of the HFN  (High Frequency Noise)
//  tag: Hfn
//
//  HighFqNoise(c) = E_e[ sqrt{ E_s[ (A(c ,s*,e*) - E_s[A(c ,s*,e*)] )^2 ] } ]
//  A(c,s,e) : ADC value for channel c, sample s, event e
//  E_e : average over the events
//  E_s : average over the samples
//
//---------------------------------------------------------------------------------
void TEcnaRun::HighFrequencyNoise()
{
// Calculation, for each channel, of the mean (over the events)
// of the ADC sigmas (over the samples 0 to fNbSampForCalc-1)

  //................... Allocation fT1d_evevt_of_sigsamp + init to zero (mandatory)
  if( fT1d_evevt_of_sigsamp == nullptr ){
    fT1d_evevt_of_sigsamp = new Double_t[fEcal->MaxCrysEcnaInStex()];              fCnew++;  
  }
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( fT1d_evevt_of_sigsamp[i0StexEcha] != (Double_t)0 )
    {fMiscDiag[14]++; fT1d_evevt_of_sigsamp[i0StexEcha] = (Double_t)0;}
    }

  //................... Allocations mean_over_samples, sigma_over_sample
  TVectorD  mean_over_samples(fNumberOfEvents);
  for(Int_t i=0; i<fNumberOfEvents; i++){mean_over_samples(i)=(Double_t)0.;} 
  TVectorD  sigma_over_samples(fNumberOfEvents);
  for(Int_t i=0; i<fNumberOfEvents; i++){sigma_over_samples(i)=(Double_t)0.;}

  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::HighFrequencyNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for each channel, of the sigma (over the events)" << std::endl
     << "           of the ADC expectation values (over the samples 1 to "
     <<  fNbSampForCalc << ")." << std::endl;}

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      //..................... Calculation of the sigma over samples
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
    {
      // Calculation, for each event, of the mean over the samples  
      mean_over_samples(n_event) = (Double_t)0.;
      for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
        {mean_over_samples(n_event) += fT3d_AdcValues[i0StexEcha][i0Sample][n_event];}
      mean_over_samples(n_event) /= (Double_t)fNbSampForCalc;
      
      // Calculation, for each event, of the sigma over the samples
      Double_t var_over_samples =  (Double_t)0;
      for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
        {Double_t deviation = fT3d_AdcValues[i0StexEcha][i0Sample][n_event] - mean_over_samples(n_event);
        var_over_samples += deviation*deviation;}
      var_over_samples /= (Double_t)fNbSampForCalc;
      
      if( var_over_samples < 0)
        {std::cout << "!TEcnaRun::HighFrequencyNoise() *** ERROR ***> Negative variance! " << fTTBELL << std::endl;}
      else
        {sigma_over_samples(n_event) = sqrt(var_over_samples);}
    }

      //....... Calculation of the mean over the events of the sigma over samples
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
    {fT1d_evevt_of_sigsamp[i0StexEcha] += sigma_over_samples(n_event);}

      fT1d_evevt_of_sigsamp[i0StexEcha] /= (Double_t)fNumberOfEvents;
    }
  fTagHfn[0] = 1;                      fFileHeader->fHfnCalc++;
}
//------------------------ (end of HighFrequencyNoise) ----------------------------

//-------------------------------------------------------------------------
//
//  Calculation of the expectation values of (sample,sample)
//  correlations for all the channels (mean cor(s,s))
//  tag: MeanCorss
//
//  MeanCorss(c)   = E_s,s'[ Cor(c;s,s') ]
//  E_s,s': average  over couples of samples (half correlation matrix)
//
//-------------------------------------------------------------------------
void TEcnaRun::MeanCorrelationsBetweenSamples()
{
  // Calculation, for all the channels, of the expectation values
  // of the correlations between the first fNbSampForCalc samples
  
  //... preliminary calculation of the correlations if not done yet
  //    (test only the first element since the cor are computed globaly)
  if ( fTagCorCss[0] != 1 ){CorrelationsBetweenSamples(); fTagCorCss[0]=0;}

  //................... Allocations ev_cor_ss + init to zero (mandatory)
  if( fT1d_ev_cor_ss == nullptr ){
    Int_t n_StexEcha =  fEcal->MaxCrysEcnaInStex();
    fT1d_ev_cor_ss = new Double_t[n_StexEcha];               fCnew++;  
  }
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( fT1d_ev_cor_ss[i0StexEcha] != (Double_t)0 )
    {fMiscDiag[15]++; fT1d_ev_cor_ss[i0StexEcha] = (Double_t)0;}
    }
  
  //.......... 1D array half_cor_ss[N(N-1)/2] to put the N (sample,sample) correlations
  //           ( half of (them minus the diagonal) ) 
  Int_t ndim = (Int_t)(fNbSampForCalc*(fNbSampForCalc - 1)/2);

  TVectorD  half_cor_ss(ndim); for(Int_t i=0; i<ndim; i++){half_cor_ss(i)=(Double_t)0.;}

  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::MeanCorrelationsBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the channels, of the expectation values of the" << std::endl
     << "           correlations between the first " << fNbSampForCalc << " samples." << std::endl;}

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      //..................... half_cor_ss() array filling
      Int_t i_count = 0;
      for (Int_t i0Sample = 0; i0Sample < fNbSampForCalc; i0Sample++)
    {
      for (Int_t j0Sample = 0; j0Sample < i0Sample; j0Sample++)
        {
          half_cor_ss(i_count) = fT3d_cor_ss[i0StexEcha][i0Sample][j0Sample];
          i_count++;
        }
    }
      //...................... mean cor(s,s') calculation
      fT1d_ev_cor_ss[i0StexEcha] = (Double_t)0;
      for(Int_t i_rcor = 0; i_rcor < ndim; i_rcor++)
    {
      fT1d_ev_cor_ss[i0StexEcha] += half_cor_ss(i_rcor);
    }
      fT1d_ev_cor_ss[i0StexEcha] /= (Double_t)ndim;
    }
  fTagMeanCorss[0] = 1;               fFileHeader->fMeanCorssCalc++;
}
//--------------- (end of MeanCorrelationsBetweenSamples) -----------

//-------------------------------------------------------------------------
//
// Calculation of the sigmas of the (sample,sample) correlations
// for all the channels (sigma of cor(s,s))
// tag: SigCorss
//
// SigmaCorss(c)  = E_s,s'[ Cor(c;s,s') - E_s,s'[ Cor(c;s,s') ] ]
// E_s,s': average  over couples of samples (half correlation matrix)
//
//--------------------------------------------------------------------------
void  TEcnaRun::SigmaOfCorrelationsBetweenSamples()
{
  //Calculation of the sigmas of the (sample,sample) correlations for all the StexEchas
 
  //... preliminary calculation of the mean cor(s,s') if not done yet
  //    (test only the first element since the cor are computed globaly)
  //    Results available in array fT1d_ev_cor_ss[i0StexEcha]
  if ( fTagMeanCorss[0] != 1 ){MeanCorrelationsBetweenSamples(); fTagMeanCorss[0]=0;}

  //................... Allocations sig_cor_ss + init to zero
  if( fT1d_sig_cor_ss == nullptr ){
    Int_t n_StexEcha =  fEcal->MaxCrysEcnaInStex();
    fT1d_sig_cor_ss = new Double_t[n_StexEcha];                fCnew++;  
  }
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( fT1d_sig_cor_ss[i0StexEcha] != (Double_t)0 )
    {fMiscDiag[16]++; fT1d_sig_cor_ss[i0StexEcha] = (Double_t)0;}
    }

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::SigmasOfCorrelationsBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the sigmas of the (sample,sample)" << std::endl
     << "           correlations for all the channels." << std::endl;}

  //.......... 1D array half_cor_ss[N(N-1)/2] to put the N (sample,sample) correlations
  //           (half of them minus the diagonal)
  Int_t ndim = (Int_t)(fNbSampForCalc*(fNbSampForCalc - 1)/2);

  TVectorD  half_cor_ss(ndim); for(Int_t i=0; i<ndim; i++){half_cor_ss(i)=(Double_t)0.;}

  //.................. Calculation
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      //..................... half_cor_ss() array filling
      Int_t i_count = 0;
      for (Int_t i0Sample = 0; i0Sample < fNbSampForCalc; i0Sample++)
    {
      for (Int_t j0Sample = 0; j0Sample < i0Sample; j0Sample++)
        {
          half_cor_ss(i_count) = fT3d_cor_ss[i0StexEcha][i0Sample][j0Sample];
          i_count++;
        }
    }

      //...................... sigma of cor(s,s') calculation
      Double_t var = (Double_t)0;
      for(Int_t i_rcor = 0; i_rcor < ndim; i_rcor++)
    {
      Double_t ecart = half_cor_ss(i_rcor) - fT1d_ev_cor_ss[i0StexEcha];
      var += ecart*ecart;
    }
      var /= (Double_t)ndim;
      fT1d_sig_cor_ss[i0StexEcha] = sqrt(var);
    }
  fTagSigCorss[0] = 1;              fFileHeader->fSigCorssCalc++;
}
//--------------- (end of SigmaOfCorrelationsBetweenSamples) -----------

//-----------------------------------------------------------------------------
//
//  Calculation of the average Pedestals for each Stin in Stex
//  tag: AvPed
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AveragePedestals()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the Pedestals

  //... preliminary calculation of the Pedestals if not done yet
  if ( fTagPed[0] != 1 ){Pedestals(); fTagPed[0]=0;}
  //................... Allocation av_mped + init to zero (mandatory)
  if( fT1d_av_mped == nullptr ){fT1d_av_mped = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_mped[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[41]++; fT1d_av_mped[i0StexStinEcna] = (Double_t)0;}}
 
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AveragePedestals()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average Pedestals" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_mped[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);

      if( fStexName == "SM " )
        {fT1d_av_mped[i0StexStinEcna] += fT1d_ev_ev[i0StexEcha];}

      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 && n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 && n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_mped[i0StexStinEcna] += fT1d_ev_ev[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}
      
      fT1d_av_mped[i0StexStinEcna] = fT1d_av_mped[i0StexStinEcna]/xdivis;
    }
  
  fTagAvPed[0] = 1;                      fFileHeader->fAvPedCalc++;
}
//-----------------------------------------------------------------------------
//
// Calculation of the average total noise for each Stin in Stex
// tag: AvTno
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AverageTotalNoise()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the Total Noise
  
  //... preliminary calculation of the averaged Total Noise if not done yet
  if ( fTagTno[0] != 1 ){TotalNoise(); fTagTno[0]=0;}
  //................... Allocation av_totn + init to zero (mandatory)
  if( fT1d_av_totn == nullptr ){fT1d_av_totn = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_totn[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[42]++; fT1d_av_totn[i0StexStinEcna] = (Double_t)0;}}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AverageTotalNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average total Noise" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_totn[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);
      
      if( fStexName == "SM " )
        {fT1d_av_totn[i0StexStinEcna] += fT1d_evsamp_of_sigevt[i0StexEcha];}
      
      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 && n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 && n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_totn[i0StexStinEcna] += fT1d_evsamp_of_sigevt[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}
      
      fT1d_av_totn[i0StexStinEcna] = fT1d_av_totn[i0StexStinEcna]/xdivis;
    }
  fTagAvTno[0] = 1;                      fFileHeader->fAvTnoCalc++;
}
//-----------------------------------------------------------------------------
//
// Calculation of the average Low Frequency noise for each Stin in Stex
// tag: AvLfn
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AverageLowFrequencyNoise()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the Low Frequency Noise
  
  //... preliminary calculation of the Low Frequency Noise if not done yet
  if ( fTagLfn[0] != 1 ){LowFrequencyNoise(); fTagLfn[0]=0;}
  //................... Allocation av_lofn + init to zero (mandatory)
  if( fT1d_av_lofn == nullptr ){fT1d_av_lofn = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_lofn[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[43]++; fT1d_av_lofn[i0StexStinEcna] = (Double_t)0;}}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AverageLowFrequencyNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average Low Frequency Noise" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_lofn[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);

      if( fStexName == "SM " )
        {fT1d_av_lofn[i0StexStinEcna] += fT1d_sigevt_of_evsamp[i0StexEcha];}
      
      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_lofn[i0StexStinEcna] += fT1d_sigevt_of_evsamp[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}
      
      fT1d_av_lofn[i0StexStinEcna] = fT1d_av_lofn[i0StexStinEcna]/xdivis;
    }
  fTagAvLfn[0] = 1;                      fFileHeader->fAvLfnCalc++; 
}
//-----------------------------------------------------------------------------
//
// Calculation of the average high frequency noise for each Stin in Stex
// tag: AvHfn
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AverageHighFrequencyNoise()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the High Frequency Noise
  
  //... preliminary calculation of the High Frequency Noise if not done yet
  if ( fTagHfn[0] != 1 ){HighFrequencyNoise(); fTagHfn[0]=0;}
  //................... Allocation av_hifn + init to zero (mandatory)
  if( fT1d_av_hifn == nullptr ){fT1d_av_hifn = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_hifn[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[44]++; fT1d_av_hifn[i0StexStinEcna] = (Double_t)0;}}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AverageHighFrequencyNoise()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average High Frequency Noise" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {   
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_hifn[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);

      if( fStexName == "SM " )
        {fT1d_av_hifn[i0StexStinEcna] += fT1d_evevt_of_sigsamp[i0StexEcha];}
      
      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_hifn[i0StexStinEcna] += fT1d_evevt_of_sigsamp[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}

      fT1d_av_hifn[i0StexStinEcna] = fT1d_av_hifn[i0StexStinEcna]/xdivis;
    }
  fTagAvHfn[0] = 1;                      fFileHeader->fAvHfnCalc++;  
}
//-----------------------------------------------------------------------------
//
// Calculation of the average mean cor(s,s) for each Stin in Stex
// tag: AvMeanCorss
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AverageMeanCorrelationsBetweenSamples()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the mean cor(s,s)
  
  //... preliminary calculation of the mean cor(s,s) if not done yet
  if ( fTagMeanCorss[0] != 1 ){MeanCorrelationsBetweenSamples(); fTagMeanCorss[0]=0;}
  //................... Allocation av_ev_corss + init to zero (mandatory)
  if( fT1d_av_ev_corss == nullptr ){fT1d_av_ev_corss = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_ev_corss[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[45]++; fT1d_av_ev_corss[i0StexStinEcna] = (Double_t)0;}}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AverageMeanCorrelationsBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average mean cor(s,s)" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_ev_corss[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);

      if( fStexName == "SM " )
        {fT1d_av_ev_corss[i0StexStinEcna] += fT1d_ev_cor_ss[i0StexEcha];}
      
      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_ev_corss[i0StexStinEcna] += fT1d_ev_cor_ss[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}

      fT1d_av_ev_corss[i0StexStinEcna] = fT1d_av_ev_corss[i0StexStinEcna]/xdivis;
    }
  fTagAvMeanCorss[0] = 1;                      fFileHeader->fAvMeanCorssCalc++;    
}
//-----------------------------------------------------------------------------
//
// Calculation of the average sigma of cor(s,s) for each Stin in Stex
// tag: AvSigCorss
//      
//-----------------------------------------------------------------------------
void TEcnaRun::AverageSigmaOfCorrelationsBetweenSamples()
{
// Calculation of the average 
// (over the Stin's 0 to fEcal->MaxStinInStex()) of the sigma of cor(s,s)
  
  //... preliminary calculation of the sigma of cor(s,s) if not done yet
  if ( fTagSigCorss[0] != 1 ){SigmaOfCorrelationsBetweenSamples(); fTagSigCorss[0]=0;}
  //................... Allocation av_sig_corss + init to zero (mandatory)
  if( fT1d_av_sig_corss == nullptr ){fT1d_av_sig_corss = new Double_t[fEcal->MaxStinEcnaInStex()]; fCnew++;}
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {if( fT1d_av_sig_corss[i0StexStinEcna] != (Double_t)0 )
      {fMiscDiag[46]++; fT1d_av_sig_corss[i0StexStinEcna] = (Double_t)0;}}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::AverageSigmaOfCorrelationsBetweenSamples()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for all the "
     << fStinName.Data() << "s, of the average sigma of cor(s,s)" << std::endl;}

  //................... Calculation
  for(Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      Int_t n1StexStinEcna = i0StexStinEcna+1;
      fT1d_av_sig_corss[i0StexStinEcna] = (Double_t)0;
      for(Int_t i0StinEcha = 0; i0StinEcha < fEcal->MaxCrysInStin(); i0StinEcha++)
    {
      Int_t i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(n1StexStinEcna, i0StinEcha);

      if( fStexName == "SM " )
        {fT1d_av_sig_corss[i0StexStinEcna] += fT1d_sig_cor_ss[i0StexEcha];}
      
      if( fStexName == "Dee" )
        {
          //---------------- Special translation for mixed SCEcna (29 and 32)
          //                 Xtal 11 of SCEcna 29 -> Xtal 11 of SCEcna 10
          //                 Xtal 11 of SCEcna 32 -> Xtal 11 of SCEcna 11
          Int_t n1StinEcha = i0StinEcha+1;
          if( n1StexStinEcna == 10 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(29, i0StinEcha);}
          if( n1StexStinEcna == 11 &&  n1StinEcha == 11 )
        {i0StexEcha = fEcalNumbering->Get0StexEchaFrom1StexStinAnd0StinEcha(32, i0StinEcha);}
          if( !( (n1StexStinEcna == 29 || n1StexStinEcna == 32) && n1StinEcha == 11 )  )
        {fT1d_av_sig_corss[i0StexStinEcna] += fT1d_sig_cor_ss[i0StexEcha];}
        }
    }
      Double_t xdivis = (Double_t)0.;
      if( fStexName == "SM "  )
    {xdivis = (Double_t)fEcal->MaxCrysInStin();}
      if( fStexName == "Dee" )
    {xdivis = (Double_t)fEcalNumbering->MaxCrysInStinEcna(fFileHeader->fStex, n1StexStinEcna, "TEcnaRun");}

      fT1d_av_sig_corss[i0StexStinEcna] = fT1d_av_sig_corss[i0StexStinEcna]/xdivis;
    }
  fTagAvSigCorss[0] = 1;                      fFileHeader->fAvSigCorssCalc++;       
}

//======================================================================
//
//       C O V A R I A N C E S   &   C O R R E L A T I O N S
//
//                 B E T W E E N   C H A N N E L S
//
//======================================================================
//----------------------------------------------------------------------
//
//  Calculation of the Low Frequency Covariances between channels
//
//  LFCov(Ci,Cj) = Cov_e[ E_s[A(Ci,s*,e*)] , E_s[A(Cj,s*,e*) ]
//
//               = E_e[ ( E_s[A(Ci,s*,e*)] - E_e[ E_s[A(Ci,s*,e*)] ] )*
//                      ( E_s[A(Cj,s*,e*)] - E_e[ E_s[A(Cj,s*,e*)] ] ) ]
//
//   A(Ci,s,e) : ADC value for channel Ci, sample s, event e
//
//   E_e , Cov_e : average, covariance over the events
//   E_s :         average over the samples
//
//   e* : random variable associated to events
//   s* : random variable associated to samples
//
//----------------------------------------------------------------------
void TEcnaRun::LowFrequencyCovariancesBetweenChannels()
{
//Calculation of the Low Frequency Covariances between channels
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::LowFrequencyCovariancesBetweenChannels()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the Low Frequency Covariances between channels" << std::endl;}

  //................. allocation fT2d_lf_cov + init to zero (mandatory)
  if( fT2d_lf_cov == nullptr ){
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT2d_lf_cov  = new Double_t*[n_StexEcha];                   fCnew++;
    fT2d1_lf_cov = new  Double_t[n_StexEcha*n_StexEcha];          fCnew++;
    for(Int_t i0StexEcha = 0 ; i0StexEcha < n_StexEcha ; i0StexEcha++){
    fT2d_lf_cov[i0StexEcha] = &fT2d1_lf_cov[0] + i0StexEcha*n_StexEcha;}}

  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {
      if( fT2d_lf_cov[i0StexEcha][j0StexEcha] != (Double_t)0 )
        {fMiscDiag[21]++; fT2d_lf_cov[i0StexEcha][j0StexEcha] = (Double_t)0;}
    }
    }
  //........................................... Calculation  (LowFrequencyCovariancesBetweenChannels)
  //................... Allocation mean_over_samples(i0StexEcha, n_event)
  TMatrixD mean_over_samples(fEcal->MaxCrysEcnaInStex(), fNumberOfEvents);
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {for(Int_t n_event=0; n_event<fNumberOfEvents; n_event++)
    {mean_over_samples(i0StexEcha, n_event)=(Double_t)0.;}}
  //................... Allocation MoeOfMos(i0StexEcha)
  TVectorD MoeOfMos(fEcal->MaxCrysEcnaInStex());
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    { MoeOfMos(i0StexEcha)=(Double_t)0.;}

  //................... Calculation  
  if(fFlagPrint != fCodePrintNoComment){
    std::cout << "          Calculation, for each pair of channels, of the covariance (over the events)" << std::endl
     << "          between the ADC expectation values (over the samples 1 to "
     <<  fNbSampForCalc << ")." << std::endl;}

  std::cout << " Please, wait (end at i= " << fEcal->MaxCrysEcnaInStex() << "): " << std::endl;

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      MoeOfMos(i0StexEcha) = (Double_t)0;

      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
        {
          // Calculation, for each event, of the mean over the samples  ( = E_s[A(c_i,s*,e_n] )
          mean_over_samples(i0StexEcha, n_event) = (Double_t)0.;
          for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
        {
          mean_over_samples(i0StexEcha, n_event) += fT3d_AdcValues[i0StexEcha][i0Sample][n_event];
        }
          mean_over_samples(i0StexEcha, n_event) /= (Double_t)fNbSampForCalc;
        }
      //Calculation of the mean over the events of E_s[A(c_i,s*,e_n] ( = E_e[E_s[A(c_i,s*,e*]] )
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
        {MoeOfMos(i0StexEcha) += mean_over_samples(i0StexEcha, n_event);}
      MoeOfMos(i0StexEcha) /= (Double_t)fNumberOfEvents;
    }
    }

  //... Calculation of half of the matrix, diagonal included        (LowFrequencyCovariancesBetweenChannels)
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for(Int_t j0StexEcha = 0; j0StexEcha <= i0StexEcha; j0StexEcha++)
        {
          if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, j0StexEcha) > 0)
          || ( fFlagSubDet == "EB") )
        {
          fT2d_lf_cov[i0StexEcha][j0StexEcha] = (Double_t)0;
          for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
            {
              fT2d_lf_cov[i0StexEcha][j0StexEcha] +=
            ( mean_over_samples(i0StexEcha, n_event) - MoeOfMos(i0StexEcha) )
            *( mean_over_samples(j0StexEcha, n_event) - MoeOfMos(j0StexEcha) );    
            }
          fT2d_lf_cov[i0StexEcha][j0StexEcha] /= (Double_t)fNumberOfEvents;
          
          fT2d_lf_cov[j0StexEcha][i0StexEcha] = fT2d_lf_cov[i0StexEcha][j0StexEcha];
        }
        }
      if( i0StexEcha%100 == 0 ){std::cout << i0StexEcha << "[LFN Cov], ";}
    }
    }
  std::cout << std::endl;
  fTagLfCov[0] = 1;    fFileHeader->fLfCovCalc++;
}
//---------- (end of LowFrequencyCovariancesBetweenChannels ) --------------------

//------------------------------------------------------------------
//
//  Calculation of the Low Frequency Correlations between channels
//
//  LFCor(Ci,Cj) = LFCov(Ci,Cj)/sqrt(LFCov(Ci,Ci)*LFCov(Cj,Cj))
//
//------------------------------------------------------------------
void TEcnaRun::LowFrequencyCorrelationsBetweenChannels()
{
//Calculation of the Low Frequency Correlations between channels

  //... preliminary calculation of the covariances if not done yet.
  if ( fTagLfCov[0] != 1 )
    {LowFrequencyCovariancesBetweenChannels(); fTagLfCov[0] = 0;}
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::LowFrequencyCorrelationsBetweenChannels()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "          Calculation of the Low Frequency Correlations between channels" << std::endl
     << "          Starting allocation. "
     << std::endl;}

  //................. allocation fT2d_lf_cor + init to zero (mandatory)
  if( fT2d_lf_cor == nullptr ){
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT2d_lf_cor  = new Double_t*[n_StexEcha];                   fCnew++;
    fT2d1_lf_cor = new  Double_t[n_StexEcha*n_StexEcha];        fCnew++;
    for(Int_t i0StexEcha = 0 ; i0StexEcha < n_StexEcha ; i0StexEcha++){
    fT2d_lf_cor[i0StexEcha] = &fT2d1_lf_cor[0] + i0StexEcha*n_StexEcha;} 
  }
  
  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {
      if( fT2d_lf_cor[i0StexEcha][j0StexEcha] != (Double_t)0 )
        {fMiscDiag[22]++; fT2d_lf_cor[i0StexEcha][j0StexEcha] = (Double_t)0;}
    }
    }
  
  //................. calculation  
  //........................... computation of half of the matrix, diagonal included
  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha <= i0StexEcha ; j0StexEcha++)
        {
          if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, j0StexEcha) > 0)
          || ( fFlagSubDet == "EB") )
        {
          if(  fT2d_lf_cov[i0StexEcha][i0StexEcha] > 0
               && fT2d_lf_cov[j0StexEcha][j0StexEcha] > 0 )
            {
              fT2d_lf_cor[i0StexEcha][j0StexEcha] =
            fT2d_lf_cov[i0StexEcha][j0StexEcha]/
            ( (Double_t)sqrt( fT2d_lf_cov[i0StexEcha][i0StexEcha]*
                      fT2d_lf_cov[j0StexEcha][j0StexEcha] ) );
            }
          else
            {
              fT2d_lf_cor[i0StexEcha][j0StexEcha] = (Double_t)0.;
            }
          fT2d_lf_cor[j0StexEcha][i0StexEcha] = fT2d_lf_cor[i0StexEcha][j0StexEcha];
        }
        }
    }
      if( i0StexEcha%100 == 0 ){std::cout << i0StexEcha << "[LFN Cor], ";}
    }
  std::cout << std::endl;
  
  fTagLfCor[0] = 1;    fFileHeader->fLfCorCalc++;
}
//--------------- (end of LowFrequencyCorrelationsBetweenChannels) --------------------

//------------------------------------------------------------------
//
//  Calculation of the High Frequency Covariances between channels
//
//  HFCov(Ci,Cj) = E_e[ Cov_s[ A(Ci,s*,e*) , A(Cj,s*,e*) ] ] 
//
//               = E_e[ E_s[ ( A(Ci,s*,e*) - E_s[A(Ci,s*,e*)] )*
//                           ( A(Cj,s*,e*) - E_s[A(Cj,s*,e*)] ) ] ]
//   
//   A(Ci,s,e) : ADC value for channel Ci, sample s, event e
//
//   E_e         : average over the events
//   E_s , Cov_s : average, covariance over the samples
//
//------------------------------------------------------------------
void TEcnaRun::HighFrequencyCovariancesBetweenChannels()
{
//Calculation of the High Frequency Covariances between channels
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::HighFrequencyCovariancesBetweenChannels()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the High Frequency Covariances between channels" << std::endl;}

  //................. allocation fT2d_hf_cov + init to zero (mandatory)
  if( fT2d_hf_cov == nullptr ){
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT2d_hf_cov   = new Double_t*[n_StexEcha];                   fCnew++;
    fT2d1_hf_cov  = new  Double_t[n_StexEcha*n_StexEcha];        fCnew++;
    for(Int_t i0StexEcha = 0 ; i0StexEcha < n_StexEcha ; i0StexEcha++){
    fT2d_hf_cov[i0StexEcha] = &fT2d1_hf_cov[0] + i0StexEcha*n_StexEcha;} 
  }

  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {
      if( fT2d_hf_cov[i0StexEcha][j0StexEcha] != (Double_t)0 )
        {fMiscDiag[23]++; fT2d_hf_cov[i0StexEcha][j0StexEcha] = (Double_t)0;}
    }
    }

  //................... Allocation mean_over_samples(i0StexEcha, n_event) 
  TMatrixD mean_over_samples(fEcal->MaxCrysEcnaInStex(), fNumberOfEvents);
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {for(Int_t n_event=0; n_event<fNumberOfEvents; n_event++)
    {mean_over_samples(i0StexEcha, n_event)=(Double_t)0.;}}
  //................... Allocation cov_over_samp(i0StexEcha,j0StexEcha)
  TMatrixD cov_over_samp(fEcal->MaxCrysEcnaInStex(), fEcal->MaxCrysEcnaInStex());
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {for(Int_t j0StexEcha=0; j0StexEcha<fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {cov_over_samp(i0StexEcha, j0StexEcha)=(Double_t)0.;}}

  //........................................... Calculation    (HighFrequencyCovariancesBetweenChannels)
  if(fFlagPrint != fCodePrintNoComment){
    std::cout << "          Calculation of the mean (over the events)" << std::endl
     << "          of the covariances between the channels (over the samples 1 to "
     <<  fNbSampForCalc << ")." << std::endl;}

  std::cout << " Please, wait (end at i= " << fEcal->MaxCrysEcnaInStex() << "): " << std::endl;

  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
        {
          // Calculation, for each event, of the mean over the samples  ( = E_s[A(c_i,s*,e_n] )
          mean_over_samples(i0StexEcha, n_event) = (Double_t)0.;
          for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
        {mean_over_samples(i0StexEcha, n_event) += fT3d_AdcValues[i0StexEcha][i0Sample][n_event];}
          mean_over_samples(i0StexEcha, n_event) /= (Double_t)fNbSampForCalc;
        }
    }
      if( i0StexEcha%100 == 0 ){std::cout << i0StexEcha << "[HFNa Cov], ";}
    }
  std::cout << std::endl;

  std::cout << " Please, wait (end at i= " << fEcal->MaxCrysEcnaInStex() << "): " << std::endl;

  //... Calculation of half of the matrix, diagonal included    (HighFrequencyCovariancesBetweenChannels)
  for(Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for(Int_t j0StexEcha = 0; j0StexEcha <= i0StexEcha; j0StexEcha++)
        {
          if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, j0StexEcha) > 0)
          || ( fFlagSubDet == "EB") )
        {
          for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
            {
              // Calculation, for each event, of the covariance over the samples
              cov_over_samp(i0StexEcha,j0StexEcha) = (Double_t)0;
              for(Int_t i0Sample=0; i0Sample<fNbSampForCalc; i0Sample++)
            {
              cov_over_samp(i0StexEcha,j0StexEcha) +=
                ( fT3d_AdcValues[i0StexEcha][i0Sample][n_event] - mean_over_samples(i0StexEcha, n_event) )*
                ( fT3d_AdcValues[j0StexEcha][i0Sample][n_event] - mean_over_samples(j0StexEcha, n_event) );
            }
              cov_over_samp(i0StexEcha,j0StexEcha) /= (Double_t)fNbSampForCalc;

            }
          //....... Calculation of the mean over the events of Cov_s[A(c_i,s*,e*),A(c_j,s*,e*)]
          //......... Calculation of half of the matrix, diagonal included
          fT2d_hf_cov[i0StexEcha][j0StexEcha] = (Double_t)0;
          for (Int_t n_event = 0; n_event < fNumberOfEvents; n_event++)
            {
              fT2d_hf_cov[i0StexEcha][j0StexEcha] += cov_over_samp(i0StexEcha,j0StexEcha);
            }
            fT2d_hf_cov[i0StexEcha][j0StexEcha] /= (Double_t)fNumberOfEvents;
          
          fT2d_hf_cov[j0StexEcha][i0StexEcha] = fT2d_hf_cov[i0StexEcha][j0StexEcha];
        }
        }
    }
      if( i0StexEcha%100 == 0 ){std::cout << i0StexEcha << "[HFNb Cov], ";}
    }
  std::cout << std::endl;

  fTagHfCov[0] = 1;    fFileHeader->fHfCovCalc++;
}
//---------- (end of HighFrequencyCovariancesBetweenChannels ) --------------------

//------------------------------------------------------------------
//
//  Calculation of the High Frequency Correlations between channels
//
//  HFCor(Ci,Cj) = HFCov(Ci,Cj)/sqrt(HFCov(Ci,Ci)*HFCov(Cj,Cj))
//
//------------------------------------------------------------------
void TEcnaRun::HighFrequencyCorrelationsBetweenChannels()
{
//Calculation of the High Frequency Correlations between channels

  //... preliminary calculation of the covariances if not done yet.
  if ( fTagHfCov[0] != 1 )
    {HighFrequencyCovariancesBetweenChannels(); fTagHfCov[0] = 0;}
  
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::HighFrequencyCorrelationsBetweenChannels()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the High Frequency Correlations between channels" << std::endl
     << "          Starting allocation. "
     << std::endl;}

  //................. allocation fT2d_hf_cor + init to zero (mandatory)
  if( fT2d_hf_cor == nullptr ){
    const Int_t n_StexEcha = fEcal->MaxCrysEcnaInStex();
    fT2d_hf_cor  = new Double_t*[n_StexEcha];                   fCnew++;
    fT2d1_hf_cor = new  Double_t[n_StexEcha*n_StexEcha];        fCnew++;
    for(Int_t i0StexEcha = 0 ; i0StexEcha < n_StexEcha ; i0StexEcha++){
    fT2d_hf_cor[i0StexEcha] = &fT2d1_hf_cor[0] + i0StexEcha*n_StexEcha;} 
  }

  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha < fEcal->MaxCrysEcnaInStex() ; j0StexEcha++)
    {
      if( fT2d_hf_cor[i0StexEcha][j0StexEcha] != (Double_t)0 )
        {fMiscDiag[24]++; fT2d_hf_cor[i0StexEcha][j0StexEcha] = (Double_t)0;}
    }
    }
  
  //................. calculation  
  //........................... computation of half of the matrix, diagonal included
  
  for (Int_t i0StexEcha = 0 ; i0StexEcha < fEcal->MaxCrysEcnaInStex() ; i0StexEcha++)
    {
      if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, i0StexEcha) > 0)
      || ( fFlagSubDet == "EB") )
    {
      for (Int_t j0StexEcha = 0 ; j0StexEcha <= i0StexEcha ; j0StexEcha++)
        {
          if( (fFlagSubDet == "EE" && fEcalNumbering->StexEchaForCons(fFileHeader->fStex, j0StexEcha) > 0)
          || ( fFlagSubDet == "EB") )
        {
          if(  fT2d_hf_cov[i0StexEcha][i0StexEcha] > 0
               && fT2d_hf_cov[j0StexEcha][j0StexEcha] > 0 )
            {
              fT2d_hf_cor[i0StexEcha][j0StexEcha] =
            fT2d_hf_cov[i0StexEcha][j0StexEcha]/
            (  (Double_t)sqrt(fT2d_hf_cov[i0StexEcha][i0StexEcha]) *
               (Double_t)sqrt(fT2d_hf_cov[j0StexEcha][j0StexEcha])  );
            }
          else
            {
              fT2d_hf_cor[i0StexEcha][j0StexEcha] = (Double_t)0.;
            }
          
          fT2d_hf_cor[j0StexEcha][i0StexEcha] = fT2d_hf_cor[i0StexEcha][j0StexEcha];
        }
        }
    }
      if( i0StexEcha%100 == 0 ){std::cout << i0StexEcha << "[HFN Cor], ";}
    }
  std::cout << std::endl;  
  
  fTagHfCor[0] = 1;    fFileHeader->fHfCorCalc++;
}
//------- (end of HighFrequencyCorrelationsBetweenChannels) ----------

//=================================================================================
//
//          L O W  &  H I G H    F R E Q U E N C Y    C O R R E L A T I O N S
//
//       B E T W E E N   T O W E R S  ( E B )  O R   S C s   ( E E )
//
//=================================================================================
//-----------------------------------------------------------------------------
//      Calculation of the mean Low Frequency Correlations
//      between channels for each Stin
//-----------------------------------------------------------------------------
void TEcnaRun::LowFrequencyMeanCorrelationsBetweenTowers()
{LowFrequencyMeanCorrelationsBetweenStins();}
void TEcnaRun::LowFrequencyMeanCorrelationsBetweenSCs()
{LowFrequencyMeanCorrelationsBetweenStins();}

void TEcnaRun::LowFrequencyMeanCorrelationsBetweenStins()
{
//Calculation of the mean Low Frequency Correlations
//between channels for each Stin

  //... preliminary calculation of the Low Frequency Cor(c,c) if not done yet
  //    Only one tag (dim=1) to set to 0 (no write in the result ROOT file)
  if(fTagLfCor[0] != 1){LowFrequencyCorrelationsBetweenChannels(); fTagLfCor[0]=0;}

  //..... mean fT2d_lfcc_mostins for each pair (Stin_X,Stin_Y)
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::LowFrequencyMeanCorrelationsBetweenStins()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the mean, for each "
     << fStinName.Data() << ", of the" << std::endl
     << "           Low Frequency Correlations between channels." << std::endl;}
  
  //................. allocation fT2d_lfcc_mostins + init to zero (mandatory)
  if( fT2d_lfcc_mostins == nullptr ){
    const Int_t n_Stin = fEcal->MaxStinEcnaInStex();
    fT2d_lfcc_mostins  = new Double_t*[n_Stin];                 fCnew++;
    fT2d1_lfcc_mostins = new  Double_t[n_Stin*n_Stin];          fCnew++;  
    for(Int_t i0StexStinEcna = 0 ; i0StexStinEcna < n_Stin ; i0StexStinEcna++){
    fT2d_lfcc_mostins[i0StexStinEcna] = &fT2d1_lfcc_mostins[0] + i0StexStinEcna*n_Stin;} 
  }
  
  for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for(Int_t j0StexStinEcna=0; j0StexStinEcna<fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      if( fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna] != (Double_t)0 )
        {fMiscDiag[31]++; fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna] = (Double_t)0;}
    }
    }

  //..... Calculation of the mean LF Cor(c,c) for each pair (Stin_X,Stin_Y)
  //
  //           ! => Warning: this matrix is NOT symmetric => take N*N elements
  //                Only (Stin,Stin) matrix is symmetric.
  //                (StinEcha,StinEcha) matrix inside a (Stin,Stin) element is NOT symmetric
  //                (except for the (Stin,Stin) DIAGONAL elements)
  //      Then:
  //            1D array half_LFccMos[N*N] to put the (channel,channel) correlations

  Int_t ndim = (Int_t)(fEcal->MaxCrysInStin()*fEcal->MaxCrysInStin());

  TVectorD  half_LFccMos(ndim); for(Int_t i=0; i<ndim; i++){half_LFccMos(i)=(Double_t)0.;}

  //..................... Calculation
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::LowFrequencyMeanCorrelationsBetweenStins()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for each "
     << fStinName.Data() << ", of the mean " << std::endl
     << "           Low Frequency cor(c,c)." << std::endl;}

  for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for(Int_t j0StexStinEcna=0; j0StexStinEcna<fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      //................... .Copy the Mean Correlations(c,c') in 1D vector half_LFccMos() 
      Int_t i_count = 0;
      for(Int_t i0StinCrys=0; i0StinCrys<fEcal->MaxCrysInStin(); i0StinCrys++)
        {
          Int_t i0StexEcha = i0StexStinEcna*fEcal->MaxCrysInStin() + i0StinCrys;
          for(Int_t j0StinCrys=0; j0StinCrys<fEcal->MaxCrysInStin(); j0StinCrys++)
        {
          Int_t j0StexEcha = j0StexStinEcna*fEcal->MaxCrysInStin() + j0StinCrys;
          if(  (i0StexEcha>= 0 && i0StexEcha < fEcal->MaxCrysEcnaInStex()) &&
               (j0StexEcha>= 0 && j0StexEcha < fEcal->MaxCrysEcnaInStex())  )
          {half_LFccMos(i_count) = fT2d_lf_cor[i0StexEcha][j0StexEcha]; i_count++;}
          else
            {std::cout << "!TEcnaRun::LowFrequencyMeanCorrelationsBetweenStins()> Channel number out of range."
              << "i0StexEcha = " << i0StexEcha <<", j0StexEcha = " << j0StexEcha << fTTBELL << std::endl; }
        }
        }
      //...... Calculation of the mean absolute values of the LF mean Correlations(c,c')
      fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna] = (Double_t)0;
      for(Int_t i_rcor = 0; i_rcor < ndim; i_rcor++)
        {
          fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna] += fabs(half_LFccMos(i_rcor));
        }
      fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna] /= (Double_t)ndim;
    }
      if( i0StexStinEcna%10 == 0 ){std::cout << i0StexStinEcna << "[LFN MCtt], ";}
    }
  std::cout << std::endl;

  fTagLFccMoStins[0] = 1;           fFileHeader->fLFccMoStinsCalc++;
} // ------- end of LowFrequencyMeanCorrelationsBetweenStins() -------

//-----------------------------------------------------------------------------
//      Calculation of the mean High Frequency Correlations
//      between channels for each Stin
//-----------------------------------------------------------------------------
void TEcnaRun::HighFrequencyMeanCorrelationsBetweenTowers()
{HighFrequencyMeanCorrelationsBetweenStins();}
void TEcnaRun::HighFrequencyMeanCorrelationsBetweenSCs()
{HighFrequencyMeanCorrelationsBetweenStins();}

void TEcnaRun::HighFrequencyMeanCorrelationsBetweenStins()
{
//Calculation of the mean High Frequency Correlations
//between channels for each Stin

  //... preliminary calculation of the High Frequency Cor(c,c) if not done yet
  //    Only one tag (dim=1) to set to 0 (no write in the result ROOT file)
  if(fTagHfCor[0] != 1){HighFrequencyCorrelationsBetweenChannels();fTagHfCor[0]=0;}

  //..... mean fT2d_hfcc_mostins for each pair (Stin_X,Stin_Y)
  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::HighFrequencyMeanCorrelationsBetweenStins()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation of the mean, for each "
     << fFlagSubDet.Data() << ", of the" << std::endl
         << "           High Frequency Correlations between channels." << std::endl;}

  //................. allocation fT2d_hfcc_mostins + init to zero (mandatory)
  if( fT2d_hfcc_mostins == nullptr ){
    const Int_t n_Stin = fEcal->MaxStinEcnaInStex();
    fT2d_hfcc_mostins  = new Double_t*[n_Stin];                 fCnew++;
    fT2d1_hfcc_mostins = new  Double_t[n_Stin*n_Stin];           fCnew++;  
    for(Int_t i0StexStinEcna = 0 ; i0StexStinEcna < n_Stin ; i0StexStinEcna++){
    fT2d_hfcc_mostins[i0StexStinEcna] = &fT2d1_hfcc_mostins[0] + i0StexStinEcna*n_Stin;} 
  }
  
  for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for(Int_t j0StexStinEcna=0; j0StexStinEcna<fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      if( fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna] != (Double_t)0 )
        {fMiscDiag[32]++; fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna] = (Double_t)0;}
    }
    }

  //..... Calculation of the mean HF Cor(c,c) for each pair (Stin_X,Stin_Y)
  //
  //           ! => Warning: this matrix is NOT symmetric => take N*N elements
  //                Only (Stin,Stin) matrix is symmetric.
  //                (StinEcha,StinEcha) matrix inside a (Stin,Stin) element is NOT symmetric
  //                (except for the (Stin,Stin) DIAGONAL elements)
  //      Then:
  //            1D array half_LFccMos[N*N] to put the (channel,channel) correlations

  Int_t ndim = (Int_t)(fEcal->MaxCrysInStin()*fEcal->MaxCrysInStin());

  TVectorD half_HFccMos(ndim); for(Int_t i=0; i<ndim; i++){half_HFccMos(i)=(Double_t)0.;}

  if(fFlagPrint != fCodePrintNoComment){std::cout << "*TEcnaRun::HighFrequencyMeanCorrelationsBetweenStins()" << std::endl;}
  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "           Calculation, for each "
     << fFlagSubDet.Data()  << ", of the mean " << std::endl
     << "           High Frequency cor(c,c)." << std::endl;}

  //..................... Calculation
  for(Int_t i0StexStinEcna=0; i0StexStinEcna<fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for(Int_t j0StexStinEcna=0; j0StexStinEcna<fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      //.................... Copy the relevant Mean Correlations(c,c') in 1D vector half_HFccMos()
      Int_t i_count = 0;
      for(Int_t i0StinCrys=0; i0StinCrys<fEcal->MaxCrysInStin(); i0StinCrys++)
        {
          Int_t i0StexEcha = i0StexStinEcna*fEcal->MaxCrysInStin() + i0StinCrys;
          for(Int_t j0StinCrys=0; j0StinCrys<fEcal->MaxCrysInStin(); j0StinCrys++)
        {
          Int_t j0StexEcha = j0StexStinEcna*fEcal->MaxCrysInStin() + j0StinCrys;
          if(  (i0StexEcha>= 0 && i0StexEcha < fEcal->MaxCrysEcnaInStex()) &&
               (j0StexEcha>= 0 && j0StexEcha < fEcal->MaxCrysEcnaInStex())  )
            {half_HFccMos(i_count) = fT2d_hf_cor[i0StexEcha][j0StexEcha]; i_count++;}
          else
            {std::cout << "!TEcnaRun::HighFrequencyMeanCorrelationsBetweenStins()> Channel number out of range."
              << "i0StexEcha = " << i0StexEcha <<", j0StexEcha = " << j0StexEcha << fTTBELL << std::endl; }
        }
        }
      //..... Calculation of the mean absolute values of the HF mean Correlations(c,c')
      fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna] = (Double_t)0;
      for(Int_t i_rcor = 0; i_rcor < ndim; i_rcor++)
        {
          fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna] += fabs(half_HFccMos(i_rcor));
        }
      fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna] /= (Double_t)ndim;
    }
      if( i0StexStinEcna%10 == 0 ){std::cout << i0StexStinEcna << "[HFN MCtt], ";}
    }
  std::cout << std::endl;

  fTagHFccMoStins[0] = 1;                   fFileHeader->fHFccMoStinsCalc++;
} // ------- end of HighFrequencyMeanCorrelationsBetweenStins() -------

//=========================================================================
//
//                  W R I T I N G     M E T H O D S
//
//=========================================================================

//=========================================================================
//
//         W R I T I N G   M E T H O D S :    R O O T    F I L E S
//
//=========================================================================
//-------------------------------------------------------------
//
//                      OpenRootFile
//
//-------------------------------------------------------------
Bool_t TEcnaRun::OpenRootFile(const Text_t *name, const TString& status) {
//Open the Root file

  Bool_t ok_open = kFALSE;

  TString s_name;
  s_name = fCnaParPaths->ResultsRootFilePath();
  s_name.Append('/');
  s_name.Append(name);
  
  //gCnaRootFile = new TEcnaRootFile(fObjectManager, s_name.Data(), status);     fCnew++;

      Long_t iCnaRootFile = fObjectManager->GetPointerValue("TEcnaRootFile");
      if( iCnaRootFile == 0 )
    {
      gCnaRootFile = new TEcnaRootFile(fObjectManager, s_name.Data(), status); /* Anew("gCnaRootFile");*/
    }
      else
    {
      gCnaRootFile = (TEcnaRootFile*)iCnaRootFile;
      gCnaRootFile->ReStart(s_name.Data(), status);
    }

  if ( gCnaRootFile->fRootFileStatus == "RECREATE" ){ok_open = gCnaRootFile->OpenW();}
  if ( gCnaRootFile->fRootFileStatus == "READ"     ){ok_open = gCnaRootFile->OpenR();}

  if (!ok_open) // unable to open file
    {
      std::cout << "TEcnaRun::OpenRootFile> Cannot open file " << s_name.Data() << std::endl;
    }
  else
    {
      if(fFlagPrint == fCodePrintAllComments)
    {std::cout << "*TEcnaRun::OpenRootFile> Open ROOT file OK for file " << s_name.Data() << std::endl;}  
      fOpenRootFile  = kTRUE;
    }
  return ok_open;
}
//-------------------------------------------------------------
//
//                      CloseRootFile
//
//-------------------------------------------------------------
Bool_t TEcnaRun::CloseRootFile(const Text_t *name) {
//Close the Root file
 
  TString s_name;
  s_name = fCnaParPaths->ResultsRootFilePath();
  s_name.Append('/');
  s_name.Append(name);

  Bool_t ok_close = kFALSE;

  if (fOpenRootFile == kTRUE ) 
    {
      gCnaRootFile->CloseFile();

      if(fFlagPrint != fCodePrintAllComments){
    std::cout << "*TEcnaRun::CloseRootFile> ROOT file " << s_name.Data() << " closed." << std::endl;}

      //     delete gCnaRootFile;     gCnaRootFile = 0;          fCdelete++;

      fOpenRootFile = kFALSE;
      ok_close      = kTRUE;
    }
  else
    {
      std::cout << "*TEcnaRun::CloseRootFile(...)> No close since no file is open."
       << fTTBELL << std::endl;
    }
  return ok_close;
}
//-------------------------------------------------------------
//
//   WriteRootFile without arguments.
//   Call WriteRootFile WITH argument (file name)
//   after an automatic generation of the file name.
//
//   Codification for the file name:
//            see comment at the beginning of this file
//
//-------------------------------------------------------------

//=================================================================================
//
//         WriteRootFile()   ====>  D O N ' T    S U P P R E S S  ! ! !
//                                  Called by the analyzer in package: "Modules"
//
//=================================================================================
Bool_t TEcnaRun::WriteRootFile(){
//Write the Root file.
//File name automatically generated by fCnaWrite->fMakeResultsFileName()
//previously called in GetReadyToCompute().

  Bool_t ok_write = kFALSE;

  //============================= check number of found events
  Int_t nCountEvts = 0;
  
  for(Int_t i0StexEcha=0; i0StexEcha<fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {for(Int_t i0Sample=0; i0Sample<fFileHeader->fNbOfSamples; i0Sample++)
      {nCountEvts += fT2d_NbOfEvts[i0StexEcha][i0Sample];}}
  
  if ( nCountEvts <= 0 )
    {
      //============== no write if no event found
      std::cout << "!TEcnaRun::WriteRootFile()> No event found for file " << fCnaWrite->GetRootFileNameShort().Data()
       << ". File will not be written." << std::endl;
      ok_write = kTRUE;
    }
  else
    {    
      if(fFlagPrint == fCodePrintAllComments){
    std::cout << "*TEcnaRun::WriteRootFile()> Results are going to be written in the ROOT file: " << std::endl
         << "                           " << fCnaWrite->GetRootFileName().Data() << std::endl;}

      const Text_t *FileShortName = (const Text_t *)fCnaWrite->GetRootFileNameShort().Data();
      ok_write = WriteRootFile(FileShortName, fFileHeader->fNbOfSamples);

      if( ok_write == kTRUE )
    {
      if(fFlagPrint != fCodePrintNoComment)
        {std::cout << "*TEcnaRun::WriteRootFile()> Writing OK for file " << fCnaWrite->GetRootFileName().Data()
          << std::endl;}
    }
      else
    {
      std::cout << "!TEcnaRun::WriteRootFile()> Writing FAILLED for file " << fCnaWrite->GetRootFileName().Data()
           << fTTBELL << std::endl;
    }
    }
  return ok_write;
} // end of WriteRootFile()

//--------------------------------------------------------------------
//
//               WriteNewRootFile with argument
//    Called by TEcnaGui for results file of the Calculations method
//    analysis type and nb of samples changed, other arguments kept
//
//--------------------------------------------------------------------
Bool_t TEcnaRun::WriteNewRootFile(const TString& TypAna){
//Write a new Root file. File name automatically generated by fCnaWrite->fMakeResultsFileName()
//called here.

  Bool_t ok_write = kFALSE;

  fCnaWrite->RegisterFileParameters(TypAna.Data(),   fNbSampForCalc, fFileHeader->fRunNumber,
                    fFileHeader->fFirstReqEvtNumber, fFileHeader->fLastReqEvtNumber,
                    fFileHeader->fReqNbOfEvts,       fFileHeader->fStex,
                    fFileHeader->fStartDate, fFileHeader->fStopDate,
                    fFileHeader->fStartTime, fFileHeader->fStopTime);

  fCnaWrite->fMakeResultsFileName();  // set fRootFileName, fRootFileNameShort
  fNewRootFileName      = fCnaWrite->GetRootFileName();
  fNewRootFileNameShort = fCnaWrite->GetRootFileNameShort();

  const Text_t *FileShortName = (const Text_t *)fNewRootFileNameShort.Data();

  if(fFlagPrint == fCodePrintAllComments){
    std::cout << "*TEcnaRun::WriteNewRootFile()> Results are going to be written in the ROOT file: " << std::endl
     << "                              " << fNewRootFileNameShort.Data() << std::endl;}

  ok_write = WriteRootFile(FileShortName, fNbSampForCalc);

  return ok_write;
}

//-------------------------------------------------------------------------
//
//    Get the new ROOT file name (long and short)
//   (called by TEcnaGui in Calculations method)
//
//-------------------------------------------------------------------------
TString TEcnaRun::GetNewRootFileName(){return fNewRootFileName;}
TString TEcnaRun::GetNewRootFileNameShort(){return fNewRootFileNameShort;}

//--------------------------------------------------------------------
//
//               WriteRootFile with argument
//
//--------------------------------------------------------------------
Bool_t TEcnaRun::WriteRootFile(const Text_t* name, Int_t& argNbSampWrite) 
{
//Write the Root file

  const Text_t* file_name = name;

  Bool_t ok_write = kFALSE;

  if ( fOpenRootFile )
    {
      std::cout << "!TEcnaRun::WriteRootFile(...) *** ERROR ***> Writing on file already open."
       << fTTBELL << std::endl;
    }
  else
    {
      // List of the different element types and associated parameters as ordered in the ROOT file (smaller -> larger)
      //                                                                 ==========
      //
      //         WARNING  *** HERE SIZES ARE THESE FOR THE BARREL (1700 Xtals) and for 10 samples ***
      //
      //   Nb of   Type of element            Type      Type                                    Size    Comment
      // elements                             Number    Name
      //
      //        1  fMatHis(1,StexStin)         ( 0)  cTypNumbers             1*(   1,  68) =         68

      //        1  fMatHis(1,StexStin)         (12)  cTypAvPed               1*(   1,  68) =         68
      //        1  fMatHis(1,StexStin)         ( 3)  cTypAvTno               1*(   1,  68) =         68
      //        1  fMatHis(1,StexStin)         ( 4)  cTypAvLfn               1*(   1,  68) =         68
      //        1  fMatHis(1,StexStin)         ( 5)  cTypAvHfn               1*(   1,  68) =         68
      //        1  fMatHis(1,StexStin)         (13)  cTypAvMeanCorss         1*(   1,  68) =         68
      //        1  fMatHis(1,StexStin)         (14)  cTypAvSigCorss          1*(   1,  68) =         68

      //        1  fMatHis(1,StexEcha)         (16)  cTypPed                 1*(   1,1700) =      1 700
      //        1  fMatHis(1,StexEcha)         (17)  cTypTno                 1*(   1,1700) =      1 700
      //        1  fMatHis(1,StexEcha)         (10)  cTypMeanCorss           1*(   1,1700) =      1 700
      //        1  fMatHis(1,StexEcha)         (18)  cTypLfn                 1*(   1,1700) =      1 700
      //        1  fMatHis(1,StexEcha)         (19)  cTypHfn                 1*(   1,1700) =      1 700
      //        1  fMatHis(1,StexEcha)         (11)  cTypSigCorss            1*(   1,1700) =      1 700

      //        1  fMatMat(Stin,Stin)          (23)  cTypLFccMoStins         1*(  68,  68) =      4 624
      //        1  fMatMat(Stin,Stin)          (24)  cTypHFccMoStins         1*(  68,  68) =      4 624

      //        1  fMatHis(StexEcha, sample)   (15)  cTypNbOfEvts            1*(1700,  10) =     17 000
      //        1  fMatHis(StexEcha, sample)   ( 1)  cTypMSp                 1*(1700,  10) =     17 000
      //        1  fMatHis(StexEcha, sample)   ( 2)  cTypSSp                 1*(1700,  10) =     17 000
 
      //   StexEcha  fMatMat(sample, sample)   ( 8)  cTypCovCss           1700*(  10,  10) =    170 000
      //   StexEcha  fMatMat(sample, sample    ( 9)  cTypCorCss           1700*(  10,  10) =    170 000

      //   StexEcha  fMatHis(sample, bin_evt)  (20)  cTypAdcEvt,          1700*(  10, 150) =  2 550 000

      //        1  fMatMat(StexEcha, StexEcha) (21)  cTypLfCov               1*(1700,1700) =  2 890 000
      //        1  fMatMat(StexEcha, StexEcha) (22)  cTypLfCor               1*(1700,1700) =  2 890 000

      //        1  fMatMat(StexEcha, StexEcha) ( 6)  cTypHfCov               1*(1700,1700) =  2 890 000 // (06/05/08)
      //        1  fMatMat(StexEcha, StexEcha) ( 7)  cTypHfCor               1*(1700,1700) =  2 890 000 // (06/05/08)


      //......................................................................................................

      OpenRootFile(file_name, "RECREATE");

      TString typ_name = "?";
      Int_t v_nb_times = 0;
      Int_t v_dim_one  = 0;
      Int_t v_dim_two  = 0;
      Int_t v_size     = 0;
      Int_t v_tot      = 0;
      Int_t v_tot_writ = 0;

      //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      //       
      //  ===> no general method and no translation to TEcnaWrite
      //       because the fT1d.. and fT2d... arrays
      //       are attributes of TEcnaRun (calls to the "TRootXXXX" methods)
      //
      //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

      //-------------------------- Stin numbers 
      //       1   fMatHis(1,Stin)           ( 0)  cTypNumbers        1*(   1,  68) =         68

      Int_t MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "StinNumbers";
      v_nb_times = fFileHeader->fStinNumbersCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagStinNumbers[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypNumbers;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootStinNumbers();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Average Pedestals (1 value per Stin)
      //       1   fMatHis(1, StexStin)   (12)  cTypAvPed      1*(1,  68) =     68

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AvPed";
      v_nb_times = fFileHeader->fAvPedCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagAvPed[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvPed;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvPed();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Average Total noise
      // StexEcha   fMatHis(1, StexStin)     ( 3)  cTypAvTno      1*(1,  68) =     68

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AvTno";
      v_nb_times = fFileHeader->fAvTnoCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagAvTno[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvTno;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvTno();
          gCnaRootFile->fCnaResultsTree->Fill(); 
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Average Low frequency noise
      //       1   fMatHis(1, StexStin)   ( 4)  cTypAvLfn      1*(1,  68) =     68

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AvLfn";
      v_nb_times = fFileHeader->fAvLfnCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}
      
      for (Int_t i = 0; i < v_nb_times; i++)
        {
          if ( fTagAvLfn[0] == 1 )
            {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvLfn;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvLfn();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
            }
        }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Average High frequency noise
      //       1   fMatHis(1, StexStin)   ( 5)  cTypAvHfn      1*(1,  68) =     68

      MaxCar = fgMaxCar;      
      typ_name.Resize(MaxCar);
      typ_name   = "AvHfn";
      v_nb_times = fFileHeader->fAvHfnCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}
      
      for (Int_t i = 0; i < v_nb_times; i++)
        {
          if ( fTagAvHfn[0] == 1 )
            {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvHfn;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvHfn();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
            }
        }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Average mean cor(s,s)
      //       1   fMatHis(1, StexStin)   (13)  cTypAvMeanCorss      1*(1,  68) =     68

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AvMeanCorss";
      v_nb_times = fFileHeader->fAvMeanCorssCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagAvMeanCorss[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvMeanCorss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvEvCorss();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //--------------------------  Average sigma of cor(s,s)
      //       1   fMatHis(1, StexStin)    (14)  cTypAvSigCorss      1*(1,  68) =     68

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AvSigCorss";
      v_nb_times = fFileHeader->fAvSigCorssCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagAvSigCorss[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAvSigCorss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAvSigCorss();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Expectation values of the expectation values of the samples (pedestals)
      //       1   fMatHis(1,StexEcha)         (16)  cTypPed                1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "Ped";
      v_nb_times = fFileHeader->fPedCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagPed[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypPed;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootPed();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
                  
      //-------------------------- Expectation values of the sigmas the samples
      //       1   fMatHis(1,StexEcha)         (17)  cTypTno               1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "Tno";
      v_nb_times = fFileHeader->fTnoCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagTno[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypTno;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootTno();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
                                     
      //-------------------------- Expectation values of the correlations between the samples
      //       1   fMatHis(1,StexEcha)         (10)  cTypMeanCorss            1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "MeanCorss";
      v_nb_times = fFileHeader->fMeanCorssCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagMeanCorss[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypMeanCorss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootMeanCorss();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    } 
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Sigmas of the expectation values of the samples  
      //       1   fMatHis(1,StexEcha)         (18)  cTypLfn               1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "Lfn";
      v_nb_times = fFileHeader->fLfnCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagLfn[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypLfn;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootLfn();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    } 
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Sigmas of the sigmas of the samples  
      //       1   fMatHis(1,StexEcha)         (19)  cTypHfn              1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "Hfn";
      v_nb_times = fFileHeader->fHfnCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}
 
      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagHfn[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypHfn;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootHfn();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
                    
      //-------------------------- Sigmas of the correlations between the samples  
      //       1   fMatHis(1,StexEcha)         (11)  cTypSigCorss           1*(   1,1700) =      1 700

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "SigCorss";
      v_nb_times = fFileHeader->fSigCorssCalc;
      v_dim_one  = 1;
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}
   
      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagSigCorss[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypSigCorss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootSigCorss();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //----- Mean Covariances between StexEchas (averaged over samples) for all (Stin_X,Stin_Y)
      //       1   fMatMat(Stin,Stin)       (23)  cTypLFccMoStins         1*(  68,  68) =      4 624

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "LFccMoStins";
      v_nb_times = fFileHeader->fLFccMoStinsCalc;
      v_dim_one  = fEcal->MaxStinEcnaInStex();
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagLFccMoStins[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypLFccMoStins;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootLFccMoStins();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    } 
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
            
      //----- Mean Correlations between StexEchas (averaged over samples) for all (Stin_X,Stin_Y)
      //       1   fMatMat(Stin,Stin)       (24)  cTypHFccMoStins         1*(  68,  68) =      4 624

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "HFccMoStins";
      v_nb_times = fFileHeader->fHFccMoStinsCalc;
      v_dim_one  = fEcal->MaxStinEcnaInStex();
      v_dim_two  = fEcal->MaxStinEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagHFccMoStins[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypHFccMoStins;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootHFccMoStins();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Numbers of found events (NbOfEvts)
      //       1   fMatHis(StexEcha, sample)   (15)  cTypNbOfEvts       1*(1700,  10) =     17 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "NbOfEvts";
      v_nb_times = fFileHeader->fNbOfEvtsCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = argNbSampWrite;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagNbOfEvts[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypNbOfEvts;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootNbOfEvts(argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
                  
      //-------------------------- Expectation values of the samples
      //       1   fMatHis(StexEcha, sample)   ( 1)  cTypMSp                  1*(1700,  10) =     17 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "MSp";
      v_nb_times = fFileHeader->fMSpCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = argNbSampWrite;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagMSp[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypMSp;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootMSp(argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Sigmas of the samples     
      //       1   fMatHis(StexEcha, sample)   ( 2)  cTypSSp                 1*(1700,  10) =     17 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "SSp";
      v_nb_times = fFileHeader->fSSpCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = argNbSampWrite;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagSSp[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypSSp;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootSSp(argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Covariances between samples

      // StexEcha   fMatMat(sample,  sample)   ( 8)  cTypCovCss           1700*(  10,  10) =    170 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "CovCss";
      v_nb_times = fFileHeader->fCovCssCalc;
      v_dim_one  = argNbSampWrite;
      v_dim_two  = argNbSampWrite;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i0StexEcha = 0; i0StexEcha < v_nb_times; i0StexEcha++)
    {
      if ( fTagCovCss[i0StexEcha] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypCovCss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i0StexEcha;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootCovCss(i0StexEcha, argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i0StexEcha == 0  && fFlagPrint == fCodePrintAllComments)
        {std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Correlations between samples   
      // StexEcha   fMatMat(sample,  sample)   ( 9)  cTypCorCss           1700*(  10,  10) =    170 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "CorCss";
      v_nb_times = fFileHeader->fCorCssCalc;
      v_dim_one  = argNbSampWrite;
      v_dim_two  = argNbSampWrite;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}
 
      for (Int_t i0StexEcha = 0; i0StexEcha < v_nb_times; i0StexEcha++)
    {
      if ( fTagCorCss[i0StexEcha] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypCorCss;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i0StexEcha;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootCorCss(i0StexEcha, argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i0StexEcha == 0  && fFlagPrint == fCodePrintAllComments)
        {std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }     
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Samples as a function of event = events distributions
      // StexEcha   fMatHis(sample,  bins)     (20)  cTypAdcEvt,        1700*(  10, 150) =  2 550 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "AdcEvt";
      v_nb_times = fFileHeader->fAdcEvtCalc;
      v_dim_one  = argNbSampWrite;
      v_dim_two  = fFileHeader->fReqNbOfEvts;
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i0StexEcha = 0; i0StexEcha < v_nb_times; i0StexEcha++)
    {
      if ( fTagAdcEvt[i0StexEcha] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypAdcEvt;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i0StexEcha;
          gCnaRootFile->fCnaIndivResult->SetSizeHis(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatMat.ReSet(1,1);
          TRootAdcEvt(i0StexEcha, argNbSampWrite);
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i0StexEcha == 0  && fFlagPrint == fCodePrintAllComments )
        {std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- Low Frequency Covariances between StexEchas
      //  sample   fMatMat(StexEcha, StexEcha)  (21)  cTypLfCov           1*(1700,1700) =  2 890 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "LfCov";
      v_nb_times = fFileHeader->fLfCovCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {      //=================================== Record type EB
      if ( fTagLfCov[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypLfCov;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootLfCov();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- Low Frequency Correlations between StexEchas
      //  sample   fMatMat(StexEcha, StexEcha)  (22)  cTypLfCor           1*(1700,1700) =  2 890 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "LfCor";
      v_nb_times = fFileHeader->fLfCorCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagLfCor[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypLfCor;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootLfCor();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //-------------------------- High Frequency Covariances between StexEchas
      //  sample   fMatMat(StexEcha, StexEcha)  (6)  cTypHfCov           1*(1700,1700) =  2 890 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "HfCov";
      v_nb_times = fFileHeader->fHfCovCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;
 
      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagHfCov[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypHfCov;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootHfCov();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}
      
      //-------------------------- High Frequency Correlations between StexEchas
      //  sample   fMatMat(StexEcha, StexEcha)  (7)  cTypHfCor           1*(1700,1700) =  2 890 000

      MaxCar = fgMaxCar;
      typ_name.Resize(MaxCar);
      typ_name   = "HfCor";
      v_nb_times = fFileHeader->fHfCorCalc;
      v_dim_one  = fEcal->MaxCrysEcnaInStex();
      v_dim_two  = fEcal->MaxCrysEcnaInStex();
      v_size     = v_nb_times*v_dim_one*v_dim_two;
      v_tot     += v_size;

      if(fFlagPrint == fCodePrintAllComments){
      std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(18) << typ_name << ": " << std::setw(4) << v_nb_times
       << " * ("  << std::setw(4) << v_dim_one << ","  << std::setw(4) << v_dim_two  << ") = "
       << std::setw(9) << v_size;}

      for (Int_t i = 0; i < v_nb_times; i++)
    {
      if ( fTagHfCor[0] == 1 )
        {
          gCnaRootFile->fCnaIndivResult->fTypOfCnaResult = cTypHfCor;
          gCnaRootFile->fCnaIndivResult->fIthElement     = i;
          gCnaRootFile->fCnaIndivResult->SetSizeMat(v_dim_one,v_dim_two);
          gCnaRootFile->fCnaIndivResult->fMatHis.ReSet(1,1);
          TRootHfCor();
          gCnaRootFile->fCnaResultsTree->Fill();
          if( i == 0 && fFlagPrint == fCodePrintAllComments ){std::cout << " => WRITTEN ON FILE "; v_tot_writ += v_size;}
        }
    }
      if(fFlagPrint == fCodePrintAllComments){std::cout << std::endl;}

      //---------------------------------------------- WRITING 
      //...................................... file 
      gCnaRootFile->fRootFile->Write();
      //...................................... header
      fFileHeader->Write();

      //...................................... status message
      if(fFlagPrint == fCodePrintAllComments){
    std::cout << "*TEcnaRun::WriteRootFile(...)> " << std::setw(20) << "TOTAL: "
         << std::setw(21) << "CALCULATED = " << std::setw(9) <<  v_tot
         << " => WRITTEN ON FILE = "    << std::setw(9) << v_tot_writ << std::endl;}

      if(fFlagPrint == fCodePrintAllComments){
    std::cout << "*TEcnaRun::WriteRootFile(...)> Write OK in file " << file_name << " in directory:" << std::endl
         << "                           " << fCnaParPaths->ResultsRootFilePath().Data()
         << std::endl;}

      ok_write = kTRUE;

      //...................................... close
      CloseRootFile(file_name);
    }
  return ok_write;
}  //-------------- End of WriteRootFile(...) -----------------------

//======================== "PREPA FILL" METHODS ===========================

//-------------------------------------------------------------------------
//
//  Prepa Fill Stin numbers as a function of the Stin index 
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootStinNumbers()
{
  if (fTagStinNumbers[0] == 1 )
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_StexStinFromIndex[j0StexStinEcna];
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill last evt numbers for all the (StexEcha,sample)
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootNbOfEvts(const Int_t& argNbSampWrite)
{
  if (fTagNbOfEvts[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      for (Int_t i0Sample = 0; i0Sample < argNbSampWrite; i0Sample++)
        {
          gCnaRootFile->fCnaIndivResult->fMatHis(j0StexEcha, i0Sample) =
            fT2d_NbOfEvts[j0StexEcha][i0Sample];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill histogram of samples as a function of event
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAdcEvt(const Int_t& user_StexEcha, const Int_t& argNbSampWrite)
{
  if (fTagAdcEvt[user_StexEcha] == 1 )
    {
      for (Int_t i0Sample = 0; i0Sample < argNbSampWrite; i0Sample++)
    {
      //...................... all the bins set to zero
      for (Int_t j_bin = 0; j_bin < fFileHeader->fReqNbOfEvts; j_bin++)
        {
          gCnaRootFile->fCnaIndivResult->fMatHis(i0Sample, j_bin) = (Double_t)0.;
        }
      //...................... fill the non-zero bins 
      for (Int_t j_bin = 0; j_bin < fFileHeader->fReqNbOfEvts; j_bin++)
        {
          gCnaRootFile->fCnaIndivResult->fMatHis(i0Sample, j_bin) =
        fT3d_AdcValues[user_StexEcha][i0Sample][j_bin]; 
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill expectation values of the samples for all the StexEchas
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootMSp(const Int_t& argNbSampWrite)
{
  if (fTagMSp[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      for (Int_t i0Sample = 0; i0Sample < argNbSampWrite; i0Sample++)
        {
          gCnaRootFile->fCnaIndivResult->fMatHis( j0StexEcha, i0Sample) =
            fT2d_ev[j0StexEcha][i0Sample];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill sigmas of the samples for all the StexEchas
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootSSp(const Int_t& argNbSampWrite)
{
  if (fTagSSp[0] == 1 )
    {
      for (Int_t j0StexEcha=0; j0StexEcha<fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      for (Int_t i0Sample=0; i0Sample<argNbSampWrite; i0Sample++)
        {
          gCnaRootFile->fCnaIndivResult->fMatHis(j0StexEcha, i0Sample) =
            fT2d_sig[j0StexEcha][i0Sample];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill mean covariances between StexEchas, mean over samples
//  for all (Stin_X, Stin_Y)
//                           (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootLFccMoStins()
{
  if (fTagLFccMoStins[0] == 1 )
    {
      for (Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexStinEcna, j0StexStinEcna) =
            fT2d_lfcc_mostins[i0StexStinEcna][j0StexStinEcna];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill mean correlations between StexEchas, mean over samples
//  for all (Stin_X, Stin_Y)
//                           (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootHFccMoStins()
{
  if (fTagHFccMoStins[0] == 1 )
    {
      for (Int_t i0StexStinEcna = 0; i0StexStinEcna < fEcal->MaxStinEcnaInStex(); i0StexStinEcna++)
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexStinEcna, j0StexStinEcna) =
            fT2d_hfcc_mostins[i0StexStinEcna][j0StexStinEcna];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill ADC distributions of the samples for all the StexEchas
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvTno()
{
  if (fTagAvTno[0] == 1 )
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_totn[j0StexStinEcna];
    }
    }
}
//-------------------------------------------------------------------------
//
//  Prepa Fill ADC distributions xmin of the samples for all the StexEchas
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvLfn()
{
  if (fTagAvLfn[0] == 1 )
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_lofn[j0StexStinEcna];  
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill ADC distributions xmax of the samples for all the StexEchas
//                       (for writing in the ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvHfn()
{
  if (fTagAvHfn[0] == 1 )
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_hifn[j0StexStinEcna];  
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill Low Frequency covariances between StexEchas
//                           (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootLfCov()
{
  if (fTagLfCov[0] == 1 )
    {
      for (Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexEcha, j0StexEcha) =
            fT2d_lf_cov[i0StexEcha][j0StexEcha];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill Low Frequency correlations between StexEchas
//                         (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootLfCor()
{
  if (fTagLfCor[0] == 1 )
    {
      for (Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexEcha, j0StexEcha) =
            fT2d_lf_cor[i0StexEcha][j0StexEcha];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill High Frequency covariances between StexEchas
//                           (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootHfCov()
{
  if (fTagHfCov[0] == 1 )
    {
      for (Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexEcha, j0StexEcha) =
            fT2d_hf_cov[i0StexEcha][j0StexEcha];
        }
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill High Frequency correlations between StexEchas
//                         (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootHfCor()
{
  if (fTagHfCor[0] == 1 )
    {
      for (Int_t i0StexEcha = 0; i0StexEcha < fEcal->MaxCrysEcnaInStex(); i0StexEcha++)
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
        {
          gCnaRootFile->fCnaIndivResult->fMatMat(i0StexEcha, j0StexEcha) =
            fT2d_hf_cor[i0StexEcha][j0StexEcha];
        }
    }
    }
}


//-------------------------------------------------------------------------
//
//  Prepa Fill covariances between samples for a given StexEcha
//                      (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootCovCss(const Int_t& user_StexEcha, const Int_t& argNbSampWrite)
{
   if (fTagCovCss[user_StexEcha] == 1 )
     {
       for (Int_t i0Sample = 0; i0Sample < argNbSampWrite; i0Sample++)
     {
       for (Int_t j0Sample = 0; j0Sample < argNbSampWrite; j0Sample++)
         {
           gCnaRootFile->fCnaIndivResult->fMatMat(i0Sample, j0Sample) =
             fT3d_cov_ss[user_StexEcha][i0Sample][j0Sample];
         }
     }
     }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill correlations between samples for a given StexEcha
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootCorCss(const Int_t& user_StexEcha, const Int_t& argNbSampWrite)
{
   if (fTagCorCss[user_StexEcha] == 1 )
     {
       for (Int_t i0Sample = 0; i0Sample < argNbSampWrite; i0Sample++)
     {
       for (Int_t j0Sample = 0; j0Sample < argNbSampWrite; j0Sample++)
         {
           gCnaRootFile->fCnaIndivResult->fMatMat(i0Sample, j0Sample) =
             fT3d_cor_ss[user_StexEcha][i0Sample][j0Sample];
         }
     }
     }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill expectation values of the expectation values of the samples
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootPed()
{
  if (fTagPed[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_ev_ev[j0StexEcha];
    }      
    }
}
//-------------------------------------------------------------------------
//
//  Prepa Fill expectation values of the sigmas of the samples
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootTno()
{
  if (fTagTno[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_evsamp_of_sigevt[j0StexEcha];
    }      
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill expectation values of the (sample,sample) correlations
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootMeanCorss()
{
  if (fTagMeanCorss[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_ev_cor_ss[j0StexEcha];
    }      
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill sigmas of the expectation values of the samples
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootLfn()
{
  if (fTagLfn[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_sigevt_of_evsamp[j0StexEcha];
    }      
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill sigmas of the expectation values of the sigmas
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootHfn()
{
  if (fTagHfn[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_evevt_of_sigsamp[j0StexEcha];
    }      
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill sigmas of the (sample,sample) correlations
//  for all the StexEchas
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootSigCorss()
{
  if (fTagSigCorss[0] == 1 )
    {
      for (Int_t j0StexEcha = 0; j0StexEcha < fEcal->MaxCrysEcnaInStex(); j0StexEcha++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexEcha) =
        fT1d_sig_cor_ss[j0StexEcha];
    }      
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill Average Pedestals
//  for all the StexStins
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvPed()
{
  if (fTagAvPed[0] == 1 )
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_mped[j0StexStinEcna];
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill 
//  
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvEvCorss()
{
  if (fTagAvMeanCorss[0] == 1 )   // test 1st elt only since global calc
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_ev_corss[j0StexStinEcna];
    }
    }
}

//-------------------------------------------------------------------------
//
//  Prepa Fill 
//  
//                        (for writing in ROOT file)
//
//-------------------------------------------------------------------------
void TEcnaRun::TRootAvSigCorss()
{
  if (fTagAvSigCorss[0] == 1 )   // test 1st elt only since global calc
    {
      for (Int_t j0StexStinEcna = 0; j0StexStinEcna < fEcal->MaxStinEcnaInStex(); j0StexStinEcna++)
    {
      gCnaRootFile->fCnaIndivResult->fMatHis(0, j0StexStinEcna) =
        fT1d_av_sig_corss[j0StexStinEcna];
    } 
    }
}

//=========================================================================
//
//         METHODS TO SET FLAGS TO PRINT (OR NOT) COMMENTS (DEBUG)
//
//=========================================================================

void TEcnaRun::PrintComments()
{
// Set flags to authorize printing of some comments concerning initialisations (default)

  fFlagPrint = fCodePrintComments;
  std::cout << "*TEcnaRun::PrintComments()> Warnings and some comments on init will be printed" << std::endl;
}

void TEcnaRun::PrintWarnings()
{
// Set flags to authorize printing of warnings

  fFlagPrint = fCodePrintWarnings;
  std::cout << "*TEcnaRun::PrintWarnings()> Warnings will be printed" << std::endl;
}

void TEcnaRun::PrintAllComments()
{
// Set flags to authorize printing of the comments of all the methods

  fFlagPrint = fCodePrintAllComments;
  std::cout << "*TEcnaRun::PrintAllComments()> All the comments will be printed" << std::endl;
}

void  TEcnaRun::PrintNoComment()
{
// Set flags to forbid the printing of all the comments

  fFlagPrint = fCodePrintNoComment;
}
//=========================== E N D ======================================