[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / PiKaPr / TestAOD / runAODProof.C

void runAODProof(Int_t c=1, const char * proofMode = "full")
{ //1 data AOD049
  //2 MC AOD048
  //3 data AOD086
  //4 MC AOD090
  
  //  gEnv->SetValue("XSec.GSI.DelegProxy", "2");
  
  gSystem->Load("libTree.so");
  gSystem->Load("libGeom.so");
   gSystem->Load("libVMC.so");
   gSystem->Load("libPhysics.so");
   gSystem->Load("libSTEERBase.so");
   gSystem->Load("libESD.so");
   gSystem->Load("libAOD.so");
   gSystem->Load("libANALYSIS.so");
   gSystem->Load("libOADB.so");
   gSystem->Load("libANALYSISalice.so");
   gSystem->AddIncludePath("-I$ALICE_ROOT/include");

   AliAnalysisAlien * handler = new AliAnalysisAlien("test");
   handler->SetOverwriteMode();
   handler->SetRunMode(proofMode);
   handler->SetProofReset(0);
   //handler->SetROOTVersion("v5-33-02a");
   //handler->SetAliROOTVersion("v5-03-11-AN");
   handler->SetAliROOTVersion("v5-03-17-AN");
   
   //handler->SetNproofWorkers(28);
   //handler->SetNproofWorkersPerSlave(1);
   handler->SetProofCluster(Form("%s@alice-caf.cern.ch", gSystem->Getenv("CAFUSER")));
   //handler->SetProofCluster(Form("%s@skaf.saske.sk",gSystem->Getenv("CAFUSER")));
   // Set handler for Real DATA:
   if (c == 1){
     //handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree");
     handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree|/alice/data/LHC10h_000138662_p2_AOD049#aodTree|/alice/data/LHC10h_000138666_p2_AOD049#aodTree|/alice/data/LHC10h_000138795_p2_AOD049#aodTree|/alice/data/LHC10h_000139107_p2_AOD049#aodTree|/alice/data/LHC10h_000139110_p2_AOD049#aodTree");
   }
   if (c == 2){
     handler->SetProofDataSet("/alice/sim/LHC11a10a_000138653_AOD048#aodTree|/alice/sim/LHC11a10a_000138662_AOD048#aodTree|/alice/sim/LHC11a10a_000138666_AOD048#aodTree|/alice/sim/LHC11a10a_000138795_AOD048#aodTree|/alice/sim/LHC11a10a_000139107_AOD048#aodTree|/alice/sim/LHC11a10a_000139110_AOD048#aodTree");      
   }
   if (c == 3){
     handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138662_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138666_AOD086_p2#aodTree|/default/lmilano/LHC10h_000139107_AOD086_p2#aodTree");      
   }
   if (c == 4){
     handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138662_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138666_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_139107_AOD090#aodTree");      
   }
   
   gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
   gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g");
   gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
   gROOT->LoadMacro("AliSpectraAODPID.cxx+g");
   gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");

   handler->SetAliRootMode("default");
   handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODPID.cxx AliSpectraAODPID.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
   handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODPID.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
   //   handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
   //  handler->SetAliRootMode("");
   handler->SetClearPackages();
   
   AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
   mgr->SetGridHandler(handler);
   AliAODInputHandler* aodH = new AliAODInputHandler();
   mgr->SetInputEventHandler(aodH);
   
   // Add PID task
   gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
   Bool_t isMC = kFALSE;
   if (c == 2 || c == 4) isMC = kTRUE;   
   AliAnalysisTask * taskPID = AddTaskPIDResponse(isMC);
   mgr->AddTask(taskPID);
   
   //LOOP OVER SELECTION
   Double_t CentCutMin[4]={0,0,20,20};
   Double_t CentCutMax[4]={100,5,40,40};
   Double_t QvecCutMin[4]={0,0,0,3};
   Double_t QvecCutMax[4]={100,100,2,100};
   // Double_t CentCutMin[4]={5,10,20,30};
   // Double_t CentCutMax[4]={10,20,30,40};
   // Double_t QvecCutMin[4]={0,0,0,0};
   // Double_t QvecCutMax[4]={100,100,100,100};
   using namespace AliSpectraNameSpace;
   AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF); 
   pid->SetNSigmaCut(5.);
   
   for(Int_t iCut=1;iCut<2;iCut++){
     AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
     mgr->AddTask(task);
     //physics selection
     task->SelectCollisionCandidates();     
     // set pid object
     task->SetPID(pid);
     // Set the cuts
     AliSpectraAODEventCuts * vcuts = new AliSpectraAODEventCuts("Event Cuts");
     AliSpectraAODTrackCuts  * tcuts = new AliSpectraAODTrackCuts("Track Cuts");
     //if(c==1 || c==2)tcuts->SetTrackType(5); //AOD 046 & 047. Standard Cut with loose DCA
     if(c==1 || c==2)tcuts->SetTrackType(6); //AOD 046 & 047. Standard Cut with tight DCA
     if(c==3 || c==4)tcuts->SetTrackType(10); //AOD 086 & 090. Standard Raa cut
     tcuts->SetEta(.8);
     //tcuts->SetDCA(.1);
     tcuts->SetPt(5);
     tcuts->SetY(.5);
     tcuts->SetPtTOFMatching(0.6);   
     tcuts->SetQvecMin(QvecCutMin[iCut]);   
     tcuts->SetQvecMax(QvecCutMax[iCut]);    
     vcuts->SetCentralityCutMax(CentCutMax[iCut]);  
     vcuts->SetCentralityCutMin(CentCutMin[iCut]);
     task->SetEventCuts(vcuts);
     task->SetTrackCuts(tcuts);
     vcuts->PrintCuts();
     tcuts->PrintCuts();
     
     // check for MC or real data
     if (c == 2 || c == 4)
       {
	 task->SetIsMC(kTRUE);
	 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	 AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
       }
     if (c == 1 || c==3)
       {
	 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	 AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
								     Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	 
       }
     mgr->ConnectInput(task, 0, cinput);
     mgr->ConnectOutput(task, 1, coutputpt1);
     mgr->ConnectOutput(task, 2, coutputpt2);
     mgr->ConnectOutput(task, 3, coutputpt3);
     mgr->ConnectOutput(task, 4, coutputpt4);
   }
   mgr->SetDebugLevel(2);
   
   if (!mgr->InitAnalysis()) return;
   mgr->PrintStatus();
   mgr->StartAnalysis("proof");
}
Commit	Line	Data
97b01bb0	1	void runAODProof(Int_t c=1, const char * proofMode = "full")
ed06106e	2	{ //1 data AOD049
	3	//2 MC AOD048
	4	//3 data AOD086
	5	//4 MC AOD090
	6
0e268798	7	// gEnv->SetValue("XSec.GSI.DelegProxy", "2");
ed06106e	8
ed06106e	9	gSystem->Load("libTree.so");
0e268798	10	gSystem->Load("libGeom.so");
c88234ad	11	gSystem->Load("libVMC.so");
	12	gSystem->Load("libPhysics.so");
	13	gSystem->Load("libSTEERBase.so");
	14	gSystem->Load("libESD.so");
	15	gSystem->Load("libAOD.so");
	16	gSystem->Load("libANALYSIS.so");
	17	gSystem->Load("libOADB.so");
	18	gSystem->Load("libANALYSISalice.so");
	19	gSystem->AddIncludePath("-I$ALICE_ROOT/include");
	20
c88234ad	21	AliAnalysisAlien * handler = new AliAnalysisAlien("test");
	22	handler->SetOverwriteMode();
	23	handler->SetRunMode(proofMode);
	24	handler->SetProofReset(0);
823864bf	25	//handler->SetROOTVersion("v5-33-02a");
823864bf	26	//handler->SetAliROOTVersion("v5-03-11-AN");
97b01bb0	27	handler->SetAliROOTVersion("v5-03-17-AN");
6522a8cb	28
742d45e0	29	//handler->SetNproofWorkers(28);
742d45e0	30	//handler->SetNproofWorkersPerSlave(1);
c88234ad	31	handler->SetProofCluster(Form("%s@alice-caf.cern.ch", gSystem->Getenv("CAFUSER")));
742d45e0	32	//handler->SetProofCluster(Form("%s@skaf.saske.sk",gSystem->Getenv("CAFUSER")));
823864bf	33	// Set handler for Real DATA:
ed06106e	34	if (c == 1){
742d45e0	35	//handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree");
742d45e0	36	handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree\|/alice/data/LHC10h_000138662_p2_AOD049#aodTree\|/alice/data/LHC10h_000138666_p2_AOD049#aodTree\|/alice/data/LHC10h_000138795_p2_AOD049#aodTree\|/alice/data/LHC10h_000139107_p2_AOD049#aodTree\|/alice/data/LHC10h_000139110_p2_AOD049#aodTree");
ed06106e	37	}
	38	if (c == 2){
	39	handler->SetProofDataSet("/alice/sim/LHC11a10a_000138653_AOD048#aodTree\|/alice/sim/LHC11a10a_000138662_AOD048#aodTree\|/alice/sim/LHC11a10a_000138666_AOD048#aodTree\|/alice/sim/LHC11a10a_000138795_AOD048#aodTree\|/alice/sim/LHC11a10a_000139107_AOD048#aodTree\|/alice/sim/LHC11a10a_000139110_AOD048#aodTree");
	40	}
	41	if (c == 3){
	42	handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138662_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138666_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000139107_AOD086_p2#aodTree");
	43	}
	44	if (c == 4){
	45	handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138662_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138666_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_139107_AOD090#aodTree");
	46	}
	47
97b01bb0	48	gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
	49	gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g");
	50	gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
	51	gROOT->LoadMacro("AliSpectraAODPID.cxx+g");
	52	gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
0e268798	53
c88234ad	54	handler->SetAliRootMode("default");
97b01bb0	55	handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODPID.cxx AliSpectraAODPID.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
97b01bb0	56	handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODPID.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
0e268798	57	// handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
c88234ad	58	// handler->SetAliRootMode("");
c88234ad	59	handler->SetClearPackages();
ed06106e	60
c88234ad	61	AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
	62	mgr->SetGridHandler(handler);
	63	AliAODInputHandler* aodH = new AliAODInputHandler();
	64	mgr->SetInputEventHandler(aodH);
ed06106e	65
c88234ad	66	// Add PID task
c88234ad	67	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
823864bf	68	Bool_t isMC = kFALSE;
ed06106e	69	if (c == 2 \|\| c == 4) isMC = kTRUE;
823864bf	70	AliAnalysisTask * taskPID = AddTaskPIDResponse(isMC);
c88234ad	71	mgr->AddTask(taskPID);
ed06106e	72
823864bf	73	//LOOP OVER SELECTION
	74	Double_t CentCutMin[4]={0,0,20,20};
	75	Double_t CentCutMax[4]={100,5,40,40};
	76	Double_t QvecCutMin[4]={0,0,0,3};
	77	Double_t QvecCutMax[4]={100,100,2,100};
624a6bb0	78	// Double_t CentCutMin[4]={5,10,20,30};
	79	// Double_t CentCutMax[4]={10,20,30,40};
	80	// Double_t QvecCutMin[4]={0,0,0,0};
	81	// Double_t QvecCutMax[4]={100,100,100,100};
0e268798	82	using namespace AliSpectraNameSpace;
0e268798	83	AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF);
e8b0fe64	84	pid->SetNSigmaCut(5.);
97b01bb0	85
6522a8cb	86	for(Int_t iCut=1;iCut<2;iCut++){
823864bf	87	AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
	88	mgr->AddTask(task);
	89	//physics selection
	90	task->SelectCollisionCandidates();
e8b0fe64	91	// set pid object
e8b0fe64	92	task->SetPID(pid);
823864bf	93	// Set the cuts
	94	AliSpectraAODEventCuts * vcuts = new AliSpectraAODEventCuts("Event Cuts");
	95	AliSpectraAODTrackCuts * tcuts = new AliSpectraAODTrackCuts("Track Cuts");
97b01bb0	96	//if(c==1 \|\| c==2)tcuts->SetTrackType(5); //AOD 046 & 047. Standard Cut with loose DCA
742d45e0	97	if(c==1 \|\| c==2)tcuts->SetTrackType(6); //AOD 046 & 047. Standard Cut with tight DCA
ed06106e	98	if(c==3 \|\| c==4)tcuts->SetTrackType(10); //AOD 086 & 090. Standard Raa cut
823864bf	99	tcuts->SetEta(.8);
823864bf	100	//tcuts->SetDCA(.1);
624a6bb0	101	tcuts->SetPt(5);
97b01bb0	102	tcuts->SetY(.5);
823864bf	103	tcuts->SetPtTOFMatching(0.6);
	104	tcuts->SetQvecMin(QvecCutMin[iCut]);
	105	tcuts->SetQvecMax(QvecCutMax[iCut]);
	106	vcuts->SetCentralityCutMax(CentCutMax[iCut]);
	107	vcuts->SetCentralityCutMin(CentCutMin[iCut]);
	108	task->SetEventCuts(vcuts);
	109	task->SetTrackCuts(tcuts);
823864bf	110	vcuts->PrintCuts();
	111	tcuts->PrintCuts();
	112
	113	// check for MC or real data
ed06106e	114	if (c == 2 \|\| c == 4)
823864bf	115	{
	116	task->SetIsMC(kTRUE);
	117	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	118	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
	119	Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	120	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
	121	Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	122	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
	123	Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
e8b0fe64	124	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
e8b0fe64	125	Form("Pt.AOD.1._MC_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
823864bf	126	}
ed06106e	127	if (c == 1 \|\| c==3)
823864bf	128	{
	129	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	130	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
	131	Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	132	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
	133	Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
	134	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
	135	Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
e8b0fe64	136	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
0e268798	137	Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
823864bf	138
	139	}
	140	mgr->ConnectInput(task, 0, cinput);
	141	mgr->ConnectOutput(task, 1, coutputpt1);
	142	mgr->ConnectOutput(task, 2, coutputpt2);
	143	mgr->ConnectOutput(task, 3, coutputpt3);
e8b0fe64	144	mgr->ConnectOutput(task, 4, coutputpt4);
823864bf	145	}
c88234ad	146	mgr->SetDebugLevel(2);
823864bf	147
c88234ad	148	if (!mgr->InitAnalysis()) return;
c88234ad	149	mgr->PrintStatus();
ed06106e	150	mgr->StartAnalysis("proof");
c88234ad	151	}