beamtime/tof-tdctest/go4/PLASTICS/TPlasticsParam.cxx (r4864/r4267)

Go to the documentation of this file.

/* Generated by Together */

#include "TPlasticsParam.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_PLASTICS_INDEXES  5
#define NB_DIGITS_PLASTICS_OFFSETS  9

/***** Default values definition *****/
#define DEBUG                   0
#define NB_PLASTICS             0 // Nb of strips for the beam profile monitor
#define TIME_CONV              10.0

#define PLASTICS_POSX           0
#define PLASTICS_POSY           0
#define PLASTICS_POSZ           0
#define PLASTIC_ORIENTATION     0
#define PLASTICS_WIDTH          0
#define PLASTICS_LENGTH         0
#define PLASTIC_TDC_TYPE        1
#define PLASTIC_NB_SIDES        0
#define PLASTICS_LEFT_TOT_TDC   0
#define PLASTICS_LEFT_TOT_CHA   0
#define PLASTICS_LEFT_TIME_TDC  0
#define PLASTICS_LEFT_TIME_CHA  0
#define PLASTICS_RIGHT_TIME_TDC 0
#define PLASTICS_RIGHT_TIME_CHA 0
#define PLASTICS_RIGHT_TOT_TDC  0
#define PLASTICS_RIGHT_TOT_CHA  0

#define TDC_FOR_REFERENCE     -1
#define CH_FOR_REFERENCE      -1
#define TDC_FOR_DIAMOND_VFTX  -1
#define CH_FOR_DIAMOND_VFTX   -1

#define TRIGGER_SELECTION    -1

#define PLASTICS_PROFILE_OFFS   0
#define PLASTICS_REFERENCE_OFFS 0
#define PLASTICS_OFFSET_LEFT    0
#define PLASTICS_OFFSET_RIGHT   0

#define MIN_TIME_BTWN_HITS    0
#define PLASTICS_TOT_GAIN_LEFT  1
#define PLASTICS_TOT_OFFS_LEFT  0
#define PLASTICS_TOT_GAIN_RIGHT 1
#define PLASTICS_TOT_OFFS_RIGHT 0
/*************************************/

TPlasticsParam::TPlasticsParam(const char* name) : TGo4Parameter(name)
{
   ReadOptions();
   if( 0 < uNbPlastics )
   {
      // At least one Plastic in class !
      CalculateLeftRightOffset();
      PrintOptions(); //print out all values to see, if the setting has worked correctly
   }
}

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

