]>
Commit | Line | Data |
---|---|---|
8b48b94e | 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 ConfigLStarPP_MC | |
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 = "SPARSE"; | |
41 | // if (!isPP) mode = "SPARSE"; | |
42 | ||
43 | // create output | |
44 | AliRsnMiniOutput *out = task->CreateOutput(Form("LStar_TrueMC1%s", suffix), mode.Data(), "MOTHER"); | |
45 | // selection settings | |
46 | out->SetDaughter(0, AliRsnDaughter::kProton); | |
47 | out->SetDaughter(1, AliRsnDaughter::kKaon); | |
48 | out->SetCharge(0, '+'); | |
49 | out->SetCharge(1, '-'); | |
50 | out->SetMotherPDG(3124); | |
51 | out->SetMotherMass(1.520); | |
52 | // pair cuts | |
53 | out->SetPairCuts(cutsPair); | |
54 | // binnings | |
55 | out->AddAxis(imID, 80, 1.0, 2.0); | |
56 | //S.K. out->AddAxis(ptID, 100, 0.0, 10.0); | |
57 | out->AddAxis(ptID, 100, 0.0, 10.0); | |
58 | //S.K. if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0); | |
59 | out->AddAxis(centID, 10, 0.0, 100.0); | |
60 | ||
61 | // create output | |
62 | AliRsnMiniOutput *out = task->CreateOutput(Form("LStar_TrueMC2%s", suffix), mode.Data(), "MOTHER"); | |
63 | // selection settings | |
64 | out->SetDaughter(0, AliRsnDaughter::kProton); | |
65 | out->SetDaughter(1, AliRsnDaughter::kKaon); | |
66 | out->SetCharge(0, '-'); | |
67 | out->SetCharge(1, '+'); | |
68 | out->SetMotherPDG(-3124); | |
69 | out->SetMotherMass(1.520); | |
70 | // pair cuts | |
71 | out->SetPairCuts(cutsPair); | |
72 | // binnings | |
73 | out->AddAxis(imID, 100, 1.0, 2.0); | |
74 | //S.K. out->AddAxis(ptID, 100, 0.0, 10.0); | |
75 | out->AddAxis(ptID, 100, 0.0, 10.0); | |
76 | //S.K. if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0); | |
77 | out->AddAxis(centID, 10, 0.0, 100.0); | |
78 | ||
79 | // create output | |
80 | AliRsnMiniOutput *outm = task->CreateOutput(Form("Ls_Mother"), mode.Data(), "MOTHER"); | |
81 | outm->SetDaughter(0, AliRsnDaughter::kProton); | |
82 | outm->SetDaughter(1, AliRsnDaughter::kKaon); | |
83 | outm->SetMotherPDG(3124); | |
84 | outm->SetMotherMass(1.520); | |
85 | // pair cuts | |
86 | outm->SetPairCuts(cutsPair); | |
87 | // binnings | |
88 | outm->AddAxis(imID, 100, 1.0, 2.0); | |
89 | outm->AddAxis(ptID, 100, 0.0, 10.0); | |
90 | if (!isPP){ | |
91 | outm->AddAxis(centID, 100, 0.0, 100.0); | |
92 | } else { | |
93 | outm->AddAxis(centID, 400, 0.0, 400.0); | |
94 | } | |
95 | ||
96 | AliRsnMiniOutput *outm = task->CreateOutput(Form("Ls_AntiMother"), mode.Data(), "MOTHER"); | |
97 | outm->SetDaughter(0, AliRsnDaughter::kProton); | |
98 | outm->SetDaughter(1, AliRsnDaughter::kKaon); | |
99 | outm->SetMotherPDG(-3124); | |
100 | outm->SetMotherMass(1.520); | |
101 | // pair cuts | |
102 | outm->SetPairCuts(cutsPair); | |
103 | // binnings | |
104 | outm->AddAxis(imID, 100, 1.0, 2.0); | |
105 | outm->AddAxis(ptID, 100, 0.0, 10.0); | |
106 | if (!isPP){ | |
107 | outm->AddAxis(centID, 100, 0.0, 100.0); | |
108 | } else { | |
109 | outm->AddAxis(centID, 400, 0.0, 400.0); | |
110 | } | |
111 | return kTRUE; | |
112 | } |