TEcnaWrite Class Reference

#include <TEcnaWrite.h>

Inheritance diagram for TEcnaWrite:

Public Member Functions
void	fMakeResultsFileName ()
void	fMakeResultsFileName (const Int_t &)
TString	GetAnalysisName ()
TString	GetAsciiFileName ()
Int_t	GetFirstReqEvtNumber ()
Int_t	GetNbOfSamples ()
Int_t	GetReqNbOfEvts ()
TString	GetRootFileName ()
TString	GetRootFileNameShort ()
Int_t	GetRunNumber ()
Int_t	GetStexNumber ()
void	Init ()
Int_t	NumberOfEventsAnalysis (Int_t **, const Int_t &, const Int_t &, const Int_t &)
Int_t	NumberOfEventsAnalysis (Int_t *, const Int_t &, const Int_t &, const Int_t &)
void	RegisterFileParameters (const TString &, const Int_t &, const Int_t &, const Int_t &, const Int_t &, const Int_t &, const Int_t &, const TString &, const TString &, const time_t, const time_t)
void	RegisterFileParameters (const TString &, const Int_t &, const Int_t &, const Int_t &, const Int_t &, const Int_t &, const Int_t &)
void	SetEcalSubDetector (const TString &)
void	SetEcalSubDetector (const TString &, const TEcnaParEcal *, const TEcnaNumbering *)
	TEcnaWrite ()
	TEcnaWrite (TEcnaObject *, const TString &)
	TEcnaWrite (const TString &, const TEcnaParPaths *, const TEcnaParCout *, const TEcnaParEcal *, const TEcnaNumbering *)
void	WriteAsciiCorrelationsBetweenSamples (const Int_t &, const Int_t &, const Int_t &, const TMatrixD &)
void	WriteAsciiCovariancesBetweenSamples (const Int_t &, const Int_t &, const Int_t &, const TMatrixD &)
void	WriteAsciiHisto (const TString &, const Int_t &, const TVectorD &)
virtual	~TEcnaWrite ()

Public Attributes
TString	fAnaType
TString	fAsciiFileName
TString	fAsciiFileNameShort
Int_t	fFirstReqEvtNumber
Int_t	fLastReqEvtNumber
Int_t	fNbOfSamples
TString	fPathForAsciiFiles
Int_t	fReqNbOfEvts
TString	fRootFileName
TString	fRootFileNameShort
Int_t	fRunNumber
TString	fStartDate
time_t	fStartTime
TString	fStexName
Int_t	fStexNumber
TString	fStinName
TString	fStopDate
time_t	fStopTime

Private Member Functions
void	fAsciiFileWriteHeader (const Int_t &)
void	fT2dWriteAscii (const Int_t &, const Int_t &, const Int_t &, const Int_t &, const TMatrixD &)

Private Attributes
Int_t	fCdelete
Int_t	fCdeleteRoot
Int_t	fCnaCommand
Int_t	fCnaError
TEcnaParCout *	fCnaParCout
TEcnaParPaths *	fCnaParPaths
Int_t	fCnew
Int_t	fCnewRoot
Int_t	fCodeAdcEvt
Int_t	fCodeAvMeanCorss
Int_t	fCodeAvPed
Int_t	fCodeAvSigCorss
Int_t	fCodeAvTno
Int_t	fCodeCorCss
Int_t	fCodeCovCss
Int_t	fCodeHeaderAscii
Int_t	fCodeHFccMoStins
Int_t	fCodeHfCor
Int_t	fCodeHfCov
Int_t	fCodeHfn
Int_t	fCodeLFccMoStins
Int_t	fCodeLfCor
Int_t	fCodeLfCov
Int_t	fCodeLfn
Int_t	fCodeMeanCorss
Int_t	fCodeMSp
Int_t	fCodeNbOfEvts
Int_t	fCodePed
Int_t	fCodePrintAllComments
Int_t	fCodePrintComments
Int_t	fCodePrintNoComment
Int_t	fCodePrintWarnings
Int_t	fCodeRoot
Int_t	fCodeSigCorss
Int_t	fCodeSSp
Int_t	fCodeTno
TEcnaParEcal *	fEcal
TEcnaNumbering *	fEcalNumbering
std::ofstream	fFcout_f
Int_t	fFlagPrint
TString	fFlagSubDet
Int_t	fgMaxCar
Double_t *	fjustap_1d_cc
Double_t *	fjustap_1d_ev
Double_t *	fjustap_1d_ss
Double_t *	fjustap_1d_var
Double_t **	fjustap_2d_cc
Double_t **	fjustap_2d_ev
Double_t **	fjustap_2d_ss
Double_t **	fjustap_2d_var
Int_t	fNbChanByLine
Int_t	fNbSampByLine
Int_t	fSectChanSizeX
Int_t	fSectChanSizeY
Int_t	fSectSampSizeX
Int_t	fSectSampSizeY
Int_t	fStexStinUser
Int_t	fStinEchaUser
TString	fTTBELL
Int_t	fUserSamp

Detailed Description

---

TEcnaWrite.h Update: 06/04/2011 Authors: B.Fabbro (berna.nosp@m.rd.f.nosp@m.abbro.nosp@m.@cea.nosp@m..fr), FX Gentit DSM/IRFU/SPP CEA-Saclay Copyright: Those valid for CEA sofware

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

cna_new/Correlated_Noise_Analysis/ECNA_cna_1.htm

Definition at line 34 of file TEcnaWrite.h.

Constructor & Destructor Documentation

TEcnaWrite::TEcnaWrite

(

)

Definition at line 35 of file TEcnaWrite.cc.

{

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

  Init();
}

TEcnaWrite::TEcnaWrite	(	TEcnaObject *	pObjectManager,
		const TString &	SubDet
	)

Definition at line 43 of file TEcnaWrite.cc.

References TEcnaObject::GetPointerValue(), and TEcnaObject::RegisterPointer().

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

  Init();
  Long_t i_this = (Long_t)this;
  pObjectManager->RegisterPointer("TEcnaWrite", i_this);

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

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

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

  fCnaParPaths->GetPathForResultsRootFiles();
  fCnaParPaths->GetPathForResultsAsciiFiles();

  //............................ fEcal  => should be changed in fParEcal
  fEcal = 0;
  Long_t iParEcal = pObjectManager->GetPointerValue("TEcnaParEcal");
  if( iParEcal == 0 )
    {fEcal = new TEcnaParEcal(pObjectManager, SubDet.Data()); /*fCnew++*/}
  else
    {fEcal = (TEcnaParEcal*)iParEcal;}

  //............................ fEcalNumbering
  fEcalNumbering = 0;
  Long_t iEcalNumbering = pObjectManager->GetPointerValue("TEcnaNumbering");
  if( iEcalNumbering == 0 )
    {fEcalNumbering = new TEcnaNumbering(pObjectManager, SubDet.Data()); /*fCnew++*/}
  else
    {fEcalNumbering = (TEcnaNumbering*)iEcalNumbering;}

  SetEcalSubDetector(SubDet.Data());
}

TEcnaWrite::TEcnaWrite	(	const TString &	SubDet,
		const TEcnaParPaths *	pCnaParPaths,
		const TEcnaParCout *	pCnaParCout,
		const TEcnaParEcal *	pEcal,
		const TEcnaNumbering *	pEcalNumbering
	)

Definition at line 92 of file TEcnaWrite.cc.

References TEcnaParPaths::GetPathForResultsRootFiles().

{

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

  Init();

  //----------------------------- Object management
  fCnaParPaths = 0;
  if( pCnaParPaths == 0 )
    {fCnaParPaths = new TEcnaParPaths();  /*fCnew++*/ ;}
  else
    {fCnaParPaths = (TEcnaParPaths*)pCnaParPaths;}

  //................. Get paths from ECNA directory
  fCnaParPaths->GetPathForResultsRootFiles();
  fCnaParPaths->GetPathForResultsAsciiFiles();

  fCnaParCout  = 0;
  if( pCnaParCout == 0 )
    {fCnaParCout  = new TEcnaParCout();   /*fCnew++*/;}
  else
    {fCnaParCout = (TEcnaParCout*)pCnaParCout;}

  fEcal = 0;
  if( pEcal == 0 )
    {fEcal = new TEcnaParEcal(SubDet.Data());  /*fCnew++*/;}
  else
    {fEcal = (TEcnaParEcal*)pEcal;}

  fEcalNumbering = 0;
  if( pEcalNumbering == 0 )
    {fEcalNumbering = new TEcnaNumbering(SubDet.Data(), fEcal);  /*fCnew++*/;}
  else
    {fEcalNumbering = (TEcnaNumbering*)pEcalNumbering;}

  SetEcalSubDetector(SubDet.Data(), fEcal, fEcalNumbering);
}

TEcnaWrite::~TEcnaWrite ( )

virtual

Definition at line 18 of file TEcnaWrite.cc.

{
  //destructor
  
  //if (fEcal          != 0){delete fEcal;          fCdelete++;}
  //if (fEcalNumbering != 0){delete fEcalNumbering; fCdelete++;}
  //if (fCnaParPaths   != 0){delete fCnaParPaths;   fCdelete++;}
  //if (fCnaParCout    != 0){delete fCnaParCout;    fCdelete++;}
  
  // cout << "[Info Management] CLASS: TEcnaWrite.         DESTROY OBJECT: this = " << this << endl;
}

Member Function Documentation

void TEcnaWrite::fAsciiFileWriteHeader ( const Int_t &  i_code )

private

Definition at line 1051 of file TEcnaWrite.cc.

