beamtime/gsi-aug12/hd/go4/MONITOR/TBeamMonitorParam.cxx (r4864/r3218)

Go to the documentation of this file.

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

#include "TMbsCrateEvent.h"

/***** Default values definition *****/
#define DEBUG                 0
#define CH_FOR_OFFSET_CAEN    -1
#define CH_FOR_OFFSET_VFTX    -1
#define RATE_INT_DURATION     -1
#define RATE_INT_STEP         -1

#define NB_SRPC               0
#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 TIME_CONV            10.0
/*
         //                      0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15   16   17   18   19   20   21   22   23
         //                      |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |    |
#define STRIPS_LEFT_LIST      {  23,  22,  21,  20,  19,  18,  17,  16,  15,  14,  13,  12,  11,  10,   9,   8,   7,   6,   5,   4,   3,   2,   1,   0 }
#define STRIPS_RIGHT_LIST     {  31,  30,  29,  28,  27,  26,  25,  24,  32,  33,  34,  35,  36,  37,  38,  39,  40,  41,  42,  43,  44,  45,  46,  47 }
#define STRIPS_OFFSET_LIST    { 0.0, 0.0, 0.0, 0.0, 3.0, 3.0, 3.0, 3.0,-3.6,-5.3,-2.3,-4.0,-6.0,-7.0,-5.1,-3.2, 0.0, 9.9, 5.2, 4.5, 3.5, 4.8, 7.2, 5.3 }
* */
#define STRIPS_LEFT           0
#define STRIPS_RIGHT          0
#define STRIPS_OFFSET         0
#define BEAM_PROFILE_SIZE     0

#define NB_PRPC               0
#define TDC_PAD               0
#define RPC_NAME_PAD          ""
#define RPC_POSX_PAD          0.0
#define RPC_POSY_PAD          0.0
#define RPC_POSZ_PAD          0.0
#define NB_PADS_X             0 // Nb of pads in X direction for the beam profile monitor
#define NB_PADS_Y             0 // Nb of pads in Y direction for the beam profile monitor
#define PADS_INDEX            0
/*
         //                      0    1    2    3    4    5
         //                      |    |    |    |    |    |
#define PADS_LIST            {{  57,  58,  59,  60,  61,  62 }, \
                              {  65,  66,  67,  68,  69,  71 } }
* */
/*************************************/

//***********************************************************                          
                              
TBeamMonitorParam::TBeamMonitorParam(const char* name) : TGo4Parameter(name)
{
//   Dummy = 0;
   
   uDebug             = 0;
   iOffsetFullChannelCaen = -1;
   iOffsetFullChannelVFTX = -1;
   iRateIntegrationDuration = -1;
   iRateIntegrationStep = -1;

   iNbStripRpcs = 0; // Nb of strip RPC for the beam profile monitor
   iNbPadRpcs = 0;   // Nb of pad RPC for the beam profile monitor
   
   ReadOptions();
   this->PrintOptions();
}
//***********************************************************   

Bool_t TBeamMonitorParam::UpdateFrom(TGo4Parameter *pp)
{
  if(pp->InheritsFrom("TBeamMonitorParam")) {
     TBeamMonitorParam * from = (TBeamMonitorParam *) pp;      
     cout << "**** TBeamMonitorParam " << GetName() << " updated from source parameter" << endl;
  } else
     cout << "Wrong parameter object: " << pp->ClassName() << endl;
  return kTRUE;
}
//***********************************************************  

