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 ConfigPhiPbPbNoPID
13 AliRsnMiniAnalysisTask *task,
17 AliRsnCutSet *cutsPair
21 if (strlen(suffix) > 0) suffix = Form("_%s", suffix);
24 // -- Define track cuts -------------------------------------------------------------------------
27 // BB parameterization depends on data sample (MC, data)
30 bbPar[0] = 2.15898 / 50.0;
32 bbPar[2] = 3.40030E-9;
36 bbPar[0] = 1.41543 / 50.0;
38 bbPar[2] = 5.0411E-11;
44 AliRsnCutKaonForPhi2010 *cut = new AliRsnCutKaonForPhi2010(Form("cut%s", suffix), 3.0, 3.0, 0.8);
46 // setup (set manually the TPC PID)
47 cut->SetMode(AliRsnCutKaonForPhi2010::kQuality);
50 AliRsnCutSet *cutSet = new AliRsnCutSet(Form("set%s", suffix), AliRsnTarget::kDaughter);
52 cutSet->SetCutScheme(cut->GetName());
55 Int_t icut = task->AddTrackCuts(cutSet);
56 ::Info("Config", "Cut ID = %d", icut);
59 // -- Values ------------------------------------------------------------------------------------
62 /* invariant mass */ Int_t imID = task->CreateValue(AliRsnMiniValue::kInvMass, kFALSE);
63 /* IM resolution */ Int_t resID = task->CreateValue(AliRsnMiniValue::kInvMassRes, kTRUE);
64 /* transv. momentum */ Int_t ptID = task->CreateValue(AliRsnMiniValue::kPt, kFALSE);
65 /* centrality */ Int_t centID = task->CreateValue(AliRsnMiniValue::kMult, kFALSE);
68 // -- Create all needed outputs -----------------------------------------------------------------
71 // use an array for more compact writing, which are different on mixing and charges
76 Bool_t use [6] = { 1 , 0 , 0 , 0 , isMC , isMC };
77 Bool_t useIM [6] = { 1 , 1 , 1 , 1 , 1 , 0 };
78 TString name [6] = {"Unlike", "Mixing", "LikePP", "LikeMM", "Trues" , "Res" };
79 TString comp [6] = {"PAIR" , "MIX" , "PAIR" , "PAIR" , "TRUE" , "TRUE" };
80 TString output [6] = {"SPARSE", "SPARSE", "SPARSE", "SPARSE", "SPARSE", "SPARSE" };
81 Char_t charge1 [6] = {'+' , '+' , '+' , '-' , '+' , '+' };
82 Char_t charge2 [6] = {'-' , '-' , '+' , '-' , '-' , '-' };
83 Int_t cutID [6] = { icut , icut , icut , icut , icut , icut };
85 for (Int_t i = 0; i < 6; i++) {
86 if (!use[i]) continue;
88 AliRsnMiniOutput *out = task->CreateOutput(Form("phi_%s%s", name[i].Data(), suffix), output[i].Data(), comp[i].Data());
90 out->SetCutID(0, cutID[i]);
91 out->SetCutID(1, cutID[i]);
92 out->SetDaughter(0, AliRsnDaughter::kKaon);
93 out->SetDaughter(1, AliRsnDaughter::kKaon);
94 out->SetCharge(0, charge1[i]);
95 out->SetCharge(1, charge2[i]);
96 out->SetMotherPDG(333);
97 out->SetMotherMass(1.019455);
99 out->SetPairCuts(cutsPair);
100 // axis X: invmass (or resolution)
102 out->AddAxis(imID, 500, 0.9, 1.4);
104 out->AddAxis(resID, 200, -0.02, 0.02);
105 // axis Y: transverse momentum
106 out->AddAxis(ptID, 100, 0.0, 10.0);
107 // axis Z: centrality
108 out->AddAxis(centID, 100, 0.0, 100.0);