{
//Ascii results file header writing  (private). Called by the WriteAscii...() methods
  
  //-----------------------------------------------
  //
  //     opening of the ASCII results file
  //     and writing of its header
  //
  //-----------------------------------------------

  if(fAsciiFileName.BeginsWith("$HOME"))
    {
      fAsciiFileName.Remove(0,5);
      TString EndOfAsciiFileName = fAsciiFileName;
      const Text_t *t_file_nohome = (const Text_t *)EndOfAsciiFileName.Data(); //  const Text_t* -> EndOfAsciiFileName
      // ( /scratch0/cna/... )
      TString home_path = gSystem->Getenv("HOME");
      fAsciiFileName = home_path;             // fAsciiFileName = absolute HOME path ( /afs/cern.ch/u/USER )
      fAsciiFileName.Append(t_file_nohome);   // Append  (const Text_t* -> EndOfAsciiFileName) to fAsciiFileName
      // ( /afs/cern.ch/u/USER/scratch0/cna/... )
    }

  fFcout_f.open(fAsciiFileName.Data());

  fFcout_f << "*** File: " << fAsciiFileName
       << " *** " << endl << endl;
  fFcout_f << "*Analysis name                : " << fAnaType       << endl;
  fFcout_f << "*First-Last samples           : 1 - " << fNbOfSamples   << endl; 
  fFcout_f << "*Run number                   : " << fRunNumber     << endl;
  fFcout_f << "*First requested event number : " << fFirstReqEvtNumber << endl;
  fFcout_f << "*Last  requested event number : " << fLastReqEvtNumber << endl;
  fFcout_f << "*Requested number of events   : " << fReqNbOfEvts << endl;
  if( fFlagSubDet == "EB" )
    {fFcout_f << "*SuperModule number           : " << fStexNumber   << endl;}
  if( fFlagSubDet == "EE" )
    {fFcout_f << "*Dee number                   : " << fStexNumber   << endl;}
  fFcout_f << "*Date first requested event   : " << fStartDate;
  fFcout_f << "*Date last requested event    : " << fStopDate      << endl;
  fFcout_f << endl;

  //========================================================================= 
  //   closing of the results file if i_code = fCodeHeaderAscii only.
  //   closing is done in the fT1dWriteAscii() and fT2dWriteAscii() methods
  //   except for i_code = fCodeHeaderAscii
  //=========================================================================
  if(i_code == fCodeHeaderAscii){fFcout_f.close();}
}

void TEcnaWrite::fMakeResultsFileName

(

)

Definition at line 640 of file TEcnaWrite.cc.

{fMakeResultsFileName(fCodeRoot);}

void TEcnaWrite::fMakeResultsFileName

(

const Int_t &

i_code

)

Definition at line 641 of file TEcnaWrite.cc.

References gather_cfg::cout, i, and cuy::ii.

{
//Results filename making  (private)
  
  //----------------------------------------------------------------------
  //
  //     Making of the name of the result file from the parameters set
  //     by call to RegisterFileParameters(...)
  //
  //     Put the names in the following class attributes:
  //
  //     fRootFileName,  fRootFileNameShort,
  //     fAsciiFileName, fAsciiFileNameShort
  //
  //     (Short means: without the directory path)
  //
  //     set indications (run number, type of quantity, ...)
  //     and add the extension ".ascii" or ".root"
  //     
  //     ROOT:  only one ROOT file:  i_code = fCodeRoot.
  //                                          All the types of quantities
  //
  //     ASCII: several ASCII files: i_code = code for one type of quantity
  //            each i_code which is not equal to fCodeRoot is also implicitly
  //            a code "fCodeAscii" (this last attribute is not in the class)
  //     
  //----------------------------------------------------------------------

  char* f_in       = new char[fgMaxCar];                 fCnew++;
  char* f_in_short = new char[fgMaxCar];                 fCnew++;

  Int_t MaxCar = fgMaxCar;
  fStexName.Resize(MaxCar); 
  fStexName  = "SM or Dee?";

  MaxCar = fgMaxCar;
  fStinName.Resize(MaxCar); 
  fStinName  = "tower or SC?";

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

  //  switch (i_code){  

  //===================================  R O O T  =========================  (fMakeResultsFileName)
  TString sPointInterrog = "?";
  TString sDollarHome    = "$HOME";

  if (i_code == fCodeRoot)
    {
      if( fCnaParPaths->ResultsRootFilePath().Data() == sPointInterrog.Data() )
    {
      cout << "!TEcnaWrite::fMakeResultsFileName>  * * * W A R N I N G * * * " << endl << endl
           << "    Path for results .root file not defined. Default option will be used here:" << endl
           << "    your results files will be written in your HOME directory." << endl << endl
           << "    In order to write the .root results file in a specific directory," << endl
               << "    you have to create a file named path_results_root in a subdirectory named ECNA" << endl
           << "    previously created in your home directory." << endl 
           << "    This file must have only one line containing the path of the directory" << endl 
           << "    where must be the .root result files." << endl
           << endl;

      TString home_path = gSystem->Getenv("HOME");
      fCnaParPaths->SetResultsRootFilePath(home_path.Data());  
    }

      if( fCnaParPaths->BeginningOfResultsRootFilePath().Data() == sDollarHome.Data() )
    {
      fCnaParPaths->TruncateResultsRootFilePath(0,5);
      const Text_t *t_file_nohome = (const Text_t *)fCnaParPaths->ResultsRootFilePath().Data(); //  /scratch0/cna/...
      
      TString home_path = gSystem->Getenv("HOME");
      fCnaParPaths->SetResultsRootFilePath(home_path.Data());   //  /afs/cern.ch/u/USER
      fCnaParPaths->AppendResultsRootFilePath(t_file_nohome);   //  /afs/cern.ch/u/USER/scratch0/cna/...
    }

      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d",
          fCnaParPaths->ResultsRootFilePath().Data(), fAnaType.Data(), fNbOfSamples,
          fRunNumber, fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d",
          fAnaType.Data(), fNbOfSamples,
          fRunNumber, fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }

  //===================================  A S C I I  ====================  (fMakeResultsFileName)
  //fCnaParPaths->GetPathForResultsAsciiFiles();
  if (i_code != fCodeRoot)
    {
      if (i_code == fCodeHeaderAscii)
    {
      if( fCnaParPaths->ResultsAsciiFilePath().Data() == sPointInterrog.Data() )
        {
          cout << "!TEcnaWrite::fMakeResultsFileName>  * * * W A R N I N G * * * " << endl << endl
           << "    Path for results .ascii file not defined. Default option will be used here:" << endl
           << "    your results files will be written in your HOME directory." << endl << endl
           << "    In order to write the .ascii results file in a specific directory," << endl
           << "    you have to create a file named path_results_ascii in a subdirectory named ECNA" << endl
           << "    previously created in your home directory." << endl 
           << "    This file must have only one line containing the path of the directory" << endl 
           << "    where must be the .ascii result files." << endl
           << endl;
          
          TString home_path = gSystem->Getenv("HOME");
          fCnaParPaths->SetResultsAsciiFilePath(home_path.Data());  
        } 
      
      if( fCnaParPaths->BeginningOfResultsAsciiFilePath().Data() == sDollarHome.Data() )
        {
          fCnaParPaths->TruncateResultsAsciiFilePath(0,5);
          const Text_t *t_file_nohome = (const Text_t *)fCnaParPaths->ResultsAsciiFilePath().Data(); // /scratch0/cna/...
          
          TString home_path = gSystem->Getenv("HOME");
          fCnaParPaths->SetResultsAsciiFilePath(home_path.Data());   //  /afs/cern.ch/u/USER
          fCnaParPaths->AppendResultsAsciiFilePath(t_file_nohome);   //  /afs/cern.ch/u/USER/scratch0/cna/...
        }
    }

      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_header",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_header",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }

  //--------------------------------------------------------------  (fMakeResultsFileName)
  if (i_code == fCodeNbOfEvts)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_NbOfEvents",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_NbOfEvents",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }

  if (i_code == fCodePed)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_Pedestals",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_Pedestals",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    } 

  if (i_code == fCodeTno)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_TotalNoise",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_TotalNoise",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    } 

  if (i_code == fCodeLfn)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_LFNoise",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_LFNoise",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    } 

  if (i_code == fCodeHfn)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_HFNoise",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_HFNoise",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }

  if (i_code == fCodeMeanCorss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_MeanCorss",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_MeanCorss",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }

  if (i_code == fCodeSigCorss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_SigmaCorss",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_SigmaCorss",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    } 

  if (i_code ==  fCodeCovCss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_Covss_%s%d_Channel_%d",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber,
          fStinName.Data(), fStexStinUser, fStinEchaUser);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_Covss_%s%d_Channel_%d",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber,
          fStinName.Data(), fStexStinUser, fStinEchaUser);
    }
  
  if (i_code == fCodeCorCss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_Corss_%s%d_Channel_%d",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber,
          fStinName.Data(), fStexStinUser, fStinEchaUser);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_Corss_%s%d_Channel_%d",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber,
          fStinName.Data(), fStexStinUser, fStinEchaUser);
    }

  //------- (not used yet)
