]>
Commit | Line | Data |
---|---|---|
caef2e3c | 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 | } |