beamtime/tof-tdctest/go4/MATCHING/TMatchingParam.cxx (r4864/r4745)

Go to the documentation of this file.

/* Generated by Together */

#include "TMatchingParam.h"

#include "TGo4Version.h"
#if __GO4BUILDVERSION__ > 40502
   #include "go4iostream.h"
#else
   #include "Riostream.h"
#endif
 
// define default values
#define DEBUG_OPTION             0
#define TRIGGER_SELECTION       -1
#define TRIGGER_SEL_DEBUG_HIST   0
#define NB_EV_AVERAGE_EFFICIENCY 0
#define DETECTOR_OFFSET_TO_USE   1
#define HISTOS_ENABLE            0

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

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

Bool_t TMatchingParam::ReadOptions(){
   // Loading default values
   uDebug           = DEBUG_OPTION;
   iTriggerSelection = TRIGGER_SELECTION;                 // Trigger to be required for the matching
   uTriggerSelDebugHistos = TRIGGER_SEL_DEBUG_HIST;       // Disable/Enable the debug ;ode for trigger selection: 1 matching per trigger
   uNbEventAverageEfficiency = NB_EV_AVERAGE_EFFICIENCY;  // Number of event over which the efficiency of HD counter vs Rate is evaluated
   uDetectorOffsetToUse = DETECTOR_OFFSET_TO_USE;         // Select from which detector the offset for the dt matching should be taken: 0 = HD, 1 = Buc Ref

   uEnableSingleRpcHistos = HISTOS_ENABLE;            // Generic Single RPC <- Big switch, Disable also all other single RPC histos!!!
   uEnableRpcLeftRightTot = HISTOS_ENABLE;            // Enable all Histos comparing TOT on left and right side of fired strips
   uEnableMultiRpcHistos = HISTOS_ENABLE;             // Generic RPC - RPC  <- Big switch, Disable also all other RPC-RPC histos!!!
   uEnableMul1Mul1Histos = HISTOS_ENABLE;             // Enable all RPC-RPC histos on Multiplicity 1-1 events
   uEnableMul1MulXHistos = HISTOS_ENABLE;             // Enable all RPC-RPC histos on Multiplicity 1-X, X-1 events
   uEnableMul2Mul2Histos = HISTOS_ENABLE;             // Enable all RPC-RPC histos on Multiplicity 1-1, 2-1, 1-2, 2-2 events
   uEnableMul1Mul2Histos = HISTOS_ENABLE;             // Enable all RPC-RPC histos on Multiplicity 2-1, 1-2 events
   uEnableMulXMulYHistos = HISTOS_ENABLE;             // Enable all RPC-RPC histos on Multiplicity X-Y events

   uEnableSinglePlasticsHistos = HISTOS_ENABLE;       // Generic Single Plastic <- Big switch, Disable also all other single plastic histos!!!
   uEnablePlasticLeftRightTot = HISTOS_ENABLE;        // Enable all Histos comapring TOT on left and right side of Plastic hit
   uEnableMultiPlasticsHistos = HISTOS_ENABLE;        // Generic Plastic - Plastic <- Big switch, Disable also all other plastic-plastic histos!!!

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

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

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

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

         if(line(0) == '#' || line(0) == '*')
         {
            //Comment in the text file
            ;
         }
         else if(line.Contains("debug")==kTRUE)
         {
            uDebug = ((TString)line(5,10)).Atoi();
            if(1 == uDebug)
               cout<<"Debug ON!!!!!!"<<endl;
         }
         else if(line.Contains("triggersel")==kTRUE)
         {
            iTriggerSelection = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( -1 == iTriggerSelection)
                  cout<<"Trigger used for matching selection: OFF"<<endl;
                  else cout<<"Trigger used for matching selection: "<<iTriggerSelection<<endl;
            } // if(1 == uDebug)
         }
         else if(line.Contains("trigseldeb")==kTRUE)
         {
            uTriggerSelDebugHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uTriggerSelDebugHistos)
                  cout<<"Trigger selection debug histograms:  ON"<<endl;
                  else cout<<"Trigger selection debug histograms:  OFF"<<endl;
            } // if(1 == uDebug)
         }
         else if(line.Contains("nbevaveffy")==kTRUE)
         {
            uNbEventAverageEfficiency = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               cout<<"Nb of event averaging HD efficiency: "<<uNbEventAverageEfficiency<<endl;
            }
         }
         else if(line.Contains("matoffsdet")==kTRUE)
         {
            uDetectorOffsetToUse = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uDetectorOffsetToUse)
                  cout<<"Offset for dt in matching from:      Bucarest Reference counter"<<endl;
                  else cout<<"Offset for dt in matching from:      Heidelberg counter"<<endl;
            }
         }
         else if(line.Contains("hsinglerpc")==kTRUE)
         {
            uEnableSingleRpcHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableSingleRpcHistos)
                  cout<<"Single RPC Histograms:               ON"<<endl;
                  else cout<<"Single RPC Histograms:               OFF"<<endl;
            } // if(1 == uDebug)
         }
         else if(line.Contains("htotrpclr")==kTRUE)
         {
            uEnableRpcLeftRightTot = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableRpcLeftRightTot)
                  cout<<"Rpc Left-Right Tot Histograms:       ON"<<endl;
                  else cout<<"Rpc Left-Right Tot Histograms:       OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("htotrpclr")==kTRUE)
         else if(line.Contains("hrpcrpc")==kTRUE)
         {
            uEnableMultiRpcHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMultiRpcHistos)
                  cout<<"RPC-RPC Histograms:                  ON"<<endl;
                  else cout<<"RPC-RPC Histograms:                  OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcrpc")==kTRUE)
         else if(line.Contains("hrpcm1m1")==kTRUE)
         {
            uEnableMul1Mul1Histos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMul1Mul1Histos)
                  cout<<"Rpc1-Mul1 Rpc2-Mul1 Histograms:      ON"<<endl;
                  else cout<<"Rpc1-Mul1 Rpc2-Mul1 Histograms:      OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcm1m1")==kTRUE)
         else if(line.Contains("hrpcm1mx")==kTRUE)
         {
            uEnableMul1MulXHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMul1MulXHistos)
                  cout<<"Rpc1-Mul1 Rpc2-MulX Histograms:      ON"<<endl;
                  else cout<<"Rpc1-Mul1 Rpc2-MulX Histograms:      OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcm1mx")==kTRUE)
         else if(line.Contains("hrpcm2m2")==kTRUE)
         {
            uEnableMul2Mul2Histos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMul2Mul2Histos)
                  cout<<"Rpc1-Mul<=2 Rpc2-Mul<=2 Histograms:  ON"<<endl;
                  else cout<<"Rpc1-Mul<=2 Rpc2-Mul<=2 Histograms:  OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcm2m2")==kTRUE)
         else if(line.Contains("hrpcm1m2")==kTRUE)
         {
            uEnableMul1Mul2Histos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMul1Mul2Histos)
                  cout<<"Rpc1-Mul1 Rpc2-Mul2 Histograms:      ON"<<endl;
                  else cout<<"Rpc1-Mul1 Rpc2-Mul2 Histograms:      OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcm1m2")==kTRUE)
         else if(line.Contains("hrpcmxmy")==kTRUE)
         {
            uEnableMulXMulYHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMulXMulYHistos)
                  cout<<"Rpc1-Mul X Rpc2-Mul Y Histograms:    ON"<<endl;
                  else cout<<"Rpc1-Mul X Rpc2-Mul Y Histograms:    OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hrpcmxmy")==kTRUE)
         else if(line.Contains("hsinglepla")==kTRUE)
         {
            uEnableSinglePlasticsHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableSinglePlasticsHistos)
                  cout<<"Single Plastic Histograms:           ON"<<endl;
                  else cout<<"Single Plastic Histograms:           OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hsinglepla")==kTRUE)
         else if(line.Contains("htotplalr")==kTRUE)
         {
            uEnablePlasticLeftRightTot = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnablePlasticLeftRightTot)
                  cout<<"Plastic Left-Right Tot Histograms:   ON"<<endl;
                  else cout<<"Plastic Left-Right Tot Histograms:   OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("htotplalr")==kTRUE)
         else if(line.Contains("hplapla")==kTRUE)
         {
            uEnableMultiPlasticsHistos = ((TString)line(11,10)).Atoi();
            if(1 == uDebug)
            {
               if( 1 == uEnableMultiPlasticsHistos)
                  cout<<"Plastic-Plastic Histograms:          ON"<<endl;
                  else cout<<"Plastic-Plastic Histograms:          OFF"<<endl;
            } // if(1 == uDebug)
         } // else if(line.Contains("hplapla")==kTRUE)
      } // while( !setupFile.eof() )
   } //  else of if( setupFile.fail() == true)
   
   PrintOptions(); //print out all values to see, if the setting has worked correctly
   return kTRUE;
}


