[u/mrichter/AliRoot.git] / HLT / QA / tasks / macros / compare-HLT-offline-local.C

// $Id$
/*
 * Example macro to run locally an analysis task for comparing the offline
 * with the HLT esd tree.
 *
 * The output is a root file containing the histograms defined in the
 * analysis task. There is one output file per detector.
 *
 * Run without arguments to get a few examples how to use the macro.
 *
 * Usage:
 * <pre>
 *   aliroot -b -l -q compare-HLT-offline-local.C'("/home/blabla/AliESDs.root","global","./",kTRUE,10)' 2>&1 | tee task.log
 *   aliroot -b -l -q compare-HLT-offline-local.C'("/home/blabla/AliESDs.root","phos global cb","./",kTRUE,100)' 2>&1 | tee task.log
 *   aliroot -b -l -q compare-HLT-offline-local.C'("alien:///alice/data/2010/LHC10b/000115322/ESDs/pass1/10000115322040.20/AliESDs.root","global")' 2>&1 | tee log
 * </pre>
 * 
 * If alien:// is contained in the name of the file, then the macro connects to the grid to access the file.
 * 
 * In case you want to run over many ESD files, then prepare a list of them in a .txt file and they will be chained for the analysis.
 * The .txt file takes the place of the first argument in that case.
 *
 * @ingroup alihlt_qa
 * @author Kalliopi.Kanaki@ift.uib.no, Hege.Erdal@student.uib.no
 */

