AliAnalysisTask *AddTask_jbook_JPsiQA(TString prod="", Bool_t gridconf=kFALSE){
//get the current analysis manager
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
Error("AddTask_jbook_JPsiQA", "No analysis manager found.");
return 0;
}
//Do we have an MC handler?
Bool_t hasMC = kFALSE;
TString list = gSystem->Getenv("LIST");
if( list.IsNull()) list=prod;
if( list.Contains("LHC10h") || list.Contains("LHC11h") ) hasMC=kFALSE;
if( list.Contains("LHC11a10") || list.Contains("LHC12a17") ) hasMC=kTRUE;
//Do we have an AOD handler?
Bool_t isAOD=(mgr->GetInputEventHandler()->IsA()==AliAODInputHandler::Class() ? kTRUE : kFALSE);
// set AOD debug levels
if(isAOD) {
mgr->AddClassDebug("AliAODTrack", AliLog::kFatal);
mgr->AddClassDebug("AliAODpidUtil", AliLog::kInfo);
}
//set config file name
TString configFile("");
printf("%s \n",gSystem->pwd());
TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
if (!trainRoot.IsNull())
configFile="$TRAIN_ROOT/jbook_jpsi/ConfigJpsiQA_jb_PbPb.C"; // gsi config
else if(!gSystem->Exec("alien_cp alien:///alice/cern.ch/user/j/jbook/PWGDQ/dielectron/macrosJPSI/ConfigJpsiQA_jb_PbPb.C ."))
configFile=Form("%s/ConfigJpsiQA_jb_PbPb.C",gSystem->pwd()); // alien config
else
configFile="$ALICE_ROOT/PWGDQ/dielectron/macrosJPSI/ConfigJpsiQA_jb_PbPb.C"; // aliroot config
if(!gridconf)
configFile="$ALICE_ROOT/PWGDQ/dielectron/macrosJPSI/ConfigJpsiQA_jb_PbPb.C"; // aliroot config
//create task and add it to the manager
AliAnalysisTaskMultiDielectron *task=new AliAnalysisTaskMultiDielectron("MultiDieData");
task->SetTriggerMask(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
if (!hasMC) task->UsePhysicsSelection();
mgr->AddTask(task);
//load dielectron configuration file
TString checkconfig="ConfigJpsiQA_jb_PbPb";
if (!gROOT->GetListOfGlobalFunctions()->FindObject(checkconfig.Data()))
gROOT->LoadMacro(configFile.Data());
//add dielectron analysis with different cuts to the task
for (Int_t i=0; i<nDie; ++i){ //nDie defined in config file
AliDielectron *jpsi=ConfigJpsiQA_jb_PbPb(i,list);
if (jpsi ) task->AddDielectron(jpsi);
if (jpsi ) printf("add: %s\n",jpsi->GetName());
}
//Add event filter
AliDielectronEventCuts *eventCuts=new AliDielectronEventCuts("eventCuts","Vertex Track && |vtxZ|<10 && ncontrib>0");
if (isAOD) eventCuts->SetVertexType(AliDielectronEventCuts::kVtxAny);
eventCuts->SetRequireVertex();
eventCuts->SetMinVtxContributors(1);
eventCuts->SetVertexZ(-10.,10.);
eventCuts->SetCentralityRange(0.0,80.0);
task->SetEventFilter(eventCuts);
//create output container
AliAnalysisDataContainer *coutput1 =
mgr->CreateContainer("jbookQA_tree",
TTree::Class(),
AliAnalysisManager::kExchangeContainer,
"jbook_default");
AliAnalysisDataContainer *cOutputHist1 =
mgr->CreateContainer("jbookQA",
TList::Class(),
AliAnalysisManager::kOutputContainer,
"jbookQA.root");
AliAnalysisDataContainer *cOutputHist2 =
mgr->CreateContainer("jbookQA_CF",
TList::Class(),
AliAnalysisManager::kOutputContainer,//AliAnalysisManager::kExchangeContainer,
"jbookQA.root");
AliAnalysisDataContainer *cOutputHist3 =
mgr->CreateContainer("jbookQA_EventStat",
TH1D::Class(),
AliAnalysisManager::kOutputContainer,
"jbookQA.root");
mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
mgr->ConnectOutput(task, 0, coutput1 );
mgr->ConnectOutput(task, 1, cOutputHist1);
mgr->ConnectOutput(task, 2, cOutputHist2);
mgr->ConnectOutput(task, 3, cOutputHist3);
return task;
}
AddTask_jbook_JPsiQA.C:10 AddTask_jbook_JPsiQA.C:11 AddTask_jbook_JPsiQA.C:12 AddTask_jbook_JPsiQA.C:13 AddTask_jbook_JPsiQA.C:14 AddTask_jbook_JPsiQA.C:15 AddTask_jbook_JPsiQA.C:16 AddTask_jbook_JPsiQA.C:17 AddTask_jbook_JPsiQA.C:18 AddTask_jbook_JPsiQA.C:19 AddTask_jbook_JPsiQA.C:20 AddTask_jbook_JPsiQA.C:21 AddTask_jbook_JPsiQA.C:22 AddTask_jbook_JPsiQA.C:23 AddTask_jbook_JPsiQA.C:24 AddTask_jbook_JPsiQA.C:25 AddTask_jbook_JPsiQA.C:26 AddTask_jbook_JPsiQA.C:27 AddTask_jbook_JPsiQA.C:28 AddTask_jbook_JPsiQA.C:29 AddTask_jbook_JPsiQA.C:30 AddTask_jbook_JPsiQA.C:31 AddTask_jbook_JPsiQA.C:32 AddTask_jbook_JPsiQA.C:33 AddTask_jbook_JPsiQA.C:34 AddTask_jbook_JPsiQA.C:35 AddTask_jbook_JPsiQA.C:36 AddTask_jbook_JPsiQA.C:37 AddTask_jbook_JPsiQA.C:38 AddTask_jbook_JPsiQA.C:39 AddTask_jbook_JPsiQA.C:40 AddTask_jbook_JPsiQA.C:41 AddTask_jbook_JPsiQA.C:42 AddTask_jbook_JPsiQA.C:43 AddTask_jbook_JPsiQA.C:44 AddTask_jbook_JPsiQA.C:45 AddTask_jbook_JPsiQA.C:46 AddTask_jbook_JPsiQA.C:47 AddTask_jbook_JPsiQA.C:48 AddTask_jbook_JPsiQA.C:49 AddTask_jbook_JPsiQA.C:50 AddTask_jbook_JPsiQA.C:51 AddTask_jbook_JPsiQA.C:52 AddTask_jbook_JPsiQA.C:53 AddTask_jbook_JPsiQA.C:54 AddTask_jbook_JPsiQA.C:55 AddTask_jbook_JPsiQA.C:56 AddTask_jbook_JPsiQA.C:57 AddTask_jbook_JPsiQA.C:58 AddTask_jbook_JPsiQA.C:59 AddTask_jbook_JPsiQA.C:60 AddTask_jbook_JPsiQA.C:61 AddTask_jbook_JPsiQA.C:62 AddTask_jbook_JPsiQA.C:63 AddTask_jbook_JPsiQA.C:64 AddTask_jbook_JPsiQA.C:65 AddTask_jbook_JPsiQA.C:66 AddTask_jbook_JPsiQA.C:67 AddTask_jbook_JPsiQA.C:68 AddTask_jbook_JPsiQA.C:69 AddTask_jbook_JPsiQA.C:70 AddTask_jbook_JPsiQA.C:71 AddTask_jbook_JPsiQA.C:72 AddTask_jbook_JPsiQA.C:73 AddTask_jbook_JPsiQA.C:74 AddTask_jbook_JPsiQA.C:75 AddTask_jbook_JPsiQA.C:76 AddTask_jbook_JPsiQA.C:77 AddTask_jbook_JPsiQA.C:78 AddTask_jbook_JPsiQA.C:79 AddTask_jbook_JPsiQA.C:80 AddTask_jbook_JPsiQA.C:81 AddTask_jbook_JPsiQA.C:82 AddTask_jbook_JPsiQA.C:83 AddTask_jbook_JPsiQA.C:84 AddTask_jbook_JPsiQA.C:85 AddTask_jbook_JPsiQA.C:86 AddTask_jbook_JPsiQA.C:87 AddTask_jbook_JPsiQA.C:88 AddTask_jbook_JPsiQA.C:89 AddTask_jbook_JPsiQA.C:90 AddTask_jbook_JPsiQA.C:91 AddTask_jbook_JPsiQA.C:92 AddTask_jbook_JPsiQA.C:93 AddTask_jbook_JPsiQA.C:94 AddTask_jbook_JPsiQA.C:95 AddTask_jbook_JPsiQA.C:96 AddTask_jbook_JPsiQA.C:97 AddTask_jbook_JPsiQA.C:98 AddTask_jbook_JPsiQA.C:99 AddTask_jbook_JPsiQA.C:100 AddTask_jbook_JPsiQA.C:101