1 AliAnalysisTask *AddTaskFlowEPCascade(Float_t centrMin=0, Float_t centrMax=100,
\r
2 TString fileName, TString suffixName="",
\r
3 Double_t vtxCut = 10.,
\r
4 Double_t etaCut = 0.9,
\r
6 Double_t mSigma = 0.0024)
\r
9 //-E-V-E-N-T- -c-u-t-s-----------------------------------------------------
\r
10 AliFlowEventCuts* cutsEvent
\r
11 = new AliFlowEventCuts(Form("event_cuts_%s", suffixName.Data()));
\r
12 cutsEvent->SetCentralityPercentileRange(centrMin,centrMax);
\r
13 cutsEvent->SetCentralityPercentileMethod(AliFlowEventCuts::kV0);
\r
14 cutsEvent->SetRefMultMethod(AliFlowEventCuts::kV0);
\r
15 cutsEvent->SetNContributorsRange(2);
\r
16 cutsEvent->SetPrimaryVertexZrange( -10., 10. );
\r
17 cutsEvent->SetCutSPDvertexerAnomaly();
\r
18 cutsEvent->SetCutZDCtiming();
\r
19 cutsEvent->SetCutTPCmultiplicityOutliers();
\r
21 AliFlowTrackCuts* cutsRP
\r
22 = (AliFlowTrackCuts*) AliFlowTrackCuts::GetStandardGlobalTrackCuts2010();
\r
23 cutsRP->SetPtRange(0.15, 10.); //added to extend pt range
\r
26 AliFlowTrackCuts * cutsDaughter
\r
27 = new AliFlowTrackCuts(Form("daughter_cuts_%s",suffixName.Data()));
\r
28 cutsDaughter->SetPtRange(0.15,10.0);
\r
29 cutsDaughter->SetEtaRange(-etaCut,etaCut);
\r
30 cutsDaughter->SetMinNClustersTPC(nTPCcls);
\r
31 cutsDaughter->SetMaxChi2PerClusterTPC(4.0);
\r
32 cutsDaughter->SetAODfilterBit(1); //TPC track only
\r
33 cutsDaughter->SetRequireITSRefit(kFALSE);
\r
34 cutsDaughter->SetRequireTPCRefit(kTRUE);
\r
35 cutsDaughter->SetMinNClustersITS(0);
\r
36 cutsDaughter->SetAcceptKinkDaughters(kFALSE);
\r
38 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
\r
39 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
\r
41 //____________________________________________//
\r
42 double pdg = 1.32171; double hwid = mSigma;
\r
43 double sep = 6*hwid; double XiStart = pdg-sep;
\r
44 double XiBands[3][2];
\r
45 XiBands[0][0]=(XiStart+6*hwid*0)-2*hwid;
\r
46 XiBands[0][1]=(XiStart+6*hwid*0)+2*hwid;
\r
47 XiBands[1][0]=(XiStart+6*hwid*1)-2*hwid;
\r
48 XiBands[1][1]=(XiStart+6*hwid*1)+2*hwid;
\r
49 XiBands[2][0]=(XiStart+6*hwid*2)-2*hwid;
\r
50 XiBands[2][1]=(XiStart+6*hwid*2)+2*hwid;
\r
52 double OmegaPDG = 1.67245;
\r
53 double Omegasg = mSigma; // 1 sigma?
\r
54 double OmegaStart = OmegaPDG-6*Omegasg;
\r
55 double OmegaBands[3][2];
\r
56 OmegaBands[0][0]=(OmegaStart+6*Omegasg*0)-2*Omegasg;
\r
57 OmegaBands[0][1]=(OmegaStart+6*Omegasg*0)+2*Omegasg;
\r
58 OmegaBands[1][0]=(OmegaStart+6*Omegasg*1)-2*Omegasg;
\r
59 OmegaBands[1][1]=(OmegaStart+6*Omegasg*1)+2*Omegasg;
\r
60 OmegaBands[2][0]=(OmegaStart+6*Omegasg*2)-2*Omegasg;
\r
61 OmegaBands[2][1]=(OmegaStart+6*Omegasg*2)+2*Omegasg;
\r
63 AliAnalysisTaskFlowEPCascade * task
\r
64 = new AliAnalysisTaskFlowEPCascade(Form("v2EP_%s",
\r
65 suffixName.Data()),
\r
67 XiBands, OmegaBands);
\r
68 // task->SetFlowEventCuts(cutsEvent);
\r
69 //task->SetFlowTrackCuts(cutsRP);
\r
70 task->SetVertexCut(vtxCut);
\r
71 task->SetFlowDauTrackCuts(cutsDaughter);
\r
72 task->SelectCollisionCandidates(AliVEvent::kMB);
\r
75 AliAnalysisDataContainer *coutput1 =
\r
76 mgr->CreateContainer(Form("hist_%s", suffixName.Data()),
\r
78 AliAnalysisManager::kOutputContainer,
\r
79 Form("%s.root",fileName.Data()));
\r
80 // Connect input/output
\r
81 mgr->ConnectInput(task, 0, cinput);
\r
82 //mgr->ConnectInput(task, 1, cinput1);
\r
83 mgr->ConnectOutput(task, 1, coutput1);
\r