1 //DEFINITION OF A FEW CONSTANTS
2 const Int_t chargeFrCorr = 20;
3 //----------------------------------------------------
5 AliAnalysisTaskSEDStarJets *AddTaskDStarJets(Bool_t theMCon=kTRUE)
8 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
10 ::Error("AddTaskDStarJets2", "No analysis manager to connect to.");
14 TFile* filecuts=new TFile("DStartoKpipiCuts.root");
15 if(!filecuts->IsOpen()){
16 cout<<"Input file not found: exit"<<endl;
20 AliRDHFCutsDStartoKpipi* RDHFDStartoKpipi=new AliRDHFCutsDStartoKpipi();
21 RDHFDStartoKpipi = (AliRDHFCutsDStartoKpipi*)filecuts->Get("DStartoKpipiCuts");
22 RDHFDStartoKpipi->SetName("DStartoKpipiCuts");
24 // mm let's see if everything is ok
25 if(!RDHFDStartoKpipi){
26 cout<<"Specific AliRDHFCuts not found"<<endl;
31 printf("CREATE TASK\n");
33 AliAnalysisTaskSEDStarJets *task = new AliAnalysisTaskSEDStarJets("AliAnalysisTaskSEDStarJets",RDHFDStartoKpipi);
35 task->SetChargeFractionCorrection(chargeFrCorr);
37 // Create and connect containers for input/output
38 TString outputfile = AliAnalysisManager::GetCommonFileName();
39 outputfile += ":PWG3_D2H_DStarJet";
41 // ------ input data ------
42 AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
44 // ----- output data -----
46 // output TH1I for event counting
47 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("charmJetCorr", TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
48 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("cuts",AliRDHFCutsDStartoKpipi::Class(),AliAnalysisManager::kOutputContainer, outputfile.Data()); //cuts
51 mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
52 mgr->ConnectOutput(task,1,coutput1);
53 mgr->ConnectOutput(task,2,coutput2);