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

/*************************************************************************
* Macro runProofSPDdNdEta                                                *
* To run dN/dEta reconstruction analysis                                 *
*                                                                        *
* Author:  M. Nicassio (INFN Bari)                                       *
* Contact: Maria.Nicassio@ba.infn.it, Domenico.Elia@ba.infn.it           *
**************************************************************************/

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

 gEnv->SetValue("XSec.GSI.DelegProxy","2");
 // Connecting to the PROOF cluster
 TProof::Open("alicecaf");

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

// gProof->ShowEnabledPackages(); 

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

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

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

 task->SetTrigger(MBTrigg);

 task->SetEvtGen(kpythia);
 
 // Attach input
 cInput = mgr->GetCommonInputContainer();
 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);

}
Commit	Line	Data
	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	}