]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/qaRec/run.C
last modification of the task before going into production
[u/mrichter/AliRoot.git] / TRD / qaRec / run.C
index 012d63db40eeeec1255515922bd15bc0199a6674..cc57343a60e0f8639c4855ef4f639b84400cb7f5 100644 (file)
 // Steer TRD QA train for Reconstruction (Clusterizer, Tracking and PID).
 // 
 // Usage:
-//   run.C(tasks, files, entries)
+//   run.C(tasks, files)
 //   tasks : "ALL" or one/more of the following:
 //     "EFF"  : TRD Tracking Efficiency 
 //     "EFFC" : TRD Tracking Efficiency Combined (barrel + stand alone) - only in case of simulations
 //     "RES"  : TRD tracking Resolution
+//     "CLRES": clusters Resolution
+//     "CAL"  : TRD calibration
+//     "ALGN" : TRD alignment
 //     "PID"  : TRD PID - pion efficiency 
 //     "PIDR" : TRD PID - reference data
+//     "DET"  : Basic TRD Detector checks
+//     "NOFR" : Data set does not have AliESDfriends.root 
 //     "NOMC" : Data set does not have Monte Carlo Informations (real data), so all tasks which rely
 //              on MC information are switched off
-// 
+//
+// In compiled mode : 
+// Don't forget to load first the libraries
+// gSystem->Load("libMemStat.so")
+// gSystem->Load("libMemStatGui.so")
+// gSystem->Load("libANALYSIS.so")
+// gSystem->Load("libTRDqaRec.so")
+// gSystem->Load("libNetx.so") ;
+// gSystem->Load("libRAliEn.so");
+//
 // Authors:
 //   Alex Bercuci (A.Bercuci@gsi.de) 
 //   Markus Fasel (m.Fasel@gsi.de) 
 