TBeamMonitorParam::~TBeamMonitorParam()
{
   /*
   cout<<Form("Par Strip        %p %p %p %p %p %p %p %p",
         uTdcTypeStrip, sStripRpcName, dStripRpcPosX, dStripRpcPosY, dStripRpcPosZ, uNbStrips, dStripWidth, dTimeConv
                           )<<endl;

   cout<<Form("Par Strip ch     %p %p %p %p %p %p",
         uStripLeftList, uStripRightList, uTdcStripLeftList, uTdcStripRightList, uChannelStripLeft, uChannelStripRight
                           )<<endl;
   cout<<Form("Par Strip ch tot %p %p %p %p %p %p",
         uTotStripLeftList, uTotStripRightList, uTotTdcStripLeftList, uTotTdcStripRightList, uTotChannelStripLeft, uTotChannelStripRight
                           )<<endl;

   cout<<Form("Par Strip beam   %p %p",
         dOffsetList, iBeamProfileSize
                           )<<endl;

   // Pad RPCs settings

   cout<<Form("Par Pad          %p %p %p %p %p %p %p %p %p %p %p %p",
         uTdcTypePad, sPadRpcName, dPadRpcPosX, dPadRpcPosY, dPadRpcPosZ, uNbPadsX,
         uNbPadsY, dPadSizeX, dPadSizeY, uPadsList, uPadsTdc, uPadsChan
                           )<<endl;
   */
   cout << "**** TBeamMonitorParam: Delete instance " << endl;
}
//-----------------------------------------------------------

