PWGLF/FORWARD/analysis2/AddTaskCreateRespMatr.C

   1 /**
   2  * @file   AddTaskCreateRespMatr.C
   3  * @author Valentina Zaccolo
   4  * @date   Fri Jan 11 15:34:26 2013
   5  *
   6  * @brief Script to add a task to create response matrices
   7  *
   8  *
   9  */
  10 AliAnalysisTask*
  11 AddTaskCreateRespMatr(const char* trig="V0AND",
  12                       Double_t vzMin=-4,
  13                       Double_t vzMax=4)
  14 {
  15   // Make our object.  2nd argumenent is absolute max Eta
  16   // 3rd argument is absolute max Vz
  17   AliForwardCreateResponseMatrices* task =
  18     new AliForwardCreateResponseMatrices("ResponseMatrices");
  19   // Set the Vertex Range to Use
  20   task->SetIpZRange(vzMin, vzMax);
  21   // Set the Trigger Mask to Use (INEL, NSD, ...)
  22   task->SetTriggerMask(trig);
  23   // Set the Number of Bins
  24   //  task->SetNBins(nBins);
  25
  26   //Add Full eta-ranges
  27   task->AddBin(-3.4,5.1);
  28
  29   //Add Symmetric eta bins.
  30   Double_t limits[] = { 3.4, 3.0, 2.5, 2.4, 2.0, 1.5, 1.4, 1.0, 0.5, 0. };
  31   Double_t* limit = limits;
  32   while ((*limit) > 0.1) {
  33     task->AddBin(-(*limit), +(*limit));
  34     // task->AddBin(0,+(*limit));
  35     // task->AddBin(0,-(*limit));
  36     limit++;
  37   }
  38   // task->AddBin(0,5.0);
  39   // task->AddBin(0,4.5);
  40   // task->AddBin(0,4.0);
  41   // task->AddBin(0,3.5);
  42
  43   // Add 0.5 eta intervals
  44   // for (Double_t l = -3; l < 5; l += 0.5) task->AddBin(l, l+.5);
  45
  46   // Add 0.20 eta intervals
  47   // for (Double_t l = -3; l < 5; l += 0.2) task->AddBin(l, l+.2);
  48
  49   task->Connect();
  50   return task;
  51 }
  52
  53
  54 //________________________________________________________________________
  55 //
  56 // EOF
  57 //