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 (Bool_t kpythia, Int_t MBTrigg, Bool_t kreadmc, Bool_t kppAna,
  10                         Bool_t kallStat, Int_t nEntries, Int_t firstEntry, Char_t* dataSet) {
  11
  12  gEnv->SetValue("XSec.GSI.DelegProxy","2");
  13  // Connecting to the PROOF cluster
  14  TProof::Open("alicecaf");
  15
  16  gProof->ClearPackages();
  17  // Enable the needed packages
  18  gProof->UploadPackage("AF-v4-17");
  19  gProof->EnablePackage("AF-v4-17");
  20
  21 // gProof->ShowEnabledPackages();
  22
  23  // Create the analysis manager
  24  mgr = new AliAnalysisManager("SPD dN/dEta Analysis");
  25
  26  // Create, add task
  27  gProof->Load("AliAnalysisTaskSPDdNdEta.cxx++g");
  28  task = new AliAnalysisTaskSPDdNdEta();
  29  mgr->AddTask(task);
  30
  31  // Add ESD handler
  32  AliESDInputHandler* esdH = new AliESDInputHandler;
  33  mgr->SetInputEventHandler(esdH);
  34  // Add MC handler
  35  if (kreadmc) {
  36    AliMCEventHandler*   mcH = new AliMCEventHandler();
  37    mgr->SetMCtruthEventHandler(mcH);
  38  }
  39  task->SetReadMC(kreadmc);
  40
  41  task->SetTrigger(MBTrigg);
  42
  43  task->SetEvtGen(kpythia);
  44
  45  // Attach input
  46  cInput = mgr->GetCommonInputContainer();
  47  mgr->ConnectInput(task, 0, cInput);
  48
  49  // Attach output
  50  if (kreadmc) cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaMC.root");
  51
  52  else cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaData.root");
  53
  54  mgr->ConnectOutput(task, 0, cOutput);
  55
  56  // Enable debug printouts
  57  mgr->SetDebugLevel(2);
  58
  59  // Run analysis
  60  mgr->InitAnalysis();
  61  mgr->PrintStatus();
  62
  63  if (kallStat)
  64    mgr->StartAnalysis("proof", dataSet);
  65  else
  66    mgr->StartAnalysis("proof", dataSet, nEntries, firstEntry);
  67
  68 }