PWG2/EVCHAR/runProofSPDdNdEta.C

   1 /*************************************************************************
   2 * Macro runProofSPDdNdEta                                                *
   3 * To run dN/dEta reconstruction analysis                                 *
   4 *                                                                        *
   5 * Author:  M. Nicassio (INFN Bari)                                       *
   6 * Contact: Maria.Nicassio@ba.infn.it, Domenico.Elia@ba.infn.it           *
   7 **************************************************************************/
   8
   9 void runProofSPDdNdEta (Int_t MBTrigg, Bool_t kreadmc, Bool_t kppAna,
  10                         Bool_t kallStat, Int_t nEntries, Int_t firstEntry, Char_t* dataSet) {
  11
  12  // Connecting to the PROOF cluster
  13  TProof::Open("alicecaf");
  14
  15  // Enable the needed packages
  16  gProof->UploadPackage("AF-v4-16");
  17  gProof->EnablePackage("AF-v4-16");
  18
  19 // gProof->ShowEnabledPackages();
  20
  21  // Create the analysis manager
  22  mgr = new AliAnalysisManager("SPD dN/dEta Analysis");
  23
  24  // Create, add task
  25  gProof->Load("AliAnalysisTaskSPDdNdEta.cxx++g");
  26  task = new AliAnalysisTaskSPDdNdEta();
  27
  28  // Add ESD handler
  29  AliESDInputHandler* esdH = new AliESDInputHandler;
  30  mgr->SetInputEventHandler(esdH);
  31  // Add MC handler
  32  if (kreadmc) {
  33    AliMCEventHandler*   mcH = new AliMCEventHandler();
  34    mgr->SetMCtruthEventHandler(mcH);
  35  }
  36  task->SetReadMC(kreadmc);
  37
  38  task->SetppAnalysis(kppAna);
  39
  40  task->SetTrigger(MBTrigg);
  41
  42  // Attach input
  43  cInput = mgr->CreateContainer("cInput", TChain::Class(), AliAnalysisManager::kInputContainer);
  44  mgr->ConnectInput(task, 0, cInput);
  45
  46  // Attach output
  47  if (kreadmc) cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaMC.root");
  48
  49  else cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaData.root");
  50
  51  mgr->ConnectOutput(task, 0, cOutput);
  52
  53  // Enable debug printouts
  54  mgr->SetDebugLevel(2);
  55
  56  // Run analysis
  57  mgr->InitAnalysis();
  58  mgr->PrintStatus();
  59
  60  if (kallStat)
  61    mgr->StartAnalysis("proof", dataSet);
  62  else
  63    mgr->StartAnalysis("proof", dataSet, nEntries, firstEntry);
  64
  65 }