#define OCOD
#ifndef OCOD
  if (i_code == fCodeMSp)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_SampleMeans",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_SampleMeans",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }
  
  if (i_code == fCodeSSp)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_SampleSigmas",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_SampleSigmas",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }
  
  if( i_code == fCodeAvTno)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageTotalNoise_c%d",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageTotalNoise_c%d",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
    }
  
  if(  i_code == fCodeLfCov)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_LF_cov",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_LF_cov",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }
  
  if ( i_code == fCodeLfCor)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_LF_cor",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_LF_cor",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }
  
  if (i_code == fCodeAvPed)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_AveragePedestals",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_AveragePedestals",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber);
    }
  
  if (i_code == fCodeAvMeanCorss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageMeanCorss%d",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageMeanCorss%d",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
    }

  if (i_code == fCodeAvSigCorss)
    {
      sprintf(f_in, "%s/%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageSigmaCorss%d",
          fCnaParPaths->ResultsAsciiFilePath().Data(), fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
      sprintf(f_in_short, "%s_S1_%d_R%d_%d_%d_%d_%s%d_AverageSigmaCorss%d",
          fAnaType.Data(), fNbOfSamples, fRunNumber,
          fFirstReqEvtNumber, fLastReqEvtNumber, fReqNbOfEvts, fStexName.Data(), fStexNumber, fStinEchaUser);
    } 
#endif // OCOD

  //----------------------------------------------------------- (fMakeResultsFileName)

  // default:
  //    cout << "*TEcnaWrite::fMakeResultsFileName(const Int_t&  i_code)> "
  //     << "wrong header code , i_code = " << i_code << endl; 
  //  }

  //======================================= f_name
  
  char* f_name = new char[fgMaxCar];                   fCnew++;
  
  for (Int_t i = 0 ; i < fgMaxCar ; i++){f_name[i] = '\0';}
  
  Int_t ii = 0;
  for (Int_t i = 0 ; i < fgMaxCar ; i++)
    {
      if ( f_in[i] != '\0' ){f_name[i] = f_in[i]; ii++;}
      else {break;}  // va directement a if ( ii+5 < fgMaxCar ) puis... f_name[ii] = '.';
    }
 
  if ( ii+5 < fgMaxCar )
    {
      //.......... writing of the file extension (.root or .ascii)  (fMakeResultsFileName)
      
      //------------------------------------------- extension .ascii
      if ( i_code != fCodeRoot  || i_code == fCodeNbOfEvts )
    {
      f_name[ii] = '.';   f_name[ii+1] = 'a';
      f_name[ii+2] = 's'; f_name[ii+3] = 'c';
      f_name[ii+4] = 'i'; f_name[ii+5] = 'i';
      
      fAsciiFileName = f_name;
    }
      //------------------------------------------- extension .root
      if ( i_code == fCodeRoot )
    {
      f_name[ii] = '.';   f_name[ii+1] = 'r';
      f_name[ii+2] = 'o'; f_name[ii+3] = 'o';  f_name[ii+4] = 't';
      
      fRootFileName = f_name;
    }
    }
  else
    {
      cout << "*TEcnaWrite::fMakeResultsFileName(...)> Name too long (for f_name)."
       << " No room enough for the extension. (ii = " << ii << ")"
       << fTTBELL << endl; 
    }


  //====================================== f_name_short  (fMakeResultsFileName)
  
  char* f_name_short = new char[fgMaxCar];          fCnew++;

  for (Int_t i = 0 ; i < fgMaxCar ; i++){f_name_short[i] = '\0';}
  
  ii = 0;
  for (Int_t i = 0 ; i < fgMaxCar ; i++)
    {
      if ( f_in_short[i] != '\0' ){f_name_short[i] = f_in_short[i]; ii++;}
      else {break;}  // va directement a f_name_short[ii] = '.';
    }
 
  if ( ii+5 < fgMaxCar )
    {
      //.......... writing of the file extension (.root or .ascii)
      
      //-------------------------------------------extension .ascii
      if ( i_code != fCodeRoot || i_code == fCodeNbOfEvts )
    {
      f_name_short[ii] = '.';   f_name_short[ii+1] = 'a';
      f_name_short[ii+2] = 's'; f_name_short[ii+3] = 'c';
      f_name_short[ii+4] = 'i'; f_name_short[ii+5] = 'i';
      
      fAsciiFileNameShort = f_name_short;
    }
      
      //-------------------------------------------- extension .root
      if ( i_code == fCodeRoot )
    {
      f_name_short[ii] = '.';   f_name_short[ii+1] = 'r';
      f_name_short[ii+2] = 'o'; f_name_short[ii+3] = 'o';  f_name_short[ii+4] = 't';
      
      fRootFileNameShort = f_name_short;
    }
    }
  else
    {
      cout << "*TEcnaWrite::fMakeResultsFileName(...)> Name too long (for f_name_short)."
       << " No room enough for the extension. (ii = " << ii << ")"
       << fTTBELL  << endl; 
    }
    delete [] f_name;        f_name       = 0;         fCdelete++;
    delete [] f_name_short;  f_name_short = 0;         fCdelete++;

    delete [] f_in;          f_in       = 0;           fCdelete++;
    delete [] f_in_short;    f_in_short = 0;           fCdelete++;

}  // end of fMakeResultsFileName

void TEcnaWrite::fT2dWriteAscii ( const Int_t &  i_code, const Int_t &  i_pasx_arg, const Int_t &  i_pasy_arg, const Int_t &  MatSize, const TMatrixD &  read_matrix  )

private

Definition at line 1469 of file TEcnaWrite.cc.

References gather_cfg::cout, i, and j.