-#define BIT(n)        (1 << (n))
-#define SETBIT(n,i)   ((n) |= BIT(i))
-#define TESTBIT(n,i)  ((Bool_t)(((n) & BIT(i)) != 0))
-#define CLEARBIT(n,i) ((n) &= ~BIT(i))
-
-const Int_t fknTasks = 3;
-Char_t *fTaskName[fknTasks] = {"Barrel Tracking Effiency", "Combined Tracking Efficiency", "Tracking Resolution"};
-enum AliTRDrecoTasks{
-  kTrackingEfficiency = 0
-  ,kTrackingCombinedEfficiency = 1
-  ,kTrackingResolution = 2
-};
-void run(const Char_t *files=0x0, Char_t *tasks="ALL", Int_t nmax=-1)
-{
+#ifndef __CINT__
+#include <Riostream.h>
+
+#include "TStopwatch.h"
+#include "TMemStat.h"
+#include "TMemStatViewerGUI.h"
+
+#include "TROOT.h"
+#include "TSystem.h"
+#include "TError.h"
+#include "TChain.h"
+#include "TGrid.h"
+#include "TAlienCollection.h"
+#include "TGridCollection.h"
+#include "TGridResult.h"
+#include "TGeoGlobalMagField.h"
+
+#include "AliMagF.h"
+#include "AliTracker.h"
+#include "AliLog.h"
+#include "AliCDBManager.h"
+#include "AliGeomManager.h"
+#include "AliAnalysisManager.h"
+#include "AliAnalysisDataContainer.h"
+#include "AliMCEventHandler.h"
+#include "AliESDInputHandler.h"
 
+#include "TRD/AliTRDtrackerV1.h"
+#include "TRD/AliTRDcalibDB.h"
+#include "TRD/qaRec/AliTRDtrackInfo/AliTRDeventInfo.h"
+#include "TRD/qaRec/AliTRDcheckESD.h"
+#include "TRD/qaRec/AliTRDtrackInfoGen.h"
+#include "TRD/qaRec/AliTRDtrackingEfficiency.h"
+#include "TRD/qaRec/AliTRDtrackingEfficiencyCombined.h"
+#include "TRD/qaRec/AliTRDtrackingResolution.h"
+#include "TRD/qaRec/AliTRDcalibration.h"
+#include "TRD/qaRec/AliTRDalignmentTask.h"
+#include "TRD/qaRec/AliTRDpidChecker.h"
+#include "TRD/qaRec/AliTRDpidRefMaker.h"
+#include "TRD/qaRec/AliTRDcheckDetector.h"
+#include "TRD/qaRec/AliTRDclusterResolution.h"
+#endif
+
+#include "run.h"
+
+Bool_t MEM = kFALSE;
+
+TChain* MakeChainLST(const char* filename = 0x0);
+TChain* MakeChainXML(const char* filename = 0x0);
+void run(Char_t *tasks="ALL", const Char_t *files=0x0)
+{
+  TMemStat *mem = 0x0;
+  if(MEM){ 
+    gSystem->Load("libMemStat.so");
+    gSystem->Load("libMemStatGui.so");
+    mem = new TMemStat("new, gnubuildin");
+    mem->AddStamp("Start");
+  }
 
   TStopwatch timer;
   timer.Start();
 
-  gSystem->Load("libANALYSIS.so");
-  gSystem->Load("libTRDqaRec.so");
+  if(gSystem->Load("libANALYSIS.so")<0) return;
+  if(gSystem->Load("libTRDqaRec.so")<0) return;
   
-  Int_t fSteerTask = 0; 
+  // DB INITIALIZATION
+  //TODO We should use the GRP if available similar to AliReconstruction::InitGRP()!
+  // initialize OCDB manager
+  AliCDBManager *cdbManager = AliCDBManager::Instance();
+  cdbManager->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+  cdbManager->SetRun(0);
+  cdbManager->SetCacheFlag(kFALSE);
+  // initialize magnetic field.
+  AliMagF *field = 0x0;
+  field = new AliMagF("Maps","Maps", 2, 1., 10., AliMagF::k5kG);
+  //field = new AliMagF("Maps","Maps", 2, 0., 10., AliMagF::k2kG);
+  TGeoGlobalMagField::Instance()->SetField(field);
+
+  // initialize TRD settings
+  AliTRDcalibDB *cal = AliTRDcalibDB::Instance();
+  AliTRDtrackerV1::SetNTimeBins(cal->GetNumberOfTimeBins());
+  AliGeomManager::LoadGeometry();
+
   Bool_t fHasMCdata = kTRUE;
-  TObjArray *task = TString(tasks).Tokenize(" ");
-  for(Int_t isel = 0; isel < task->GetEntriesFast(); isel++){
-    TString s = (dynamic_cast<TObjString *>(task->UncheckedAt(isel)))->String();
+  Bool_t fHasFriends = kTRUE;
+  TObjArray *tasksArray = TString(tasks).Tokenize(" ");
+
+  Int_t fSteerTask = 0;
+  for(Int_t isel = 0; isel < tasksArray->GetEntriesFast(); isel++){
+    TString s = (dynamic_cast<TObjString *>(tasksArray->UncheckedAt(isel)))->String();
     if(s.CompareTo("ALL") == 0){
-      for(Int_t itask = 0; itask < fknTasks; itask++) SETBIT(fSteerTask, itask);
-      continue;
-    } else if(s.CompareTo("EFF") == 0){
-      SETBIT(fSteerTask, kTrackingEfficiency);
-      continue;
-    } else if(s.CompareTo("EFFC") == 0){
-      SETBIT(fSteerTask, kTrackingCombinedEfficiency);
-      continue;
-    } else if(s.CompareTo("RES") == 0){
-      SETBIT(fSteerTask, kTrackingResolution);
-      continue;
-    } else if(s.CompareTo("NOMC") == 0){
-       CLEARBIT(fSteerTask, kTrackingEfficiency);
-       CLEARBIT(fSteerTask, kTrackingEfficiencyCombined);
-    } else{
-      Info("run.C", Form("Task %s not implemented (yet).", s.Data()));
+      for(Int_t itask = 0; itask < NQATASKS; itask++) SETBIT(fSteerTask, itask);
       continue;
+    } else if(s.CompareTo("NOFR") == 0){ 
+      fHasFriends = kFALSE;
+    } else if(s.CompareTo("NOMC") == 0){ 
+      fHasMCdata = kFALSE;
+    } else { 
+      Bool_t foundOpt = kFALSE;  
+      for(Int_t itask = 1; itask < NTRDTASKS; itask++){
+        if(s.CompareTo(fgkTRDtaskOpt[itask]) != 0) continue;
+        SETBIT(fSteerTask, itask); SETBIT(fSteerTask, 0);
+        foundOpt = kTRUE;
+        break;
+      }
+      if(!foundOpt) Info("run.C", Form("Task %s not implemented (yet).", s.Data()));
     }
   }
-  printf("\n\tRUNNING TRAIN FOR TASKS:\n");
-  for(itask = 0; itask < fknTasks; itask++){
-    if(TESTBIT(fSteerTask, itask)) printf("\t%s\n", fTaskName[itask]);
-  }
+  // extra rules for calibration tasks
+  if(TSTBIT(fSteerTask, kClErrParam)) SETBIT(fSteerTask, kResolution);
+  if(TSTBIT(fSteerTask, kPIDRefMaker)) SETBIT(fSteerTask, kPIDChecker);
+  if(TSTBIT(fSteerTask, kAlignment)) SETBIT(fSteerTask, kResolution);
+
+  // define task list pointers;
+  AliTRDrecoTask *taskPtr[NTRDTASKS], *task = 0x0;
+  memset(taskPtr, 0, NTRDTASKS*sizeof(AliAnalysisTask*));
 
   //____________________________________________//
-  gROOT->LoadMacro(Form("%s/TRD/qaRec/CreateESDChain.C", gSystem->ExpandPathName("$ALICE_ROOT")));
-  TChain *chain = CreateESDChain(files, nmax);
+  // DEFINE DATA CHAIN
+  TChain *chain = 0x0;
+  if(!files) chain = MakeChainLST();
+  else{
+    TString fn(files);
+    if(fn.EndsWith("xml")) chain = MakeChainXML(files);
+    else chain = MakeChainLST(files);
+  }
+  if(!chain) return;
+
   //chain->SetBranchStatus("*", 0);
   chain->SetBranchStatus("*FMD*",0);
   chain->SetBranchStatus("*Calo*",0);
@@ -83,67 +170,255 @@ void run(const Char_t *files=0x0, Char_t *tasks="ALL", Int_t nmax=-1)
 
   //____________________________________________
   // Make the analysis manager
-  AliAnalysisManager *mgr = new AliAnalysisManager("TRD Track Info Manager");
+  AliAnalysisManager *mgr = new AliAnalysisManager("TRD Reconstruction QA");
   //mgr->SetSpecialOutputLocation(source); // To Be Changed
-  AliVEventHandler* esdH = new AliESDInputHandler;
-  mgr->SetInputEventHandler(esdH);  
-  AliMCEventHandler *mc = new AliMCEventHandler();
-  mgr->SetMCtruthEventHandler(mc);
+  AliVEventHandler *esdH = 0x0, *mcH = 0x0;
+  mgr->SetInputEventHandler(esdH = new AliESDInputHandler);
+  if(fHasMCdata) mgr->SetMCtruthEventHandler(mcH = new AliMCEventHandler());
   //mgr->SetDebugLevel(10);
 
+  //____________________________________________
+  // TRD check ESD
+  AliTRDcheckESD *checkESD = new AliTRDcheckESD();
+  mgr->AddTask(checkESD);
+  checkESD->SetMC(fHasMCdata);
+  mgr->ConnectInput(checkESD, 0, mgr->GetCommonInputContainer());  mgr->ConnectOutput(checkESD, 0, mgr->CreateContainer(checkESD->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", checkESD->GetName())));
+
   //____________________________________________
   // TRD track summary generator
-  AliTRDtrackInfoGen *task1 = new AliTRDtrackInfoGen();
-  task1->SetDebugLevel(1);
-  if(!fHasMCdata) task1->SetHasMCdata(kFALSE);
-  mgr->AddTask(task1);
-  // Create containers for input/output
-  AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("data", TChain::Class(), AliAnalysisManager::kInputContainer);
-  AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("TrackInfoList", TObjArray::Class(), AliAnalysisManager::kExchangeContainer);
-  mgr->ConnectInput( task1, 0, cinput1);
-  mgr->ConnectOutput(task1, 0, coutput1);
+       if(TSTBIT(fSteerTask, kInfoGen)){
+    mgr->AddTask(task = new AliTRDtrackInfoGen());
+    taskPtr[(Int_t)kInfoGen] = task;
+    task->SetDebugLevel(0);
+    task->SetMCdata(fHasMCdata);
+    mgr->ConnectInput( task, 0, mgr->GetCommonInputContainer());
+    AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("trackInfo", TObjArray::Class(), AliAnalysisManager::kExchangeContainer);
+    AliAnalysisDataContainer *coutput1a = mgr->CreateContainer("eventInfo", AliTRDeventInfo::Class(), AliAnalysisManager::kExchangeContainer);
+    mgr->ConnectOutput(task, 0, coutput1);
+    mgr->ConnectOutput(task, 1, coutput1a);
+  }
+
+  //____________________________________________
+  // TRD detector checker
+       if(TSTBIT(fSteerTask, kCheckDetector)){
+    mgr->AddTask(task = new AliTRDcheckDetector());
+    taskPtr[(Int_t)kCheckDetector] = task;
+    task->SetDebugLevel(0);
+    task->SetMCdata(fHasMCdata);
+    
+    // Create containers for input/output
+    mgr->ConnectInput( task, 0, coutput1);
+    mgr->ConnectInput( task, 1, coutput1a);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+  }
 
   //____________________________________________
   // TRD barrel tracking efficiency
-  if(TESTBIT(fSteerTask, kTrackingEfficiency)){
-    AliTRDtrackingEfficiency *task2 = new AliTRDtrackingEfficiency();
-    task2->SetDebugLevel(1);
-    mgr->AddTask(task2);
+  if(fHasMCdata && TSTBIT(fSteerTask, kTrackingEff)){
+    mgr->AddTask(task = new AliTRDtrackingEfficiency());
+    taskPtr[(Int_t)kTrackingEff] = task;
+    task->SetDebugLevel(0);
+
     //Create containers for input/output
-    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("Efficiency", TList::Class(), AliAnalysisManager::kOutputContainer, "TRD.TrackingEfficiency.root");
-    mgr->ConnectInput( task2, 0, coutput1);
-    mgr->ConnectOutput(task2, 0, coutput2);
+    mgr->ConnectInput( task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
   }
 
   //____________________________________________
   // TRD combined tracking efficiency
-  if(TESTBIT(fSteerTask, kTrackingCombinedEfficiency)){
-    AliTRDtrackingEfficiencyCombined *task3 = new AliTRDtrackingEfficiencyCombined();
-    task3->SetDebugLevel(0);
-    mgr->AddTask(task3);
+  if(fHasMCdata && TSTBIT(fSteerTask, kTrackingEffMC)){
+    mgr->AddTask(task = new AliTRDtrackingEfficiencyCombined());
+    taskPtr[(Int_t)kTrackingEffMC] = task;
+    task->SetDebugLevel(0);
+
     // Create containers for input/output
-    AliAnalysisDataContainer *coutput3 = mgr->CreateContainer("Efficiency2", TObjArray::Class(), AliAnalysisManager::kOutputContainer, "TRD.TrackingEfficiencyCombined.root");
-    mgr->ConnectInput( task3, 0, coutput1);
-    mgr->ConnectOutput(task3, 0, coutput3);
+    mgr->ConnectInput( task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
   }
 
   //____________________________________________
   // TRD tracking resolution
-  if(TESTBIT(fSteerTask, kTrackingResolution)){
-    AliTRDtrackingResolution *task4 = new AliTRDtrackingResolution();
-    if(!fHasMCdata) task4->SetHasMCdata(kFALSE);
-    task4->SetDebugLevel(1);
-    mgr->AddTask(task4);
+  if(TSTBIT(fSteerTask, kResolution)){
+    mgr->AddTask(task = new AliTRDtrackingResolution());
+    taskPtr[(Int_t)kResolution] = task;
+    task->SetMCdata(fHasMCdata);
+    task->SetPostProcess(kFALSE);
+    task->SetDebugLevel(0);
+    
+    // Create containers for input/output
+    mgr->ConnectInput( task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+
+    // Create output containers for calibration tasks
+    const Int_t nc = 4;
+    const Char_t *cn[nc] = {"Cl", "Trklt", "MC_Cl", "MC_Trklt"}; 
+    AliAnalysisDataContainer *co[nc]; 
+    for(Int_t ic = 0; ic<nc; ic++){
+      co[ic] = mgr->CreateContainer(Form("%s%s", task->GetName(), cn[ic]), TObjArray::Class(), AliAnalysisManager::kExchangeContainer);
+      mgr->ConnectOutput(task, 1+ic, co[ic]);
+    }
+    
+    // test reconstruction calibration plugin
+    if(TSTBIT(fSteerTask, kClErrParam)){
+      mgr->AddTask(task = new AliTRDclusterResolution());
+      taskPtr[(Int_t)kClErrParam] = task;
+      ((AliTRDclusterResolution*)task)->SetExB();
+      mgr->ConnectInput(task, 0, co[0]);
+      mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+  
+      mgr->AddTask(task = new AliTRDclusterResolution("ClErrParamMC"));
+      taskPtr[(Int_t)kClErrParam+1] = task;
+      ((AliTRDclusterResolution*)task)->SetExB();
+      mgr->ConnectInput(task, 0, co[2]);
+      mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+    }
+  }
+
+  //____________________________________________
+  // TRD calibration
+  if(TSTBIT(fSteerTask, kCalibration)){
+    mgr->AddTask(task = new AliTRDcalibration());
+    taskPtr[(Int_t)kCalibration] = task;
+    ((AliTRDcalibration*)task)->SetLow(0);
+    ((AliTRDcalibration*)task)->SetHigh(30);
+    ((AliTRDcalibration*)task)->SetFillZero(kFALSE);
+    task->SetDebugLevel(0);
+
     // Create containers for input/output
-    AliAnalysisDataContainer *coutput4 = mgr->CreateContainer("Resolution", TList::Class(), AliAnalysisManager::kOutputContainer, "TRD.TrackingResolution.root");
-    mgr->ConnectInput( task4, 0, coutput1);
-    mgr->ConnectOutput(task4, 0, coutput4);
+    mgr->ConnectInput(task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
   }
+  
+  //____________________________________________
+  // TRD alignment
+  if(TSTBIT(fSteerTask, kAlignment)){
+    mgr->AddTask(task = new AliTRDalignmentTask());
+    taskPtr[(Int_t)kAlignment] = task;
+    task->SetDebugLevel(0);
+
+    // Create containers for input/output
+    mgr->ConnectInput(task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(Form("h%s", task->GetName()), TObjArray::Class(), AliAnalysisManager::kExchangeContainer));
+
+    mgr->ConnectOutput(task, 1, mgr->CreateContainer(task->GetName(), TTree::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+  }
+  
+  //____________________________________________
+  // TRD PID
+  if(TSTBIT(fSteerTask, kPIDChecker)){
+    mgr->AddTask(task = new AliTRDpidChecker());
+    taskPtr[(Int_t)kPIDChecker] = task;
+    task->SetDebugLevel(0);
+    task->SetMCdata(fHasMCdata);
+    
+    // Create containers for input/output
+    mgr->ConnectInput( task, 0, coutput1);
+    mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+
+    //____________________________________________
+    // TRD pid reference 
+    if(TSTBIT(fSteerTask, kPIDRefMaker)){
+      mgr->AddTask(task = new AliTRDpidRefMaker());
+      taskPtr[(Int_t)kPIDRefMaker] = task;
+      task->SetDebugLevel(0);
+      task->SetMCdata(fHasMCdata);
+      
+      // Create containers for input/output
+      mgr->ConnectInput( task, 0, coutput1);
+      mgr->ConnectOutput(task, 0, mgr->CreateContainer(task->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%s.root", task->GetName())));
+      mgr->ConnectOutput(task, 1, mgr->CreateContainer(Form("%sNN", task->GetName()), TTree::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%sNN.root", task->GetName())));
+      mgr->ConnectOutput(task, 2, mgr->CreateContainer(Form("%sLQ", task->GetName()), TTree::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Task%sLQ.root", task->GetName())));
+    }
+  }
+
 
   if (!mgr->InitAnalysis()) return;
-  mgr->PrintStatus();
-  mgr->StartAnalysis("local",chain);
+  printf("\n\tRUNNING TRAIN FOR TASKS:\n");
+  for(Int_t itask = 1; itask < NTRDTASKS; itask++){
+    if(TSTBIT(fSteerTask, itask)) printf("\t   %s [%s]\n",  taskPtr[itask]->GetName(), taskPtr[itask]->GetTitle());
+  }
+  printf("\n\n");
+  //mgr->PrintStatus();
+
+  mgr->StartAnalysis("local", chain);
 
   timer.Stop();
-  timer.Print();
+  timer.Print();  
+
+  cal->Terminate();
+  TGeoGlobalMagField::Instance()->SetField(NULL);
+  delete cdbManager;
+  for(Int_t it=NTRDTASKS; it--; ){ 
+    if(taskPtr[it]){ 
+      printf("Cleaning %s [%s] ...\n", fgkTRDtaskClassName[it], taskPtr[it]->GetTitle());
+      delete taskPtr[it];
+    }
+  }
+  delete checkESD;
+
+  if(mcH) delete mcH;
+  delete esdH;
+  delete mgr;
+  delete chain;
+  if(MEM) delete mem;
+  if(MEM) TMemStatViewerGUI::ShowGUI();
+}
+
+//____________________________________________
+TChain* MakeChainLST(const char* filename)
+{
+  // Create the chain
+  TChain* chain = new TChain("esdTree");
+
+  if(!filename){
+    chain->Add(Form("%s/AliESDs.root", gSystem->pwd()));
+    return chain;
+  }
+
+
+  // read ESD files from the input list.
+  ifstream in;
+  in.open(filename);
+  TString esdfile;
+  while(in.good()) {
+    in >> esdfile;
+    if (!esdfile.Contains("root")) continue; // protection
+    chain->Add(esdfile.Data());
+  }
+
+  in.close();
+
+  return chain;
+}
+
+//____________________________________________
+TChain* MakeChainXML(const char* xmlfile)
+{
+  if (!TFile::Open(xmlfile)) {
+    Error("MakeChainXML", Form("No file %s was found", xmlfile));
+    return 0x0;
+  }
+
+  if(gSystem->Load("libNetx.so")<0) return 0x0;
+  if(gSystem->Load("libRAliEn.so")<0) return 0x0;
+  TGrid::Connect("alien://") ;
+
+  TGridCollection *collection = (TGridCollection*) TAlienCollection::Open(xmlfile);
+  if (!collection) {
+    Error("MakeChainXML", Form("No collection found in %s", xmlfile)) ; 
+    return 0x0; 
+  }
+  //collection->CheckIfOnline();
+
+  TGridResult* result = collection->GetGridResult("",0 ,0);
+  if(!result->GetEntries()){
+    Error("MakeChainXML", Form("No entries found in %s", xmlfile)) ; 
+    return 0x0; 
+  }
+  // Makes the ESD chain 
+  TChain* chain = new TChain("esdTree");
+  for (Int_t idx = 0; idx < result->GetEntries(); idx++) {
+    chain->Add(result->GetKey(idx, "turl")); 
+  }
+  return chain;
 }