void compare_HLT_offline_local(TString file, 
                               const char* detectorTask="global",
			       TString taskFolder="$ALICE_ROOT/HLT/QA/tasks/", 
			       bool fUseHLTTrigger=kFALSE, 
			       Long64_t nEvents=1234567890
			      )
{

  TStopwatch timer;
  timer.Start();

  gSystem->Load("libTree.so");
  gSystem->Load("libGeom.so");
  gSystem->Load("libVMC.so");
  gSystem->Load("libPhysics.so");
 
  //----------- Loading the required libraries ---------//

  gSystem->Load("libSTEERBase.so");
  gSystem->Load("libESD.so");
  gSystem->Load("libAOD.so");
  gSystem->Load("libANALYSIS.so");
  gSystem->Load("libANALYSISalice.so");
  gSystem->Load("libHLTbase.so");
 
  gSystem->AddIncludePath("-I$(ALICE_ROOT)/HLT/BASE -I$ALICE_ROOT/PWG1/TPC -I.");
  
  gSystem->Load("libTPCcalib.so");
  gSystem->Load("libTRDbase.so");
  gSystem->Load("libTRDrec.so");
  gSystem->Load("libITSbase.so");
  gSystem->Load("libITSrec.so");
  gSystem->Load("libTENDER.so");
  gSystem->Load("libPWG1.so");
 
  gROOT->ProcessLine(".include $ALICE_ROOT/include");
  //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
  
  Bool_t bPHOS = kFALSE, bGLOBAL = kFALSE, bEMCAL = kFALSE, bPWG1 = kFALSE, bD0 = kFALSE, bCB = kFALSE;
 
  TString allArgs = detectorTask;
  TString argument;
 
  TObjArray *pTokens = allArgs.Tokenize(" ");
  if(pTokens){
    for(int i=0; i<pTokens->GetEntries(); i++){
      argument=((TObjString*)pTokens->At(i))->GetString();
      if(argument.IsNull()) continue;

      if(argument.CompareTo("phos", TString::kIgnoreCase)==0){
	bPHOS = kTRUE;
	continue;
      }         
      else if(argument.CompareTo("emcal", TString::kIgnoreCase)==0){
	bEMCAL = kTRUE;
	continue;
      }         
      if(argument.CompareTo("global", TString::kIgnoreCase)==0){
	bGLOBAL = kTRUE;
	continue;
      }      
      if(argument.CompareTo("pwg1", TString::kIgnoreCase)==0){
	bPWG1 = kTRUE;
	continue;
      }
      if(argument.CompareTo("D0", TString::kIgnoreCase)==0){
	bD0 = kTRUE;
	continue;
      }
      if(argument.CompareTo("cb", TString::kIgnoreCase)==0){
	bCB = kTRUE;
	continue;
      } 
      if(argument.CompareTo("all",TString::kIgnoreCase)==0){
	bPHOS   = kTRUE;
	bEMCAL  = kTRUE;
	bGLOBAL = kTRUE; 
	bD0     = kTRUE;  
	bCB     = kTRUE; 
	continue;
      }
      else break;
    }
  }
      
  //-------------- Compile the analysis tasks ---------- //
  
  if(bPHOS){
    //gSystem->Load("libHLTbase");
    gSystem->Load("libAliHLTUtil");
    gSystem->Load("libAliHLTGlobal");
    TString strTask1("AliAnalysisTaskHLTCalo.cxx+");
    TString strTask2("AliAnalysisTaskHLTPHOS.cxx+");    
    gROOT->LoadMacro(taskFolder+strTask1); 
    gROOT->LoadMacro(taskFolder+strTask2); 
    cout << "\n========= You are loading the following tasks --> "<< (taskFolder+strTask1).Chop()  << " and " <<  (taskFolder+strTask2).Chop() << endl;
  }
  
  if(bEMCAL){
    //gSystem->Load("libHLTbase");
    gSystem->Load("libAliHLTUtil");
    gSystem->Load("libAliHLTGlobal");
    TString strTask1("AliAnalysisTaskHLTCalo.cxx+");
    TString strTask2("AliAnalysisTaskHLTEMCAL.cxx+");
    gROOT->LoadMacro(taskFolder+strTask1); 
    gROOT->LoadMacro(taskFolder+strTask2); 
    cout << "\n========= You are loading the following tasks --> "<< (taskFolder+strTask1).Chop()  << " and " <<  (taskFolder+strTask2).Chop() << endl;
  }  
  
  if(bGLOBAL){
     TString strTask("AliAnalysisTaskHLT.cxx+");
     gROOT->LoadMacro(taskFolder+strTask);
     cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop()  << endl;
  }
  if(bD0){
     TString strTask("AliAnalysisTaskD0Trigger.cxx+");
     gROOT->LoadMacro(taskFolder+strTask); 
     cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop()  << endl;
  }
 
  if(bCB){
    TString strTask("AliAnalysisTaskHLTCentralBarrel.cxx+");
    gROOT->LoadMacro(taskFolder+strTask);
    cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop()  << endl;
  }
  
  if(bPWG1) gROOT->LoadMacro("$ALICE_ROOT/HLT/QA/tasks/macros/AddTaskPerformance.C");
   
  if(file.Contains("alien")) TGrid::Connect("alien://");
    
  if(file.Contains("AliESDs.root")){
    TChain *chain = new TChain("esdTree"); 
    chain->Add(file);
  }
  
  //Constructs chain from filenames in *.txt
  //on the form $DIR/AliESDs.root
  else if(file.Contains(".txt")){
    gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
    chain=CreateESDChain(file.Data());
  }
  
  else if(!file){
    printf("File %s does not exist or is corrupted.\n",file.Data());
    return;  
  }

  if(!chain){
    Printf("Chain is empty");
    return;
  }
   
  //-------- Make the analysis manager ---------------//
 
  AliAnalysisManager *mgr  = new AliAnalysisManager("TestManager");
  AliESDInputHandler *esdH = new AliESDInputHandler;

  //For the PWG1 task, setting HLT is handled inside AliPerformanceTask.C
  if(!bPWG1)  esdH->SetReadHLT();
  esdH->SetReadFriends(kFALSE);
  mgr->SetInputEventHandler(esdH);  
  mgr->SetNSysInfo(1000);

  //AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(kFALSE,kTRUE);
 
  //-------------- define the tasks ------------//
  
  if(bPHOS){
    AliAnalysisTaskHLTPHOS *taskPHOS = new AliAnalysisTaskHLTPHOS("offhlt_comparison_PHOS");
    taskPHOS->SetUseHLTTriggerDecision(fUseHLTTrigger);
    if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskPHOS->GetName());
    mgr->AddTask(taskPHOS);
    if(fUseHLTTrigger==kFALSE)
       AliAnalysisDataContainer *coutputPHOS =  mgr->CreateContainer("phos_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-PHOS-comparison.root");  
    else 
       AliAnalysisDataContainer *coutputPHOS =  mgr->CreateContainer("phos_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-PHOS-comparison_triggered.root");      
    mgr->ConnectInput(taskPHOS,0,mgr->GetCommonInputContainer());
    mgr->ConnectOutput(taskPHOS,1,coutputPHOS);
  }

  if(bEMCAL){
    AliAnalysisTaskHLTEMCAL *taskEMCAL = new AliAnalysisTaskHLTEMCAL("offhlt_comparison_EMCAL");
    taskEMCAL->SetUseHLTTriggerDecision(fUseHLTTrigger);
    if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskEMCAL->GetName());
    mgr->AddTask(taskEMCAL);
    if(fUseHLTTrigger==kFALSE)
       AliAnalysisDataContainer *coutputEMCAL =  mgr->CreateContainer("emcal_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-EMCAL-comparison.root");  
    else
      AliAnalysisDataContainer *coutputEMCAL =  mgr->CreateContainer("emcal_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-EMCAL-comparison_triggered.root");      
    mgr->ConnectInput(taskEMCAL,0,mgr->GetCommonInputContainer());
    mgr->ConnectOutput(taskEMCAL,1,coutputEMCAL);
  }
  
  if(bGLOBAL){
    AliAnalysisTaskHLT *taskGLOBAL = new AliAnalysisTaskHLT("offhlt_comparison_GLOBAL");
    taskGLOBAL->SetUseHLTTriggerDecision(fUseHLTTrigger);
    if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskGLOBAL->GetName());
    //taskGLOBAL->SelectCollisionCandidates();
    mgr->AddTask(taskGLOBAL);
    if(fUseHLTTrigger==kFALSE)
      AliAnalysisDataContainer *coutputGLOBAL =  mgr->CreateContainer("global_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-GLOBAL-comparison.root");  
    else 
      AliAnalysisDataContainer *coutputGLOBAL =  mgr->CreateContainer("global_histograms",TList::Class(), AliAnalysisManager::kOutputContainer,"HLT-OFFLINE-GLOBAL-comparison_triggered.root");  
    mgr->ConnectInput(taskGLOBAL,0,mgr->GetCommonInputContainer());
    mgr->ConnectOutput(taskGLOBAL,1,coutputGLOBAL);
  }

  if(bPWG1){
    Bool_t hasMC=kFALSE;  
    // -- Add Task for HLT and Offline
    AliPerformanceTask *HLTtpcQA = AddTaskPerformance(hasMC,kFALSE,kTRUE);
    AliPerformanceTask *tpcQA = AddTaskPerformance(hasMC,kFALSE); 
    if(!HLTtpcQA || !tpcQA) {
      Error("RunPerformanceTrain","AliPerformanceTask not created!");
      return;
    }
  }
  if(bD0){
    float cuts[7]={0.5,0.04,0.7,0.8,0.05,-0.00025,0.7};
    AliAnalysisTaskD0Trigger *taskD0 = new AliAnalysisTaskD0Trigger("offhlt_comparison_D0",cuts);
    mgr->AddTask(taskD0);
    AliAnalysisDataContainer *coutputD0 =  mgr->CreateContainer("D0_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-D0-comparison.root");  
    mgr->ConnectInput(taskD0,0,mgr->GetCommonInputContainer());
    mgr->ConnectOutput(taskD0,1,coutputD0);
  }  
  
  if(bCB){
     AliAnalysisTaskHLTCentralBarrel *taskCB = new AliAnalysisTaskHLTCentralBarrel("offhlt_comparison_CB");    
     AliAnalysisDataContainer *coutputCB =  mgr->CreateContainer("esd_thnsparse", TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-CentralBarrel-comparison.root");  
     mgr->ConnectInput(taskCB,0,mgr->GetCommonInputContainer());
     mgr->ConnectOutput(taskCB,1,coutputCB);
  }
  
  if (!mgr->InitAnalysis()) return;
  mgr->PrintStatus();
  mgr->StartAnalysis("local",chain, nEvents);

  timer.Stop();
  timer.Print();
}