{
//Writing of a matrix by sectors (private)

  Int_t i_pasx = i_pasx_arg;           // taille secteur en x
  Int_t i_pasy = i_pasy_arg;           // taille secteur en y
  
  //------------ formatage des nombres en faisant appel a la classe ios
  
  fFcout_f << setiosflags(ios::showpoint | ios::uppercase);
  fFcout_f.setf(ios::dec, ios::basefield);
  fFcout_f.setf(ios::fixed, ios::floatfield);
  fFcout_f.setf(ios::left, ios::adjustfield);
  fFcout_f.setf(ios::right, ios::adjustfield);
  fFcout_f << setprecision(3) << setw(6);
  
  cout << setiosflags(ios::showpoint | ios::uppercase);
  cout.setf(ios::dec, ios::basefield);
  cout.setf(ios::fixed, ios::floatfield);
  cout.setf(ios::left, ios::adjustfield);
  cout.setf(ios::right, ios::adjustfield);
  cout << setprecision(3) << setw(6);  

  //--------------------- fin du formatage standard C++ -------------------

  //-----------------------------------------------------------------------
  //  Reservation dynamique d'un array Double_t** de dimensions
  //  les multiples de 5 juste au-dessus des dimensions de l'array 2D
  //  a ecrire ( array de dimensions
  //  (fEcal->MaxSampADC(),fEcal->MaxSampADC())
  //  (fEcal->MaxCrysEcnaInStex(),fEcal->MaxCrysEcnaInStex()) )
  //-----------------------------------------------------------------------
  // Determination des tailles multiples de fSectChanSizeX ou fSectSampSizeX

#define NOUC
#ifndef NOUC

  //*************** channels *************
  Int_t justap_chan = 0;
  
  if( fEcal->MaxCrysEcnaInStex()%fSectChanSizeX == 0 )
    {justap_chan = fEcal->MaxCrysEcnaInStex();}
  else
    {justap_chan=((fEcal->MaxCrysEcnaInStex()/fSectChanSizeX)+1)*fSectChanSizeX;}
  
  //....................... Allocation fjustap_2d_cc

  if ( i_code == fCodeHfCov || i_code == fCodeHfCor||
       i_code == fCodeLfCov || i_code == fCodeLfCor ){
    if(fjustap_2d_cc == 0)
      {
    //................... Allocation
    fjustap_2d_cc = new Double_t*[justap_chan];              fCnew++;  
    fjustap_1d_cc = new  Double_t[justap_chan*justap_chan];  fCnew++;  
    for(Int_t i = 0 ; i < justap_chan ; i++){
      fjustap_2d_cc[i] = &fjustap_1d_cc[0] + i*justap_chan;}
      }

     //............................... Transfert des valeurs dans fjustap_2d_cc  (=init)
    for(Int_t i = 0; i < fEcal->MaxCrysEcnaInStex(); i++){
      for(Int_t j = 0; j < fEcal->MaxCrysEcnaInStex(); j++){
    if ( i_code == fCodeHfCov ){
      fjustap_2d_cc[i][j] = read_matrix(i,j);}
    if ( i_code == fCodeHfCor ){
      fjustap_2d_cc[i][j] = read_matrix(i,j);}
    if ( i_code == fCodeLfCov ){
      fjustap_2d_cc[i][j] = read_matrix(i,j);}
    if ( i_code == fCodeLfCor ){
      fjustap_2d_cc[i][j] = read_matrix(i,j);}
      }
    }

    //.......................... mise a zero du reste de la matrice (=init)
    for(Int_t i = fEcal->MaxCrysEcnaInStex(); i < justap_chan; i++){
      for(Int_t j = fEcal->MaxCrysEcnaInStex(); j < justap_chan; j++){
    fjustap_2d_cc[i][j] = (Double_t)0.;}}
  }

#endif //NOUC

  //************************************ Samples ***************************
  Int_t justap_samp = 0;
  
  if ( fEcal->MaxSampADC()%fSectSampSizeX == 0 ){
    justap_samp = fEcal->MaxSampADC();}
  else{
    justap_samp=((fEcal->MaxSampADC()/fSectSampSizeX)+1)*fSectSampSizeX;}

  //....................... allocation fjustap_2d_ss

  if (i_code == fCodeCovCss      || i_code == fCodeCorCss ||
      i_code == fCodeAvMeanCorss || i_code == fCodeAvSigCorss){
    if(fjustap_2d_ss == 0)
      {
    //................... Allocation
    fjustap_2d_ss = new Double_t*[justap_samp];              fCnew++;  
    fjustap_1d_ss = new  Double_t[justap_samp*justap_samp];  fCnew++;  
    for(Int_t i = 0 ; i < justap_samp ; i++){
      fjustap_2d_ss[i] = &fjustap_1d_ss[0] + i*justap_samp;}
      }

  //.............................. Transfert des valeurs dans fjustap_2d_ss (=init)
  for(Int_t i = 0; i < fEcal->MaxSampADC(); i++){
    for(Int_t j = 0; j < fEcal->MaxSampADC(); j++){
      if( i_code == fCodeCovCss ){
    fjustap_2d_ss[i][j] = read_matrix(i,j);}
      if( i_code == fCodeCorCss ){
    fjustap_2d_ss[i][j] = read_matrix(i,j);}
      if( i_code == fCodeAvMeanCorss ){
    fjustap_2d_ss[i][j] = read_matrix(i,j);}
      if( i_code == fCodeAvSigCorss ){
    fjustap_2d_ss[i][j] = read_matrix(i,j);}
    }
  }

  //.......................... mise a zero du reste de la matrice (=init)
  for(Int_t i = fEcal->MaxSampADC(); i < justap_samp; i++){
    for(Int_t j = fEcal->MaxSampADC(); j < justap_samp; j++){
      fjustap_2d_ss[i][j] = (Double_t)0.;}}
  }

  //..................... impressions + initialisations selon i_code
  
  Int_t isx_max = 0;    
  Int_t isy_max = 0;

#define COCC
#ifndef COCC
  if(i_code == fCodeHfCov)
    {
      fFcout_f << "Covariance matrix between channels "
           << "for sample number " << fUserSamp;
      isx_max = justap_chan;
      isy_max = justap_chan;
    }
  if(i_code == fCodeHfCor)
    {
      fFcout_f << "*Correlation matrix between channels "
           << "for sample number " << fUserSamp;
      isx_max = justap_chan;
      isy_max = justap_chan;
    }
  
  if(i_code == fCodeLfCov)
    {
      fFcout_f << "Covariance matrix between channels "
           << "averaged on the samples ";
      isx_max = justap_chan;
      isy_max = justap_chan;
    }
  if(i_code == fCodeLfCor)
    {
      fFcout_f << "Correlation matrix between channels "
           << "averaged on the samples ";
      isx_max = justap_chan;
      isy_max = justap_chan;
    }
#endif // COCC

  Int_t n1StexStin = 0; 
  Int_t i0StinEcha = 0;
  Int_t n1StinEcha = 0;
  
  if( fFlagSubDet == "EB")
    {
      n1StexStin = fStexStinUser;
      i0StinEcha = fEcalNumbering->Get0TowEchaFrom0SMEcha(fStinEchaUser);
    }
  if( fFlagSubDet == "EE")
    {
      n1StexStin = fStexStinUser;
      Int_t fStinEchaUser_m = fStinEchaUser-1;
      n1StinEcha = fEcalNumbering->Get1SCEchaFrom0DeeEcha(fStinEchaUser_m);
      
      TString sDeeDir = fEcalNumbering->GetDeeDirViewedFromIP(fStexNumber);
    }
  
  if(i_code == fCodeCovCss)
    {
      if( fFlagSubDet == "EB")
    {
      fFcout_f << "Covariance matrix between samples "
           << "for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << i0StinEcha << ")";
    }
      if( fFlagSubDet == "EE")
    {
      fFcout_f << "Covariance matrix between samples "
           << "for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << n1StinEcha << ")";
    }
      isx_max = justap_samp;
      isy_max = justap_samp;
    }
  if(i_code == fCodeCorCss)
    {
      if( fFlagSubDet == "EB")
    {
      fFcout_f << "Correlation matrix between samples "
           << "for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << i0StinEcha << ")"; 
    }
      if( fFlagSubDet == "EE")
    {
      fFcout_f << "Correlation matrix between samples "
           << "for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << n1StinEcha << ")"; 
    }
      isx_max = justap_samp;
      isy_max = justap_samp;
    }
 
  if(i_code == fCodeAvMeanCorss)
    {
      if( fFlagSubDet == "EB")
    {
      fFcout_f << "Correction factors to the covariances "
           << "between samples for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << i0StinEcha << ")";
    }
      if( fFlagSubDet == "EE")
    {
      fFcout_f << "Correction factors to the covariances "
           << "between samples for channel number " << fStinEchaUser
           << " (" << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << n1StinEcha << ")";
    }
      isx_max = justap_samp;
      isy_max = justap_samp;
    }
  
  if(i_code == fCodeAvSigCorss)
    {
      if( fFlagSubDet == "EB")
    {
      fFcout_f << "Correction factors to the correlations "
           << "between samples for channel number " << fStinEchaUser
           << " ( " << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << i0StinEcha << ")";
    }
      if( fFlagSubDet == "EE")
    {
      fFcout_f << "Correction factors to the correlations "
           << "between samples for channel number " << fStinEchaUser
           << " ( " << fStinName << ": " << n1StexStin
           << " , channel in " << fStinName << ": " << n1StinEcha << ")";
    }
      isx_max = justap_samp;
      isy_max = justap_samp;
    }
  
  fFcout_f << endl;

  //............... Calcul des nombres de secteurs selon x
  //                i_pasx  = taille secteur en x 
  //                isx_max = taille de la matrice en x
  //                n_sctx  = nombre de secteurs en x
  //
  if(i_pasx > isx_max){i_pasx = isx_max;}  
  Int_t n_sctx;
  Int_t max_verix; 
  n_sctx = isx_max/i_pasx;
  max_verix = n_sctx*i_pasx; 
  if(max_verix < isx_max){ n_sctx++;}

  //............... Calcul des nombres de secteurs selon y
  //                i_pasy  = taille secteur en y
  //                isy_max = taille de la matrice en y
  //                n_scty  = nombre de secteurs en x
  //
  if(i_pasy > isy_max){i_pasy = isy_max;} 
  Int_t n_scty;
  Int_t max_veriy; 
  n_scty = isy_max/i_pasy;
  max_veriy = n_scty*i_pasy; 
  if(max_veriy < isy_max){ n_scty++;}

#define NBSC
#ifndef NBSC
  //................ Ecriture de la taille et du nombre des secteurs
  if( i_code == fCodeCovCss      || i_code == fCodeCorCss ||
      i_code == fCodeAvMeanCorss || i_code == fCodeAvSigCorss)
    {
      fFcout_f << "sector size = " << fSectSampSizeX 
           << " , number of sectors = " << n_sctx << " x " << n_scty
           <<endl;
    }
 if( i_code == fCodeHfCov || i_code == fCodeHfCor ||
     i_code == fCodeLfCov || i_code == fCodeLfCor )
    {
      fFcout_f << "sector size = " << fSectChanSizeX 
           << " , number of sectors = " << n_sctx << " x " << n_scty
           << endl;
    }
#endif // NBSC

  fFcout_f << endl;

  //............... impression matrice par secteurs i_pas x i_pas  
  //........................... boucles pour display des secteurs 
  Int_t   ix_inf = -i_pasx;
  
  for(Int_t nsx = 0 ; nsx < n_sctx ; nsx++)
    { 
      //......................... calcul limites secteur
      ix_inf = ix_inf + i_pasx;
      Int_t   ix_sup = ix_inf + i_pasx; 
      
      Int_t   iy_inf = -i_pasy;
      
      for(Int_t nsy = 0 ; nsy < n_scty ; nsy++)
    {       
      iy_inf = iy_inf + i_pasy;          
      Int_t   iy_sup = iy_inf + i_pasy;
      
      //......................... display du secteur (nsx,nsy)
      
      if(i_code == fCodeHfCov || i_code == fCodeCovCss ||
         i_code == fCodeAvMeanCorss || i_code == fCodeAvSigCorss )
        {fFcout_f << "        ";}
      if(i_code == fCodeHfCor || i_code == fCodeCorCss)
        {fFcout_f << "      ";}
      
      for (Int_t iy_c = iy_inf ; iy_c < iy_sup ; iy_c++)
        {
          if(i_code == fCodeHfCov || i_code == fCodeLfCov || i_code == fCodeCovCss ||
         i_code == fCodeAvMeanCorss ||
         i_code == fCodeAvSigCorss)
        {fFcout_f.width(8);}
          if(i_code == fCodeHfCor || i_code == fCodeLfCor || i_code == fCodeCorCss)
        {fFcout_f.width(6);}
          fFcout_f << iy_c << "  ";
        }     
      fFcout_f << endl << endl;
      
      for (Int_t ix_c = ix_inf ; ix_c < ix_sup ; ix_c++)
        { 
          if(i_code == fCodeHfCov|| i_code == fCodeLfCov || i_code == fCodeCovCss ||
         i_code == fCodeAvMeanCorss ||
         i_code == fCodeAvSigCorss)
        {fFcout_f.width(8);}
          if(i_code == fCodeHfCor || i_code == fCodeLfCor || i_code == fCodeCorCss)
        {fFcout_f.width(6);}
          fFcout_f << ix_c << "   ";
          
          for (Int_t iy_c = iy_inf ; iy_c < iy_sup ; iy_c++)
        {
          if(i_code == fCodeHfCov  ||
             i_code == fCodeLfCov  ||
             i_code == fCodeCovCss ||
             i_code == fCodeAvMeanCorss || 
             i_code == fCodeAvSigCorss ){
            fFcout_f.width(8);}

          if(i_code == fCodeHfCor || i_code == fCodeLfCor ||  i_code == fCodeCorCss){
            fFcout_f.width(6);}
          
          if( i_code == fCodeHfCov ||  i_code == fCodeLfCov || i_code == fCodeHfCor){
            fFcout_f << fjustap_2d_cc[ix_c][iy_c] << "  ";}
          
          if ( i_code == fCodeCovCss ||
               i_code == fCodeCorCss ||
               i_code == fCodeAvMeanCorss ||
               i_code == fCodeAvSigCorss )
            {
              fFcout_f << fjustap_2d_ss[ix_c][iy_c] << "  ";
            }
        }
          fFcout_f << endl;
        }
      fFcout_f << endl;
    }
    }

  //........... closing of the results file
  
  fFcout_f.close();
  
  if(fFlagPrint != fCodePrintNoComment)
    {
      cout << "*TEcnaWrite::fT2dWriteAscii(....)> INFO: "
       << "matrix has been written in file: " << endl
       << "            " << fAsciiFileName.Data() << endl;
    }
  
}// end of TEcnaWrite::fT2dWriteAscii

