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

void runAODProof(Int_t c=4, 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-04-25-AN");
   
  //handler->SetNproofWorkers(1);
  //handler->SetNproofWorkersPerSlave(4);
  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");
    handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD049_p2#aodTree|/default/lmilano/LHC10h_000138662_AOD049_p2#aodTree|/default/lmilano/LHC10h_000138666_AOD049_p2#aodTree|/default/lmilano/LHC10h_000138730_AOD049_p2#aodTree|/default/lmilano/LHC10h_000138732_AOD049_p2#aodTree|/default/lmilano/LHC10h_000139507_AOD049_p2#aodTree|/default/lmilano/LHC10h_000139465_AOD049_p2#aodTree|/default/lmilano/LHC10h_000139437_AOD049_p2#aodTree|/default/lmilano/LHC10h_000139107_AOD049_p2#aodTree|/default/lmilano/LHC10h_000139510_AOD049_p2#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");      
    handler->SetProofDataSet("/default/lmilano/LHC11a10a_138653_AOD048#aodTree|/default/lmilano/LHC11a10a_138662_AOD048#aodTree|/default/lmilano/LHC11a10a_138666_AOD048#aodTree|/default/lmilano/LHC11a10a_138730_AOD048#aodTree|/default/lmilano/LHC11a10a_138732_AOD048#aodTree|/default/lmilano/LHC11a10a_139507_AOD048#aodTree|/default/lmilano/LHC11a10a_139465_AOD048#aodTree|/default/lmilano/LHC11a10a_139437_AOD048#aodTree|/default/lmilano/LHC11a10a_139107_AOD048#aodTree|/default/lmilano/LHC11a10a_139510_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|/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree");      
    handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138666_AOD086_p2#aodTree|/default/lmilano/LHC10h_000139107_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree");      
    //handler->SetProofDataSet("/default/lmilano/LHC10h_000138275_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|/default/lmilano/LHC11a10a_bis_138275_AOD090#aodTree");      
    handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138666_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_139107_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138275_AOD090#aodTree");      
    //handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_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 AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
  handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODEventCuts.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;   
  Printf("-------------------------------adding in runAOD AddTaskPIDResponse");
  AliAnalysisTask * taskPID = AddTaskPIDResponse(isMC);
  //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
  //AddTaskPIDqa();
   
  //LOOP OVER SELECTION
  Double_t CentCutMin[4]={0,0,20,20};
  Double_t CentCutMax[4]={100,5,40,40};
  Double_t QvecPosCutMin[4]={0,0,0,3};
  Double_t QvecPosCutMax[4]={100,100,2,100};
  Double_t QvecNegCutMin[4]={0,0,0,3};
  Double_t QvecNegCutMax[4]={100,100,2,100};
  Double_t EtaMin[4]={-0.8,-0.8,-0.8,-0.8};
  Double_t EtaMax[4]={0.8,0.8,0.8,0.8};
  using namespace AliSpectraNameSpace;
  AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF); 
  pid->SetNSigmaCut(3.);
   
  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
    //if(c==3 || c==4)tcuts->SetTrackType(4); //AOD 086 & 090. Jet analysis
    //if(c==3 || c==4)tcuts->SetTrackType(7); //AOD 086 & 090. TPC Only
    tcuts->SetTrackBits(1);
    Printf("\n\n\n\n\n\n\n-------------------- tcuts->GetTrackType() %d \n\n\n\n\n\n\n",tcuts->GetTrackType());
    //tcuts->SetDCA(.1);
    tcuts->SetPt(5);
    tcuts->SetY(.5);
    tcuts->SetPtTOFMatching(.6);   
    //cut on qvector and eta
    vcuts->SetQVectorPosCut(QvecPosCutMin[iCut],QvecPosCutMax[iCut]);
    vcuts->SetQVectorNegCut(QvecNegCutMin[iCut],QvecNegCutMax[iCut]);
    tcuts->SetEta(EtaMin[iCut],EtaMax[iCut]);
    if(c==1)vcuts->SetUseCentPatchAOD049(kTRUE);
    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);
	vcuts->SetIsMC(kTRUE);
	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, 
								    Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
								    Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
		AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
								    Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
								   Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
      }
    if (c == 1 || c==3)
      {
	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, 
								    Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
								  Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
								  Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
								  Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
									 vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	 
      }
    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