Bool_t TPlasticsParam::ReadOptions(){
   // Loading default values
   uDebug              = DEBUG;
   iReferenceTdc       = TDC_FOR_REFERENCE;
   iReferenceChannel   = CH_FOR_REFERENCE;
   iDiamondTdcVFTX     = TDC_FOR_DIAMOND_VFTX;
   iDiamondChannelVFTX = CH_FOR_DIAMOND_VFTX;
   uNbPlastics         = NB_PLASTICS;
   dTimeConv           = TIME_CONV;
   iTriggerSelection   = TRIGGER_SELECTION;
   dMinimalTimeBetweenHits  = MIN_TIME_BTWN_HITS;
   dReference1SignalOffset   = PLASTICS_REFERENCE_OFFS;
   dReference2SignalOffset   = PLASTICS_REFERENCE_OFFS;

   for( Int_t iPlastic = 0; iPlastic < NB_PLASTICS_MAX; iPlastic ++)
   {
      dPlasticsPosX[iPlastic]  = PLASTICS_POSX;
      dPlasticsPosY[iPlastic]  = PLASTICS_POSY;
      dPlasticsPosZ[iPlastic]  = PLASTICS_POSZ;

      uOrientation[iPlastic]   = PLASTIC_ORIENTATION;
      dPlasticWidth[iPlastic]  = PLASTICS_WIDTH;
      dPlasticLength[iPlastic] = PLASTICS_LENGTH;

      uPlasticTdcType[iPlastic]          = PLASTIC_TDC_TYPE;
      uNbSides[iPlastic]                 = PLASTIC_NB_SIDES;
      uTdcPlasticLeftTot[iPlastic]       = PLASTICS_LEFT_TOT_TDC;
      uChannelPlasticLeftTot[iPlastic]   = PLASTICS_LEFT_TOT_CHA;
      uTdcPlasticLeftTime[iPlastic]      = PLASTICS_LEFT_TIME_TDC;
      uChannelPlasticLeftTime[iPlastic]  = PLASTICS_LEFT_TIME_CHA;
      uTdcPlasticRightTime[iPlastic]     = PLASTICS_RIGHT_TIME_TDC;
      uChannelPlasticRightTime[iPlastic] = PLASTICS_RIGHT_TIME_CHA;
      uTdcPlasticRightTot[iPlastic]      = PLASTICS_RIGHT_TOT_TDC;
      uChannelPlasticRightTot[iPlastic]  = PLASTICS_RIGHT_TOT_CHA;

      dProfOffsetList[iPlastic]          = PLASTICS_PROFILE_OFFS;
      dRefOffsetList[iPlastic]           = PLASTICS_REFERENCE_OFFS;

      dOffsetListLeft[iPlastic]          = PLASTICS_OFFSET_LEFT;
      dOffsetListRight[iPlastic]         = PLASTICS_OFFSET_RIGHT;

      dToTGainListLeft[iPlastic]         = PLASTICS_TOT_GAIN_LEFT;
      dTotOffsetListLeft[iPlastic]       = PLASTICS_TOT_OFFS_LEFT;
      dToTGainListRight[iPlastic]        = PLASTICS_TOT_GAIN_RIGHT;
      dTotOffsetListRight[iPlastic]      = PLASTICS_TOT_OFFS_RIGHT;
   }
   
   // Read the option file
   TString sOptionFileName = "./Plastics_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 TPlasticsParam "<<GetName()<<" not found: using the default value ";
      cout<<sOptionFileName<<endl;
   } // if( optionFilenameFile.fail() == true)
      else
      {
         TString line;

         /*
          */
           // New
         UInt_t uIndexPlastic =((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
              ;
           }
           else if( kTRUE == line.Contains( Form("plastics%02d ", uIndexPlastic) ) )
           {
              sOptionFileName = ((TString)line(11, line.Sizeof()-11 ));
              cout<<"************************************************************************"<<endl;
              cout<<"Using "<<sOptionFileName<<" as source for the options of TPlasticsParam "<<GetName()<<" step"<<endl;
              bFilenameFound = kTRUE;
              break;
           }
         } // while( !optionFilenameFile.eof() )
         if( optionFilenameFile.eof() && kFALSE == bFilenameFound )
         {
            cout<<"************************************************************************"<<endl;
            cout<<"Field for TPlasticsParam "<<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("TPlasticsParam %s: Setup file %s can not be opened or is missing.", GetName(), sOptionFileName.Data());
          TGo4Log::Error("=> Using default values!");
   }
   else
   {
      TGo4Log::Info("TPlasticsParam %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(12,10)).Atoi();
            if(1 == uDebug)
               cout<<"Debug ON!!!!!!"<<endl;
         } // else if(line.Contains("debug")==kTRUE)
         else if(line.Contains("tdcreferenc")==kTRUE)
         {
            iReferenceTdc = ((TString)line(12,10)).Atoi();
            if( iReferenceTdc > MAX_FPGA_TDC)
            {
               cout<<"Error: VFTX tdc for reference higher than max allowed -> replaced by default value: ";
               cout<<TDC_FOR_REFERENCE<<endl;
               iReferenceTdc = TDC_FOR_REFERENCE;
            } // if( iReferenceTdc > MAX_FPGA_TDC)
         } // else if(line.Contains("tdcreferenc")==kTRUE)
         else if(line.Contains("chareferenc")==kTRUE)
         {
            iReferenceChannel = ((TString)line(12,10)).Atoi();
            if( iReferenceChannel > FPGA_TDC_NBCHAN)
            {
               cout<<"Error: VFTX channel for reference higher than max allowed -> replaced by default value: ";
               cout<<CH_FOR_REFERENCE<<endl;
               iReferenceChannel = CH_FOR_REFERENCE;
            } // if( iReferenceChannel > FPGA_TDC_NBCHAN)
         } // else if(line.Contains("chareferenc")==kTRUE)
         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;
            } // if( iDiamondTdcVFTX > MAX_FPGA_TDC)
         } // else if(line.Contains("tdcdiamvftx")==kTRUE)
         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;
            } // if( iDiamondChannelVFTX > FPGA_TDC_NBCHAN)
         } // else if(line.Contains("chadiamvftx")==kTRUE)
         else if(line.Contains("nbplastics")==kTRUE)
         {
            uNbPlastics = ((TString)line(12,10)).Atoi();
            if( uNbPlastics > NB_PLASTICS_MAX)
            {
               cout<<"Error: nb plastics higher than max allowed "<<uNbPlastics;
               cout<<" -> replaced by max value: ";
               cout<<NB_PLASTICS_MAX<<endl;
               uNbPlastics = NB_PLASTICS_MAX;
            }
         } // else if(line.Contains("nbplastics")==kTRUE)
         else if(line.Contains("timeconvcm")==kTRUE)
         {
            dTimeConv = ((TString)line(12,10)).Atof();
         } // else if(line.Contains("timeconvcm")==kTRUE)
         else if(line.Contains("triggsel")==kTRUE)
         {
            iTriggerSelection = ((TString)line(12,10)).Atoi();
         } // else if(line.Contains("triggsel")==kTRUE)
         else if(line.Contains("plastictdctype")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics tdc type array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uPlasticTdcType[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platdctype")==kTRUE)
         else if(line.Contains("plasticposx")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics pos X array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dPlasticsPosX[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atof();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plastposx")==kTRUE)
         else if(line.Contains("plasticposy")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics pos Y array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dPlasticsPosY[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atof();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plastposy")==kTRUE)
         else if(line.Contains("plasticposz")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics pos Z array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dPlasticsPosZ[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atof();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plastposz")==kTRUE)
         else if(line.Contains("plasticorient")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics orientation array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uOrientation[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plaorient")==kTRUE)
         else if(line.Contains("plasticwidth")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics width array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dPlasticWidth[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atof();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plastwidt")==kTRUE)
         else if(line.Contains("plasticlength")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics length array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dPlasticLength[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atof();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plastleng")==kTRUE)
         else if(line.Contains("planbside")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics nb read sides array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uNbSides[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("planbside")==kTRUE)
         else if(line.Contains("plalefttottdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics left tot tdc array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uTdcPlasticLeftTot[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttottdc")==kTRUE)
         else if(line.Contains("plalefttotcha")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics left tot channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uChannelPlasticLeftTot[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttotcha")==kTRUE)
         else if(line.Contains("plalefttimetdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics left time tdc array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uTdcPlasticLeftTime[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttimetdc")==kTRUE)
         else if(line.Contains("plalefttimecha")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics left time channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uChannelPlasticLeftTime[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttimecha")==kTRUE)
         else if(line.Contains("plarighttimetdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics right time tdc array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uTdcPlasticRightTime[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plarighttimetdc")==kTRUE)
         else if(line.Contains("plarighttimecha")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics right time channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uChannelPlasticRightTime[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plarighttimecha")==kTRUE)
         else if(line.Contains("plarighttottdc")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics right tot tdc array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uTdcPlasticRightTot[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttottdc")==kTRUE)
         else if(line.Contains("plarighttotcha")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: uNbPlastics is 0 and plastics right tot channel array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  uChannelPlasticRightTot[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_INDEXES)).Atoi();
                  line = line(NB_DIGITS_PLASTICS_INDEXES,NB_DIGITS_PLASTICS_INDEXES*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("plalefttotcha")==kTRUE)
         else if(line.Contains("platimeoff")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics profile offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dProfOffsetList[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platimeoff")==kTRUE)
         else if(line.Contains("reftimeoff")==kTRUE)
         {
            line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
            if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics reference offsets array not empty: skip filling"<<line<<endl;
               cout<<"(if line were inverted in the option file, it could be the error)"<<endl;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dRefOffsetList[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("reftimeoff")==kTRUE)
         else if(line.Contains("ref1signaloffs")==kTRUE)
         {
            dReference1SignalOffset = ((TString)line(16,10)).Atof();
         } // else if(line.Contains("ref1signaloffs")==kTRUE)
         else if(line.Contains("ref2signaloffs")==kTRUE)
         {
            dReference2SignalOffset = ((TString)line(16,10)).Atof();
         } // else if(line.Contains("ref2signaloffs")==kTRUE)
         else if(line.Contains("mindistbthit")==kTRUE)
         {
            dMinimalTimeBetweenHits = ((TString)line(13,10)).Atoi();
         } // else if(line.Contains("mindistbthit")==kTRUE)
         else if(line.Contains("platotgainleft")==kTRUE)
         {
           line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
           if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics 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;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dToTGainListLeft[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platotgainleft")==kTRUE)
         else if(line.Contains("platotoffsleft")==kTRUE)
         {
           line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
           if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics 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;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dTotOffsetListLeft[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platotoffsleft")==kTRUE)
         else if(line.Contains("platotgainright")==kTRUE)
         {
           line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
           if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics 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;
               }
            } // if( uNbPlastics == 0 && line.Sizeof() )
            else
            {
               for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
               {
                  dToTGainListRight[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               } // for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platotgainright")==kTRUE)
         else if(line.Contains("platotoffsright")==kTRUE)
         {
           line = line(16,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
           if( uNbPlastics == 0 && line.Sizeof() )
            {
               if( 1 == uDebug )
               {
               cout<<"Warning: nbPlastics is 0 and plastics 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 < uNbPlastics; temp_index++)
               {
                  dTotOffsetListRight[temp_index] = ((TString)line(0,NB_DIGITS_PLASTICS_OFFSETS)).Atof();
                  line = line(NB_DIGITS_PLASTICS_OFFSETS,NB_DIGITS_PLASTICS_OFFSETS*NB_PLASTICS_MAX);
               }
            } // else of if( uNbPlastics == 0 && line.Sizeof() )
         } // else if(line.Contains("platotoffsright")==kTRUE)
      } // while( !setupFile.eof() )
   } //  else of if( setupFile.fail() == true)
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
      if( 1 == uPlasticTdcType[temp_index] )
      {
         uTdcPlasticLeftTot[temp_index]     = uTdcPlasticLeftTime[temp_index];
         if( 1 == uChannelPlasticLeftTime[temp_index]%2)
            uChannelPlasticLeftTot[temp_index] = uChannelPlasticLeftTime[temp_index] - 1;
            else uChannelPlasticLeftTot[temp_index] = uChannelPlasticLeftTime[temp_index] + 1;
         uTdcPlasticRightTot[temp_index]     = uTdcPlasticRightTime[temp_index];
         if( 1 == uChannelPlasticRightTot[temp_index]%2)
            uChannelPlasticRightTime[temp_index] = uChannelPlasticRightTot[temp_index] - 1;
            else uChannelPlasticRightTime[temp_index] = uChannelPlasticRightTot[temp_index] + 1;
      }

   return kTRUE;
}


Int_t TPlasticsParam::PrintOptions(){
   cout<<"********************* PLASTICS calibration parameters **********************"<<endl;
   cout<<"------------- Options for PLASTICS --------------"<<endl;
   cout<<"-------------------------------------------------"<<endl;
   if(1 == uDebug)
      cout<<"  Debug ON!!!!!!"<<endl;
      else cout<<"  Debug OFF"<<endl;
   if( -1 < iReferenceTdc && -1 < iReferenceChannel )
   {
      cout<<"  Chan. used for diamond: "<<iReferenceChannel;
      cout<<" in TDC "<<iReferenceTdc<<" in VFTX probably"<<endl;
   }
   if( -1 < iDiamondTdcVFTX && -1 < iDiamondChannelVFTX )
   {
      cout<<"  Chan. used for diamond: "<<iDiamondChannelVFTX;
      cout<<" in TDC "<<iDiamondTdcVFTX<<" in VFTX"<<endl;
   }
   cout<<"  Nb Plastics:               "<<uNbPlastics<<endl;
   cout<<"  Time Convertion factor:    "<<dTimeConv<<" cm/ns"<<endl;
   if( -1 == iTriggerSelection )
      cout<<"  Trigger selection:         OFF"<<endl;
      else cout<<"  Trigger selection:         Trigger "<<iTriggerSelection<<endl;
   cout<<"-------------------------------------------------"<<endl;
   cout<<"  Plastic:           |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Tdc Type           |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uPlasticTdcType[temp_index]<<" ";
   }
   cout<<endl<<"  Position X         |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<dPlasticsPosX[temp_index]<<" ";
   }
   cout<<endl<<"  Position Y         |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<dPlasticsPosY[temp_index]<<" ";
   }
   cout<<endl<<"  Position Z         |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<dPlasticsPosZ[temp_index]<<" ";
   }
   cout<<endl<<"  Orientation        |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uOrientation[temp_index]<<" ";
   }
   cout<<endl<<"  Width              |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<dPlasticWidth[temp_index]<<" ";
   }
   cout<<endl<<"  Length             |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<dPlasticLength[temp_index]<<" ";
   }
   cout<<endl<<"  Nb sides           |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uNbSides[temp_index]<<" ";
   }   cout<<endl<<"  Tdc Tot Left       |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uTdcPlasticLeftTot[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Tot Left   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uChannelPlasticLeftTot[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Time Left      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uTdcPlasticLeftTime[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Time Left  |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uChannelPlasticLeftTime[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Time Right     |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uTdcPlasticRightTime[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Time Right |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uChannelPlasticRightTime[temp_index]<<" ";
   }
   cout<<endl<<"  Tdc Tot Right      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uTdcPlasticRightTot[temp_index]<<" ";
   }
   cout<<endl<<"  Channel Tot Right  |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_INDEXES - 1);
      cout<<uChannelPlasticRightTot[temp_index]<<" ";
   }
   cout<<endl;
   cout<<"-------------------------------------------------"<<endl;
   cout<<"  Plastic:           |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Beam Profile Offset|->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dProfOffsetList[temp_index]<<" ";
   }
   cout<<endl<<"  Reference Offset   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dRefOffsetList[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( 0 < dMinimalTimeBetweenHits )
      cout<<"  Minimal Dist btwn TDC hit: "<<dMinimalTimeBetweenHits<<" ps (Merged if closer)"<<endl;
      else cout<<"  Minimal Dist btwn TDC hit: OFF"<<endl;
   cout<<"  Plastic:           |-      ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<temp_index<<" ";
   }
   cout<<endl<<"  Tot Gain Left      |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dToTGainListLeft[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Offset Left    |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dTotOffsetListLeft[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Gain Right     |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dToTGainListRight[temp_index]<<" ";
   }
   cout<<endl<<"  Tot Offset Right   |->     ";
   for(UInt_t temp_index = 0; temp_index < uNbPlastics; temp_index++)
   {
      cout.width(NB_DIGITS_PLASTICS_OFFSETS - 1);
      cout<<dTotOffsetListRight[temp_index]<<" ";
   }
   cout<<endl;
   cout<<"************************************************************************"<<endl;
   
   return 0;
}

Bool_t TPlasticsParam::CalculateLeftRightOffset()
{
   for( UInt_t uPlastic = 0; uPlastic < uNbPlastics; uPlastic++)
   {
      dOffsetListLeft[uPlastic]  = 0.5*( (dProfOffsetList[uPlastic])+2.*(dRefOffsetList[uPlastic]));
      dOffsetListRight[uPlastic] = 0.5*(-(dProfOffsetList[uPlastic])+2.*(dRefOffsetList[uPlastic]));
   }
   return kTRUE;
}