TString TEcnaWrite::GetAnalysisName

(

)

Definition at line 387 of file TEcnaWrite.cc.

{return fAnaType;}

TString TEcnaWrite::GetAsciiFileName

(

)

Definition at line 384 of file TEcnaWrite.cc.

{return fAsciiFileName;}

Int_t TEcnaWrite::GetFirstReqEvtNumber

(

)

Definition at line 390 of file TEcnaWrite.cc.

{return fFirstReqEvtNumber;}

Int_t TEcnaWrite::GetNbOfSamples

(

)

Definition at line 388 of file TEcnaWrite.cc.

{return fNbOfSamples;}

Int_t TEcnaWrite::GetReqNbOfEvts

(

)

Definition at line 391 of file TEcnaWrite.cc.

{return fReqNbOfEvts;}

TString TEcnaWrite::GetRootFileName

(

)

Definition at line 385 of file TEcnaWrite.cc.

{return fRootFileName;}

TString TEcnaWrite::GetRootFileNameShort

(

)

Definition at line 386 of file TEcnaWrite.cc.

{return fRootFileNameShort;}

Int_t TEcnaWrite::GetRunNumber

(

)

Definition at line 389 of file TEcnaWrite.cc.

{return fRunNumber;}

Int_t TEcnaWrite::GetStexNumber

(

)

Definition at line 392 of file TEcnaWrite.cc.

{return fStexNumber;}

void TEcnaWrite::Init

(

void

)

Definition at line 135 of file TEcnaWrite.cc.

{
  //----------------------------- Parameters values
  fTTBELL = '\007';

  fgMaxCar = (Int_t)512;              // max number of characters in TStrings
  fCodeHeaderAscii =  0;
  fCodeRoot        =  1;

  //------------------------------------------- Codes
  //................. presently used codes
  fCodeNbOfEvts    = 101;
  fCodePed         = 102;
  fCodeTno         = 103;
  fCodeLfn         = 104;
  fCodeHfn         = 105;
  fCodeMeanCorss   = 106;
  fCodeSigCorss    = 107;

  fCodeCovCss      = 201;
  fCodeCorCss      = 202;

  //................. not yet used codes 
  //fCodeAdcEvt      =  3;
  //fCodeMSp         =  4;  
  //fCodeSSp         =  5;

  //fCodeAvPed       = 17;
  //fCodeAvTno       =  6;
  //fCodeAvMeanCorss = 18;
  //fCodeAvSigCorss  = 19;

  //fCodeLfCov       = 11;
  //fCodeLfCor       = 12;
  //fCodeHfCov       =  9;
  //fCodeHfCor       = 10;

  //fCodeLFccMoStins = 13;
  //fCodeHFccMoStins = 14;
  
  //----------------------------------------

  fUserSamp     = 0;
  fStexStinUser = 0;
  fStinEchaUser = 0;

  fjustap_2d_ev  = 0;
  fjustap_1d_ev  = 0;

  fjustap_2d_var = 0;
  fjustap_1d_var = 0;

  fjustap_2d_cc  = 0;
  fjustap_1d_cc  = 0;

  fjustap_2d_ss  = 0;
  fjustap_1d_ss  = 0;

  fCodePrintAllComments = fCnaParCout->GetCodePrint("AllComments");

}

Int_t TEcnaWrite::NumberOfEventsAnalysis	(	Int_t **	T2d_NbOfEvts,
		const Int_t &	MaxCrysEcnaInStex,
		const Int_t &	MaxNbOfSamples,
		const Int_t &	NbOfReqEvts
	)

Definition at line 520 of file TEcnaWrite.cc.

References gather_cfg::cout.

{
  //  CHECK OF THE NUMBER OF FOUND EVENTS, return rNumberOfEvents (NumberOfEvents())
  //       (number used to compute the average values over the events)
  //
  //  2D array: called by TEcnaRun object

  Int_t rNumberOfEvents = 0;
  Int_t PresentNumber   = 0;
  Int_t DifferentMinusValue  = 0;
  Int_t DifferentPlusValue  = 0;

  for(Int_t i0StexEcha = 0 ; i0StexEcha < MaxCrysEcnaInStex ; i0StexEcha++)
    {
      for(Int_t i_samp = 0 ; i_samp < MaxNbOfSamples ; i_samp++)
    {
      Int_t NbOfEvts = T2d_NbOfEvts[i0StexEcha][i_samp];

      if( NbOfEvts > 0 )
        {
          if( PresentNumber == 0 )
        {
          PresentNumber = NbOfEvts;
        }
          else
        {
          if( NbOfEvts > PresentNumber )
            { 
              PresentNumber = NbOfEvts;
              DifferentPlusValue++;
            }
          if( NbOfEvts < PresentNumber )
            {
              DifferentMinusValue++;
            }
        }
        }
    }
    }
  //.............................................................  (NumberOfEvents())
  rNumberOfEvents = PresentNumber;

  if( DifferentMinusValue > 0 || DifferentPlusValue > 0 )
    {
      cout << "!TEcnaWrite::NumberOfEventsAnalysis()> " << endl
       << "****************** W A R N I N G  :  NUMBER OF EVENTS NOT CONSTANT !  *********************************"
       << endl
       << "  Result ROOT file: " << fRootFileName << endl
       << "  The number of events is not the same for all the non-empty channels" << endl      
       << "  The maximum number (" << rNumberOfEvents << ") is considered as the number of events " << endl 
       << "  for calculations of pedestals, noises and correlations." << endl
       << "  Number of channels with 0 < nb of evts < " << rNumberOfEvents << " : " << DifferentMinusValue << endl
    // << "  Number of channels with nb of evts > " << rNumberOfEvents << " : " << DifferentPlusValue  << endl
    // << "  Number of empty channels : " << EmptyChannel << endl
       << "  Some values of pedestals, noises and correlations may be wrong for channels" << endl
       << "  with number of events different from " << rNumberOfEvents << "." << endl
       << "  Please, check the histogram 'Numbers of events'." << endl
       << "*******************************************************************************************************"
       << endl;
    }
  else
    {
      if(fFlagPrint == fCodePrintAllComments)
    {
      if( rNumberOfEvents < NbOfReqEvts )
        {
          cout << "*TEcnaWrite::NumberOfEventsAnalysis()> *** INFO *** Number of events found in data = "
           << rNumberOfEvents << ": less than number of requested events ( = " << NbOfReqEvts << ")" << endl;
        }
    }
    }
  return rNumberOfEvents;  
  
}//----- ( end of NumberOfEvents(...) ) ----------------

Int_t TEcnaWrite::NumberOfEventsAnalysis	(	Int_t *	ArrayNbOfEvts,
		const Int_t &	MaxArray,
		const Int_t &	NbOfReqEvts,
		const Int_t &	StexNumber
	)

Definition at line 408 of file TEcnaWrite.cc.

References gather_cfg::cout.

{
  //  CHECK THE NUMBER OF FOUND EVENTS, return rNumberOfEvents (NumberOfEvents())
  //  (number used to compute the average values over the events)
  //
  //  1D array: called by TEcnaRead object

  Int_t rNumberOfEvents = 0;
  Int_t PresentNumber   = 0;
  Int_t EmptyChannel    = 0;
  Int_t DifferentMinusValue = 0;
  Int_t DifferentPlusValue  = 0;

  //........................................................ i_SSoSE = StexStin or StinEcha
  for(Int_t i_SSoSE=0; i_SSoSE<MaxArray; i_SSoSE++)
    {
      Int_t NbOfEvts = ArrayNbOfEvts[i_SSoSE];

      if( NbOfEvts > 0 )
    {
      if( PresentNumber == 0 )   // first channel
        {
          PresentNumber = NbOfEvts;
        }
      else                      // current channel
        {
          if( NbOfEvts > PresentNumber )   // warning
        { 
          PresentNumber = NbOfEvts;
          DifferentPlusValue++;
        }
          if( NbOfEvts < PresentNumber )   // warning
        {
          DifferentMinusValue++;
        }
        }
    }
      else
    {
      EmptyChannel++;
    }
    }

  rNumberOfEvents = PresentNumber;

  if( EmptyChannel > 0 )
    {
      if( MaxArray == fEcal->MaxCrysInSM() )  // (EB)
    {
      cout << "!TEcnaWrite::NumberOfEventsAnalysis()> *** WARNING *** "
           << EmptyChannel << " empty channels detected in SM " << StexNumber
           << " (EB " << fEcalNumbering->PlusMinusSMNumber(StexNumber) << ")" << endl;
    }
      if( MaxArray == fEcal->MaxCrysEcnaInDee() )  // (EE)
    {
      EmptyChannel -= fEcal->EmptyChannelsInDeeMatrixIncompleteSCIncluded();
      if( EmptyChannel > 0 )
        {
          cout << "!TEcnaWrite::NumberOfEventsAnalysis()> *** WARNING *** "
           << EmptyChannel << " empty channels detected in Dee " << StexNumber << endl;
        }
    }
    }

  if( DifferentMinusValue > 0 || DifferentPlusValue > 0 )
    {
      cout << "!TEcnaWrite::NumberOfEventsAnalysis()> " << endl;

      if( MaxArray == fEcal->MaxCrysInSM() )  // (EB)
    {
      cout << "************** W A R N I N G  :  NUMBER OF EVENTS NOT CONSTANT FOR SM " << StexNumber
           << " (EB " << fEcalNumbering->PlusMinusSMNumber(StexNumber) << ")  *********************";
    }

      if( MaxArray == fEcal->MaxCrysEcnaInDee() )  // (EE)
    {
      cout << "****************** W A R N I N G  :  NUMBER OF EVENTS NOT CONSTANT FOR Dee " << StexNumber
           << " **************************";
    }

      cout << endl
       << "  Result ROOT file: " << fRootFileName << endl
       << "  The number of events is not the same for all the non-empty channels." << endl     
       << "  The maximum number (" << rNumberOfEvents << ") is considered as the number of events for calculations "
       << endl  
       << "  of pedestals, noises and correlations." << endl
       << "  Number of channels with 0 < nb of evts < " << rNumberOfEvents << " : " << DifferentMinusValue
       << endl
    // << "  Number of channels with nb of evts > " << rNumberOfEvents << " = " << DifferentPlusValue  << endl
       << "  Number of empty channels : " << EmptyChannel << endl
       << "  Some values of pedestals, noises and correlations may be wrong for channels" << endl
       << "  with number of events different from " << rNumberOfEvents << "." << endl
       << "  Please, check the histogram 'Numbers of events'." << endl
       << "*******************************************************************************************************"
       << endl;
    }
  else
    {
      if(fFlagPrint == fCodePrintAllComments)
    {
      if( rNumberOfEvents < NbOfReqEvts )
        {
          cout << "*TEcnaWrite::NumberOfEventsAnalysis()> *** INFO *** Number of events found in data = "
           << rNumberOfEvents << ": less than number of requested events ( = " << NbOfReqEvts << ")" << endl;
        }
    }
    }
  return rNumberOfEvents;  
}

