/data/refman/pasoursint/CMSSW_5_2_7_hltpatch1/src/CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaWrite.h

Go to the documentation of this file.

#ifndef ZTR_TEcnaWrite
#define ZTR_TEcnaWrite

#include <Riostream.h>

#include "TObject.h"
#include "TSystem.h"
#include "Riostream.h"
#include <math.h>
#include <time.h>

#include "TVectorD.h"
#include "TMatrixD.h"

#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaObject.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaNumbering.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParEcal.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParPaths.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaParCout.h"
#include "CalibCalorimetry/EcalCorrelatedNoiseAnalysisAlgos/interface/TEcnaHeader.h"


class TEcnaWrite : public TObject {

 private:

  //..... Attributes

  Int_t fgMaxCar;   // Max nb of caracters for char*
  Int_t fCnew,     fCdelete;
  Int_t fCnewRoot, fCdeleteRoot;

  TString fTTBELL;

  Int_t fCnaCommand, fCnaError;

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

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

  TEcnaParEcal *fEcal;
  TString fFlagSubDet;

  TEcnaNumbering *fEcalNumbering;
  TEcnaParPaths  *fCnaParPaths;
  TEcnaParCout   *fCnaParCout;

  ofstream fFcout_f;

  //...................................... Codes for file names
  Int_t fCodeHeaderAscii;
  Int_t fCodeRoot;

  Int_t fCodeNbOfEvts;
  Int_t fCodePed;
  Int_t fCodeTno;
  Int_t fCodeLfn;
  Int_t fCodeHfn;
  Int_t fCodeMeanCorss;
  Int_t fCodeSigCorss;

  Int_t fCodeCovCss;
  Int_t fCodeCorCss;
  
  Int_t fCodeAdcEvt;
  Int_t fCodeMSp;  
  Int_t fCodeSSp;
  
  Int_t fCodeAvPed;  
  Int_t fCodeAvTno;
  Int_t fCodeAvMeanCorss;
  Int_t fCodeAvSigCorss;
  
  Int_t fCodeLfCov;
  Int_t fCodeLfCor;
  Int_t fCodeHfCov;
  Int_t fCodeHfCor;
  
  Int_t fCodeLFccMoStins;
  Int_t fCodeHFccMoStins;
  
  //..........................................................................
  Int_t fSectChanSizeX, fSectChanSizeY;
  Int_t fSectSampSizeX, fSectSampSizeY;

  Int_t fNbChanByLine; // Nb channels by line (for ASCII results file)
  Int_t fNbSampByLine; // Nb samples by line  (for ASCII results file)
  Int_t fUserSamp;     // Current sample  number (for ASCII results file)
  Int_t fStexStinUser; // Current Stin number in Stex 
  Int_t fStinEchaUser; // Current electronic channel number in Stin
                       // for ASCII results file ([0,24] for EB, [1,25] for EE)


  Double_t**  fjustap_2d_ev;
  Double_t*   fjustap_1d_ev;

  Double_t**  fjustap_2d_var;
  Double_t*   fjustap_1d_var;

  Double_t**  fjustap_2d_cc;
  Double_t*   fjustap_1d_cc;

  Double_t**  fjustap_2d_ss;
  Double_t*   fjustap_1d_ss;

  //.................... Private methods

  void     fAsciiFileWriteHeader(const Int_t&);
  void     fT2dWriteAscii(const Int_t&, const Int_t&, const Int_t&, const Int_t&, const TMatrixD&);

 public:

  //..... Public attributes

  TString  fAnaType;
  Int_t    fNbOfSamples;
  Int_t    fRunNumber;
  Int_t    fFirstReqEvtNumber;
  Int_t    fLastReqEvtNumber;
  Int_t    fReqNbOfEvts;
  Int_t    fStexNumber;
  TString  fStexName;
  TString  fStinName;

  TString  fPathForAsciiFiles;
  TString  fStartDate, fStopDate;
  time_t   fStartTime, fStopTime;

  TString  fRootFileNameShort;  // name of  the results ROOT file 
  TString  fRootFileName;       // name of  the results ROOT file with its path = fPathRoot/fRootFileNameShort

  TString  fAsciiFileName;      // name of  the results ASCII file
  TString  fAsciiFileNameShort; // name of  the results ASCII file = fPathAscii/fAsciiFileNameShort

  //..... Methods
  
  TEcnaWrite();
  TEcnaWrite(TEcnaObject*, const TString);
  TEcnaWrite(const TString, const TEcnaParPaths*, const TEcnaParCout*,
             const TEcnaParEcal*, const TEcnaNumbering*);

  virtual  ~TEcnaWrite();
  
  void     Init();
  void     SetEcalSubDetector(const TString);
  void     SetEcalSubDetector(const TString, const TEcnaParEcal*, const TEcnaNumbering*);
  
  //...................................................... making file name method
  void     fMakeResultsFileName();                // => default: arg = fCodeRoot
  void     fMakeResultsFileName(const Int_t&);

  //.................................... ASCII writing file methods 
  // void    WriteAsciiSampleMeans();    // methode a remettre ?
  // void    WriteAsciiSampleSigmas();   // methode a remettre ?

  void WriteAsciiCovariancesBetweenSamples(const Int_t&, const Int_t&, const Int_t&, const TMatrixD&);
  void WriteAsciiCorrelationsBetweenSamples(const Int_t&, const Int_t&, const Int_t&, const TMatrixD&);

  void WriteAsciiHisto(const TString, const Int_t&, const TVectorD&);

  //...........................................................................
  TString GetAsciiFileName();
  TString GetRootFileName();
  TString GetRootFileNameShort();
  TString GetAnalysisName();

  Int_t GetNbOfSamples();
  Int_t GetRunNumber();
  Int_t GetFirstReqEvtNumber();
  Int_t GetReqNbOfEvts();
  Int_t GetStexNumber();

  Int_t NumberOfEventsAnalysis(Int_t**, const Int_t&, const Int_t&, const Int_t&); // Called by TEcnaRun
  Int_t NumberOfEventsAnalysis(Int_t*,  const Int_t&, const Int_t&, const Int_t&); // Called by TEcnaRead

  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&);
  
  ClassDef(TEcnaWrite,1)// Writing in file (.ascii, .root) methods for CNA (Correlated Noises Analysis)
    };
    
#endif   //    ZTR_TEcnaParameter