Removed obsolete comment
[u/mrichter/AliRoot.git] / PWGCF / Correlations / macros / AddTaskLV0Correlation.C
CommitLineData
da6ac661 1AliLeadingV0Correlation* AddTaskLV0Correlation(TString fListName = "LV0Correlation",
5b607771 2 TString fCollisiontype = "PP",
3 Bool_t fAnalysisMC = 1,
da6ac661 4 Int_t fCase = 2,
da6ac661 5 Bool_t fRemoveAutoCorr = 0,
5b607771 6 Double_t fPVzCut = 10,
7 Int_t fFilterBit = 256,
8 Double_t fRapidityCut = 0.75,
da6ac661 9 Int_t fmaxEventsinPool = 2000,
5b607771 10 Int_t fminTracksinPool = 100,
da6ac661 11 Int_t fMinEventsToMix = 5,
5b607771 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,
da6ac661 20 Double_t fSigmaPID = 3.0,
21 Double_t fCutCTK0 = 20.0,
22 Double_t fCutCTLa = 30.0,
23 Double_t fMassCutK0 = 0.0105,
5b607771 24 Double_t fMassCutLa = 0.0105,
25 Double_t fTriglow = 6.0,
26 Double_t fTrighigh = 12.0)
ff0805a7 27{
28 // Get the current analysis manager.
29 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
30 if (!mgr) {Error("AddTaskLV0Correlation.C", "No Analysis Manager ");return 0;}
31
da6ac661 32 //PVz Binning for pool PP or PbPb
bf855ed5 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};
ff0805a7 34 Int_t pvzbinnumb = sizeof(pvzbinlimits)/sizeof(Double_t) - 1;
35
da6ac661 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,
5b607771 39 200,210,220,230,240,250,260,270,280,290,300};
da6ac661 40
41 Int_t cent_mult_bin_numbPP = sizeof(cent_mult_binlimitsPP)/sizeof(Double_t) - 1;
42
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};
48
49 Int_t cent_mult_bin_numbPbPb = sizeof(cent_mult_binlimitsPbPb)/sizeof(Double_t) - 1;
ff0805a7 50 //Correlation task
51 AliLeadingV0Correlation *myTask = new AliLeadingV0Correlation(fListName.Data());
52
53 myTask->SetCollidingSystem(fCollisiontype);
54 myTask->SetMCAnalysis(fAnalysisMC);
da6ac661 55 myTask->SetCase(fCase);
da6ac661 56 myTask->SetRemoveAutoCorr(fRemoveAutoCorr);
9bcef488 57 myTask->SetMaxNEventsInPool(fmaxEventsinPool);
5b607771 58 myTask->SetMinNTracksInPool(fminTracksinPool);
da6ac661 59 myTask->SetMinEventsToMix(fMinEventsToMix);
ff0805a7 60 myTask->SetPoolPVzBinLimits(pvzbinnumb,pvzbinlimits);
da6ac661 61 if(fCollisiontype=="PP")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPP,cent_mult_binlimitsPP);
5b607771 62 if(fCollisiontype=="PbPb")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPbPb,cent_mult_binlimitsPbPb);
da6ac661 63 myTask->SetPrimeryVertexCut(fPVzCut);
da6ac661 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);
5b607771 79 myTask->SetTrigLow(fTriglow);
80 myTask->SetTrigHigh(fTrighigh);
ff0805a7 81
82 mgr->AddTask(myTask);
83
84 // Create containers for input/output
85 TString outputFileName = AliAnalysisManager::GetCommonFileName();
9bcef488 86 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
ff0805a7 87 AliAnalysisDataContainer *coutput = mgr->CreateContainer(fListName.Data(), TList::Class(),AliAnalysisManager::kOutputContainer,outputFileName);
88
89 //____________________________________________//
90 mgr->ConnectInput(myTask,0,cinput);
91 mgr->ConnectOutput(myTask,1,coutput);
92
93return myTask;
94}
95
96
5b607771 97