void TEcnaWrite::RegisterFileParameters	(	const TString &	ArgAnaType,
		const Int_t &	ArgNbOfSamples,
		const Int_t &	ArgRunNumber,
		const Int_t &	ArgFirstReqEvtNumber,
		const Int_t &	ArgLastReqEvtNumber,
		const Int_t &	ArgReqNbOfEvts,
		const Int_t &	ArgStexNumber,
		const TString &	ArgStartDate,
		const TString &	ArgStopDate,
		const time_t	ArgStartTime,
		const time_t	ArgStopTime
	)

Definition at line 610 of file TEcnaWrite.cc.

{
  fAnaType           = ArgAnaType;
  fNbOfSamples       = ArgNbOfSamples;
  fRunNumber         = ArgRunNumber;
  fFirstReqEvtNumber = ArgFirstReqEvtNumber;
  fLastReqEvtNumber  = ArgLastReqEvtNumber;
  fReqNbOfEvts       = ArgReqNbOfEvts;
  fStexNumber        = ArgStexNumber;
  fStartDate         = ArgStartDate;

  fStopDate          = ArgStopDate;
  fStartTime         = ArgStartTime;
  fStopTime          = ArgStopTime;
}

void TEcnaWrite::RegisterFileParameters	(	const TString &	ArgAnaType,
		const Int_t &	ArgNbOfSamples,
		const Int_t &	ArgRunNumber,
		const Int_t &	ArgFirstReqEvtNumber,
		const Int_t &	ArgLastReqEvtNumber,
		const Int_t &	ArgReqNbOfEvts,
		const Int_t &	ArgStexNumber
	)

Definition at line 596 of file TEcnaWrite.cc.

{
  fAnaType           = ArgAnaType;
  fNbOfSamples       = ArgNbOfSamples;
  fRunNumber         = ArgRunNumber;
  fFirstReqEvtNumber = ArgFirstReqEvtNumber;
  fLastReqEvtNumber  = ArgLastReqEvtNumber;
  fReqNbOfEvts       = ArgReqNbOfEvts;
  fStexNumber        = ArgStexNumber;
}

void TEcnaWrite::SetEcalSubDetector

(

const TString &

SubDet

)

Definition at line 204 of file TEcnaWrite.cc.

{
 // Set Subdetector (EB or EE)

  Int_t MaxCar = fgMaxCar;
  fFlagSubDet.Resize(MaxCar);
  fFlagSubDet = fEcal->GetEcalSubDetector();

  //........................................................................
  //
  //             (for ASCII files writing methods only) ...
  //
  //                DEFINITION OF THE SECTOR SIZES
  //       FOR THE CORRELATION AND COVARIANCE MATRICES DISPLAY
  //
  //            MUST BE A DIVISOR OF THE TOTAL NUMBER.
  //            ======================================
  //
  //     Examples:
  //      
  //      (1)    25 channels => size = 25 or 5 (divisors of 25)
  //
  //             25 => matrix = 1 x 1 sector  of size (25 x 25)
  //                          = (1 x 1) x (25 x 25) = 1 x 625 = 625 
  //              5 => matrix = 5 x 5 sectors of size (5 x 5)
  //                          = (5 x 5) x ( 5 x  5) = 25 x 25 = 625 
  //
  //      (2)    10 samples  => size = 10, 5 or 2 (divisors of 10)
  //
  //             10 => matrix = 1 X 1 sectors of size (10 x 10) 
  //                          = (1 x 1) x (10 x 10) =  1 x 100 = 100
  //              5 => matrix = 2 x 2 sectors of size (5 x 5) 
  //                          = (2 x 2) x ( 5 x  5) =  4 x  25 = 100
  //              2 => matrix = 5 x 5 sectors of size (2 x 2) 
  //                          = (5 x 5) x ( 2 x  2) = 25 x  4  = 100
  //
  //........................................................................
  fSectChanSizeX = fEcal->MaxCrysHocoInStin();
  fSectChanSizeY = fEcal->MaxCrysVecoInStin();
  fSectSampSizeX = fEcal->MaxSampADC(); 
  fSectSampSizeY = fEcal->MaxSampADC();

  //........................................................................
  //
  //                DEFINITION OF THE NUMBER OF VALUES BY LINE
  //                for the Expectation Values, Variances and.
  //                Event distributions by (channel,sample)
  //
  //               MUST BE A DIVISOR OF THE TOTAL NUMBER.
  //               ======================================
  //
  //     Examples: 
  //                1) For expectation values and variances:
  //
  //                25 channels => size = 5
  //                => sample sector = 5 lines of 5 values
  //                                 = 5 x 5 = 25 values 
  //
  //                10 samples  => size = 10
  //                => channel sector = 1 lines of 10 values
  //                                  = 1 x 10 = 10 values
  //
  //                2) For event distributions:
  //
  //                100 bins  => size = 10
  //                => sample sector = 10 lines of 10 values
  //                                 = 10 x 10 = 100 values
  //
  //........................................................................
  fNbChanByLine = fEcal->MaxCrysHocoInStin();
  fNbSampByLine = fEcal->MaxSampADC();  
}//---------- (end of SetEcalSubDetector) ------------------

void TEcnaWrite::SetEcalSubDetector	(	const TString &	SubDet,
		const TEcnaParEcal *	pEcal,
		const TEcnaNumbering *	pEcalNumbering
	)

Definition at line 278 of file TEcnaWrite.cc.

References TEcnaParEcal::GetEcalSubDetector().

{
 // Set Subdetector (EB or EE)


  fEcal = 0;
  if( pEcal == 0 )
    {fEcal = new TEcnaParEcal(SubDet.Data());  fCnew++;}
  else
    {fEcal = (TEcnaParEcal*)pEcal;}

  Int_t MaxCar = fgMaxCar;
  fFlagSubDet.Resize(MaxCar);
  fFlagSubDet = fEcal->GetEcalSubDetector();

  fEcalNumbering = 0;
  if( pEcalNumbering == 0 )
    {fEcalNumbering = new TEcnaNumbering(SubDet.Data(), fEcal);  fCnew++;}
  else
    {fEcalNumbering = (TEcnaNumbering*)pEcalNumbering;}

  //........................................................................
  //
  //             (for ASCII files writing methods only) ...
  //
  //                DEFINITION OF THE SECTOR SIZES
  //       FOR THE CORRELATION AND COVARIANCE MATRICES DISPLAY
  //
  //            MUST BE A DIVISOR OF THE TOTAL NUMBER.
  //            ======================================
  //
  //     Examples:
  //      
  //      (1)    25 channels => size = 25 or 5 (divisors of 25)
  //
  //             25 => matrix = 1 x 1 sector  of size (25 x 25)
  //                          = (1 x 1) x (25 x 25) = 1 x 625 = 625 
  //              5 => matrix = 5 x 5 sectors of size (5 x 5)
  //                          = (5 x 5) x ( 5 x  5) = 25 x 25 = 625 
  //
  //      (2)    10 samples  => size = 10, 5 or 2 (divisors of 10)
  //
  //             10 => matrix = 1 X 1 sectors of size (10 x 10) 
  //                          = (1 x 1) x (10 x 10) =  1 x 100 = 100
  //              5 => matrix = 2 x 2 sectors of size (5 x 5) 
  //                          = (2 x 2) x ( 5 x  5) =  4 x  25 = 100
  //              2 => matrix = 5 x 5 sectors of size (2 x 2) 
  //                          = (5 x 5) x ( 2 x  2) = 25 x  4  = 100
  //
  //........................................................................
  fSectChanSizeX = fEcal->MaxCrysHocoInStin();
  fSectChanSizeY = fEcal->MaxCrysVecoInStin();
  fSectSampSizeX = fEcal->MaxSampADC(); 
  fSectSampSizeY = fEcal->MaxSampADC();

  //........................................................................
  //
  //                DEFINITION OF THE NUMBER OF VALUES BY LINE
  //                for the Expectation Values, Variances and.
  //                Event distributions by (channel,sample)
  //
  //               MUST BE A DIVISOR OF THE TOTAL NUMBER.
  //               ======================================
  //
  //     Examples: 
  //                1) For expectation values and variances:
  //
  //                25 channels => size = 5
  //                => sample sector = 5 lines of 5 values
  //                                 = 5 x 5 = 25 values 
  //
  //                10 samples  => size = 10
  //                => channel sector = 1 lines of 10 values
  //                                  = 1 x 10 = 10 values
  //
  //                2) For event distributions:
  //
  //                100 bins  => size = 10
  //                => sample sector = 10 lines of 10 values
  //                                 = 10 x 10 = 100 values
  //
  //........................................................................
  fNbChanByLine = fEcal->MaxCrysHocoInStin();
  fNbSampByLine = fEcal->MaxSampADC();  
}//---------- (end of SetEcalSubDetector) ------------------

