1 AliAnalysisTaskKinkResonance *AddTaskKinkResonanceL1520ESD(Short_t lCollidingSystems=0 /*0 = pp, 1 = AA*/)
3 // Creates, configures and attaches to the train a kink resonance task.
4 // Get the pointer to the existing analysis manager via the static access method.
5 //==============================================================================
6 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
8 ::Error("AddTaskKinkResonanceL1520ESD", "No analysis manager to connect to.");
12 // Check the analysis type using the event handlers connected to the analysis manager.
13 //==============================================================================
14 if (!mgr->GetInputEventHandler()) {
15 ::Error("AddTaskKinkResonanceL1520ESD", "This task requires an input event handler");
18 TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
20 ::Error("AddTaskKinkResonanceL1520ESD", "This task needs an ESD input handler");
23 if (!mgr->GetMCtruthEventHandler()) {
24 ::Error("AddTaskKinkResonanceL1520ESD", "This task needs an MC handler");
28 // Setup the analysis object
30 AliResonanceKink *kinkResonanceObjectLESD=new AliResonanceKink();
31 kinkResonanceObjectLESD->InitOutputHistograms(100,1.4,1.8, 36, -0.9, 0.9, 100, 0.0, 10.0);
32 kinkResonanceObjectLESD->SetPDGCodes(kProton, kKPlus, AliResonanceKink::kLambda1520);
33 kinkResonanceObjectLESD->SetAnalysisType("ESD"); // "ESD" or "MC" or "DATA"
34 kinkResonanceObjectLESD->SetMaxNsigmaToVertex(4.0);
35 kinkResonanceObjectLESD->SetMaxDCAxy(3.0);
36 kinkResonanceObjectLESD->SetMaxDCAzaxis(3.0);
37 kinkResonanceObjectLESD->SetPtTrackCut(0.25);
38 kinkResonanceObjectLESD->SetMinTPCclusters(50);
39 kinkResonanceObjectLESD->SetMaxChi2PerTPCcluster(3.5);
40 kinkResonanceObjectLESD->SetMaxCov0(2.0);
41 kinkResonanceObjectLESD->SetMaxCov2(2.0);
42 kinkResonanceObjectLESD->SetMaxCov5(0.5);
43 kinkResonanceObjectLESD->SetMaxCov9(0.5);
44 kinkResonanceObjectLESD->SetMaxCov14(2.0);
45 kinkResonanceObjectLESD->SetMinKinkRadius(120.);
46 kinkResonanceObjectLESD->SetMaxKinkRadius(220.);
47 kinkResonanceObjectLESD->SetQtLimits(0.05, 0.5);
48 kinkResonanceObjectLESD->SetUpperAbsEtaCut(0.9);
50 // Create and configure the task
51 AliAnalysisTaskKinkResonance *taskresonanceL1520ESD = new AliAnalysisTaskKinkResonance("TaskResL1520ESDKinkPID");
52 taskresonanceL1520ESD->SetAnalysisKinkObject(kinkResonanceObjectLESD);
53 mgr->AddTask(taskresonanceL1520ESD);
55 // Create ONLY the output containers for the data produced by the task.
56 // Get and connect other common input/output containers via the manager as below
57 //==============================================================================
58 TString outputFileName = AliAnalysisManager::GetCommonFileName();
59 outputFileName += ":PWG2KINKResonanceL1520ESD";
60 if (lCollidingSystems) outputFileName += "_AA";
61 else outputFileName += "_PP";
62 if (mgr->GetMCtruthEventHandler()) outputFileName += "_MC";
64 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("KinkResL1520ESD",
66 AliAnalysisManager::kOutputContainer,
69 mgr->ConnectInput(taskresonanceL1520ESD, 0, mgr->GetCommonInputContainer());
70 mgr->ConnectOutput(taskresonanceL1520ESD, 1, coutput1);
71 return taskresonanceL1520ESD;