Bool_t TBeamMonitorParam::ReadOptions()
{
   // Loading default values for safety
   uDebug             = DEBUG;
   iOffsetFullChannelCaen = CH_FOR_OFFSET_CAEN;
   iOffsetFullChannelVFTX = CH_FOR_OFFSET_VFTX;
   iNbStripRpcs = NB_SRPC;
   iNbPadRpcs   = NB_PRPC;
   iRateIntegrationDuration = RATE_INT_DURATION;
   iRateIntegrationStep     = RATE_INT_STEP;

   // Reading out parameters from the option file
   fstream setupFile("./Monitor_options.txt", ios::in);
   if( setupFile.fail() == true)
   {
      cout<<"Error opening setup file rocget4_options.txt"<<endl;
      cout<<"=> Using default values!"<<endl;
      
      for(UInt_t temp_index = 0; temp_index < NB_STRIP_RPC_MAX; temp_index++)
      {
         iBeamProfileSize[temp_index]     = BEAM_PROFILE_SIZE;
      }

   }
   else
   {
      TString line;
      Int_t iStripRpcIndex = -1;
      Int_t iPadRpcIndex = -1;
      line.ReadLine(setupFile, kFALSE);
      while( !setupFile.eof() )
      {
         if(line(0) == '#' || line(0) == '*')
         {
            //Comment
            ;
         }
         else if(line.Contains("debug")==kTRUE)
         {
            uDebug = ((TString)line(6,10)).Atoi();
         }
         else if(line.Contains("choffscaen")==kTRUE)
         {
            iOffsetFullChannelCaen = ((TString)line(11,10)).Atoi();
            if( iOffsetFullChannelCaen > MAX_1290*T1290Data::NumChan)
            {
               cout<<"Error: CAEN channel for offset higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_OFFSET_CAEN<<endl;
               iOffsetFullChannelCaen = CH_FOR_OFFSET_CAEN;
            }
         }
         else if(line.Contains("choffsvftx")==kTRUE)
         {
            iOffsetFullChannelVFTX = ((TString)line(11,10)).Atoi();
            if( iOffsetFullChannelCaen > MAX_FPGA_TDC*FPGA_TDC_NBCHAN)
            {
               cout<<"Error: VFTX channel for offset higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_OFFSET_VFTX<<endl;
               iOffsetFullChannelVFTX = CH_FOR_OFFSET_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("nbstriprpc")==kTRUE)
         {
            iNbStripRpcs = ((TString)line(11,10)).Atoi();
            if( iNbStripRpcs > NB_STRIP_RPC_MAX)
            {
                cout<<"Error: nb strip RPC higher than max allowed -> replaced by max value: ";
                cout<<NB_STRIP_RPC_MAX<<endl;
                iNbStripRpcs = NB_STRIP_RPC_MAX;
            }
         }
         else if(line.Contains("tdctypestrip")==kTRUE)
         {
            iStripRpcIndex ++;
            if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               uTdcTypeStrip[iStripRpcIndex] = ((TString)line(13,10)).Atoi();
            }
         }
         else if(line.Contains("striprpcname")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               sStripRpcName[iStripRpcIndex] = ((TString)line(13,20));
            }
         }
         else if(line.Contains("striprpcposx")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dStripRpcPosX[iStripRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("striprpcposy")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dStripRpcPosY[iStripRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("striprpcposz")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dStripRpcPosZ[iStripRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("nbstrips")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               uNbStrips[iStripRpcIndex] = ((TString)line(9,10)).Atoi();
               if( uNbStrips[iStripRpcIndex] > NB_STRIPS_MAX)
               {
                  cout<<"Error: nb strips higher than max allowed for RPC "<<uNbStrips[iStripRpcIndex];
                  cout<<" -> replaced by max value: ";
                  cout<<NB_STRIPS_MAX<<endl;
                  uNbStrips[iStripRpcIndex] = NB_STRIPS_MAX;
               }
            }
         }
         else if(line.Contains("stripwidth")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dStripWidth[iStripRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("timeconvcm")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dTimeConv[iStripRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("stripsleft")==kTRUE)
         {
            line = line(12,5*NB_STRIPS_MAX);
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else if( uNbStrips[iStripRpcIndex] == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left 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[iStripRpcIndex]; temp_index++)
               {
                  uStripLeftList[iStripRpcIndex][temp_index] = ((TString)line(0,5)).Atoi();
                  line = line(5,5*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("stripsright")==kTRUE)
         {
            line = line(12,5*NB_STRIPS_MAX);
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channel index!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else if( uNbStrips[iStripRpcIndex] == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right 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[iStripRpcIndex]; temp_index++)
               {
                  uStripRightList[iStripRpcIndex][temp_index] = ((TString)line(0,5)).Atoi();
                  line = line(5,5*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("stripsoffs")==kTRUE)
         {
#define OFFSET_DIGITS 7
            line = line(12,OFFSET_DIGITS*NB_STRIPS_MAX);
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channel index!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else if( uNbStrips[iStripRpcIndex] == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and strips 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[iStripRpcIndex]; temp_index++)
               {
                  dOffsetList[iStripRpcIndex][temp_index] = ((TString)line(0,OFFSET_DIGITS)).Atof();
                  line = line(OFFSET_DIGITS,OFFSET_DIGITS*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strtotleft")==kTRUE)
         {
            line = line(12,5*NB_STRIPS_MAX);
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channels!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else if( uNbStrips[iStripRpcIndex] == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and left 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[iStripRpcIndex]; temp_index++)
               {
                  uTotStripLeftList[iStripRpcIndex][temp_index] = ((TString)line(0,5)).Atoi();
                  line = line(5,5*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("strtotright")==kTRUE)
         {
            line = line(12,5*NB_STRIPS_MAX);
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channel index!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else if( uNbStrips[iStripRpcIndex] == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbStrips is 0 and right 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[iStripRpcIndex]; temp_index++)
               {
                  uTotStripRightList[iStripRpcIndex][temp_index] = ((TString)line(0,5)).Atoi();
                  line = line(5,5*NB_STRIPS_MAX);
               }
            }
         }
         else if(line.Contains("beamsize")==kTRUE)
         {
            if( -1 == iStripRpcIndex )
            {
               cout<<" Index for the strip RPCs is -1 => please declare number of strips in this RPC before channel index!"<<endl;
            }
            else if( iNbStripRpcs <= iStripRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the strip RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
               iBeamProfileSize[iStripRpcIndex] = ((TString)line(9,10)).Atoi();
         }
         
         else if(line.Contains("nbpadsrpc")==kTRUE)
         {
            iNbPadRpcs = ((TString)line(11,10)).Atoi();
            if( iNbPadRpcs > NB_PAD_RPC_MAX)
            {
                cout<<"Error: nb strip RPC higher than max allowed -> replaced by max value: ";
                cout<<NB_PAD_RPC_MAX<<endl;
                iNbPadRpcs = NB_PAD_RPC_MAX;
            }
         }
         else if(line.Contains("tdctypepad")==kTRUE)
         {
            iPadRpcIndex ++;
            if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               uTdcTypePad[iPadRpcIndex] = ((TString)line(12,10)).Atoi();
            }
         }
         else if(line.Contains("padrpcname")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declare tdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               sPadRpcName[iPadRpcIndex] = ((TString)line(11,20));
            }
         }
         else if(line.Contains("padrpcposx")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dPadRpcPosX[iPadRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("padrpcposy")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dPadRpcPosY[iPadRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("padrpcposz")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dPadRpcPosZ[iPadRpcIndex] = ((TString)line(12,10)).Atof();
            }
         }
         else if(line.Contains("nbpadsx")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               uNbPadsX[iPadRpcIndex] = ((TString)line(9,10)).Atoi();
               if( uNbPadsX[iPadRpcIndex] > NB_PADS_X_MAX)
               {
                  cout<<"Error: nb pads on X axis higher than max allowed -> replaced by max value ";
                  cout<<NB_PADS_X_MAX<<endl;
                  uNbPadsX[iPadRpcIndex] = NB_PADS_X_MAX;
               }
            }
         }
         else if(line.Contains("nbpadsy")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of pad RPC declared"<<endl;
            }
            else 
            {
               uNbPadsY[iPadRpcIndex] = ((TString)line(9,10)).Atoi();
               if( uNbPadsY[iPadRpcIndex] > NB_PADS_Y_MAX)
               {
                  cout<<"Error: nb pads on Y axis higher than max allowed -> replaced by max value ";
                  cout<<NB_PADS_Y_MAX<<endl;
                  uNbPadsY[iPadRpcIndex] = NB_PADS_Y_MAX;
               }
            }
         }
         else if(line.Contains("sizepadx")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dPadSizeX[iPadRpcIndex] = ((TString)line(9,10)).Atof();
            }
         }
         else if(line.Contains("sizepady")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of strip RPC declared"<<endl;
            }
            else
            {
               dPadSizeY[iPadRpcIndex] = ((TString)line(9,10)).Atof();
            }
         }
         else if(line.Contains("padsrow")==kTRUE)
         {
            if( -1 == iPadRpcIndex )
            {
               cout<<" Index for the pad RPCs is -1 => please declaretdc type for this RPC before other caracs!"<<endl;
            }
            else if( iNbPadRpcs <= iPadRpcIndex )
            {
               if( 1 == uDebug )
                  cout<<" Index for the pad RPCs out of bound => Please check the number of pad RPC declared"<<endl;
            }
            else 
            {
               for(UInt_t temp_indexY = 0; temp_indexY < uNbPadsY[iPadRpcIndex]; temp_indexY++)
               {
                  if( line.Contains( Form("padsrow%u", temp_indexY+1)) ==kTRUE)
                  {
                     line = line(12,280);
                     for(UInt_t temp_indexX = 0; temp_indexX < uNbPadsX[iPadRpcIndex]; temp_indexX++)
                     {
                        uPadsList[iPadRpcIndex][temp_indexY][temp_indexX] = ((TString)line(0,5)).Atoi();
                        line = line(5,280);
                     }
                  }
                  line.ReadLine(setupFile, kFALSE);
               }
            }
         }
         
         line.ReadLine(setupFile, kFALSE);
      }
   }
   
   setupFile.close();
   
   return kTRUE;
}
//-----------------------------------------------------------

Int_t TBeamMonitorParam::PrintOptions()
{   
   cout<<"*********************** Beam monitor Parameters ************************"<<endl;
   if( 1 == uDebug) 
      cout<<"Debug output:               ON"<<endl;
   if( -1 < iOffsetChannelCaen && -1 < iOffsetChannelVFTX )
   {
      cout<<"Chan. used for time offset: "<<iOffsetChannelCaen;
      cout<<" in CAEN and "<<iOffsetChannelVFTX<<" in VFTX"<<endl;
   }
   if( -1 < iRateIntegrationDuration && -1 < iRateIntegrationStep)
   {
      cout<<"Rate hist. reset spacing:   "<<iRateIntegrationDuration<<" s"<<endl;
      cout<<"Rate hist. update spacing:  "<<iRateIntegrationStep<<" s"<<endl;
   }
   cout<<"-------------------------------------------------"<<endl;
   cout<<"Nb of strip RPC monitored: "<<iNbStripRpcs<<endl;
   for( Int_t iStripRpcIndex = 0; iStripRpcIndex < iNbStripRpcs; iStripRpcIndex++)
   {
      cout<<"-> Options for strip RPC "<<sStripRpcName[iStripRpcIndex]<<endl;
      if( 1 == uTdcTypeStrip[iStripRpcIndex] )
         cout<<"   TDC type:                  VFTX"<<endl;
         else cout<<"   TDC type:                  CAEN"<<endl;
      cout<<"   Position:                  X= "<<dStripRpcPosX[iStripRpcIndex];
      cout<<" Y= "<<dStripRpcPosY[iStripRpcIndex]<<" Z= "<<dStripRpcPosZ[iStripRpcIndex];
      cout<<endl;
      cout<<"   Nb strips in beam profile: "<<uNbStrips[iStripRpcIndex]<<endl;
      cout<<"   Strip width (incl gap.):   "<<dStripWidth[iStripRpcIndex]<<" cm"<<endl;
      cout<<"   Time Convertion factor:    "<<dTimeConv[iStripRpcIndex]<<" cm/ns"<<endl;
      cout<<"   Strip channel:    |-      ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips[iStripRpcIndex]; temp_index++)
      {
         cout.width(3);
         cout<<temp_index<<" ";
      }
      cout<<endl<<"   Left              |->     ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips[iStripRpcIndex]; temp_index++)
      {
         cout.width(3);
         cout<<uStripLeftList[iStripRpcIndex][temp_index]<<" ";
      }
      cout<<endl<<"   Right             |->     ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips[iStripRpcIndex]; temp_index++)
      {
         cout.width(3);
         cout<<uStripRightList[iStripRpcIndex][temp_index]<<" ";
      }
      cout<<endl;
      cout<<"   Strip offset:     |-      ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips[iStripRpcIndex]; temp_index++)
      {
         cout.width(OFFSET_DIGITS -1);
         cout<<temp_index<<" ";
      }
      cout<<endl<<"                     |->     ";
      for(UInt_t temp_index = 0; temp_index < uNbStrips[iStripRpcIndex]; temp_index++)
      {
         cout.width(OFFSET_DIGITS -1);
         cout<<dOffsetList[iStripRpcIndex][temp_index]<<" ";
      }
      cout<<endl;
      cout<<"   Size of beam profile in X: "<<iBeamProfileSize[iStripRpcIndex]<<" cm "<<endl;
      if( iStripRpcIndex + 1 < iNbStripRpcs )
         cout<<"              ...................                "<<endl;
   } // for( Int_t iStripRpcIndex = 0; iNbStripRpcs < iStripRpcIndex; iStripRpcIndex++)
   cout<<"-------------------------------------------------"<<endl;
   cout<<"Nb of pad RPC monitored:   "<<iNbPadRpcs<<endl;
   for( Int_t iPadRpcIndex = 0; iPadRpcIndex < iNbPadRpcs; iPadRpcIndex++)
   {
      cout<<"-> Options for pad RPC "<<sPadRpcName[iPadRpcIndex]<<endl;
      if( 1 == uTdcTypePad[iPadRpcIndex] )
         cout<<"   TDC type:                  VFTX"<<endl;
         else cout<<"   TDC type:                  CAEN"<<endl;
      cout<<"   Position:                  X= "<<dPadRpcPosX[iPadRpcIndex];
      cout<<" Y= "<<dPadRpcPosY[iPadRpcIndex]<<" Z= "<<dPadRpcPosZ[iPadRpcIndex];
      cout<<endl;
      cout<<"   Nb pads in X axis in beam profile: "<<uNbPadsX[iPadRpcIndex]<<endl;
      cout<<"   Nb pads in Y axis in beam profile: "<<uNbPadsY[iPadRpcIndex]<<endl;
      cout<<"   Pad size in X dir. (incl gap.):    "<<dPadSizeX[iPadRpcIndex]<<endl;
      cout<<"   Pad size in Y dir. (incl gap.):    "<<dPadSizeY[iPadRpcIndex]<<endl;
      cout<<"   Pad channel:      |-      ";
      for(UInt_t temp_index = 0; temp_index < uNbPadsX[iPadRpcIndex]; temp_index++)
      {
         cout.width(2);
         cout<<temp_index<<" ";
      }
      for(UInt_t temp_indexY = 0; temp_indexY < uNbPadsY[iPadRpcIndex]; temp_indexY++)
      {
         cout<<endl<<"   Row ";
         cout.width(2);
         cout<<temp_indexY<<"            |->     ";
         for(UInt_t temp_indexX = 0; temp_indexX < uNbPadsX[iPadRpcIndex]; temp_indexX++)
         {
            cout.width(2);
            cout<<uPadsList[iPadRpcIndex][temp_indexY][temp_indexX]<<" ";
         }
      }
      cout<<endl;
      if( iPadRpcIndex + 1 < iNbPadRpcs )
         cout<<"              ...................                "<<endl;
   } // for( Int_t iPadRpcIndex = 0; iNbPadRpcs < iPadRpcIndex; iPadRpcIndex++)
   cout<<"************************************************************************"<<endl;
   
   return 0;
}

Int_t   TBeamMonitorParam::ConvertChToTdcAndCh( Int_t iChannelNumberCaenTdc, Int_t iChannelNumberVftxTdc )
{
   // Time offset between CAEN and VFTX
   iOffsetTdcCaen     = iOffsetFullChannelCaen / iChannelNumberCaenTdc;
   iOffsetChannelCaen = iOffsetFullChannelCaen % iChannelNumberCaenTdc;
   iOffsetTdcVFTX     = iOffsetFullChannelVFTX / iChannelNumberVftxTdc;
   iOffsetChannelVFTX = iOffsetFullChannelVFTX % iChannelNumberVftxTdc;

   // Strip RPCs
   for( Int_t iStripRpcIndex = 0; iStripRpcIndex < iNbStripRpcs; iStripRpcIndex++)
   {
      if( 1 == uTdcTypeStrip[iStripRpcIndex] )
      {
         for( UInt_t uStripIndex = 0; uStripIndex < uNbStrips[iStripRpcIndex]; uStripIndex++)
         {
            uTdcStripLeftList[iStripRpcIndex][uStripIndex]  = uStripLeftList[iStripRpcIndex][uStripIndex]  / iChannelNumberVftxTdc;
            uChannelStripLeft[iStripRpcIndex][uStripIndex]  = uStripLeftList[iStripRpcIndex][uStripIndex]  % iChannelNumberVftxTdc;
            uTdcStripRightList[iStripRpcIndex][uStripIndex] = uStripRightList[iStripRpcIndex][uStripIndex] / iChannelNumberVftxTdc;
            uChannelStripRight[iStripRpcIndex][uStripIndex] = uStripRightList[iStripRpcIndex][uStripIndex] % iChannelNumberVftxTdc;
         } // for( UInt_t uStripIndex = 0; uStripIndex < fPar->uNbStrips[iStripRpcIndex]; uStripIndex++)
      } // if( 1 == uTdcTypeStrip[iStripRpcIndex] )
         else
         {
            for( UInt_t uStripIndex = 0; uStripIndex < uNbStrips[iStripRpcIndex]; uStripIndex++)
            {
               uTdcStripLeftList[iStripRpcIndex][uStripIndex]  = uStripLeftList[iStripRpcIndex][uStripIndex]  / iChannelNumberCaenTdc;
               uChannelStripLeft[iStripRpcIndex][uStripIndex]  = uStripLeftList[iStripRpcIndex][uStripIndex]  % iChannelNumberCaenTdc;
               uTdcStripRightList[iStripRpcIndex][uStripIndex] = uStripRightList[iStripRpcIndex][uStripIndex] / iChannelNumberCaenTdc;
               uChannelStripRight[iStripRpcIndex][uStripIndex] = uStripRightList[iStripRpcIndex][uStripIndex] % iChannelNumberCaenTdc;

               uTotTdcStripLeftList[iStripRpcIndex][uStripIndex]  = uTotStripLeftList[iStripRpcIndex][uStripIndex]  / iChannelNumberCaenTdc;
               uTotChannelStripLeft[iStripRpcIndex][uStripIndex]  = uTotStripLeftList[iStripRpcIndex][uStripIndex]  % iChannelNumberCaenTdc;
               uTotTdcStripRightList[iStripRpcIndex][uStripIndex] = uTotStripRightList[iStripRpcIndex][uStripIndex] / iChannelNumberCaenTdc;
               uTotChannelStripRight[iStripRpcIndex][uStripIndex] = uTotStripRightList[iStripRpcIndex][uStripIndex] % iChannelNumberCaenTdc;
            } // for( UInt_t uStripIndex = 0; uStripIndex < fPar->uNbStrips[iStripRpcIndex]; uStripIndex++)
         } // else of if( 1 == uTdcTypeStrip[iStripRpcIndex] )
   } // for( Int_t iStripRpcIndex = 0; iStripRpcIndex < iNbStripRpcs; iStripRpcIndex++)
   
   // Pad RPCs
   for( Int_t iPadRpcIndex = 0; iPadRpcIndex < iNbPadRpcs; iPadRpcIndex++)
   {
      if( 1 == uTdcTypePad[iPadRpcIndex] )
      {
         for(UInt_t uPadIndexY = 0; uPadIndexY < uNbPadsY[iPadRpcIndex]; uPadIndexY++)
            for(UInt_t uPadIndexX = 0; uPadIndexX < uNbPadsX[iPadRpcIndex]; uPadIndexX++)
         {      
            uPadsTdc[iPadRpcIndex][uPadIndexY][uPadIndexX]  = uPadsList[iPadRpcIndex][uPadIndexY][uPadIndexX] / iChannelNumberVftxTdc;
            uPadsChan[iPadRpcIndex][uPadIndexY][uPadIndexX] = uPadsList[iPadRpcIndex][uPadIndexY][uPadIndexX] % iChannelNumberVftxTdc;
         } // for uPadIndexY uPadIndexX
      } // if( 1 == uTdcTypePad[iPadRpcIndex] )
         else
         {
            for(UInt_t uPadIndexY = 0; uPadIndexY < uNbPadsY[iPadRpcIndex]; uPadIndexY++)
               for(UInt_t uPadIndexX = 0; uPadIndexX < uNbPadsX[iPadRpcIndex]; uPadIndexX++)
            { 
               uPadsTdc[iPadRpcIndex][uPadIndexY][uPadIndexX]  = uPadsList[iPadRpcIndex][uPadIndexY][uPadIndexX] / iChannelNumberCaenTdc;
               uPadsChan[iPadRpcIndex][uPadIndexY][uPadIndexX] = uPadsList[iPadRpcIndex][uPadIndexY][uPadIndexX] % iChannelNumberCaenTdc;
            } // for uPadIndexY uPadIndexX
         }// else of if( 1 == uTdcTypePad[iPadRpcIndex] )
   } //for( Int_t iPadRpcIndex = 0; iPadRpcIndex < iNbPadRpcs; iPadRpcIndex++)
   
   return 0;
}