void TEcnaWrite::WriteAsciiCorrelationsBetweenSamples	(	const Int_t &	user_StexStin,
		const Int_t &	user_StinEcha,
		const Int_t &	MatSize,
		const TMatrixD &	read_matrix
	)

Definition at line 1435 of file TEcnaWrite.cc.

{
//Writing of the correlations between samples for a given StexEcha in an ASCII file


  if( fFlagSubDet == "EB" ){fStexStinUser = user_StexStin;}
  if( fFlagSubDet == "EE" ){fStexStinUser =
                  fEcalNumbering->GetDeeSCConsFrom1DeeSCEcna(fStexNumber, user_StexStin);}

  if( fFlagSubDet == "EB" ){fStinEchaUser = user_StinEcha;}
  if( fFlagSubDet == "EE" ){fStinEchaUser = user_StinEcha+1;}

  Int_t i_code = fCodeCorCss; // code for correlations between samples
  fMakeResultsFileName(i_code); 
  fAsciiFileWriteHeader(i_code);

  Int_t i_pasx = fSectSampSizeX;
  Int_t i_pasy = fSectSampSizeY;

  fT2dWriteAscii(i_code, i_pasx, i_pasy, MatSize, read_matrix);
}

void TEcnaWrite::WriteAsciiCovariancesBetweenSamples	(	const Int_t &	user_StexStin,
		const Int_t &	user_StinEcha,
		const Int_t &	MatSize,
		const TMatrixD &	read_matrix
	)

Definition at line 1405 of file TEcnaWrite.cc.

{
//Writing of the covariances between samples for a given StexEcha in an ASCII file

  if( fFlagSubDet == "EB" ){fStexStinUser = user_StexStin;}
  if( fFlagSubDet == "EE" ){fStexStinUser =
                  fEcalNumbering->GetDeeSCConsFrom1DeeSCEcna(fStexNumber, user_StexStin);}

  if( fFlagSubDet == "EB" ){fStinEchaUser = user_StinEcha;}
  if( fFlagSubDet == "EE" ){fStinEchaUser = user_StinEcha+1;}

  Int_t i_code = fCodeCovCss;  // code for covariances between samples
  fMakeResultsFileName(i_code);  
  fAsciiFileWriteHeader(i_code);

  Int_t i_pasx = fSectSampSizeX;
  Int_t i_pasy = fSectSampSizeY;

  fT2dWriteAscii(i_code, i_pasx, i_pasy, MatSize, read_matrix);
}

void TEcnaWrite::WriteAsciiHisto	(	const TString &	HistoCode,
		const Int_t &	HisSize,
		const TVectorD &	read_histo
	)

Definition at line 1188 of file TEcnaWrite.cc.

References gather_cfg::cout, and relativeConstraints::value.

{
//Write histo with correspondance CNA-channel <-> xtal number in SM or Dee

  // BuildCrysTable() is called in the method Init() which is called by the constructor

  Int_t i_code = fCodeNbOfEvts;

  //--------------------------------------------------------------------------------
  if( HistoCode == "D_NOE_ChNb" ){i_code = fCodeNbOfEvts;}
  if( HistoCode == "D_Ped_ChNb" ){i_code = fCodePed;}
  if( HistoCode == "D_TNo_ChNb" ){i_code = fCodeTno;}
  if( HistoCode == "D_LFN_ChNb" ){i_code = fCodeLfn;}
  if( HistoCode == "D_HFN_ChNb" ){i_code = fCodeHfn;}
  if( HistoCode == "D_MCs_ChNb" ){i_code = fCodeMeanCorss;}
  if( HistoCode == "D_SCs_ChNb" ){i_code = fCodeSigCorss;}

  fMakeResultsFileName(i_code);    // => Making of the results .ascii file name
  fAsciiFileWriteHeader(i_code);   // => Open of the file associated with stream fFcout_f

  //..................................... format numerical values
  fFcout_f << setiosflags(ios::showpoint | ios::uppercase);
  fFcout_f << setprecision(3) << setw(6);
  fFcout_f.setf(ios::dec, ios::basefield);
  fFcout_f.setf(ios::fixed, ios::floatfield);
  fFcout_f.setf(ios::left, ios::adjustfield);
  fFcout_f.setf(ios::right, ios::adjustfield);
  
  cout << setiosflags(ios::showpoint | ios::uppercase);
  cout << setprecision(3) << setw(6);
  cout.setf(ios::dec, ios::basefield);
  cout.setf(ios::fixed, ios::floatfield);
  cout.setf(ios::left, ios::adjustfield);
  cout.setf(ios::right, ios::adjustfield);

  //........................................................ WriteAsciiHisto
  TString aStexName;
  Int_t MaxCar = fgMaxCar;
  aStexName.Resize(MaxCar); 
  aStexName  = "SM or Dee?";
  
  TString aStinName;     
  MaxCar = fgMaxCar;
  aStinName.Resize(MaxCar); 
  aStinName  = "Tower or SC?";
  
  TString aHoco;
  MaxCar = fgMaxCar;
  aHoco.Resize(MaxCar); 
  aHoco  = "Eta or IX?";
  
  TString aVeco;
  MaxCar = fgMaxCar;
  aVeco.Resize(MaxCar); 
  aVeco  = "Phi or IY?";
  
  TString aSpecifa;     
  MaxCar = fgMaxCar;
  aSpecifa.Resize(MaxCar); 
  aSpecifa  = " ";
  
  TString aSpecifc;     
  MaxCar = fgMaxCar;
  aSpecifc.Resize(MaxCar); 
  aSpecifc  = " ";

  TString aSpecifd;     
  MaxCar = fgMaxCar;
  aSpecifd.Resize(MaxCar); 
  aSpecifd  = " ";

  TString aSpecife;     
  MaxCar = fgMaxCar;
  aSpecife.Resize(MaxCar); 
  aSpecife  = " ";
    
  TString aSpecif1;     
  MaxCar = fgMaxCar;
  aSpecif1.Resize(MaxCar); 
  aSpecif1  = " ";
  
  TString aSpecif2;
  MaxCar = fgMaxCar;
  aSpecif2.Resize(MaxCar); 
  aSpecif2  = " ";
  
  if( fFlagSubDet == "EB" )
    {
      aStexName = "SM   "     ; aStinName = "tower"     ; aSpecifa = " channel# ";
      aHoco     = "   Eta    "; aVeco     = "   Phi    "; aSpecifc = " channel# ";
      aSpecifd  = " crystal# "; aSpecife  = "SM    ";
    }
  if( fFlagSubDet == "EE" )
    {
      aStexName = "Dee  "     ; aStinName = " SC  "     ; aSpecifa = "  Sector# ";
      aHoco     = "   IX     "; aVeco     = "   IY     "; aSpecifc = " crystal# ";
      aSpecifd  = "   SC #   "; aSpecife  = "Sector";
    }
  
  //.............................................................. WriteAsciiHisto
  for (Int_t i0StexEcha=0; i0StexEcha<HisSize; i0StexEcha++)
    {
      Int_t n1StexStin   = 0;
      Int_t StexStinEcna = 0;
      Int_t i0StinEcha   = 0;
      Int_t n1StinEcha   = 0;
      Int_t n1StexCrys   = 0;
      Int_t n1DataSector = 0;
      Int_t n1SCinDS     = 0;

      if( fFlagSubDet == "EB")
    {
      n1StexStin   = fEcalNumbering->Get1SMTowFrom0SMEcha(i0StexEcha);
      StexStinEcna = n1StexStin;
      i0StinEcha   = fEcalNumbering->Get0TowEchaFrom0SMEcha(i0StexEcha);
      n1StexCrys   = fEcalNumbering->Get1SMCrysFrom1SMTowAnd0TowEcha(n1StexStin, i0StinEcha);
    }
      if( fFlagSubDet == "EE")
    {
      StexStinEcna = fEcalNumbering->Get1DeeSCEcnaFrom0DeeEcha(i0StexEcha);
      n1DataSector = fEcalNumbering->GetDSFrom1DeeSCEcna(fStexNumber, StexStinEcna);
      n1SCinDS     = fEcalNumbering->GetDSSCFrom1DeeSCEcna(fStexNumber, StexStinEcna);
      n1StexStin   = fEcalNumbering->GetDeeSCConsFrom1DeeSCEcna(fStexNumber, StexStinEcna);
      n1StinEcha   = fEcalNumbering->Get1SCEchaFrom0DeeEcha(i0StexEcha);
    }
      
      if( n1StexStin > 0 )
    {
      if( (fFlagSubDet == "EB" && i0StinEcha == 0) || (fFlagSubDet == "EE" && n1StinEcha == 1) )
        {
          if( HistoCode == "D_NOE_ChNb" ){aSpecif1 = "Number of"; aSpecif2 = "     events (requested)";}
          if( HistoCode == "D_Ped_ChNb" ){aSpecif1 = "Pedestals"; aSpecif2 = "             ";}
          if( HistoCode == "D_TNo_ChNb" ){aSpecif1 = "   Total "; aSpecif2 = "   noise     ";}
          if( HistoCode == "D_MCs_ChNb" ){aSpecif1 = "    Mean "; aSpecif2 = "   cor(s,s)  ";}
          if( HistoCode == "D_LFN_ChNb" ){aSpecif1 = "   Low Fq"; aSpecif2 = "   noise     ";}
          if( HistoCode == "D_HFN_ChNb" ){aSpecif1 = "  High Fq"; aSpecif2 = "   noise     ";}
          if( HistoCode == "D_SCs_ChNb" ){aSpecif1 = " Sigma of"; aSpecif2 = "   cor(s,s)  ";}
          
          fFcout_f << endl;
      
          fFcout_f << aSpecifa.Data()
               << "  " << aStinName.Data() << "#   "
               << aSpecifc.Data()
               << aSpecifd.Data()
               << aHoco.Data()
               << aVeco.Data()
               << aSpecif1.Data() << endl;
      
          fFcout_f << "  in " << aStexName.Data()
               << "  in " << aStexName.Data()
               << "  in " << aStinName.Data()
               << "  in " << aSpecife.Data()
               << "  in " << aStexName.Data()
               << "  in " << aStexName.Data()
               << aSpecif2.Data()  << endl << endl;
        }
 
      Double_t value = read_histo(i0StexEcha);

      if( fFlagSubDet == "EB" )
        {
          fFcout_f  << setw(7)  << i0StexEcha
            << setw(8)  << n1StexStin
            << setw(11) << i0StinEcha   // (Electronic channel number in tower)
            << setw(10) << n1StexCrys
            << setw(10) << (Int_t)fEcalNumbering->GetEta(fStexNumber, StexStinEcna, i0StinEcha)
            << setw(10) << (Int_t)fEcalNumbering->GetPhiInSM(fStexNumber, StexStinEcna, i0StinEcha);
        }
      if( fFlagSubDet == "EE" )
        {
          Int_t n1StinEcha_m = n1StinEcha-1;
          fFcout_f  << setw(7)  << n1DataSector
            << setw(8)  << n1StexStin
            << setw(11) << n1StinEcha   // (Xtal number for construction in SC)
            << setw(10) << n1SCinDS
            << setw(10) << fEcalNumbering->GetIXCrysInDee(fStexNumber, StexStinEcna, n1StinEcha_m)
            << setw(10) << fEcalNumbering->GetJYCrysInDee(fStexNumber, StexStinEcna, n1StinEcha_m);
        }

      if( HistoCode == "D_NOE_ChNb")
        {
          Int_t ivalue = (Int_t)value;
          fFcout_f  << setw(13) << ivalue;
          fFcout_f  << setw(4) << "(" << setw(6) << fReqNbOfEvts << ")";
        }
      else
        {
          fFcout_f  << setw(13) << value;
        }

      fFcout_f  << endl;
    }
    } // end of loop:  for (Int_t i0StexEcha=0; i0StexEcha<HisSize; i0StexEcha++)

  fFcout_f.close();
  
  // if(fFlagPrint != fCodePrintNoComment)
  // {
      cout << "*TEcnaWrite::WriteAsciiHisto(...)> INFO: "
       << "histo has been written in file: " << endl
       << "            " << fAsciiFileName.Data() << endl;
      // }
} // end of TEcnaWrite::WriteAsciiHisto

