SPD multiplicity analysis (Domenico, Mariella)

[u/mrichter/AliRoot.git] / PWG2 / EVCHAR / runProofSPDdNdEta.C

void runProofSPDdNdEta (Int_t MBTrigg, Bool_t kreadmc, Bool_t kallStat, Int_t nEntries, Int_t firstEntry, Char_t* dataSet) {

 // Connecting to the PROOF cluster
 TProof::Open("alicecaf");

 // Enable the needed packages
 gProof->UploadPackage("AF-v4-16");
 gProof->EnablePackage("AF-v4-16");

// gProof->ShowEnabledPackages(); // show a list of enabled packages

 // Create the analysis manager
 mgr = new AliAnalysisManager("SPD dN/dEta Analysis");

 // Create, add task
 gProof->Load("AliAnalysisTaskSPDdNdEta.cxx++g");
 task = new AliAnalysisTaskSPDdNdEta();

 // Add ESD handler
 AliESDInputHandler* esdH = new AliESDInputHandler;
 mgr->SetInputEventHandler(esdH);
 // Add MC handler
 if (kreadmc) {
 AliMCEventHandler*   mcH = new AliMCEventHandler();
 mgr->SetMCtruthEventHandler(mcH);
 task->SetReadMC(kTRUE);  
 } else task->SetReadMC(kFALSE);

 task->SetTrigger(MBTrigg);

 // Attach input
 cInput = mgr->CreateContainer("cInput", TChain::Class(), AliAnalysisManager::kInputContainer);
 mgr->ConnectInput(task, 0, cInput);

 // Attach output
 if (kreadmc) cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaMC.root");

 else cOutput= mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer, "SPDdNdEtaData.root");

 mgr->ConnectOutput(task, 0, cOutput);

 // Enable debug printouts
 mgr->SetDebugLevel(2);

 // Run analysis
 mgr->InitAnalysis();
 mgr->PrintStatus();
 
 if (kallStat)
   mgr->StartAnalysis("proof", dataSet);
 else 
   mgr->StartAnalysis("proof", dataSet, nEntries, firstEntry);

}