void AddTask_GammaConvV1AOD(TString mode){

	gSystem->Load("libCore.so");  
	gSystem->Load("libTree.so");
	gSystem->Load("libGeom.so");
	gSystem->Load("libVMC.so");
	gSystem->Load("libPhysics.so");
	gSystem->Load("libMinuit");
	gSystem->Load("libSTEERBase");
	gSystem->Load("libESD");
	gSystem->Load("libAOD");
	gSystem->Load("libANALYSIS");
	gSystem->Load("libANALYSISalice");  
	gSystem->Load("libEve.so");   
	gSystem->Load("libCDB.so");
	gSystem->Load("libProof.so");
	gSystem->Load("libRAWDatabase.so");
	gSystem->Load("libSTEER.so");
	gSystem->Load("libSTEERBase.so");
	gSystem->Load("libTENDER.so");
	gSystem->Load("libTRDbase.so");
	gSystem->Load("libVZERObase.so");
	gSystem->Load("libVZEROrec.so");
	gSystem->Load("libTENDERSupplies.so");
   	gSystem->Load("libPWGflowBase.so");
	gSystem->Load("libPWGflowTasks.so");
	gSystem->Load("libPWGGAGammaConv.so");

   Bool_t isMC	= kTRUE;
   Bool_t readTR = kTRUE;

//   TString Pattern="*/AliESDs.root";
   TString Pattern="AOD/*/AliAOD.root";
   TString Prefix ="";



     Int_t run_numbers[100] = {137161, 137162, 137165, 137230, 137231, 137232, 137235, 137236, 137243, 137366, 137430, 137431, 137432, 137434, 137439, 137440, 137441, 137443, 137530, 137531, 137539, 137541, 137544, 137546, 137549, 137595, 137608, 137638, 137639, 137685, 137686, 137691, 137692, 137693, 137704, 137718, 137722, 137724, 137751, 137752, 137844, 137848, 138190, 138192, 138197, 138201, 138225, 138275, 138364, 138396, 138438, 138439, 138442, 138469, 138534, 138578, 138579, 138582, 138583, 138621,  138638, 138652, 138653, 138662, 138666, 138730, 138732, 138837, 138870, 138871, 138872, 139028, 139029, 139036, 139037, 139038, 139104, 139105, 139107, 139173, 139308, 139309, 139310, 139311, 139314, 139328, 139329, 139360, 139437, 139438, 139439, 139465, 139503, 139504, 139505, 139507, 139510};
  // run has no AODs's: 138624 (70440 events), 
  //Int_t  run_numbers[100] = {138638, 138652, 138653, 138662, 138666, 138730, 138732, 138837, 138870, 138871, 138872, 139028, 139029, 139036, 139037, 139038, 139104, 139105, 139107, 139173, 139308, 139309, 139310, 139311, 139314, 139328, 139329, 139360, 139437, 139438, 139439, 139465, 139503, 139504, 139505, 139507, 139510};
   //Int_t run_numbers[10] = {137161, 137162, 137165, 137230, 137231, 137232, 137235, 137236, 137243, 137366};
   //Int_t run_numbers[10] = {137430, 137431, 137432, 137434, 137439, 137440, 137441, 137443, 137530, 137531};
   //Int_t run_numbers[10] = {137539, 137541, 137544, 137546, 137549, 137595, 137608, 137638, 137639, 137685};
   //Int_t run_numbers[10] = {137686, 137691, 137692, 137693, 137704, 137718, 137722, 137724, 137751, 137752};
   //Int_t run_numbers[10] = {137844, 137848, 138190, 138192, 138197, 138201, 138225, 138275, 138364, 138396};
   //Int_t run_numbers[10] = {138438, 138439, 138442, 138469, 138534, 138578, 138579, 138582, 138583, 138621};
   //Int_t run_numbers[10] = {138638, 138652, 138653, 138662, 138666, 138730, 138732, 138837, 138870};
   //Int_t run_numbers[10] = {138871, 138872, 139028, 139029, 139036, 139037, 139038, 139104, 139105, 139107};
   //Int_t run_numbers[10] = {139173, 139308, 139309, 139310, 139311, 139314, 139328, 139329, 139360, 139437};
   //Int_t run_numbers[10] = {139438, 139439, 139465, 139503, 139504, 139505, 139507, 139510};
   gSystem->Setenv("alien_CLOSE_SE","ALICE::Grenoble::SE");

   AliAnalysisManager *mgr = new AliAnalysisManager("Analysis");
   if (!mgr) {
      Error("AddTask_GammaConvV1", "No analysis manager found.");
      return 0;
   }

   gROOT->LoadMacro("CreateAlienHandlerAOD.C");
   AliAnalysisAlien *alienHandler = CreateAlienHandler(mode);
   if (!alienHandler) return;

   TString WorkingDir = "MCPbPb_LHC13d2_v3";
   TString data = "/alice/sim/2013/LHC13d2/";
   
   mgr->SetGridHandler(alienHandler);
   alienHandler->SetGridWorkingDir(WorkingDir.Data());
   alienHandler->SetGridOutputDir(Form("OUTPUT_%s",WorkingDir.Data()));

   alienHandler->SetGridDataDir(data);
   alienHandler->SetDataPattern(Pattern.Data());
   alienHandler->SetRunPrefix(Prefix.Data());
   //alienHandler->SetRunRange(137161,139510);
   for (Int_t i=0; i<100; i++) {
      if (run_numbers[i]==0) break;
      alienHandler->AddRunNumber(run_numbers[i]);
   }
      
   AliAODInputHandler * handler;
   handler = new AliAODInputHandler(); mgr->SetInputEventHandler((AliAODInputHandler*)handler);


   //========= Add PID Reponse to ANALYSIS manager ====
   if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
      gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
      AddTaskPIDResponse(isMC);
   }//========= Add PID Reponse to ANALYSIS manager ====
      
   gROOT->LoadMacro("$ALICE_ROOT/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C");
   
   cout<<"Connecting to Alien..."<<endl;
   TGrid::Connect("alien://");
   cout<<"==============================="<<endl;
   cout<<"Successfully connected to Alien"<<endl;
   cout<<"==============================="<<endl;

   AliAnalysisTask *taskA = AddTask_GammaConvV1_PbPb(1,  //change different set of cuts
                                                     isMC, //run MC 
                                                     kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
                                                     kTRUE, // enable additional Photon QA task
                                                     "alien:///alice/cern.ch/user/f/fbock/MCSpectraInput.root", // path to file for weigting input
                                                     kTRUE,  //enable Weighting
                                                     "1000000060084000001500000");

   AliAnalysisTask *taskB = AddTask_GammaConvV1_PbPb(2,  //change different set of cuts
                                                     isMC, //run MC 
                                                     kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
                                                     kTRUE, // enable additional Photon QA task
                                                     "alien:///alice/cern.ch/user/f/fbock/MCSpectraInput.root", // path to file for weigting input
                                                     kTRUE,  //enable Weighting
                                                     "1000000060084000001500000");

   AliAnalysisTask *taskC = AddTask_GammaConvV1_PbPb(5,  //change different set of cuts
                                                     isMC, //run MC 
                                                     kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
                                                     kTRUE, // enable additional Photon QA task
                                                     "alien:///alice/cern.ch/user/f/fbock/MCSpectraInput.root", // path to file for weigting input
                                                     kTRUE,  //enable Weighting
                                                     "1000000060084000001500000");

   AliAnalysisTask *taskD = AddTask_GammaConvV1_PbPb(6,  //change different set of cuts
                                                     isMC, //run MC 
                                                     kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
                                                     kTRUE, // enable additional Photon QA task
                                                     "alien:///alice/cern.ch/user/f/fbock/MCSpectraInput.root", // path to file for weigting input
                                                     kTRUE,  //enable Weighting
                                                     "1000000060084000001500000");
   

   if (!mgr->InitAnalysis())return;
   mgr->StartAnalysis("grid");
}
Int_t CheckLoadLibrary(const char* library)                                                                                                        
{                                                                                                                                                  
   // checks if a library is already loaded, if not loads the library                                                                               
                                                                                                                                                   
   if (strlen(gSystem->GetLibraries(Form("%s.so", library), "", kFALSE)) > 0)                                                                       
      return 1;                                                                                                                                      
                                                                                                                                                   
   return gSystem->Load(library);                                                                                                                   
}