decf69d9	1	void runAODProof(Int_t c=4, 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");
ae0fdd7d	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");
c88234ad	20
ae0fdd7d	21	AliAnalysisAlien * handler = new AliAnalysisAlien("test");
	22	handler->SetOverwriteMode();
	23	handler->SetRunMode(proofMode);
	24	handler->SetProofReset(0);
	25	//handler->SetROOTVersion("v5-33-02a");
	26	//handler->SetAliROOTVersion("v5-03-11-AN");
	27	handler->SetAliROOTVersion("v5-04-25-AN");
6522a8cb	28
972a21ad	29	//handler->SetNproofWorkers(1);
ae0fdd7d	30	//handler->SetNproofWorkersPerSlave(4);
	31	handler->SetProofCluster(Form("%s@alice-caf.cern.ch", gSystem->Getenv("CAFUSER")));
	32	//handler->SetProofCluster(Form("%s@skaf.saske.sk",gSystem->Getenv("CAFUSER")));
	33	// Set handler for Real DATA:
	34	if (c == 1){
	35	//handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree");
	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");
	37	handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000138662_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000138666_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000138730_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000138732_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000139507_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000139465_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000139437_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000139107_AOD049_p2#aodTree\|/default/lmilano/LHC10h_000139510_AOD049_p2#aodTree");
	38	}
	39	if (c == 2){
	40	//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");
decf69d9	41	handler->SetProofDataSet("/default/lmilano/LHC11a10a_138653_AOD048#aodTree\|/default/lmilano/LHC11a10a_138662_AOD048#aodTree\|/default/lmilano/LHC11a10a_138666_AOD048#aodTree\|/default/lmilano/LHC11a10a_138730_AOD048#aodTree\|/default/lmilano/LHC11a10a_138732_AOD048#aodTree\|/default/lmilano/LHC11a10a_139507_AOD048#aodTree\|/default/lmilano/LHC11a10a_139465_AOD048#aodTree\|/default/lmilano/LHC11a10a_139437_AOD048#aodTree\|/default/lmilano/LHC11a10a_139107_AOD048#aodTree\|/default/lmilano/LHC11a10a_139510_AOD048#aodTree");
ae0fdd7d	42	}
ae0fdd7d	43	if (c == 3){
936e2842	44	//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\|/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree");
8961bb9a	45	handler->SetProofDataSet("/default/lmilano/LHC10h_000138653_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138666_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000139107_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree");
8961bb9a	46	//handler->SetProofDataSet("/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree");
ae0fdd7d	47	}
ae0fdd7d	48	if (c == 4){
8961bb9a	49	//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\|/default/lmilano/LHC11a10a_bis_138275_AOD090#aodTree");
	50	handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138666_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_139107_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138275_AOD090#aodTree");
	51	//handler->SetProofDataSet("/default/lmilano/LHC11a10a_bis_138653_AOD090#aodTree");
ae0fdd7d	52	}
ed06106e	53
ae0fdd7d	54	gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
	55	gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g");
	56	gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
	57	gROOT->LoadMacro("AliSpectraAODPID.cxx+g");
	58	gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
0e268798	59
ae0fdd7d	60	handler->SetAliRootMode("default");
	61	handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODPID.cxx AliSpectraAODPID.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
	62	handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODPID.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
	63	// handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
	64	// handler->SetAliRootMode("");
	65	handler->SetClearPackages();
ed06106e	66
ae0fdd7d	67	AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
	68	mgr->SetGridHandler(handler);
	69	AliAODInputHandler* aodH = new AliAODInputHandler();
	70	mgr->SetInputEventHandler(aodH);
ed06106e	71
ae0fdd7d	72	// Add PID task
	73	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
	74	Bool_t isMC = kFALSE;
	75	if (c == 2 \|\| c == 4) isMC = kTRUE;
972a21ad	76	Printf("-------------------------------adding in runAOD AddTaskPIDResponse");
ae0fdd7d	77	AliAnalysisTask * taskPID = AddTaskPIDResponse(isMC);
972a21ad	78	//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
972a21ad	79	//AddTaskPIDqa();
ed06106e	80
ae0fdd7d	81	//LOOP OVER SELECTION
	82	Double_t CentCutMin[4]={0,0,20,20};
	83	Double_t CentCutMax[4]={100,5,40,40};
decf69d9	84	Double_t QvecPosCutMin[4]={0,0,0,3};
	85	Double_t QvecPosCutMax[4]={100,100,2,100};
	86	Double_t QvecNegCutMin[4]={0,0,0,3};
	87	Double_t QvecNegCutMax[4]={100,100,2,100};
	88	Double_t EtaMin[4]={-0.8,-0.8,-0.8,-0.8};
	89	Double_t EtaMax[4]={0.8,0.8,0.8,0.8};
ae0fdd7d	90	using namespace AliSpectraNameSpace;
	91	AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF);
	92	pid->SetNSigmaCut(3.);
