beamtime/tof-tdctest/go4/RPC/TRpcParam.cxx (r4864/r4267)

Go to the documentation of this file.

/* Generated by Together */

#include "TRpcParam.h"
#include "TGo4Log.h"

#include "TGo4Version.h"
#if __GO4BUILDVERSION__ > 40502
   #include "go4iostream.h"
#else
   #include "Riostream.h"
#endif

#include "MBSUNPACK/TMbsCrateEvent.h"

#define NB_DIGITS_STRIP_INDEXES  5
#define NB_DIGITS_STRIP_OFFSETS  9

/***** Default values definition *****/
#define DEBUG                 0
#define TDC_FOR_OFFSET_CAEN   -1
#define CH_FOR_OFFSET_CAEN    -1
#define TDC_FOR_OFFSET_VFTX   -1
#define CH_FOR_OFFSET_VFTX    -1

#define TDC_FOR_DIAMOND_CAEN  -1
#define CH_FOR_DIAMOND_CAEN   -1
#define TDC_FOR_DIAMOND_VFTX  -1
#define CH_FOR_DIAMOND_VFTX   -1

#define RATE_INT_DURATION     -1
#define RATE_INT_STEP         -1

#define TDC_STRIP             0
#define RPC_NAME_STRIP        ""
#define RPC_POSX_STRIP        0.0
#define RPC_POSY_STRIP        0.0
#define RPC_POSZ_STRIP        0.0
#define NB_STRIPS             0 // Nb of strips for the beam profile monitor
#define STRIP_WIDTH           0
#define STRIP_LENGTH          0
#define TIME_CONV            10.0
#define MAX_DIST_HITS         5.0
#define REJECT_MULTI_HITS     0
#define ORIENTATION           0
#define BEAM_PROFILE_SIZE     0
#define TRIGGER_SELECTION    -1

#define NOISY_STRIP           0
#define DEAD_STRIP            0
#define STRIPS_LEFT_TOT_TDC   0
#define STRIPS_LEFT_TOT_CHA   0
#define STRIPS_LEFT_TIME_TDC  0
#define STRIPS_LEFT_TIME_CHA  0
#define STRIPS_RIGHT_TIME_TDC 0
#define STRIPS_RIGHT_TIME_CHA 0
#define STRIPS_RIGHT_TOT_TDC  0
#define STRIPS_RIGHT_TOT_CHA  0

#define STRIPS_PROFILE_OFFS   0
#define STRIPS_REFERENCE_OFFS 0
#define REFERENCE_SIGNAL_OFFS 0
#define STRIPS_OFFSET_LEFT    0
#define STRIPS_OFFSET_RIGHT   0

#define USE_CENTER_OF_GRAVITY 0
#define MIN_TIME_BTWN_HITS    0
#define STRIPS_TOT_GAIN_LEFT  1
#define STRIPS_TOT_OFFS_LEFT  0
#define STRIPS_TOT_GAIN_RIGHT 1
#define STRIPS_TOT_OFFS_RIGHT 0

#define INITIAL_CALIB_OFFSET  0

#define MATCHING_OFFSET       0
/*************************************/

TRpcParam::TRpcParam(const char* name) : TGo4Parameter(name)
{
   ReadOptions();
   if( 0 < uNbStrips )
   {
      // not en empty RPC !
      CalculateLeftRightOffset();
      PrintOptions(); //print out all values to see, if the setting has worked correctly
   }
}

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

