1 AliAnalysisTaskFlavourJetCorrelations *AddTaskFlavourJetCorrelations(
2 AliAnalysisTaskFlavourJetCorrelations::ECandidateType cand = AliAnalysisTaskFlavourJetCorrelations::kDstartoKpipi,
3 TString filename = "DStartoKpipiCuts.root",
4 Bool_t theMCon = kFALSE,
5 Bool_t reco = kTRUE /*must be true if theMCon is false*/,
6 TString jetArrname = "",
8 Bool_t triggerOnLeadingJet = kFALSE,
9 Int_t leadingHadType = 0 /*0 = charged, 1 = neutral, 2 = both*/,
12 const char *cutType = "TPC",
13 Double_t percjetareacut = 1.)
15 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
17 ::Error("AddTaskFlavourJetCorrelations::AddTaskFlavourJetCorrelations", "No analysis manager to connect to.");
21 Bool_t useStdC = kFALSE;
22 TFile* filecuts=TFile::Open(filename);
23 if (!filecuts || (filecuts && !filecuts->IsOpen())) {
24 cout<<"Input file not found: use std cuts"<<endl;
28 AliRDHFCuts *analysiscuts = 0x0;
32 analysiscuts = new AliRDHFCutsD0toKpi();
33 analysiscuts->SetStandardCutsPP2010();
35 analysiscuts = (AliRDHFCutsD0toKpi*)filecuts->Get("D0toKpiCuts");
39 analysiscuts = new AliRDHFCutsDStartoKpipi();
40 analysiscuts->SetStandardCutsPP2010();
42 analysiscuts = (AliRDHFCutsDStartoKpipi*)filecuts->Get("DStartoKpipiCuts");
43 analysiscuts->SetName("DStartoKpipiCuts");
47 if (!analysiscuts) { // mm let's see if everything is ok
48 AliFatal("Specific AliRDHFCuts not found");
52 printf("CREATE TASK\n"); //CREATE THE TASK
55 AliAnalysisTaskFlavourJetCorrelations *task = new AliAnalysisTaskFlavourJetCorrelations("AnaTaskFlavourJetCorrelations",
57 task->SetJetsName(jetArrname);
59 task->SetUseReco(reco);
60 task->SetTriggerOnLeadingJet(triggerOnLeadingJet);
61 task->SetJetAcceptanceType(cutType);
67 if(reco) suffix+="rec";
70 // Create and connect containers for input/output
71 TString outputfile = AliAnalysisManager::GetCommonFileName();
72 outputfile += ":PWG3_D2H_DEmcalJet";
75 TString candname="DStar";
76 if(cand==0) candname="D0";
78 TString nameContainer1="hCor";
79 TString nameContainer2="cutsJ";
81 nameContainer1 += candname;
82 nameContainer2 += candname;
84 nameContainer1 += suffix;
85 nameContainer2 += suffix;
86 // ------ input data ------
87 AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
89 // ----- output data -----
90 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(nameContainer1, TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
91 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer(nameContainer2, AliRDHFCuts::Class(),AliAnalysisManager::kOutputContainer, outputfile.Data());
93 mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
94 mgr->ConnectOutput(task,1,coutput1);
95 mgr->ConnectOutput(task,2,coutput2);
98 Printf("Input and Output connected to the manager");