1 AliLeadingV0Correlation* AddTaskLV0Correlation(TString fListName = "LV0Correlation",
2 TString fCollisiontype = "PP",
3 Bool_t fAnalysisMC = 1,
5 Bool_t fRemoveAutoCorr = 0,
7 Int_t fFilterBit = 256,
8 Double_t fRapidityCut = 0.75,
9 Int_t fmaxEventsinPool = 2000,
10 Int_t fminTracksinPool = 100,
11 Int_t fMinEventsToMix = 5,
12 Double_t fV0radius = 0.3,
13 Double_t fV0PostoPVz = 0.060,
14 Double_t fV0NegtoPVz = 0.060,
15 Double_t fDCAV0Daughters = 1.00,
16 Double_t fCPAK0 = 0.960,
17 Double_t fCPALam = 0.997,
18 Double_t fRejectLamK0 = 0.005,
19 Double_t fRejectK0Lam = 0.010,
20 Double_t fSigmaPID = 3.0,
21 Double_t fCutCTK0 = 20.0,
22 Double_t fCutCTLa = 30.0,
23 Double_t fMassCutK0 = 0.0105,
24 Double_t fMassCutLa = 0.0105,
25 Double_t fTriglow = 6.0,
26 Double_t fTrighigh = 12.0)
28 // Get the current analysis manager.
29 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
30 if (!mgr) {Error("AddTaskLV0Correlation.C", "No Analysis Manager ");return 0;}
32 //PVz Binning for pool PP or PbPb
33 Double_t pvzbinlimits[] = {-12,-11,-10,-9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12};
34 Int_t pvzbinnumb = sizeof(pvzbinlimits)/sizeof(Double_t) - 1;
36 //Mult Binning for pool pp
37 Double_t cent_mult_binlimitsPP[] = { 0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
38 100,110,120,130,140,150,160,170,180,190,
39 200,210,220,230,240,250,260,270,280,290,300};
41 Int_t cent_mult_bin_numbPP = sizeof(cent_mult_binlimitsPP)/sizeof(Double_t) - 1;
43 //Cent Binning for pool PbPb
44 Double_t cent_mult_binlimitsPbPb[] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,
45 14,15,16,17,18,19,20,21,22,23,24,25,26,27,
46 28,29,30,31,32,33,34,35,36,37,38,39,40,42,
47 44,46,48,50,52,54,56,58,60,65,70,75,80,90};
49 Int_t cent_mult_bin_numbPbPb = sizeof(cent_mult_binlimitsPbPb)/sizeof(Double_t) - 1;
51 AliLeadingV0Correlation *myTask = new AliLeadingV0Correlation(fListName.Data());
53 myTask->SetCollidingSystem(fCollisiontype);
54 myTask->SetMCAnalysis(fAnalysisMC);
55 myTask->SetCase(fCase);
56 myTask->SetRemoveAutoCorr(fRemoveAutoCorr);
57 myTask->SetMaxNEventsInPool(fmaxEventsinPool);
58 myTask->SetMinNTracksInPool(fminTracksinPool);
59 myTask->SetMinEventsToMix(fMinEventsToMix);
60 myTask->SetPoolPVzBinLimits(pvzbinnumb,pvzbinlimits);
61 if(fCollisiontype=="PP")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPP,cent_mult_binlimitsPP);
62 if(fCollisiontype=="PbPb")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPbPb,cent_mult_binlimitsPbPb);
63 myTask->SetPrimeryVertexCut(fPVzCut);
64 myTask->SetFilterBit(fFilterBit);
65 myTask->SetCutRap(fRapidityCut);
66 myTask->SetV0Radius(fV0radius);
67 myTask->SetV0PostoPVz(fV0PostoPVz);
68 myTask->SetV0NegtoPVz(fV0NegtoPVz);
69 myTask->SetDCAV0Daughters(fDCAV0Daughters);
70 myTask->SetCPAK0(fCPAK0);
71 myTask->SetCPALam(fCPALam);
72 myTask->SetRejectLamK0(fRejectLamK0);
73 myTask->SetRejectK0Lam(fRejectK0Lam);
74 myTask->SetSigmaPID(fSigmaPID);
75 myTask->SetCTK0(fCutCTK0);
76 myTask->SetCTLa(fCutCTLa);
77 myTask->SetMassCutK0(fMassCutK0);
78 myTask->SetMassCutLa(fMassCutLa);
79 myTask->SetTrigLow(fTriglow);
80 myTask->SetTrigHigh(fTrighigh);
84 // Create containers for input/output
85 TString outputFileName = AliAnalysisManager::GetCommonFileName();
86 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
87 AliAnalysisDataContainer *coutput = mgr->CreateContainer(fListName.Data(), TList::Class(),AliAnalysisManager::kOutputContainer,outputFileName);
89 //____________________________________________//
90 mgr->ConnectInput(myTask,0,cinput);
91 mgr->ConnectOutput(myTask,1,coutput);