Member Data Documentation

TString TEcnaWrite::fAnaType

Definition at line 127 of file TEcnaWrite.h.

TString TEcnaWrite::fAsciiFileName

Definition at line 144 of file TEcnaWrite.h.

TString TEcnaWrite::fAsciiFileNameShort

Definition at line 145 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCdelete

private

Definition at line 41 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCdeleteRoot

private

Definition at line 42 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCnaCommand

private

Definition at line 46 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCnaError

private

Definition at line 46 of file TEcnaWrite.h.

TEcnaParCout* TEcnaWrite::fCnaParCout

private

Definition at line 58 of file TEcnaWrite.h.

TEcnaParPaths* TEcnaWrite::fCnaParPaths

private

Definition at line 57 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCnew

private

Definition at line 41 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCnewRoot

private

Definition at line 42 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeAdcEvt

private

Definition at line 77 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeAvMeanCorss

private

Definition at line 83 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeAvPed

private

Definition at line 81 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeAvSigCorss

private

Definition at line 84 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeAvTno

private

Definition at line 82 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeCorCss

private

Definition at line 75 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeCovCss

private

Definition at line 74 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeHeaderAscii

private

Definition at line 63 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeHFccMoStins

private

Definition at line 92 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeHfCor

private

Definition at line 89 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeHfCov

private

Definition at line 88 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeHfn

private

Definition at line 70 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeLFccMoStins

private

Definition at line 91 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeLfCor

private

Definition at line 87 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeLfCov

private

Definition at line 86 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeLfn

private

Definition at line 69 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeMeanCorss

private

Definition at line 71 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeMSp

private

Definition at line 78 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeNbOfEvts

private

Definition at line 66 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodePed

private

Definition at line 67 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodePrintAllComments

private

Definition at line 51 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodePrintComments

private

Definition at line 51 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodePrintNoComment

private

Definition at line 51 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodePrintWarnings

private

Definition at line 51 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeRoot

private

Definition at line 64 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeSigCorss

private

Definition at line 72 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeSSp

private

Definition at line 79 of file TEcnaWrite.h.

Int_t TEcnaWrite::fCodeTno

private

Definition at line 68 of file TEcnaWrite.h.

TEcnaParEcal* TEcnaWrite::fEcal

private

Definition at line 53 of file TEcnaWrite.h.

TEcnaNumbering* TEcnaWrite::fEcalNumbering

private

Definition at line 56 of file TEcnaWrite.h.

std::ofstream TEcnaWrite::fFcout_f

private

Definition at line 60 of file TEcnaWrite.h.

Int_t TEcnaWrite::fFirstReqEvtNumber

Definition at line 130 of file TEcnaWrite.h.

Int_t TEcnaWrite::fFlagPrint

private

Definition at line 50 of file TEcnaWrite.h.

TString TEcnaWrite::fFlagSubDet

private

Definition at line 54 of file TEcnaWrite.h.

Int_t TEcnaWrite::fgMaxCar

private

Definition at line 40 of file TEcnaWrite.h.

Double_t* TEcnaWrite::fjustap_1d_cc

private

Definition at line 113 of file TEcnaWrite.h.

Double_t* TEcnaWrite::fjustap_1d_ev

private

Definition at line 107 of file TEcnaWrite.h.

Double_t* TEcnaWrite::fjustap_1d_ss

private

Definition at line 116 of file TEcnaWrite.h.

Double_t* TEcnaWrite::fjustap_1d_var

private

Definition at line 110 of file TEcnaWrite.h.

Double_t** TEcnaWrite::fjustap_2d_cc

private

Definition at line 112 of file TEcnaWrite.h.

Double_t** TEcnaWrite::fjustap_2d_ev

private

Definition at line 106 of file TEcnaWrite.h.

Double_t** TEcnaWrite::fjustap_2d_ss

private

Definition at line 115 of file TEcnaWrite.h.

Double_t** TEcnaWrite::fjustap_2d_var

private

Definition at line 109 of file TEcnaWrite.h.

Int_t TEcnaWrite::fLastReqEvtNumber

Definition at line 131 of file TEcnaWrite.h.

Int_t TEcnaWrite::fNbChanByLine

private

Definition at line 98 of file TEcnaWrite.h.

Int_t TEcnaWrite::fNbOfSamples

Definition at line 128 of file TEcnaWrite.h.

Int_t TEcnaWrite::fNbSampByLine

private

Definition at line 99 of file TEcnaWrite.h.

TString TEcnaWrite::fPathForAsciiFiles

Definition at line 137 of file TEcnaWrite.h.

Int_t TEcnaWrite::fReqNbOfEvts

Definition at line 132 of file TEcnaWrite.h.

TString TEcnaWrite::fRootFileName

Definition at line 142 of file TEcnaWrite.h.

TString TEcnaWrite::fRootFileNameShort

Definition at line 141 of file TEcnaWrite.h.

Int_t TEcnaWrite::fRunNumber

Definition at line 129 of file TEcnaWrite.h.

Int_t TEcnaWrite::fSectChanSizeX

private

Definition at line 95 of file TEcnaWrite.h.

Int_t TEcnaWrite::fSectChanSizeY

private

Definition at line 95 of file TEcnaWrite.h.

Int_t TEcnaWrite::fSectSampSizeX

private

Definition at line 96 of file TEcnaWrite.h.

Int_t TEcnaWrite::fSectSampSizeY

private

Definition at line 96 of file TEcnaWrite.h.

TString TEcnaWrite::fStartDate

Definition at line 138 of file TEcnaWrite.h.

time_t TEcnaWrite::fStartTime

Definition at line 139 of file TEcnaWrite.h.

TString TEcnaWrite::fStexName

Definition at line 134 of file TEcnaWrite.h.

Int_t TEcnaWrite::fStexNumber

Definition at line 133 of file TEcnaWrite.h.

Int_t TEcnaWrite::fStexStinUser

private

Definition at line 101 of file TEcnaWrite.h.

Int_t TEcnaWrite::fStinEchaUser

private

Definition at line 102 of file TEcnaWrite.h.

TString TEcnaWrite::fStinName

Definition at line 135 of file TEcnaWrite.h.

TString TEcnaWrite::fStopDate

Definition at line 138 of file TEcnaWrite.h.

time_t TEcnaWrite::fStopTime

Definition at line 139 of file TEcnaWrite.h.

TString TEcnaWrite::fTTBELL

private

Definition at line 44 of file TEcnaWrite.h.

Int_t TEcnaWrite::fUserSamp

private

Definition at line 100 of file TEcnaWrite.h.