PWG2/FORWARD/analysis2/Pass2.C

   1 /**
   2  * Read in AOD and generate @f$ dN/d\eta@f$ for the selected
   3  * trigger classes and vertex ranges
   4  *
   5  * @param file     Input file (AOD)
   6  * @param triggers Triggers to investigate
   7  * @param energy   Energy (only used for comparisons)
   8  * @param vzMin    Minimum interaction point z coordinate
   9  * @param vzMax    Maximum interaction point z coordinate
  10  * @param rebin    How many bins to group
  11  * @param title    Title to put on the plot
  12  * @param hhd      Whether to do HHD comparison
  13  * @param comp     Whether to do comparisons
  14  *
  15  * @ingroup pwg2_forward_analysis_scripts
  16  */
  17 void
  18 Pass2(const char* file="AliAODs.root",
  19       const char* triggers="INEL",
  20       Int_t       energy=900,
  21       Double_t    vzMin=-10,
  22       Double_t    vzMax=10,
  23       Int_t       rebin=5,
  24       const char* title="",
  25       bool        hhd=true,
  26       bool        comp=true)
  27 {
  28   gROOT->LoadMacro("$ALICE_ROOT/PWG2/FORWARD/analysis2/scripts/Compile.C");
  29   Compile("$ALICE_ROOT/PWG2/FORWARD/analysis2/DrawRes.C","g");
  30
  31   Int_t trgMask;
  32   TString     trgs(triggers);
  33   trgs.ToUpper();
  34   TObjString* trg;
  35   TIter       next(trgs.Tokenize(" ,|"));
  36   while ((trg = static_cast<TObjString*>(next()))) {
  37     TString s(trg->GetString());
  38     if      (s.IsNull()) continue;
  39     if      (s.CompareTo("INEL")  == 0) trgMask = AliAODForwardMult::kInel;
  40     else if (s.CompareTo("INEL>0")== 0) trgMask = AliAODForwardMult::kInelGt0;
  41     else if (s.CompareTo("NSD")   == 0) trgMask = AliAODForwardMult::kNSD;
  42     else
  43       Warning("Pass2", "Unknown trigger %s", s.Data());
  44   }
  45   if (trgMask == 0) {
  46     trgMask = 1;
  47     trgs.Append("INEL");
  48   }
  49   TString tit(title);
  50   tit.ReplaceAll("@", " ");
  51
  52   printf("--------------------------------------\n"
  53          "Settings for this:\n"
  54          "  Input AOD:    %s\n"
  55          "  Vertex range: %+4.1f -> %+4.1f cm\n"
  56          "  Rebinning:    %d\n"
  57          "  Trigger mask: 0x%02x (%s)\n"
  58          "  Energy:       %dGeV\n"
  59          "  Title:        %s\n"
  60          "  HHD comp.:    %s\n"
  61          "  Other comp.:  %s\n"
  62          "--------------------------------------\n",
  63          file, vzMin, vzMax, rebin, trgMask, trgs.Data(), energy, tit.Data(),
  64          hhd ? "yes" : "no", comp ? "yes" : "no");
  65
  66   DrawRes dr;
  67   TStopwatch t;
  68   t.Start();
  69   dr.Run(file, vzMin, vzMax, rebin, trgMask, energy, tit.Data(), hhd, comp);
  70   t.Stop();
  71   t.Print();
  72 }
  73 //
  74 // EOF
  75 //
  76
  77