beamtime/gsi-aug12/prod/go4/VFTX/TVftxParam.cxx (r4864/r4011)

Go to the documentation of this file.

/* Generated by Together */

#include "TVftxParam.h"

#include "TGo4Version.h"
#if __GO4BUILDVERSION__ > 40502
   #include "go4iostream.h"
#else
   #include "Riostream.h"
#endif
 
// define default values
#define DEBUG_OPTION              0
#define CALIB_HISTOS_ENABLE       0
#define REFERENCE_CHANNEL        -1
#define MAIN_REFERENCE_TDC       -1
#define DIAMOND_CHANNEL          -1
#define COINCIDENCE_MAPPING       0
#define BIT_CONTROL               0
#define FINE_TIME_HISTO_ENABLE    0
#define COARSE_TIME_HIST_ENABLE   0
#define CALIBRATION_ENABLE        0
#define MIN_NB_HITS_CALIB      5000
#define INITIAL_CAL_FILENAME      ""
#define CALIBRATION_OUT_ENABLE    0
#define TIME_HISTO_ENABLE         0
#define TIME_HISTO_SINGLES_ENABLE 0
#define TOT_HISTO_ENABLE          0
#define MULTI_DIST_HISTO_ENABLE   0
#define AUTO_TDC_OFFSET_DET       0

#define SINGLE_CHAN_CALIB_INPUT   0
#define SINGLE_CHAN_CALIB_OUTPUT  0

TVftxParam::TVftxParam(const char* name) : TGo4Parameter(name)
{
   ReadOptions();
}

//***********************************************************
TVftxParam::~TVftxParam()
{
}

