2 // *** Configuration script for phi->KK analysis with 2010 runs ***
4 // A configuration script for RSN package needs to define the followings:
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
11 Bool_t ConfigKStarSimple
13 AliRsnMiniAnalysisTask *task,
18 const char *computationType,
23 AliRsnCutSet *cutsPair
27 // -- Define track cuts -------------------------------------------------------------------------
30 // integrated pion cut
31 AliRsnCutDaughterKStar2010PP *cutPi = new AliRsnCutDaughterKStar2010PP("cutPionForKStar", AliPID::kPion);
33 AliRsnCutSet *cutSetPi = new AliRsnCutSet("setPionForKStar", AliRsnTarget::kDaughter);
34 cutSetPi->AddCut(cutPi);
35 cutSetPi->SetCutScheme(cutPi->GetName());
37 Int_t iCutPi = task->AddTrackCuts(cutSetPi);
39 // integrated kaon cut
40 AliRsnCutDaughterKStar2010PP *cutK = new AliRsnCutDaughterKStar2010PP("cutKaonForKStar", AliPID::kKaon);
42 AliRsnCutSet *cutSetK = new AliRsnCutSet("setKaonForKStar", AliRsnTarget::kDaughter);
43 cutSetK->AddCut(cutK);
44 cutSetK->SetCutScheme(cutK->GetName());
46 Int_t iCutK = task->AddTrackCuts(cutSetK);
49 // -- Values ------------------------------------------------------------------------------------
52 /* invariant mass */ Int_t imID = task->CreateValue(AliRsnMiniValue::kInvMass, kFALSE);
53 /* IM resolution */ Int_t resID = task->CreateValue(AliRsnMiniValue::kInvMassRes, kTRUE);
54 /* transv. momentum */ Int_t ptID = task->CreateValue(AliRsnMiniValue::kPt, kFALSE);
55 /* centrality */ Int_t centID = task->CreateValue(AliRsnMiniValue::kMult, kFALSE);
58 // -- Create all needed outputs -----------------------------------------------------------------
62 AliRsnMiniOutput *out = task->CreateOutput(name, outType, computationType);
65 out->SetCutID(0, iCutK);
66 out->SetCutID(1, iCutPi);
69 out->SetDaughter(0, AliRsnDaughter::kKaon);
70 out->SetDaughter(1, AliRsnDaughter::kPion);
71 out->SetCharge(0, charge1);
72 out->SetCharge(1, charge2);
74 // resonance properties
75 out->SetMotherPDG(313);
76 out->SetMotherMass(0.896);
79 out->SetPairCuts(cutsPair);
81 // axis X: invmass (or resolution)
83 out->AddAxis(imID, 90, 0.6, 1.5);
85 out->AddAxis(resID, 200, -0.02, 0.02);
87 // axis Y: transverse momentum
88 out->AddAxis(ptID, 100, 0.0, 10.0);