PWG0/genLevelSimulation/runTaskProof.C

   1 void runTaskProof(const char * dataset, const char * datasetpath="/COMMON/COMMON/", const char * outdir="") {
   2   gEnv->SetValue("XSec.GSI.DelegProxy","2");
   3   TProof::Open("alice-caf","workers=50");
   4
   5   //  gSystem->AddIncludePath("-I${ALICE_ROOT}/include/ -I${ALICE_ROOT}/PWG0/ -I${ALICE_ROOT}/PWG0/dNdEta/");
   6   gSystem->AddIncludePath("-I${ALICE_ROOT}/include/");
   7   gProof->UploadPackage("$ALICE_ROOT/STEERBase");
   8   gProof->EnablePackage("$ALICE_ROOT/STEERBase");
   9   gProof->UploadPackage("$ALICE_ROOT/ESD");
  10   gProof->EnablePackage("$ALICE_ROOT/ESD");
  11   gProof->UploadPackage("$ALICE_ROOT/AOD");
  12   gProof->EnablePackage("$ALICE_ROOT/AOD");
  13   gProof->UploadPackage("$ALICE_ROOT/ANALYSIS");
  14   gProof->EnablePackage("$ALICE_ROOT/ANALYSIS");
  15   gProof->UploadPackage("$ALICE_ROOT/ANALYSISalice");
  16   gProof->EnablePackage("$ALICE_ROOT/ANALYSISalice");
  17
  18
  19     // Make the analysis manager
  20   AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
  21   //AliVEventHandler *esdH = new AliESDInputHandler();
  22   AliESDInputHandler* esdH = new AliESDInputHandler();
  23   //esdH->SetActiveBranches("ESDfriend");
  24   mgr->SetInputEventHandler(esdH);
  25
  26   AliMCEventHandler *mc = new AliMCEventHandler();
  27   mc->SetReadTR(kFALSE);
  28   mgr->SetMCtruthEventHandler(mc);
  29
  30   // assign simple task
  31   gProof->Load(gSystem->ExpandPathName("$(ALICE_ROOT)/PWG0/genLevelSimulation/AliAnalysisTaskdNdetaMC.cxx+g"));
  32   AliAnalysisTask *task = new AliAnalysisTaskdNdetaMC("TaskdNdeta");
  33   mgr->AddTask(task);
  34
  35   char outfname[1000];
  36   sprintf(outfname, "dndeta_%s.root",dataset);
  37
  38   // Create containers for input/output
  39   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
  40
  41   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("clist1",TList::Class(),
  42                                                             AliAnalysisManager::kOutputContainer,outfname);
  43
  44   mgr->ConnectInput(task,0,cinput1);
  45   mgr->ConnectOutput(task,1,coutput1);
  46
  47   if (!mgr->InitAnalysis()) return;
  48
  49   mgr->PrintStatus();
  50   mgr->StartAnalysis("proof",Form("%s%s#TE",datasetpath,dataset),5000);
  51   //  mgr->StartAnalysis("proof","/COMMON/COMMON/LHC09b14_7TeV_0.5T_Phojet#esdTree");
  52
  53   if (!strcmp(outdir,"")){
  54     gSystem->Exec(Form("mv %s %s", outfname, outdir));
  55   }
  56
  57 }