void compare_HLT_offline_local(){
  cout << " " << endl;
  cout << " Usage examples:" << endl;
  cout << "    compare-HLT-offline-local.C'(file, taskOption, taskFolder, fUseHLTTrigger, nEvents)' 2>&1 | tee log" << endl;
  cout << "    compare-HLT-offline-local.C'(\"AliESDs.root\",\"global\")' 2>&1 | tee log" << endl;
  cout << "    compare-HLT-offline-local.C'(\"AliESDs.root\",\"global\",\"./\",kFALSE,nEvents)' 2>&1 | tee log" << endl;
  cout << "    compare-HLT-offline-local.C'(\"AliESDs.root\",\"global phos cb D0\", \"./\", kTRUE, nEvents)' 2>&1 | tee log" << endl;
  cout << "    compare-HLT-offline-local.C'(\"files.txt\",\"all\")' 2>&1 | tee log" << endl;
  cout << "    compare-HLT-offline-local.C'(\"alien:///alice/data/2010/LHC10b/000115322/ESDs/pass1/10000115322040.20/AliESDs.root\",\"global\")' 2>&1 | tee log" << endl;
  cout << " " << endl;
}
Commit	Line	Data
be1fa04a	1	// $Id$
	2	/*
	3	* Example macro to run locally an analysis task for comparing the offline
	4	* with the HLT esd tree.
	5	*
7742dbd5	6	* The output is a root file containing the histograms defined in the
7742dbd5	7	* analysis task. There is one output file per detector.
be1fa04a	8	*
4b621d9c	9	* Run without arguments to get a few examples how to use the macro.
4b621d9c	10	*
be1fa04a	11	* Usage:
be1fa04a	12	* <pre>
ff84363f	13	* aliroot -b -l -q compare-HLT-offline-local.C'("/home/blabla/AliESDs.root","global","./",kTRUE,10)' 2>&1 \| tee task.log
3d9c49a9	14	* aliroot -b -l -q compare-HLT-offline-local.C'("/home/blabla/AliESDs.root","phos global cb","./",kTRUE,100)' 2>&1 \| tee task.log
ff84363f	15	* aliroot -b -l -q compare-HLT-offline-local.C'("alien:///alice/data/2010/LHC10b/000115322/ESDs/pass1/10000115322040.20/AliESDs.root","global")' 2>&1 \| tee log
be1fa04a	16	* </pre>
b3406717	17	*
b3406717	18	* If alien:// is contained in the name of the file, then the macro connects to the grid to access the file.
1e81c36d	19	*
	20	* In case you want to run over many ESD files, then prepare a list of them in a .txt file and they will be chained for the analysis.
	21	* The .txt file takes the place of the first argument in that case.
7742dbd5	22	*
b71630ca	23	* @ingroup alihlt_qa
4b621d9c	24	* @author Kalliopi.Kanaki@ift.uib.no, Hege.Erdal@student.uib.no
be1fa04a	25	*/
be1fa04a	26
4b621d9c	27	void compare_HLT_offline_local(TString file,
	28	const char* detectorTask="global",
	29	TString taskFolder="$ALICE_ROOT/HLT/QA/tasks/",
	30	bool fUseHLTTrigger=kFALSE,
	31	Long64_t nEvents=1234567890
	32	)
	33	{
fc0eba5b	34
be1fa04a	35	TStopwatch timer;
	36	timer.Start();
	37
	38	gSystem->Load("libTree.so");
	39	gSystem->Load("libGeom.so");
	40	gSystem->Load("libVMC.so");
	41	gSystem->Load("libPhysics.so");
	42
	43	//----------- Loading the required libraries ---------//
	44
	45	gSystem->Load("libSTEERBase.so");
	46	gSystem->Load("libESD.so");
	47	gSystem->Load("libAOD.so");
	48	gSystem->Load("libANALYSIS.so");
	49	gSystem->Load("libANALYSISalice.so");
	50	gSystem->Load("libHLTbase.so");
fc0eba5b	51
394b9d00	52	gSystem->AddIncludePath("-I$(ALICE_ROOT)/HLT/BASE -I$ALICE_ROOT/PWG1/TPC -I.");
da6744ee	53
	54	gSystem->Load("libTPCcalib.so");
	55	gSystem->Load("libTRDbase.so");
	56	gSystem->Load("libTRDrec.so");
	57	gSystem->Load("libITSbase.so");
	58	gSystem->Load("libITSrec.so");
	59	gSystem->Load("libTENDER.so");
	60	gSystem->Load("libPWG1.so");
	61
be1fa04a	62	gROOT->ProcessLine(".include $ALICE_ROOT/include");
4b621d9c	63	//gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
7742dbd5	64
ff84363f	65	Bool_t bPHOS = kFALSE, bGLOBAL = kFALSE, bEMCAL = kFALSE, bPWG1 = kFALSE, bD0 = kFALSE, bCB = kFALSE;
7742dbd5	66
	67	TString allArgs = detectorTask;
	68	TString argument;
	69
	70	TObjArray *pTokens = allArgs.Tokenize(" ");
	71	if(pTokens){
da6744ee	72	for(int i=0; i<pTokens->GetEntries(); i++){
	73	argument=((TObjString*)pTokens->At(i))->GetString();
	74	if(argument.IsNull()) continue;
be1fa04a	75
da6744ee	76	if(argument.CompareTo("phos", TString::kIgnoreCase)==0){
	77	bPHOS = kTRUE;
	78	continue;
	79	}
	80	else if(argument.CompareTo("emcal", TString::kIgnoreCase)==0){
	81	bEMCAL = kTRUE;
	82	continue;
	83	}
da6744ee	84	if(argument.CompareTo("global", TString::kIgnoreCase)==0){
	85	bGLOBAL = kTRUE;
	86	continue;
	87	}
	88	if(argument.CompareTo("pwg1", TString::kIgnoreCase)==0){
	89	bPWG1 = kTRUE;
	90	continue;
	91	}
cd785861	92	if(argument.CompareTo("D0", TString::kIgnoreCase)==0){
	93	bD0 = kTRUE;
	94	continue;
	95	}
ff84363f	96	if(argument.CompareTo("cb", TString::kIgnoreCase)==0){
	97	bCB = kTRUE;
	98	continue;
	99	}
da6744ee	100	if(argument.CompareTo("all",TString::kIgnoreCase)==0){
da6744ee	101	bPHOS = kTRUE;
da6744ee	102	bEMCAL = kTRUE;
ea46fe48	103	bGLOBAL = kTRUE;
ff84363f	104	bD0 = kTRUE;
ff84363f	105	bCB = kTRUE;
da6744ee	106	continue;
	107	}
	108	else break;
7742dbd5	109	}
7742dbd5	110	}
ea46fe48	111
7742dbd5	112	//-------------- Compile the analysis tasks ---------- //
fc0eba5b	113
ea46fe48	114	if(bPHOS){
7e9b72dc	115	//gSystem->Load("libHLTbase");
ea46fe48	116	gSystem->Load("libAliHLTUtil");
ea46fe48	117	gSystem->Load("libAliHLTGlobal");
7e9b72dc	118	TString strTask1("AliAnalysisTaskHLTCalo.cxx+");
7e9b72dc	119	TString strTask2("AliAnalysisTaskHLTPHOS.cxx+");
4b621d9c	120	gROOT->LoadMacro(taskFolder+strTask1);
4b621d9c	121	gROOT->LoadMacro(taskFolder+strTask2);
394b9d00	122	cout << "\n========= You are loading the following tasks --> "<< (taskFolder+strTask1).Chop() << " and " << (taskFolder+strTask2).Chop() << endl;
fc0eba5b	123	}
fc0eba5b	124
ea46fe48	125	if(bEMCAL){
7e9b72dc	126	//gSystem->Load("libHLTbase");
ea46fe48	127	gSystem->Load("libAliHLTUtil");
ea46fe48	128	gSystem->Load("libAliHLTGlobal");
7e9b72dc	129	TString strTask1("AliAnalysisTaskHLTCalo.cxx+");
7e9b72dc	130	TString strTask2("AliAnalysisTaskHLTEMCAL.cxx+");
4b621d9c	131	gROOT->LoadMacro(taskFolder+strTask1);
4b621d9c	132	gROOT->LoadMacro(taskFolder+strTask2);
394b9d00	133	cout << "\n========= You are loading the following tasks --> "<< (taskFolder+strTask1).Chop() << " and " << (taskFolder+strTask2).Chop() << endl;
fc0eba5b	134	}
fc0eba5b	135
4b621d9c	136	if(bGLOBAL){
	137	TString strTask("AliAnalysisTaskHLT.cxx+");
	138	gROOT->LoadMacro(taskFolder+strTask);
394b9d00	139	cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop() << endl;
4b621d9c	140	}
	141	if(bD0){
	142	TString strTask("AliAnalysisTaskD0Trigger.cxx+");
	143	gROOT->LoadMacro(taskFolder+strTask);
394b9d00	144	cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop() << endl;
4b621d9c	145	}
ff84363f	146
	147	if(bCB){
	148	TString strTask("AliAnalysisTaskHLTCentralBarrel.cxx+");
	149	gROOT->LoadMacro(taskFolder+strTask);
	150	cout << "\n========= You are loading the following task --> "<< (taskFolder+strTask).Chop() << endl;
	151	}
4b621d9c	152
4b621d9c	153	if(bPWG1) gROOT->LoadMacro("$ALICE_ROOT/HLT/QA/tasks/macros/AddTaskPerformance.C");
ea46fe48	154
b3406717	155	if(file.Contains("alien")) TGrid::Connect("alien://");
1e81c36d	156
	157	if(file.Contains("AliESDs.root")){
	158	TChain *chain = new TChain("esdTree");
	159	chain->Add(file);
	160	}
	161
	162	//Constructs chain from filenames in *.txt
	163	//on the form $DIR/AliESDs.root
	164	else if(file.Contains(".txt")){
	165	gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
	166	chain=CreateESDChain(file.Data());
	167	}
3d9c49a9	168
d8cdf034	169	else if(!file){
3d9c49a9	170	printf("File %s does not exist or is corrupted.\n",file.Data());
	171	return;
	172	}
1e81c36d	173
	174	if(!chain){
	175	Printf("Chain is empty");
	176	return;
	177	}
d4598f43	178
be1fa04a	179	//-------- Make the analysis manager ---------------//
	180
	181	AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
	182	AliESDInputHandler *esdH = new AliESDInputHandler;
ea46fe48	183
	184	//For the PWG1 task, setting HLT is handled inside AliPerformanceTask.C
	185	if(!bPWG1) esdH->SetReadHLT();
fcf71d87	186	esdH->SetReadFriends(kFALSE);
be1fa04a	187	mgr->SetInputEventHandler(esdH);
be1fa04a	188	mgr->SetNSysInfo(1000);
da6744ee	189
ea46fe48	190	//AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(kFALSE,kTRUE);
be1fa04a	191
	192	//-------------- define the tasks ------------//
	193
7742dbd5	194	if(bPHOS){
da6744ee	195	AliAnalysisTaskHLTPHOS *taskPHOS = new AliAnalysisTaskHLTPHOS("offhlt_comparison_PHOS");
7e9b72dc	196	taskPHOS->SetUseHLTTriggerDecision(fUseHLTTrigger);
7e9b72dc	197	if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskPHOS->GetName());
da6744ee	198	mgr->AddTask(taskPHOS);
7e9b72dc	199	if(fUseHLTTrigger==kFALSE)
	200	AliAnalysisDataContainer *coutputPHOS = mgr->CreateContainer("phos_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-PHOS-comparison.root");
	201	else
	202	AliAnalysisDataContainer *coutputPHOS = mgr->CreateContainer("phos_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-PHOS-comparison_triggered.root");
da6744ee	203	mgr->ConnectInput(taskPHOS,0,mgr->GetCommonInputContainer());
4b621d9c	204	mgr->ConnectOutput(taskPHOS,1,coutputPHOS);
7742dbd5	205	}
fc0eba5b	206
fc0eba5b	207	if(bEMCAL){
da6744ee	208	AliAnalysisTaskHLTEMCAL *taskEMCAL = new AliAnalysisTaskHLTEMCAL("offhlt_comparison_EMCAL");
7e9b72dc	209	taskEMCAL->SetUseHLTTriggerDecision(fUseHLTTrigger);
7e9b72dc	210	if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskEMCAL->GetName());
da6744ee	211	mgr->AddTask(taskEMCAL);
7e9b72dc	212	if(fUseHLTTrigger==kFALSE)
	213	AliAnalysisDataContainer *coutputEMCAL = mgr->CreateContainer("emcal_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-EMCAL-comparison.root");
	214	else
	215	AliAnalysisDataContainer *coutputEMCAL = mgr->CreateContainer("emcal_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-EMCAL-comparison_triggered.root");
