2 AliAnalysisTaskFragFuncBB *AddTaskFragFuncBB(UInt_t filterMask=16, Int_t iPhysicsSelection=1){
\r
3 AddTaskFragFuncBB("jets", "", filterMask, iPhysicsSelection);
\r
6 AliAnalysisTaskFragFuncBB *AddTaskFragFuncBB(UInt_t filterMask=16, Bool_t kUseAODMC=kFALSE, Int_t iPhysicsSelection=1, UInt_t iFlag){
\r
7 AliAnalysisTaskFragFuncBB *ff=0;
\r
9 if(iFlag&(1<<0)) ff = AddTaskFragFuncBB("jets", "jetsAODMC2b_UA104", filterMask, iPhysicsSelection);
\r
10 if(iFlag&(1<<1)) ff = AddTaskFragFuncBB("jetsAOD_UA107", "jetsAODMC2b_UA107", filterMask, iPhysicsSelection);
\r
12 if(iFlag&(1<<2)) ff = AddTaskFragFuncBB("jets", "", filterMask, iPhysicsSelection);
\r
17 // _______________________________________________________________________________________
\r
19 AliAnalysisTaskFragFuncBB *AddTaskFragFuncBB(
\r
20 const char* bRecJets,
\r
21 const char* bGenJets,
\r
23 Int_t iPhysicsSelection)
\r
25 // Creates a fragmentation function task,
\r
26 // configures it and adds it to the analysis manager.
\r
28 // Get the pointer to the existing analysis manager via the static access method.
\r
29 //==============================================================================
\r
30 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
\r
32 ::Error("AddTaskFragFuncBB", "No analysis manager to connect to.");
\r
36 // Check the analysis type using the event handlers connected to the analysis manager.
\r
37 //==============================================================================
\r
38 if (!mgr->GetInputEventHandler()) {
\r
39 ::Error("AddTaskFragFunc", "This task requires an input event handler");
\r
43 TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
\r
44 TString typeRec(bRecJets);
\r
45 TString typeGen(bGenJets);
\r
49 // Create the task and configure it.
\r
50 //===========================================================================
\r
52 AliAnalysisTaskFragFuncBB *fragfunc = new AliAnalysisTaskFragFuncBB(Form("Fragmenation Function %s %s", bRecJets, bGenJets));
\r
54 // pwg4fragfunc->SetAnalysisType(AliAnalysisTaskFragFunc::kAnaMC);
\r
55 // if(iAODanalysis) pwg4fragfunc->SetAODInput(kTRUE);
\r
56 // pwg4fragfunc->SetDebugLevel(11);
\r
58 Printf("Rec Jets %s", bRecJets);
\r
59 Printf("Gen Jets %s", bGenJets);
\r
61 fragfunc->SetBranchGenJets(bGenJets);
\r
62 fragfunc->SetBranchRecJets(bRecJets);
\r
64 fragfunc->SetFilterMask(filterMask);
\r
65 //fragfunc->SetUseGlobalSelection(kTRUE);
\r
68 // Assume all jets are produced already
\r
69 fragfunc->SetAODJetInput(kTRUE);
\r
70 fragfunc->SetAODTrackInput(kTRUE);
\r
71 fragfunc->SetAODMCInput(kTRUE);
\r
74 if(typeRec.Contains("AODMC2b")){// work down from the top AODMC2b -> AODMC2 -> AODMC -> AOD
\r
75 fragfunc->SetTrackTypeRec(AliAnalysisTaskFragFunc::kTrackAODMCChargedAcceptance);
\r
77 else if (typeRec.Contains("AODMC2")){
\r
78 fragfunc->SetTrackTypeRec(AliAnalysisTaskFragFunc::kTrackAODMCCharged);
\r
80 else if (typeRec.Contains("AODMC")){
\r
81 fragfunc->SetTrackTypeRec(AliAnalysisTaskFragFunc::kTrackAODMCAll);
\r
83 else { // catch akk use AOD
\r
84 fragfunc->SetTrackTypeRec(AliAnalysisTaskFragFunc::kTrackAOD);
\r
87 if(typeGen.Contains("AODMC2b")){// work down from the top AODMC2b -> AODMC2 -> AODMC -> AOD
\r
88 fragfunc->SetTrackTypeGen(AliAnalysisTaskFragFunc::kTrackAODMCChargedAcceptance);
\r
90 else if (typeGen.Contains("AODMC2")){
\r
91 fragfunc->SetTrackTypeGen(AliAnalysisTaskFragFunc::kTrackAODMCCharged);
\r
93 else if (typeGen.Contains("AODMC")){
\r
94 fragfunc->SetTrackTypeGen(AliAnalysisTaskFragFunc::kTrackAODMCAll);
\r
96 else if (typeGen.Length()>0){ // catch all use AOD
\r
97 fragfunc->SetTrackTypeGen(AliAnalysisTaskFragFunc::kTrackAOD);
\r
100 fragfunc->SetTrackTypeGen(AliAnalysisTaskFragFunc::kTrackKineCharged);
\r
103 //if(iPhysicsSelection)fragfunc->SelectCollisionCandidates();
\r
105 mgr->AddTask(fragfunc);
\r
107 // Create ONLY the output containers for the data produced by the task.
\r
108 // Get and connect other common input/output containers via the manager as below
\r
109 //==============================================================================
\r
110 AliAnalysisDataContainer *coutput1_FragFunc = mgr->CreateContainer(Form("fracfuncBB_%s_%s",bRecJets,bGenJets), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_fragfuncBB_%s_%s",AliAnalysisManager::GetCommonFileName(),bRecJets,bGenJets));
\r
112 mgr->ConnectInput (fragfunc, 0, mgr->GetCommonInputContainer());
\r
113 mgr->ConnectOutput (fragfunc, 0, mgr->GetCommonOutputContainer());
\r
114 mgr->ConnectOutput (fragfunc, 1, coutput1_FragFunc);
\r