hrich700taskset.cc

Go to the documentation of this file.

//////////////////////////////////////////////////////////////////////////////
//
// $Id: $
//
//*-- Author  : Dan Magestro
//*-- Revised : Martin Jurkovic <martin.jurkovic@ph.tum.de> 2010
//
//_HADES_CLASS_DESCRIPTION
//////////////////////////////////////////////////////////////////////////////
//
//  HRich700TaskSet
//
//  This HTaskSet contains the task sets for the Rich detector.
//  It can be distinguished between "real" and "simulation" task set by
//  passing the corresponding string to the first argument of the make()
//  function. Embedding is automatically recognized by return value of
//  gHades->getEmbeddingMode(). For embedding ONLY the "real" task set
//  has to be selected, otherwise the program will finish with an error.
//
//  To control the operation of the RICH tasks, HRich700TaskSet has following
//  switches with preset default values:
//    - inoiseOn   = kTRUE;
//    - iSkipEvent = kFALSE;
//
//  The simulation of noise in HRichDigitizer behaves in pure
//  simulation differently as compared to the embedding modus.
//     * pure simulation: simulation of both correlated and random noise is ON
//     * embedding:       only simulation of correlated noise is ON
//
//  Default values can be changed by passing corresponding strings
//  to the second argument of make() function. For more details,
//  please read the documentation for this function below.
//
//  BEGIN_HTML For more information, see
//  <a target="_top" href="http://webdb.gsi.de/pls/hades_webdbs/hanal.hal_posting_query.show_posting_text?p_id=436">
//  this analysis logbook entry<a>. END_HTML
//
//////////////////////////////////////////////////////////////////////////////


#include "TObjString.h"

#include "hades.h"
#include "hrich700taskset.h"
#include "hrich700digitizer.h"
#include "hrich700calibrater.h"
#include "hrich700ringfinderhough.h"

using namespace std;

ClassImp(HRich700TaskSet)


HRich700TaskSet::HRich700TaskSet(void)
   : HTaskSet()
{
   // Default constructor
    fNoRingFinder = kFALSE;
    fNoCalibrater = kFALSE;

}

HRich700TaskSet::HRich700TaskSet(const Text_t name[], const Text_t title[])
   : HTaskSet(name, title)
{
    // Constructor
    fNoRingFinder = kFALSE;
    fNoCalibrater = kFALSE;
}

HRich700TaskSet::~HRich700TaskSet(void)
{
   // Destructor.
}


HTask*
HRich700TaskSet::make(const Char_t *select, const Option_t *option)
{
// Returns a pointer to the Rich task or taskset specified by 'select'.
// There are two possibilities to be passed:
//   - "real" for real data analysis and embedding
//   - "simulation" for simulation only
//
// Following options can be passed to the HRich700TaskSet to control the
// behavior of the RICH analysis and simulation:
//  NORINGFINDER   - HRich700RingFinderHough will be disabled
//  NOCALIBRATER   - HRich700calibrater will be disabled
//

   HTaskSet* tasks      = new HTaskSet("Rich", "List of Rich tasks");
   TString   sel        = select;
   TString   opt        = option;

   parseArguments(opt);

   sel.ToLower();
   if (0 == sel.CompareTo("simulation") &&
       0 != gHades->getEmbeddingMode()) {
      Error("make", "embeddingMode > 0 in simulation NOT ALLOWED!");
      return NULL;
   }

   if (0 == sel.CompareTo("real"))    // real + embedding
   {
       if(!fNoCalibrater){
           tasks->add(new HRich700Calibrater("rich.calibrater",  "Rich Calibrater"));
       }
       if( 0 != gHades->getEmbeddingMode()){ //  embedding
           tasks->add(new HRich700Digitizer("rich.digi",  "Rich digitizer"));
       }
   }

   if (0 == sel.CompareTo("simulation"))
   {
       tasks->add(new HRich700Digitizer("rich.digi",  "Rich digitizer"));
   }

   if(!fNoRingFinder)// hit finder is the same in all cases
   {
       tasks->add(new HRich700RingFinderHough("rich.ringfinder",  "Rich ring finder"));
   }

   return tasks;
}

void HRich700TaskSet::parseArguments(TString s1)
{

   s1.ToLower();
   s1.ReplaceAll(" ", "");


   if (0 != s1.Length()) {

      // iterate over the list of arguments and compare it
       // to the known key words.

       TObjArray* myarguments = s1.Tokenize(",");
       TIterator*   myiter     = NULL;
       TObjString*  stemp      = NULL;
       TString      argument;

       myiter = myarguments->MakeIterator();
       myiter->Reset();
       while (0 != (stemp = static_cast<TObjString*>(myiter->Next()))) {
           argument = stemp->GetString();
           Info("parseArguments()","option: %s", argument.Data() );
           if (0 == argument.CompareTo("noringfinder")) {
               fNoRingFinder = kTRUE;
           }
           else if (0 == argument.CompareTo("nocalibrater")) {
               fNoCalibrater = kTRUE;
           } else {
               Error("parseArguments","Unknown Option %s found!",argument.Data());
           }
       }
       delete myiter;
       myiter = NULL;

       myarguments->Delete();
       delete myarguments;

   }
}