da6744ee	216	mgr->ConnectInput(taskEMCAL,0,mgr->GetCommonInputContainer());
4b621d9c	217	mgr->ConnectOutput(taskEMCAL,1,coutputEMCAL);
fc0eba5b	218	}
a63c69d8	219
db575e59	220	if(bGLOBAL){
da6744ee	221	AliAnalysisTaskHLT *taskGLOBAL = new AliAnalysisTaskHLT("offhlt_comparison_GLOBAL");
4b621d9c	222	taskGLOBAL->SetUseHLTTriggerDecision(fUseHLTTrigger);
4b621d9c	223	if(fUseHLTTrigger==kTRUE) printf("\n\nOnly HLT triggered events will be used to fill the distributions for task %s.\n\n", taskGLOBAL->GetName());
ea46fe48	224	//taskGLOBAL->SelectCollisionCandidates();
da6744ee	225	mgr->AddTask(taskGLOBAL);
4b621d9c	226	if(fUseHLTTrigger==kFALSE)
	227	AliAnalysisDataContainer *coutputGLOBAL = mgr->CreateContainer("global_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-GLOBAL-comparison.root");
	228	else
	229	AliAnalysisDataContainer *coutputGLOBAL = mgr->CreateContainer("global_histograms",TList::Class(), AliAnalysisManager::kOutputContainer,"HLT-OFFLINE-GLOBAL-comparison_triggered.root");
da6744ee	230	mgr->ConnectInput(taskGLOBAL,0,mgr->GetCommonInputContainer());
4b621d9c	231	mgr->ConnectOutput(taskGLOBAL,1,coutputGLOBAL);
db575e59	232	}
da6744ee	233
	234	if(bPWG1){
	235	Bool_t hasMC=kFALSE;
	236	// -- Add Task for HLT and Offline
	237	AliPerformanceTask *HLTtpcQA = AddTaskPerformance(hasMC,kFALSE,kTRUE);
	238	AliPerformanceTask *tpcQA = AddTaskPerformance(hasMC,kFALSE);
	239	if(!HLTtpcQA \|\| !tpcQA) {
	240	Error("RunPerformanceTrain","AliPerformanceTask not created!");
	241	return;
	242	}
	243	}
