[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(80);
  //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_000138666_AOD086_p2#aodTree|/default/lmilano/LHC10h_000139107_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138275_AOD086_p2#aodTree|/default/lmilano/LHC10h_000139465_AOD086_p2#aodTree|/default/lmilano/LHC10h_000139437_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138442_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138396_AOD086_p2#aodTree|/default/lmilano/LHC10h_000138364_AOD086_p2#aodTree");      
  }
  if (c == 4){
    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|/default/lmilano/LHC11a10a_bis_139465_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_139437_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138442_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138396_AOD090#aodTree|/default/lmilano/LHC11a10a_bis_138364_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);
  //Printf("\n\n\n\n\n\n\n\nsetting specific OADB for TPC Response");
  //taskPID->SetOADBPath("/alice/cern.ch/user/a/akalweit/ForLeornado/OADB");
  
  AliAnalysisTaskPIDResponse *taskPID=AddTaskPIDResponse(isMC);
  //taskPID->SetOADBPath("alien:///alice/cern.ch/user/a/akalweit/ForLeornado/OADB");

  //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
  //AddTaskPIDqa();
   
  gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskVZEROEPSelection.C");
  AliVZEROEPSelectionTask *selTask = AddTaskVZEROEPSelection();

  //LOOP OVER SELECTION
  //                            0    1    2    3    4    5
  Double_t CentCutMin[4]= {     0,  20,  20,  20};
  Double_t CentCutMax[4]= {     5,  50,  50,  50};
  Double_t QvecCutMin[4]=    {  0,   0,   0, 1.5};
  Double_t QvecCutMax[4]=   { 100, 100, 0.4, 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};
  if(c==2 || c==4){//for the MC the efficiency doesn't depend on q, we increase the stat
    for(Int_t i=0;i<4;i++){
      QvecCutMin[i]=0;
      QvecCutMax[i]=100;
    }
  }
  AliAnalysisTaskSpectraAOD* taskAOD=AddTaskSpectraAOD();

  // using namespace AliSpectraNameSpace;
  // AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF); 
  // pid->SetNSigmaCut(3.);
  
  // for(Int_t iCut=0;iCut<4;iCut++){
  //   if(iCut!=1)continue; //uncomemnt if you want to run only on 1 subsample
  //   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(1024);
  //   Printf("\n-------------------- tcuts->GetTrackType() %d \n",tcuts->GetTrackType());
  //   //tcuts->SetDCA(.1);
  //   tcuts->SetPt(5);
  //   tcuts->SetY(.5);
  //   tcuts->SetPtTOFMatching(.6);   
  //   //cut on qvector and eta
  //   vcuts->SetQVectorCut(QvecCutMin[iCut],QvecCutMax[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_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
  // 								    Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
  // 								   Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
  // 								    Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),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_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, 
  // 								    Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
  // 								    Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
  // 	AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
  // 								    Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),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
16ee3540	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
829b5a81	29	//handler->SetNproofWorkers(80);
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){
00493191	44	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\|/default/lmilano/LHC10h_000139465_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000139437_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138442_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138396_AOD086_p2#aodTree\|/default/lmilano/LHC10h_000138364_AOD086_p2#aodTree");
ae0fdd7d	45	}
ae0fdd7d	46	if (c == 4){
00493191	47	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\|/default/lmilano/LHC11a10a_bis_139465_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_139437_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138442_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138396_AOD090#aodTree\|/default/lmilano/LHC11a10a_bis_138364_AOD090#aodTree"); }
ed06106e	48
ae0fdd7d	49	gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
	50	gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g");
	51	gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
	52	gROOT->LoadMacro("AliSpectraAODPID.cxx+g");
	53	gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
0e268798	54
ae0fdd7d	55	handler->SetAliRootMode("default");
	56	handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODPID.cxx AliSpectraAODPID.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
	57	handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODPID.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
	58	// handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
	59	// handler->SetAliRootMode("");
	60	handler->SetClearPackages();
ed06106e	61
ae0fdd7d	62	AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
	63	mgr->SetGridHandler(handler);
	64	AliAODInputHandler* aodH = new AliAODInputHandler();
	65	mgr->SetInputEventHandler(aodH);
ed06106e	66
ae0fdd7d	67	// Add PID task
	68	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
	69	Bool_t isMC = kFALSE;
	70	if (c == 2 \|\| c == 4) isMC = kTRUE;
972a21ad	71	Printf("-------------------------------adding in runAOD AddTaskPIDResponse");
16ee3540	72	//AliAnalysisTask * taskPID = AddTaskPIDResponse(isMC);
	73	//Printf("\n\n\n\n\n\n\n\nsetting specific OADB for TPC Response");
	74	//taskPID->SetOADBPath("/alice/cern.ch/user/a/akalweit/ForLeornado/OADB");
	75
	76	AliAnalysisTaskPIDResponse *taskPID=AddTaskPIDResponse(isMC);
	77	//taskPID->SetOADBPath("alien:///alice/cern.ch/user/a/akalweit/ForLeornado/OADB");
	78
972a21ad	79	//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
972a21ad	80	//AddTaskPIDqa();
ed06106e	81
16ee3540	82	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskVZEROEPSelection.C");
	83	AliVZEROEPSelectionTask *selTask = AddTaskVZEROEPSelection();
	84
ae0fdd7d	85	//LOOP OVER SELECTION
829b5a81	86	// 0 1 2 3 4 5
16ee3540	87	Double_t CentCutMin[4]= { 0, 20, 20, 20};
	88	Double_t CentCutMax[4]= { 5, 50, 50, 50};
	89	Double_t QvecCutMin[4]= { 0, 0, 0, 1.5};
	90	Double_t QvecCutMax[4]= { 100, 100, 0.4, 100};
	91	Double_t EtaMin[4]={ -0.8,-0.8,-0.8,-0.8};
	92	Double_t EtaMax[4]={ 0.8, 0.8, 0.8, 0.8};
829b5a81	93	if(c==2 \|\| c==4){//for the MC the efficiency doesn't depend on q, we increase the stat
16ee3540	94	for(Int_t i=0;i<4;i++){
	95	QvecCutMin[i]=0;
	96	QvecCutMax[i]=100;
829b5a81	97	}
829b5a81	98	}
f76247d9	99	AliAnalysisTaskSpectraAOD* taskAOD=AddTaskSpectraAOD();
	100
	101	// using namespace AliSpectraNameSpace;
	102	// AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF);
	103	// pid->SetNSigmaCut(3.);
829b5a81	104
f76247d9	105	// for(Int_t iCut=0;iCut<4;iCut++){
	106	// if(iCut!=1)continue; //uncomemnt if you want to run only on 1 subsample
	107	// AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
	108	// mgr->AddTask(task);
	109	// //physics selection
	110	// task->SelectCollisionCandidates();
	111	// // set pid object
	112	// task->SetPID(pid);
	113	// // Set the cuts
	114	// AliSpectraAODEventCuts * vcuts = new AliSpectraAODEventCuts("Event Cuts");
	115	// AliSpectraAODTrackCuts * tcuts = new AliSpectraAODTrackCuts("Track Cuts");
	116	// //if(c==1 \|\| c==2)tcuts->SetTrackType(5); //AOD 046 & 047. Standard Cut with loose DCA
	117	// //if(c==1 \|\| c==2)tcuts->SetTrackType(6); //AOD 046 & 047. Standard Cut with tight DCA
	118	// //if(c==3 \|\| c==4)tcuts->SetTrackType(10); //AOD 086 & 090. Standard Raa cut
	119	// //if(c==3 \|\| c==4)tcuts->SetTrackType(4); //AOD 086 & 090. Jet analysis
	120	// //if(c==3 \|\| c==4)tcuts->SetTrackType(7); //AOD 086 & 090. TPC Only
	121	// tcuts->SetTrackBits(1024);
	122	// Printf("\n-------------------- tcuts->GetTrackType() %d \n",tcuts->GetTrackType());
	123	// //tcuts->SetDCA(.1);
	124	// tcuts->SetPt(5);
	125	// tcuts->SetY(.5);
	126	// tcuts->SetPtTOFMatching(.6);
	127	// //cut on qvector and eta
	128	// vcuts->SetQVectorCut(QvecCutMin[iCut],QvecCutMax[iCut]);
	129	// tcuts->SetEta(EtaMin[iCut],EtaMax[iCut]);
	130	// if(c==1)vcuts->SetUseCentPatchAOD049(kTRUE);
	131	// vcuts->SetCentralityCutMax(CentCutMax[iCut]);
	132	// vcuts->SetCentralityCutMin(CentCutMin[iCut]);
	133	// task->SetEventCuts(vcuts);
	134	// task->SetTrackCuts(tcuts);
	135	// vcuts->PrintCuts();
	136	// tcuts->PrintCuts();
decf69d9	137
f76247d9	138	// // check for MC or real data
	139	// if (c == 2 \|\| c == 4)
	140	// {
	141	// task->SetIsMC(kTRUE);
	142	// vcuts->SetIsMC(kTRUE);
	143	// AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	144	// AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
	145	// Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	146	// AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
	147	// Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	148	// AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
	149	// Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	150	// AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
	151	// Form("OutputAODSpectraTask_mc_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	152	// }
	153	// if (c == 1 \|\| c==3)
	154	// {
	155	// AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	156	// AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer(Form("chistpt%d",iCut), AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer,
	157	// Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	158	// AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer(Form("cvcutpt%d",iCut), AliSpectraAODEventCuts::Class(), AliAnalysisManager::kOutputContainer,
	159	// Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	160	// AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer,
	161	// Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
	162	// AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut), AliSpectraAODPID::Class(), AliAnalysisManager::kOutputContainer,
	163	// Form("OutputAODSpectraTask_data_Cent%.0fto%.0f_QVec%.1fto%.1f_Eta%.1fto%.1f.root",vcuts->GetCentralityMin(),vcuts->GetCentralityMax(),vcuts->GetQVectorCutMin(), vcuts->GetQVectorCutMax(),tcuts->GetEtaMin(),tcuts->GetEtaMax()));
16ee3540	164
f76247d9	165	// }
	166	// mgr->ConnectInput(task, 0, cinput);
	167	// mgr->ConnectOutput(task, 1, coutputpt1);
	168	// mgr->ConnectOutput(task, 2, coutputpt2);
	169	// mgr->ConnectOutput(task, 3, coutputpt3);
	170	// mgr->ConnectOutput(task, 4, coutputpt4);
	171	// }
	172	// mgr->SetDebugLevel(2);
16ee3540	173
ae0fdd7d	174	if (!mgr->InitAnalysis()) return;
	175	mgr->PrintStatus();
	176	mgr->StartAnalysis("proof");
c88234ad	177	}