Bool_t TRpcParam::ReadOptions(){
   // Loading default values
   uDebug                   = DEBUG;
   iOffsetTdcCaen           = TDC_FOR_OFFSET_CAEN;
   iOffsetChannelCaen       = CH_FOR_OFFSET_CAEN;
   iOffsetTdcVFTX           = TDC_FOR_OFFSET_VFTX;
   iOffsetChannelVFTX       = CH_FOR_OFFSET_VFTX;

   iDiamondTdcCaen          = TDC_FOR_DIAMOND_VFTX;
   iDiamondChannelCaen      = CH_FOR_DIAMOND_CAEN;
   iDiamondTdcVFTX          = TDC_FOR_DIAMOND_VFTX;
   iDiamondChannelVFTX      = CH_FOR_DIAMOND_VFTX;

   iRateIntegrationDuration = RATE_INT_DURATION;
   iRateIntegrationStep     = RATE_INT_STEP;

   // Strip RPCs settings
   uTdcTypeStrip            = TDC_STRIP;
   sStripRpcName            = RPC_NAME_STRIP;
   dStripRpcPosX            = RPC_POSX_STRIP;
   dStripRpcPosY            = RPC_POSY_STRIP;
   dStripRpcPosZ            = RPC_POSZ_STRIP;
   uNbStrips                = NB_STRIPS;
   dStripWidth              = STRIP_WIDTH;
   dStripLength             = STRIP_LENGTH;
   dTimeConv                = TIME_CONV;
   dStripHitMaxDistance     = MAX_DIST_HITS;
   uRejectMultipleHits      = REJECT_MULTI_HITS;
   uOrientation             = ORIENTATION;
   iBeamProfileSize         = BEAM_PROFILE_SIZE;
   iTriggerSelection        = TRIGGER_SELECTION;
   uUseChargeGravityCenter  = USE_CENTER_OF_GRAVITY;
   dMinimalTimeBetweenHits  = MIN_TIME_BTWN_HITS;
   dReference1SignalOffset   = REFERENCE_SIGNAL_OFFS;
   dReference2SignalOffset   = REFERENCE_SIGNAL_OFFS;

   for( Int_t iStrip = 0; iStrip < NB_STRIPS_MAX; iStrip ++)
   {
      uTdcStripType[iStrip]          = TDC_STRIP;
      uTdcStripNoisy[iStrip]         = NOISY_STRIP;
      uTdcStripDead[iStrip]          = DEAD_STRIP;
      uTdcStripLeftTot[iStrip]       = STRIPS_LEFT_TOT_TDC;
      uChannelStripLeftTot[iStrip]   = STRIPS_LEFT_TOT_CHA;
      uTdcStripLeftTime[iStrip]      = STRIPS_LEFT_TIME_TDC;
      uChannelStripLeftTime[iStrip]  = STRIPS_LEFT_TIME_CHA;
      uTdcStripRightTime[iStrip]     = STRIPS_RIGHT_TIME_TDC;
      uChannelStripRightTime[iStrip] = STRIPS_RIGHT_TIME_CHA;
      uTdcStripRightTot[iStrip]      = STRIPS_RIGHT_TOT_TDC;
      uChannelStripRightTot[iStrip]  = STRIPS_RIGHT_TOT_CHA;

      dProfOffsetList[iStrip]        = STRIPS_REFERENCE_OFFS;
      dRefOffsetList[iStrip]         = STRIPS_REFERENCE_OFFS;
      dOffsetListLeft[iStrip]        = STRIPS_OFFSET_LEFT;
      dOffsetListRight[iStrip]       = STRIPS_OFFSET_RIGHT;

      dToTGainListLeft[iStrip]       = STRIPS_TOT_GAIN_LEFT;
      dTotOffsetListLeft[iStrip]     = STRIPS_TOT_OFFS_LEFT;
      dToTGainListRight[iStrip]      = STRIPS_TOT_GAIN_RIGHT;
      dTotOffsetListRight[iStrip]    = STRIPS_TOT_OFFS_RIGHT;

      dOffsetMatching[iStrip] = MATCHING_OFFSET;
   }
   uInitialOffsetCalib               = INITIAL_CALIB_OFFSET;
   dLowerBoundFitBeamProfile         = 0;
   dUpperBoundFitBeamProfile         = 0;
   dLowerBoundFitRefProfile          = 0;
   dUpperBoundFitRefProfile          = 0;
   for(Int_t iFitParam = 0; iFitParam < 4 ; iFitParam++)
   {
      dInitialValueFitParBeamProfile[iFitParam] = 0;
      dInitialValueFitParRefProfile[iFitParam]  = 0;
   }
   
   // Read the option file
   TString sOptionFileName = "./Rpc_options.txt";
   TString sName = GetName();

   // Try to read the option filename from the general option files text file
   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 TRpcParam "<<GetName()<<" not found: using the default value ";
      cout<<sOptionFileName<<endl;
   } // if( optionFilenameFile.fail() == true)
      else
      {
         TString line;

         // New
         UInt_t uIndexRpc =((TString)sName( sName.Length()-2,2)).Atoi();

         // 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
              ;
           }
           // New
           else if( kTRUE == line.Contains( Form("rpc%02d ", uIndexRpc) ) )
           {
              sOptionFileName = ((TString)line(11, line.Sizeof()-11 ));
              cout<<"************************************************************************"<<endl;
              cout<<"Using "<<sOptionFileName<<" as source for the options of TRpcParam "<<GetName()<<" step"<<endl;
              bFilenameFound = kTRUE;
              break;
           }
         } // while( !optionFilenameFile.eof() )
         if( optionFilenameFile.eof() && kFALSE == bFilenameFound )
         {
            cout<<"************************************************************************"<<endl;
            cout<<"Field for TRpcParam "<<GetName()<<" step not found in File with option filename => Using "<<sOptionFileName<<endl;
         } // if( optionFilenameFile.eof() && kFALSE == bFilenameFound )
      } // else of if( optionFilenameFile.fail() == true)

   fstream setupFile(sOptionFileName, ios::in);
   
   // Testing if the option does exist. If not I use default values.
   if( setupFile.fail() == true)
   {
          TGo4Log::Error("TRpcParam %s: Setup file %s can not be opened or is missing.", GetName(), sOptionFileName.Data());
          TGo4Log::Error("=> Using default values!");
   }
   else
   {
      TGo4Log::Info("TRpcParam %s: Load Setup file %s.", GetName(), sOptionFileName.Data());
      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("tdcoffscaen")==kTRUE)
         {
            iOffsetTdcCaen = ((TString)line(12,10)).Atoi();
            if( iOffsetTdcCaen > MAX_1290)
            {
               cout<<"Error: CAEN tdc for offset higher than max allowed -> replaced by default value: ";
               cout<<TDC_FOR_OFFSET_CAEN<<endl;
               iOffsetTdcCaen = TDC_FOR_OFFSET_CAEN;
            }
         }
         else if(line.Contains("chaoffscaen")==kTRUE)
         {
            iOffsetChannelCaen = ((TString)line(12,10)).Atoi();
            if( iOffsetChannelCaen > T1290Data::NumChan)
            {
               cout<<"Error: CAEN channel for offset higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_OFFSET_CAEN<<endl;
               iOffsetChannelCaen = CH_FOR_OFFSET_CAEN;
            }
         }
         else if(line.Contains("tdcoffsvftx")==kTRUE)
         {
            iOffsetTdcVFTX = ((TString)line(12,10)).Atoi();
            if( iOffsetTdcVFTX > MAX_FPGA_TDC)
            {
               cout<<"Error: VFTX tdc for offset higher than max allowed -> replaced by default value: ";
               cout<<TDC_FOR_OFFSET_VFTX<<endl;
               iOffsetTdcVFTX = TDC_FOR_OFFSET_VFTX;
            }
         }
         else if(line.Contains("chaoffsvftx")==kTRUE)
         {
            iOffsetChannelVFTX = ((TString)line(12,10)).Atoi();
            if( iOffsetChannelVFTX > FPGA_TDC_NBCHAN)
            {
               cout<<"Error: VFTX channel for offset higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_OFFSET_VFTX<<endl;
               iOffsetChannelVFTX = CH_FOR_OFFSET_VFTX;
            }
         }
         else if(line.Contains("tdcdiamcaen")==kTRUE)
         {
            iDiamondTdcCaen = ((TString)line(12,10)).Atoi();
            if( iDiamondTdcCaen > MAX_1290)
            {
               cout<<"Error: CAEN tdc for diamond higher than max allowed -> replaced by default value: ";
               cout<<TDC_FOR_DIAMOND_CAEN<<endl;
               iDiamondTdcCaen = TDC_FOR_DIAMOND_CAEN;
            }
         }
         else if(line.Contains("chadiamcaen")==kTRUE)
         {
            iDiamondChannelCaen = ((TString)line(12,10)).Atoi();
            if( iDiamondChannelCaen > T1290Data::NumChan)
            {
               cout<<"Error: CAEN channel for diamond higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_DIAMOND_CAEN<<endl;
               iDiamondChannelCaen = CH_FOR_DIAMOND_CAEN;
            }
         }
         else if(line.Contains("tdcdiamvftx")==kTRUE)
         {
            iDiamondTdcVFTX = ((TString)line(12,10)).Atoi();
            if( iDiamondTdcVFTX > MAX_FPGA_TDC)
            {
               cout<<"Error: VFTX tdc for diamond higher than max allowed -> replaced by default value: ";
               cout<<TDC_FOR_DIAMOND_VFTX<<endl;
               iDiamondTdcVFTX = TDC_FOR_DIAMOND_VFTX;
            }
         }
         else if(line.Contains("chadiamvftx")==kTRUE)
         {
            iDiamondChannelVFTX = ((TString)line(12,10)).Atoi();
            if( iDiamondChannelVFTX > FPGA_TDC_NBCHAN)
            {
               cout<<"Error: VFTX channel for diamond higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_DIAMOND_VFTX<<endl;
               iDiamondChannelVFTX = CH_FOR_DIAMOND_VFTX;
            }
         }
         else if(line.Contains("rateintdur")==kTRUE)
         {
            iRateIntegrationDuration = ((TString)line(11,10)).Atoi();
         }
         else if(line.Contains("ratestep")==kTRUE)
         {
            iRateIntegrationStep = ((TString)line(11,10)).Atoi();
            if( iRateIntegrationStep < 2 )
            {
               cout<<"Rate step size smaller than Go4 minimal Hist. update step of 2s, set step size to 2s"<<endl;
               iRateIntegrationStep = 2;
            }
            if( iRateIntegrationDuration < iRateIntegrationStep )
            {
               cout<<"Rate step size smaller than rate integration size, set step size to integration size"<<endl;
               iRateIntegrationStep = iRateIntegrationDuration;
            }
         }
         else if(line.Contains("tdctypestrip")==kTRUE)
         {
            uTdcTypeStrip = ((TString)line(13,10)).Atoi();
         }
         else if(line.Contains("striprpcname")==kTRUE)
         {
            sStripRpcName = ((TString)line(13,20));
         }
         else if(line.Contains("striprpcposx")==kTRUE)
         {
            dStripRpcPosX = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("striprpcposy")==kTRUE)
         {
            dStripRpcPosY = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("striprpcposz")==kTRUE)
         {
            dStripRpcPosZ = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("nbstrips")==kTRUE)
         {
            uNbStrips = ((TString)line(9,10)).Atoi();
            if( uNbStrips > NB_STRIPS_MAX)
            {
               cout<<"Error: nb strips higher than max allowed for RPC "<<uNbStrips;
               cout<<" -> replaced by max value: ";
               cout<<NB_STRIPS_MAX<<endl;
               uNbStrips = NB_STRIPS_MAX;
            }
         }
         else if(line.Contains("stripwidth")==kTRUE)
         {
            dStripWidth = ((TString)line(12,10)).Atof();
         }
         else if(line.Contains("striplengt")==kTRUE)
         {
            dStripLength = ((TString)line(12,10)).Atof();
         }
         else if(line.Contains("timeconvcm")==kTRUE)
         {
            dTimeConv = ((TString)line(12,10)).Atof();
         }
         else if(line.Contains("hitsmaxdist")==kTRUE)
         {
            dStripHitMaxDistance = ((TString)line(12,10)).Atof();
         }
         else if(line.Contains("rejectmulti")==kTRUE)
         {
            uRejectMultipleHits = ((TString)line(13,10)).Atoi();
         }
         else if(line.Contains("orientation")==kTRUE)
         {
            uOrientation = ((TString)line(13,10)).Atoi();
         }
         else if(line.Contains("beamsize")==kTRUE)
         {
            iBeamProfileSize = ((TString)line(9,10)).Atoi();
         }
         else if(line.Contains("triggsel")==kTRUE)
         {
            iTriggerSelection = ((TString)line(9,10)).Atoi();
         }
         else if(line.Contains("strtdctype")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips tdc type array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripType[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("noisystrip")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and noisy strips array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripNoisy[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("deadstrip")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and dead strips array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripDead[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strlefttottdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left strips Tot TDC array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripLeftTot[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strlefttotcha")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left strips Tot Channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uChannelStripLeftTot[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strlefttimetdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left strips Time TDC array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripLeftTime[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strlefttimecha")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left strips Time channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uChannelStripLeftTime[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strrighttimetdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right strips Time TDC array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripRightTime[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strrighttimecha")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right strips Time channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uChannelStripRightTime[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strrighttottdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right strips Tot TDC array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uTdcStripRightTot[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strrighttotcha")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
            if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right strips Tot channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  uChannelStripRightTot[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_INDEXES)).Atoi();
                  line = line(NB_DIGITS_STRIP_INDEXES,NB_DIGITS_STRIP_INDEXES*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("stripsoffs")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips profile offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dProfOffsetList[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("reftimeoff")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips reference offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dRefOffsetList[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("ref1signaloffs")==kTRUE)
         {
            dReference1SignalOffset = ((TString)line(16,10)).Atof();
         }
         else if(line.Contains("ref2signaloffs")==kTRUE)
         {
            dReference2SignalOffset = ((TString)line(16,10)).Atof();
         }
         else if(line.Contains("stripsoffsleft")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips left offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dOffsetListLeft[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("stripsoffsright")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips right offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dOffsetListRight[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("totgravcentr")==kTRUE)
         {
            uUseChargeGravityCenter = ((TString)line(13,10)).Atoi();
         }
         else if(line.Contains("mindistbthit")==kTRUE)
         {
            dMinimalTimeBetweenHits = ((TString)line(13,10)).Atoi();
         }
         else if(line.Contains("strtotgainleft")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips left tot gain array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dToTGainListLeft[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strtotoffsleft")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips left tot offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dTotOffsetListLeft[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strtotgainright")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips right tot gain array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dToTGainListRight[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strtotoffsright")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips right tot offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dTotOffsetListRight[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("initoffscal")==kTRUE)
         {
            uInitialOffsetCalib = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("lowlimfitbea")==kTRUE)
         {
            dLowerBoundFitBeamProfile = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("upplimfitbea")==kTRUE)
         {
            dUpperBoundFitBeamProfile = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("lowlimfitref")==kTRUE)
         {
            dLowerBoundFitRefProfile = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("upplimfitref")==kTRUE)
         {
            dUpperBoundFitRefProfile = ((TString)line(13,10)).Atof();
         }
         else if(line.Contains("initvalfitparbe")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_OFFSETS*4);
            for(Int_t iFitParam = 0; iFitParam < 4 ; iFitParam++)
            {
               if( 1 == uDebug )
                  cout<<"initvalfitparbe "<<iFitParam<<" "<<line<<endl;
               dInitialValueFitParBeamProfile[iFitParam] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
               line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*4);
            }
         }
         else if(line.Contains("initvalfitparre")==kTRUE)
         {
            line = line(16,NB_DIGITS_STRIP_OFFSETS*4);
            for(Int_t iFitParam = 0; iFitParam < 4 ; iFitParam++)
            {
               dInitialValueFitParRefProfile[iFitParam] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
               line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*4);
            }
         }
         else if(line.Contains("strmatchingoffs")==kTRUE)
         {
           line = line(16,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
           if( uNbStrips == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips matching offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            }
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
               {
                  dOffsetMatching[temp_index] = ((TString)line(0,NB_DIGITS_STRIP_OFFSETS)).Atof();
                  line = line(NB_DIGITS_STRIP_OFFSETS,NB_DIGITS_STRIP_OFFSETS*NB_STRIPS_MAX);
               }
            }
         }
      } // while( !setupFile.eof() )
   } //  else of if( setupFile.fail() == true)
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
      if( 1 == uTdcStripType[temp_index] )
      {
         uTdcStripLeftTot[temp_index]     = uTdcStripLeftTime[temp_index];
         if( 1 == uChannelStripLeftTime[temp_index]%2)
            uChannelStripLeftTot[temp_index] = uChannelStripLeftTime[temp_index] - 1;
            else uChannelStripLeftTot[temp_index] = uChannelStripLeftTime[temp_index] + 1;
         uTdcStripRightTot[temp_index]     = uTdcStripRightTime[temp_index];
         if( 1 == uChannelStripRightTime[temp_index]%2)
            uChannelStripRightTot[temp_index] = uChannelStripRightTime[temp_index] - 1;
            else uChannelStripRightTot[temp_index] = uChannelStripRightTime[temp_index] + 1;
      }
   
   return kTRUE;
}


Int_t TRpcParam::PrintOptions(){
   cout<<"********************* RPC calibration parameters **********************"<<endl;
   cout<<"------------- Options for strip RPC ";
   cout.width(20);
   cout<<sStripRpcName<<" --------------"<<endl;
   cout<<"-------------------------------------------------"<<endl;
   if(1 == uDebug)
      cout<<"  Debug ON!!!!!!"<<endl;
      else cout<<"  Debug OFF"<<endl;

   if( -1 < iOffsetTdcCaen && -1 < iOffsetChannelCaen &&
       -1 < iOffsetTdcVFTX && -1 < iOffsetChannelVFTX )
   {
      cout<<"  Chan. used for time offset: "<<iOffsetChannelCaen;
      cout<<" in TDC "<<iOffsetTdcCaen<<" in CAEN"<<endl;
      cout<<"  Chan. used for time offset: "<<iOffsetChannelVFTX;
      cout<<" in TDC "<<iOffsetTdcVFTX<<" in VFTX"<<endl;
   }
   if( -1 < iRateIntegrationDuration && -1 < iRateIntegrationStep)
   {
      cout<<"  Rate hist. reset spacing:   "<<iRateIntegrationDuration<<" s"<<endl;
      cout<<"  Rate hist. update spacing:  "<<iRateIntegrationStep<<" s"<<endl;
   }
   if( 0 == uTdcTypeStrip )
   {
      cout<<"  TDC type:                  CAEN"<<endl;
   }
   else if( 1 == uTdcTypeStrip )
   {
      cout<<"  TDC type:                  VFTX"<<endl;
   }
   else if( 2 == uTdcTypeStrip )
   {
      cout<<"  TDC type:                  MIXED"<<endl;
   }
   cout<<"  Position:                  X= "<<dStripRpcPosX;
   cout<<" Y= "<<dStripRpcPosY<<" Z= "<<dStripRpcPosZ;
   cout<<endl;
   cout<<"  Nb strips in beam profile: "<<uNbStrips<<endl;
   cout<<"  Strip width (incl gap.):   "<<dStripWidth<<" cm"<<endl;
   cout<<"  Strip length:              "<<dStripLength<<" cm"<<endl;
   cout<<"  Time Convertion factor:    "<<dTimeConv<<" cm/ns"<<endl;
   cout<<"  Max Hit dist for cluster:  "<<dStripHitMaxDistance<<" cm"<<endl;
   if( 1 == uRejectMultipleHits )
      cout<<"  Multiple TDC hits reject.: ON"<<endl;
      else cout<<"  Multiple TDC hits reject.: OFF"<<endl;
   if( 0 == uOrientation )
   {
      cout<<"  Strips orientation:        Horizontal"<<endl;
   }
   else if( 1 == uOrientation )
   {
      cout<<"  Strips orientation:        Vertical"<<endl;
   }
   cout<<"  Size of beam profile in X: "<<iBeamProfileSize<<" cm "<<endl;
   if( -1 == iTriggerSelection )
      cout<<"  Trigger selection:         OFF"<<endl;
      else cout<<"  Trigger selection:         Trigger "<<iTriggerSelection<<endl;
   cout<<"-------------------------------------------------"<<endl;
   cout<<"  Strip channel:     |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<temp_index<<" ";
   }
   if( 2 == uTdcTypeStrip )
   {
      cout<<endl<<"  Tdc Type           |->     ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
      {
         cout.width(NB_DIGITS_STRIP_INDEXES - 1);
         cout<<uTdcStripType[temp_index]<<" ";
      }
   }
   cout<<endl<<"  Noisy Strip?       |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripNoisy[temp_index]<<" ";
   }
   cout<<endl<<"  Dead Strip?        |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripDead[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Tot Left       |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripLeftTot[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Tot Left   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uChannelStripLeftTot[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Time Left      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripLeftTime[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Time Left  |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uChannelStripLeftTime[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Time Right     |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripRightTime[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Time Right |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uChannelStripRightTime[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Tot Right      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uTdcStripRightTot[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Tot Right  |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_INDEXES - 1);
      cout<<uChannelStripRightTot[temp_index]<<" ";
   }
   cout<<endl;
   cout<<"-------------------------------------------------"<<endl;
   cout<<"  Strip channel:     |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Beam Profile Offset|->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dRefOffsetList[temp_index]<<" ";
   }
   cout<<endl<<"  Reference Offset   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dRefOffsetList[temp_index]<<" ";
   }
   cout<<endl<<"  Time Offset Left   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dOffsetListLeft[temp_index]<<" ";
   }
   cout<<endl<<"  Time Offset Right  |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dOffsetListRight[temp_index]<<" ";
   }
   cout<<endl;
   if( 0 != dReference1SignalOffset )
      cout<<"  Reference 1 signal offset: "<<dReference1SignalOffset<<" ps"<<endl;
   if( 0 != dReference2SignalOffset )
      cout<<"  Reference 2 signal offset: "<<dReference2SignalOffset<<" ps"<<endl;
   cout<<"-------------------------------------------------"<<endl;
   if( 1 == uUseChargeGravityCenter )
      cout<<"  Cluster center of gravity: ON"<<endl;
      else cout<<"  Cluster center of gravity: OFF"<<endl;
   if( 0 < dMinimalTimeBetweenHits )
      cout<<"  Minimal Dist btwn TDC hit: "<<dMinimalTimeBetweenHits<<" ps (Merged if closer)"<<endl;
      else cout<<"  Minimal Dist btwn TDC hit: OFF"<<endl;
   cout<<"  Strip channel:     |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Tot Gain Left      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dToTGainListLeft[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Offset Left    |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dTotOffsetListLeft[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Gain Right     |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dToTGainListRight[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Offset Right   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dTotOffsetListRight[temp_index]<<" ";
   }
   cout<<endl;
   cout<<"-------------------------------------------------"<<endl;
   if( 0 <= uInitialOffsetCalib )
      cout<<"  Initial offset calibration: "<<uInitialOffsetCalib<<endl;
      else cout<<"  Initial offset calibration: OFF"<<endl;
   cout<<"  Lower lim. Fit funct. beam:"<<dLowerBoundFitBeamProfile<<endl;
   cout<<"  Upper lim. Fit funct. beam:"<<dUpperBoundFitBeamProfile<<endl;
   cout<<"  Lower lim. Fit funct. Ref: "<<dLowerBoundFitRefProfile<<endl;
   cout<<"  Upper lim. Fit funct. Ref: "<<dUpperBoundFitRefProfile<<endl;

   cout<<"  Fit Parameter:     |-      ";
   for( Int_t iInitialParam = 0; iInitialParam < 4; iInitialParam++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<iInitialParam<<" ";
   }
   cout<<endl<<"  Initial value Beam |->     ";
   for( Int_t iInitialParam = 0; iInitialParam < 4; iInitialParam++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dInitialValueFitParBeamProfile[iInitialParam]<<" ";
   }
   cout<<endl<<"  Initial value Ref  |->     ";
   for( Int_t iInitialParam = 0; iInitialParam < 4; iInitialParam++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dInitialValueFitParRefProfile[iInitialParam]<<" ";
   }
   cout<<endl;
   cout<<"-------------------------------------------------"<<endl;
   cout<<"  Strip channel:     |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Matching Offset    |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbStrips; temp_index++)
   {
      cout.width(NB_DIGITS_STRIP_OFFSETS - 1);
      cout<<dOffsetMatching[temp_index]<<" ";
   }
   cout<<endl;
   cout<<"************************************************************************"<<endl;
   
   return 0;
}

/*
 * Compute the time offsets in ps for the left and right end of the strips using
 * the full offset used to center the beam profile and reference profile histograms
 */
Bool_t TRpcParam::CalculateLeftRightOffset()
{
   for( UInt_t uStrip = 0; uStrip < uNbStrips; uStrip++)
   {
      dOffsetListLeft[uStrip]  = 0.5*( (dProfOffsetList[uStrip])+2.*(dRefOffsetList[uStrip]));
      dOffsetListRight[uStrip] = 0.5*(-(dProfOffsetList[uStrip])+2.*(dRefOffsetList[uStrip]));
   }
   return kTRUE;
}