Int_t TMatchingParam::PrintOptions(){
   cout<<"********************* MATCHING calibration parameters **********************"<<endl;
   if(1 == uDebug)
      cout<<"Debug ON!!!!!!"<<endl;
      else cout<<"Debug OFF"<<endl;
   if( -1 == iTriggerSelection)
      cout<<"Trigger used for matching selection: OFF"<<endl;
      else cout<<"Trigger used for matching selection: "<<iTriggerSelection<<endl;
   if( 1 == uTriggerSelDebugHistos)
      cout<<"Trigger selection debug histograms:  ON"<<endl;
      else cout<<"Trigger selection debug histograms:  OFF"<<endl;
   cout<<"Nb of event averaging HD efficiency: "<<uNbEventAverageEfficiency<<endl;
   if( 1 == uDetectorOffsetToUse)
      cout<<"Offset for dt in matching from:      Bucarest Reference counter"<<endl;
      else cout<<"Offset for dt in matching from:      Heidelberg counter"<<endl;
   if( 1 == uEnableSingleRpcHistos )
   {
      cout<<"Single RPC Histograms:               ON"<<endl;
      if( 1 == uEnableRpcLeftRightTot)
         cout<<"Rpc Left-Right Tot Histograms:       ON"<<endl;
   }
      else cout<<"Single RPC Histograms:               OFF"<<endl;
   if( 1 == uEnableMultiRpcHistos )
   {
      cout<<"RPC-RPC Histograms:                  ON"<<endl;
      if( 1 == uEnableMul1Mul1Histos)
         cout<<"Rpc1-Mul1 Rpc2-Mul1 Histograms:      ON"<<endl;
      if( 1 == uEnableMul1MulXHistos)
         cout<<"Rpc1-Mul1 Rpc2-MulX Histograms:      ON"<<endl;
      if( 1 == uEnableMul2Mul2Histos)
         cout<<"Rpc1-Mul<=2 Rpc2-Mul<=2 Histograms:  ON"<<endl;
      if( 1 == uEnableMul1Mul2Histos)
         cout<<"Rpc1-Mul1 Rpc2-Mul2 Histograms:      ON"<<endl;
      if( 1 == uEnableMulXMulYHistos)
         cout<<"Rpc1-Mul X Rpc2-Mul Y Histograms:    ON"<<endl;
   }
      else cout<<"RPC-RPC Histograms:                  OFF"<<endl;
   if( 1 == uEnableSinglePlasticsHistos )
   {
      cout<<"Single Plastic Histograms:           ON"<<endl;
      if( 1 == uEnablePlasticLeftRightTot)
         cout<<"Plastic Left-Right Tot Histograms:   ON"<<endl;
   }
      else cout<<"Single Plastic Histograms:           OFF"<<endl;
   if( 1 == uEnableMultiPlasticsHistos )
      cout<<"Plastic-Plastic Histograms:          ON"<<endl;
      else cout<<"Plastic-Plastic Histograms:          OFF"<<endl;
   cout<<"************************************************************************"<<endl;
   
   return 0;
}