Bool_t TVftxParam::ReadOptions(){
   // Loading default values
   uDebug           = DEBUG_OPTION;
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
   {
      uVftxCalHistosEnabled[iIndexVftx] = CALIB_HISTOS_ENABLE;
      iVftxReference1Channel[iIndexVftx] = REFERENCE_CHANNEL;
      iVftxReference2Channel[iIndexVftx] = DIAMOND_CHANNEL;
   }
   iMainReferenceTdc = MAIN_REFERENCE_TDC;
   uVftxCoincidenceMap = COINCIDENCE_MAPPING;
   uBitControl      = BIT_CONTROL;
   uFineHistEnable  = FINE_TIME_HISTO_ENABLE;
   uCoarseTimeHisto = COARSE_TIME_HIST_ENABLE; 
   uEnableCalib     = CALIBRATION_ENABLE;
   iMinNbHitsCalib  = MIN_NB_HITS_CALIB;
   sInitialCalibrationFilename = INITIAL_CAL_FILENAME;
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
      uVftxInitalCalIndex[iIndexVftx] = iIndexVftx;
   uEnableCalibOutput = CALIBRATION_OUT_ENABLE;
   uTimeHistEnable  = TIME_HISTO_ENABLE;
   uTimeHistSinglesEnable = TIME_HISTO_SINGLES_ENABLE;
   uTotHistoEnable  = TOT_HISTO_ENABLE;
   uMultiHitsDistanceHistoEnable = MULTI_DIST_HISTO_ENABLE;

   uAutomaticOffsetEnable = AUTO_TDC_OFFSET_DET;
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
      iAutomaticTdcOffset[iIndexVftx] = 0;
   
   uSingleChannelCalibFilesEnable = SINGLE_CHAN_CALIB_INPUT;
   uSingleChannelCalibFilesOutput = SINGLE_CHAN_CALIB_OUTPUT;

   // Try to read the option filename from the general option files text file
   TString sOptionFileName = "./Vftx_options.txt";
   Bool_t bFilenameFound = kFALSE;
   fstream optionFilenameFile("./Option_files.txt", ios::in);
   // Testing if the file does exist. If not I use default values
   if( optionFilenameFile.fail() == true)
   {
      cout<<"************************************************************************"<<endl;
      cout<<" File with option filename for VFTX not found: using the default value ";
      cout<<sOptionFileName<<endl;
   } // if( optionFilenameFile.fail() == true)
      else
      {
         TString line;

         // Goes through all lines
         while( !optionFilenameFile.eof() )
         {
           line.ReadLine(optionFilenameFile, kFALSE); // Read a line of the option file and set the iterator automatically to the next line

           if(line(0) == '#' || line(0) == '*'||
               (line(0) == '/' && line(1) == '/') )
           {
              //Comment in the text file
              ;
           }
           else if(line.Contains("vftxcalib")==kTRUE)
           {
              sOptionFileName = ((TString)line(11, line.Sizeof()-11 ));
              cout<<"************************************************************************"<<endl;
              cout<<"Using "<<sOptionFileName<<" as source for the options of Vftx step"<<endl;
              bFilenameFound = kTRUE;
              break;
           }
         } // while( !optionFilenameFile.eof() )
         if( optionFilenameFile.eof() && kFALSE == bFilenameFound )
         {
            cout<<"************************************************************************"<<endl;
            cout<<"Field for Vftx step not found in File with option filename => Using "<<sOptionFileName<<endl;
         } // if( optionFilenameFile.eof() && kFALSE == bFilenameFound )
      } // else of if( optionFilenameFile.fail() == true)

   // Read the option file
   fstream setupFile(sOptionFileName, ios::in);
   
   // Testing if the option does exist. If not I use default values.
   if( setupFile.fail() == true)
   {
      cout<<"Error for TVftxParam: Setup file "<<sOptionFileName<<" can not be opened or is missing."<<endl;
      cout<<"=> Using default values!"<<endl;
   }
   else
   {
      TString line;
      
      // Goes through all lines
      while( !setupFile.eof() )
      {
         line.ReadLine(setupFile, kFALSE); // Read a line of the option file and set the iterator automatically to the next line
         if(1 == uDebug)
            cout<<line<<endl;

         if(line(0) == '#' || line(0) == '*')
         {
            //Comment in the text file
            ;
         }
         else if(line.Contains("debug")==kTRUE)
         {
            uDebug = ((TString)line(5,10)).Atoi();
            if(1 == uDebug)
               cout<<"Debug ON!!!!!!"<<endl;
         }
         else if(line.Contains("vftxcalena")==kTRUE)
         {
            line = line(11,3*MAX_FPGA_TDC);
            if( 0 != line.Sizeof() )
            {
               for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
               {
                  uVftxCalHistosEnabled[temp_index] = ((TString)line(0,3)).Atoi();
                  line = line(3,3*MAX_FPGA_TDC);
               } // for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
            } // if( 0 != line.Sizeof() )
         }
         else if(line.Contains("refonechan")==kTRUE)
         {
            line = line(11,3*MAX_FPGA_TDC);
            if( 0 != line.Sizeof() )
            {
               for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
               {
                  iVftxReference1Channel[temp_index] = ((TString)line(0,3)).Atoi();
                  line = line(3,3*MAX_FPGA_TDC);
               } // for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
            } // if( 0 != line.Sizeof() )
         }
         else if(line.Contains("reftwochan")==kTRUE)
         {
            line = line(11,3*MAX_FPGA_TDC);
            if( 0 != line.Sizeof() )
            {
               for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
               {
                  iVftxReference2Channel[temp_index] = ((TString)line(0,3)).Atoi();
                  line = line(3,3*MAX_FPGA_TDC);
               } // for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
            } // if( 0 != line.Sizeof() )
         }
         else if(line.Contains("mainreftdc")==kTRUE)
         {
            iMainReferenceTdc = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test iMainReferenceTdc: "<<iMainReferenceTdc<<endl;
         }
         else if(line.Contains("chcoimap")==kTRUE)
         {
            uVftxCoincidenceMap = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uVftxCoincidenceMap: "<<uVftxCoincidenceMap<<endl;
         }
         else if(line.Contains("bitcontrol")==kTRUE)
         {
            uBitControl = ((TString)line(10,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uBitControl: "<<uBitControl<<endl;
         }
         else if(line.Contains("coarsetimehisto")==kTRUE)
         {
            uCoarseTimeHisto = ((TString)line(15,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uCoarseTimeHisto: "<<uCoarseTimeHisto<<endl;
         }
         else if(line.Contains("finetime")==kTRUE)
         {
            uFineHistEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uFineHistEnable: "<<uFineHistEnable<<endl;
         }
         else if(line.Contains("finetevo")==kTRUE)
         {
            uFineEvolutionHistEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uFineEvolutionHistEnable: "<<uFineEvolutionHistEnable<<endl;
         }
         else if(line.Contains("calibena")==kTRUE)
         {
            uEnableCalib = ((TString)line(9,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uEnableCalib: "<<uEnableCalib<<endl;
         }
         else if(line.Contains("calibnb")==kTRUE)
         {
            iMinNbHitsCalib = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test iMinNbHitsCalib: "<<iMinNbHitsCalib<<endl;
         }
         else if(line.Contains("initialcalibfile")==kTRUE)
         {
            sInitialCalibrationFilename = (TString)line(17,128);
            if(1 == uDebug)
               cout<<"Test sInitialCalibrationFilename: "<<sInitialCalibrationFilename<<endl;
         } 
         else if(line.Contains("inicalindex")==kTRUE)
         {
            line = line(12,3*MAX_FPGA_TDC);
            if( 0 != line.Sizeof() )
            {
               for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
               {
                  uVftxInitalCalIndex[temp_index] = ((TString)line(0,3)).Atoi();
                  line = line(3,3*MAX_FPGA_TDC);
               } // for(UInt_t temp_index = 0; temp_index < MAX_FPGA_TDC && 0 != line.Sizeof(); temp_index++)
            } // if( 0 != line.Sizeof() )
         }
         else if(line.Contains("caliboutena")==kTRUE)
         {
            uEnableCalibOutput = ((TString)line(12,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uEnableCalibOutput: "<<uEnableCalibOutput<<endl;
         }
         else if(line.Contains("timeplot")==kTRUE)
         {
            uTimeHistEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uTimeHistEnable: "<<uTimeHistEnable<<endl;
         }
         else if(line.Contains("timesing")==kTRUE)
         {
            uTimeHistSinglesEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uTimeHistSinglesEnable: "<<uTimeHistSinglesEnable<<endl;
         }
         else if(line.Contains("totplot")==kTRUE)
         {
            uTotHistoEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uTotHistoEnable: "<<uTotHistoEnable<<endl;
         }
         else if(line.Contains("muldist")==kTRUE)
         {
            uMultiHitsDistanceHistoEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uMultiHitsDistanceHistoEnable: "<<uMultiHitsDistanceHistoEnable<<endl;
         }
         else if(line.Contains("autooffs")==kTRUE)
         {
            uAutomaticOffsetEnable = ((TString)line(8,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uAutomaticOffsetEnable: "<<uAutomaticOffsetEnable<<endl;
         }
         else if(line.Contains("singchcalib")==kTRUE)
         {
            uSingleChannelCalibFilesEnable = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uSingleChannelCalibFilesEnable: "<<uSingleChannelCalibFilesEnable<<endl;
         }
         else if(line.Contains("schcalibout")==kTRUE)
         {
            uSingleChannelCalibFilesOutput = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
               cout<<"Test uSingleChannelCalibFilesOutput: "<<uSingleChannelCalibFilesOutput<<endl;
         }
      } // while( !setupFile.eof() )
   } //  else of if( setupFile.fail() == true)
   if( uEnableCalib )
   {
      cout<<" Fine time Histos needed for corrections => Enable!!"<<endl;
      uFineHistEnable = 1;

      if( -1 == iMinNbHitsCalib &&
          ( kTRUE == sInitialCalibrationFilename.EqualTo("") ||
            kTRUE == sInitialCalibrationFilename.EqualTo("-") ) )
      {
         cout<<" Nb of hits per channel needed for corrections disabled while calibration file disabled ";
         cout<<"=> Enable iMinNbHitsCalib and set it to 400!!"<<endl;
         iMinNbHitsCalib = 400;
      }
   }
   
   this->PrintOptions(); //print out all values to see, if the setting has worked correctly
   return kTRUE;
}


Int_t TVftxParam::PrintOptions(){
   cout<<"********************* VFTX calibration parameters **********************"<<endl;
   if(1 == uDebug)
      cout<<"Debug ON!!!!!!"<<endl;
      else cout<<"Debug OFF"<<endl;
   cout<<"VFTX FPGA TDC index:       |-      ";
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
   {
      cout.width(2);
      cout<<iIndexVftx<<" ";
   }
   cout<<endl<<"Calibration histo. enabled |->     ";
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
   {
      cout.width(2);
      cout<<uVftxCalHistosEnabled[iIndexVftx] <<" ";
   }
   cout<<endl<<"1st reference channel      |->     ";
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
   {
      cout.width(2);
      cout<<iVftxReference1Channel[iIndexVftx] <<" ";
   }
   cout<<endl<<"2nd reference channel      |->     ";
   for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
   {
      cout.width(2);
      cout<<iVftxReference2Channel[iIndexVftx] <<" ";
   }
   cout<<endl;
   if( -1 < iMainReferenceTdc )
      cout<<"Main reference tdc            :     "<<iMainReferenceTdc<<endl;
      else cout<<"Main reference tdc            :     OFF"<<endl;

   cout<<"Channel coincidence mapping   :     "<<uVftxCoincidenceMap<<endl;
   cout<<"Control bit Histograms        :     "<<uBitControl<<endl; 
   cout<<"Fine Time Histograms          :     "<<uFineHistEnable<<endl;
   cout<<"Fine Time Evolution Histograms:     "<<uFineEvolutionHistEnable<<endl;
   cout<<"Coarse time Histograms        :     "<<uCoarseTimeHisto<<endl;
   if(1 == uEnableCalib)
      cout<<"Calibration                   :     ON"<<endl;
      else cout<<"Calibration                   :     OFF"<<endl;
   if( -1 == iMinNbHitsCalib &&
         kTRUE != sInitialCalibrationFilename.EqualTo("") && kTRUE != sInitialCalibrationFilename.EqualTo("-") )
   {
      cout<<"Filename for calib hist:            "<<sInitialCalibrationFilename<<endl;
      cout<<"VFTX FPGA TDC index:       |-      ";
      for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
      {
         cout.width(2);
         cout<<iIndexVftx<<" ";
      }
      cout<<endl<<"Tdc initial calib. index:  |->     ";
      for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
      {
         cout.width(2);
         cout<<uVftxInitalCalIndex[iIndexVftx] <<" ";
      }
      cout<<endl;
   }
      else
      {
         cout<<"Min nb hits for calib         :     "<<iMinNbHitsCalib<<endl;
         if( kTRUE != sInitialCalibrationFilename.EqualTo("") && kTRUE != sInitialCalibrationFilename.EqualTo("-") )
         {
            cout<<"Filename for intial calib hist:     "<<sInitialCalibrationFilename<<endl;
            cout<<"VFTX FPGA TDC index:       |-      ";
            for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
            {
               cout.width(2);
               cout<<iIndexVftx<<" ";
            }
            cout<<endl<<"Tdc initial calib. index:  |->     ";
            for( Int_t iIndexVftx = 0; iIndexVftx < MAX_FPGA_TDC; iIndexVftx++)
            {
               cout.width(2);
               cout<<uVftxInitalCalIndex[iIndexVftx] <<" ";
            }
            cout<<endl;
         }
      }
   if(1 == uEnableCalibOutput)
      cout<<"Calibration saved in root file:     ON"<<endl;
      else cout<<"Calibration saved in root file:     OFF"<<endl;
   if( 1 == uTimeHistEnable )
      cout<<"Time Histograms               :     ON"<<endl;
      else cout<<"Time Histograms               :     OFF"<<endl;
   if( 1 == uTimeHistSinglesEnable )
      cout<<"Time Histograms Multiplicity 1:     ON"<<endl;
      else cout<<"Time Histograms Multiplicity 1:     OFF"<<endl;
   if( 1 == uTotHistoEnable )
      cout<<"Tot  Histograms               :     ON"<<endl;
      else cout<<"Tot  Histograms               :     OFF"<<endl;
   if( 1 == uMultiHitsDistanceHistoEnable )
      cout<<"Histograms Distance Multi Hits:     ON"<<endl;
      else cout<<"Histograms Distance Multi Hits:     OFF"<<endl;
   if(1 == uAutomaticOffsetEnable)
      cout<<"Automatic tdc offset detection:     ON"<<endl;
      else cout<<"Automatic tdc offset detection:     OFF"<<endl;
   if(1 == uSingleChannelCalibFilesEnable)
      cout<<"Use single channel calib files:     ON"<<endl;
      else cout<<"Use single channel calib files:     OFF"<<endl;
   switch( uSingleChannelCalibFilesOutput )
   {
      case 0:
         cout<<"Single channel calib files output:  OFF"<<endl;
         break;
      case 1:
         cout<<"Single channel calib files output:  New files from data"<<endl;
         break;
      case 2:
         cout<<"Single channel calib files output:  Update files with data"<<endl;
         break;
      default:
         break;
   }

   cout<<"************************************************************************"<<endl;
   
   return 0;
}