97b01bb0	93
ae0fdd7d	94	for(Int_t iCut=1;iCut<2;iCut++){
	95	AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
	96	mgr->AddTask(task);
	97	//physics selection
	98	task->SelectCollisionCandidates();
	99	// set pid object
	100	task->SetPID(pid);
	101	// Set the cuts
	102	AliSpectraAODEventCuts * vcuts = new AliSpectraAODEventCuts("Event Cuts");
	103	AliSpectraAODTrackCuts * tcuts = new AliSpectraAODTrackCuts("Track Cuts");
936e2842	104	//if(c==1 \|\| c==2)tcuts->SetTrackType(5); //AOD 046 & 047. Standard Cut with loose DCA
ae0fdd7d	105	//if(c==1 \|\| c==2)tcuts->SetTrackType(6); //AOD 046 & 047. Standard Cut with tight DCA
972a21ad	106	//if(c==3 \|\| c==4)tcuts->SetTrackType(10); //AOD 086 & 090. Standard Raa cut
972a21ad	107	//if(c==3 \|\| c==4)tcuts->SetTrackType(4); //AOD 086 & 090. Jet analysis
936e2842	108	//if(c==3 \|\| c==4)tcuts->SetTrackType(7); //AOD 086 & 090. TPC Only
8961bb9a	109	tcuts->SetTrackBits(1);
8961bb9a	110	Printf("\n\n\n\n\n\n\n-------------------- tcuts->GetTrackType() %d \n\n\n\n\n\n\n",tcuts->GetTrackType());
ae0fdd7d	111	//tcuts->SetDCA(.1);
	112	tcuts->SetPt(5);
	113	tcuts->SetY(.5);
	114	tcuts->SetPtTOFMatching(.6);
decf69d9	115	//cut on qvector and eta
	116	vcuts->SetQVectorPosCut(QvecPosCutMin[iCut],QvecPosCutMax[iCut]);
	117	vcuts->SetQVectorNegCut(QvecNegCutMin[iCut],QvecNegCutMax[iCut]);
	118	tcuts->SetEta(EtaMin[iCut],EtaMax[iCut]);
936e2842	119	if(c==1)vcuts->SetUseCentPatchAOD049(kTRUE);
ae0fdd7d	120	vcuts->SetCentralityCutMax(CentCutMax[iCut]);
	121	vcuts->SetCentralityCutMin(CentCutMin[iCut]);
	122	task->SetEventCuts(vcuts);
	123	task->SetTrackCuts(tcuts);
	124	vcuts->PrintCuts();
	125	tcuts->PrintCuts();
decf69d9	126
ae0fdd7d	127	// check for MC or real data
	128	if (c == 2 \|\| c == 4)
	129	{
	130	task->SetIsMC(kTRUE);
972a21ad	131	vcuts->SetIsMC(kTRUE);
ae0fdd7d	132	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
ae0fdd7d	133	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	134	Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	135	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	136	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	137	Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
	138	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	139	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
	140	Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
	141	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	142	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	143	Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	144	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	145	}
	146	if (c == 1 \|\| c==3)
	147	{
	148	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	149	AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	150	Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	151	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	152	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	153	Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	154	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	155	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	156	Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	157	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
ae0fdd7d	158	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
decf69d9	159	Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVecPos%.1fto%.1f_QVecNeg%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),
decf69d9	160	vcuts->GetQVectorPosCutMin(), vcuts->GetQVectorPosCutMax(),vcuts->GetQVectorNegCutMin(), vcuts->GetQVectorNegCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
823864bf	161
ae0fdd7d	162	}
	163	mgr->ConnectInput(task, 0, cinput);
	164	mgr->ConnectOutput(task, 1, coutputpt1);
	165	mgr->ConnectOutput(task, 2, coutputpt2);
	166	mgr->ConnectOutput(task, 3, coutputpt3);
	167	mgr->ConnectOutput(task, 4, coutputpt4);
	168	}
	169	mgr->SetDebugLevel(2);
823864bf	170
ae0fdd7d	171	if (!mgr->InitAnalysis()) return;
	172	mgr->PrintStatus();
	173	mgr->StartAnalysis("proof");
c88234ad	174	}