PWGLF/RESONANCES/macros/mini/ConfigKStarMC.C

   1 //
   2 // *** Configuration script for phi->KK analysis with 2010 runs ***
   3 //
   4 // A configuration script for RSN package needs to define the followings:
   5 //
   6 // (1) decay tree of each resonance to be studied, which is needed to select
   7 //     true pairs and to assign the right mass to all candidate daughters
   8 // (2) cuts at all levels: single daughters, tracks, events
   9 // (3) output objects: histograms or trees
  10 //
  11 Bool_t ConfigKStarMC
  12 (
  13    AliRsnMiniAnalysisTask *task,
  14    Bool_t                  isPP,
  15    const char             *suffix,
  16    AliRsnCutSet           *cutsPair
  17 )
  18 {
  19    // manage suffix
  20    if (strlen(suffix) > 0) suffix = Form("_%s", suffix);
  21
  22    //
  23    // -- Define track cuts -------------------------------------------------------------------------
  24    //
  25
  26    /*** EMPTY FOR TRUE PAIRS COMPUTATION ***/
  27
  28    //
  29    // -- Values ------------------------------------------------------------------------------------
  30    //
  31
  32    /* invariant mass   */ Int_t imID   = task->CreateValue(AliRsnMiniValue::kInvMass, kFALSE);
  33    /* transv. momentum */ Int_t ptID   = task->CreateValue(AliRsnMiniValue::kPt, kFALSE);
  34    /* centrality       */ Int_t centID = task->CreateValue(AliRsnMiniValue::kMult, kFALSE);
  35
  36    //
  37    // -- Create all needed outputs -----------------------------------------------------------------
  38    //
  39
  40    TString mode = "HIST";
  41    if (!isPP) mode = "SPARSE";
  42
  43    // create output
  44    AliRsnMiniOutput *out = task->CreateOutput(Form("kstar_TrueMC%s", suffix), mode.Data(), "MOTHER");
  45    // selection settings
  46    out->SetDaughter(0, AliRsnDaughter::kKaon);
  47    out->SetDaughter(1, AliRsnDaughter::kPion);
  48    out->SetMotherPDG(313);
  49    out->SetMotherMass(0.896);
  50    // pair cuts
  51    out->SetPairCuts(cutsPair);
  52    // binnings
  53    out->AddAxis(imID,  90, 0.6,  1.5);
  54    out->AddAxis(ptID, 100, 0.0, 10.0);
  55    if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0);
  56
  57    return kTRUE;
  58 }