cd785861	244	if(bD0){
	245	float cuts[7]={0.5,0.04,0.7,0.8,0.05,-0.00025,0.7};
	246	AliAnalysisTaskD0Trigger *taskD0 = new AliAnalysisTaskD0Trigger("offhlt_comparison_D0",cuts);
	247	mgr->AddTask(taskD0);
4b621d9c	248	AliAnalysisDataContainer *coutputD0 = mgr->CreateContainer("D0_histograms",TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-D0-comparison.root");
cd785861	249	mgr->ConnectInput(taskD0,0,mgr->GetCommonInputContainer());
4b621d9c	250	mgr->ConnectOutput(taskD0,1,coutputD0);
ff84363f	251	}
	252
	253	if(bCB){
	254	AliAnalysisTaskHLTCentralBarrel *taskCB = new AliAnalysisTaskHLTCentralBarrel("offhlt_comparison_CB");
	255	AliAnalysisDataContainer *coutputCB = mgr->CreateContainer("esd_thnsparse", TList::Class(), AliAnalysisManager::kOutputContainer, "HLT-OFFLINE-CentralBarrel-comparison.root");
	256	mgr->ConnectInput(taskCB,0,mgr->GetCommonInputContainer());
	257	mgr->ConnectOutput(taskCB,1,coutputCB);
cd785861	258	}
89ad7eb5	259
be1fa04a	260	if (!mgr->InitAnalysis()) return;
be1fa04a	261	mgr->PrintStatus();
ea46fe48	262	mgr->StartAnalysis("local",chain, nEvents);
be1fa04a	263
	264	timer.Stop();
	265	timer.Print();
	266	}
4b621d9c	267
	268	void compare_HLT_offline_local(){
	269	cout << " " << endl;
	270	cout << " Usage examples:" << endl;
	271	cout << " compare-HLT-offline-local.C'(file, taskOption, taskFolder, fUseHLTTrigger, nEvents)' 2>&1 \| tee log" << endl;
	272	cout << " compare-HLT-offline-local.C'(\"AliESDs.root\",\"global\")' 2>&1 \| tee log" << endl;
	273	cout << " compare-HLT-offline-local.C'(\"AliESDs.root\",\"global\",\"./\",kFALSE,nEvents)' 2>&1 \| tee log" << endl;
ff84363f	274	cout << " compare-HLT-offline-local.C'(\"AliESDs.root\",\"global phos cb D0\", \"./\", kTRUE, nEvents)' 2>&1 \| tee log" << endl;
4b621d9c	275	cout << " compare-HLT-offline-local.C'(\"files.txt\",\"all\")' 2>&1 \| tee log" << endl;
	276	cout << " compare-HLT-offline-local.C'(\"alien:///alice/data/2010/LHC10b/000115322/ESDs/pass1/10000115322040.20/AliESDs.root\",\"global\")' 2>&1 \| tee log" << endl;
	277	cout << " " << endl;
	278	}