]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ANALYSIS/AliAnalysisAlien.cxx
remove left over print statement
[u/mrichter/AliRoot.git] / ANALYSIS / AliAnalysisAlien.cxx
index bae04c34993ec1c6af653f4d11df96f0dba6d6a4..6003707254862cda024d71cbc3ccf4e3ea469e31 100644 (file)
@@ -27,6 +27,8 @@
 #include "TROOT.h"
 #include "TSystem.h"
 #include "TFile.h"
+#include "TFileCollection.h"
+#include "TChain.h"
 #include "TObjString.h"
 #include "TObjArray.h"
 #include "TGrid.h"
@@ -63,6 +65,9 @@ AliAnalysisAlien::AliAnalysisAlien()
                   fFastReadOption(0),
                   fOverwriteMode(1),
                   fNreplicas(2),
+                  fNproofWorkers(0),
+                  fNproofWorkersPerSlave(0),
+                  fProofReset(0),
                   fRunNumbers(),
                   fExecutable(),
                   fExecutableCommand(),
@@ -94,6 +99,11 @@ AliAnalysisAlien::AliAnalysisAlien()
                   fJobTag(),
                   fOutputSingle(),
                   fRunPrefix(),
+                  fProofCluster(),
+                  fProofDataSet(),
+                  fFileForTestMode(),
+                  fRootVersionForProof(),
+                  fAliRootMode(),
                   fInputFiles(0),
                   fPackages(0)
 {
@@ -121,6 +131,9 @@ AliAnalysisAlien::AliAnalysisAlien(const char *name)
                   fFastReadOption(0),
                   fOverwriteMode(1),
                   fNreplicas(2),
+                  fNproofWorkers(0),
+                  fNproofWorkersPerSlave(0),
+                  fProofReset(0),
                   fRunNumbers(),
                   fExecutable(),
                   fExecutableCommand(),
@@ -152,6 +165,11 @@ AliAnalysisAlien::AliAnalysisAlien(const char *name)
                   fJobTag(),
                   fOutputSingle(),
                   fRunPrefix(),
+                  fProofCluster(),
+                  fProofDataSet(),
+                  fFileForTestMode(),
+                  fRootVersionForProof(),
+                  fAliRootMode(),
                   fInputFiles(0),
                   fPackages(0)
 {
@@ -179,6 +197,9 @@ AliAnalysisAlien::AliAnalysisAlien(const AliAnalysisAlien& other)
                   fFastReadOption(other.fFastReadOption),
                   fOverwriteMode(other.fOverwriteMode),
                   fNreplicas(other.fNreplicas),
+                  fNproofWorkers(other.fNproofWorkers),
+                  fNproofWorkersPerSlave(other.fNproofWorkersPerSlave),
+                  fProofReset(other.fProofReset),
                   fRunNumbers(other.fRunNumbers),
                   fExecutable(other.fExecutable),
                   fExecutableCommand(other.fExecutableCommand),
@@ -210,6 +231,11 @@ AliAnalysisAlien::AliAnalysisAlien(const AliAnalysisAlien& other)
                   fJobTag(other.fJobTag),
                   fOutputSingle(other.fOutputSingle),
                   fRunPrefix(other.fRunPrefix),
+                  fProofCluster(other.fProofCluster),
+                  fProofDataSet(other.fProofDataSet),
+                  fFileForTestMode(other.fFileForTestMode),
+                  fRootVersionForProof(other.fRootVersionForProof),
+                  fAliRootMode(other.fAliRootMode),
                   fInputFiles(0),
                   fPackages(0)
 {
@@ -267,6 +293,9 @@ AliAnalysisAlien &AliAnalysisAlien::operator=(const AliAnalysisAlien& other)
       fFastReadOption          = other.fFastReadOption;
       fOverwriteMode           = other.fOverwriteMode;
       fNreplicas               = other.fNreplicas;
+      fNproofWorkers           = other.fNproofWorkers;
+      fNproofWorkersPerSlave   = other.fNproofWorkersPerSlave;
+      fProofReset              = other.fProofReset;
       fRunNumbers              = other.fRunNumbers;
       fExecutable              = other.fExecutable;
       fExecutableCommand       = other.fExecutableCommand;
@@ -298,6 +327,11 @@ AliAnalysisAlien &AliAnalysisAlien::operator=(const AliAnalysisAlien& other)
       fJobTag                  = other.fJobTag;
       fOutputSingle            = other.fOutputSingle;
       fRunPrefix               = other.fRunPrefix;
+      fProofCluster            = other.fProofCluster;
+      fProofDataSet            = other.fProofDataSet;
+      fFileForTestMode         = other.fFileForTestMode;
+      fRootVersionForProof     = other.fRootVersionForProof;
+      fAliRootMode             = other.fAliRootMode;
       if (other.fInputFiles) {
          fInputFiles = new TObjArray();
          TIter next(other.fInputFiles);
@@ -392,7 +426,7 @@ void AliAnalysisAlien::CdWork()
    }   
    // Work directory not existing - create it
    gGrid->Cd(homedir);
-   if (gGrid->Mkdir(workdir)) {
+   if (gGrid->Mkdir(workdir, "-p")) {
       gGrid->Cd(fGridWorkingDir);
       Info("CreateJDL", "\n#####   Created alien working directory %s", fGridWorkingDir.Data());
    } else {
@@ -460,6 +494,8 @@ Bool_t AliAnalysisAlien::CheckInputData()
          return kFALSE;
       }
       Info("CheckInputData", "Analysis will make a single xml for base data directory %s",fGridDataDir.Data());
+      if (fDataPattern.Contains("tag") && TestBit(AliAnalysisGrid::kTest))
+         TObject::SetBit(AliAnalysisGrid::kUseTags, kTRUE); // ADDED (fix problem in determining the tag usage in test mode) 
       return kTRUE;
    }
    // Process declared files
@@ -555,7 +591,7 @@ Bool_t AliAnalysisAlien::CheckInputData()
          msg += " type: xml_collection;";
          if (useTags) msg += " using_tags: Yes";
          else          msg += " using_tags: No";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          if (fNrunsPerMaster<2) {
             AddDataFile(Form("%s.xml", os->GetString().Data()));
          } else {
@@ -583,7 +619,7 @@ Bool_t AliAnalysisAlien::CheckInputData()
          msg += " type: xml_collection;";
          if (useTags) msg += " using_tags: Yes";
          else          msg += " using_tags: No";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          if (fNrunsPerMaster<2) {
             AddDataFile(Form("%s%d.xml",fRunPrefix.Data(),irun));
          } else {
@@ -917,7 +953,7 @@ Bool_t AliAnalysisAlien::CreateJDL()
       } else {
          if (!fGridOutputDir.Contains("/")) fGridOutputDir = Form("%s/%s", workdir.Data(), fGridOutputDir.Data());
          if (!DirectoryExists(fGridOutputDir)) {
-            if (gGrid->Mkdir(fGridOutputDir)) {
+            if (gGrid->Mkdir(fGridOutputDir,"-p")) {
                Info("CreateJDL", "\n#####   Created alien output directory %s", fGridOutputDir.Data());
             } else {
                Error("CreateJDL", "Could not create alien output directory %s", fGridOutputDir.Data());
@@ -1042,10 +1078,13 @@ Bool_t AliAnalysisAlien::CreateJDL()
          first = kTRUE;
          while ((os=(TObjString*)next2())) {
             if (!first) comment = NULL;
-            if (!os->GetString().Contains("@") && fCloseSE.Length())
-               fMergingJDL->AddToOutputArchive(Form("%s@%s",os->GetString().Data(), fCloseSE.Data()), comment); 
+            TString currentfile = os->GetString();
+            currentfile.ReplaceAll(".root", "*.root");
+            currentfile.ReplaceAll(".zip", "-Stage$2_$3.zip");
+            if (!currentfile.Contains("@") && fCloseSE.Length())
+               fMergingJDL->AddToOutputArchive(Form("%s@%s",currentfile.Data(), fCloseSE.Data()), comment); 
             else
-               fMergingJDL->AddToOutputArchive(os->GetString(), comment);
+               fMergingJDL->AddToOutputArchive(currentfile, comment);
             first = kFALSE;   
          }      
          delete arr;         
@@ -1093,6 +1132,7 @@ Bool_t AliAnalysisAlien::CreateJDL()
                fMergingJDL->AddToOutputSandbox(Form("%s@%s",os->GetString().Data(), fCloseSE.Data()), comment); 
             else
                fMergingJDL->AddToOutputSandbox(os->GetString(), comment);
+            first = kFALSE;   
          }   
          delete arr;
       }
@@ -1120,7 +1160,7 @@ Bool_t AliAnalysisAlien::CreateJDL()
       } else {
          if (!fGridOutputDir.Contains("/")) fGridOutputDir = Form("%s/%s", workdir.Data(), fGridOutputDir.Data());
          if (!fProductionMode && !DirectoryExists(fGridOutputDir)) {
-            if (gGrid->Mkdir(fGridOutputDir)) {
+            if (gGrid->Mkdir(fGridOutputDir,"-p")) {
                Info("CreateJDL", "\n#####   Created alien output directory %s", fGridOutputDir.Data());
             } else {
                Error("CreateJDL", "Could not create alien output directory %s", fGridOutputDir.Data());
@@ -1202,7 +1242,7 @@ Bool_t AliAnalysisAlien::WriteJDL(Bool_t copy)
          else fGridJDL->SetOutputDirectory(Form("%s/$2",fGridOutputDir.Data()), "Output directory");
       } else {   
          fGridJDL->SetOutputDirectory(Form("%s/$2/#alien_counter_03i#", fGridOutputDir.Data()), "Output directory");
-         fMergingJDL->SetOutputDirectory(Form("%s/$1", fGridOutputDir.Data()), "Output directory");
+         fMergingJDL->SetOutputDirectory("$1", "Output directory");
       }   
    }
       
@@ -1231,11 +1271,13 @@ Bool_t AliAnalysisAlien::WriteJDL(Bool_t copy)
    sjdl.Prepend(Form("Jobtag = {\n   \"comment:%s\"\n};\n", fJobTag.Data()));
    index = sjdl.Index("JDLVariables");
    if (index >= 0) sjdl.Insert(index, "\n# JDL variables\n");
+   sjdl += "Workdirectorysize = {\"5000MB\"};";
    sjdl1 += "JDLVariables = \n{\n   \"Packages\",\n   \"OutputDir\"\n};\n";
    sjdl1.Prepend(Form("Jobtag = {\n   \"comment:%s_Merging\"\n};\n", fJobTag.Data()));
    sjdl1.Prepend("# Generated merging jdl\n# $1 = full alien path to output directory to be merged\n# $2 = merging stage\n# $3 = merged chunk\n");
    index = sjdl1.Index("JDLVariables");
    if (index >= 0) sjdl1.Insert(index, "\n# JDL variables\n");
+   sjdl1 += "Workdirectorysize = {\"5000MB\"};";
    // Write jdl to file
    ofstream out;
    out.open(fJDLName.Data(), ios::out);
@@ -1350,20 +1392,20 @@ void AliAnalysisAlien::CheckDataType(const char *lfn, Bool_t &isCollection, Bool
       TGridResult *res = gGrid->Command(Form("listFilesFromCollection -z -v %s",lfn), kFALSE);
       if (!res) {
          msg += " using_tags: No (unknown)";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          return;
       }   
       const char* typeStr = res->GetKey(0, "origLFN");
       if (!typeStr || !strlen(typeStr)) {
          msg += " using_tags: No (unknown)";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          return;
       }   
       TString file = typeStr;
       useTags = file.Contains(".tag");
       if (useTags) msg += " using_tags: Yes";
       else          msg += " using_tags: No";
-      Info("CheckDataType", msg.Data());
+      Info("CheckDataType", "%s", msg.Data());
       return;
    }
    TString slfn(lfn);
@@ -1375,13 +1417,13 @@ void AliAnalysisAlien::CheckDataType(const char *lfn, Bool_t &isCollection, Bool
       TGridCollection *coll = (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"alien://%s\",1);",lfn));
       if (!coll) {
          msg += " using_tags: No (unknown)";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          return;
       }   
       TMap *map = coll->Next();
       if (!map) {
          msg += " using_tags: No (unknown)";
-         Info("CheckDataType", msg.Data());
+         Info("CheckDataType", "%s", msg.Data());
          return;
       }   
       map = (TMap*)map->GetValue("");
@@ -1391,7 +1433,7 @@ void AliAnalysisAlien::CheckDataType(const char *lfn, Bool_t &isCollection, Bool
       delete coll;
       if (useTags) msg += " using_tags: Yes";
       else          msg += " using_tags: No";
-      Info("CheckDataType", msg.Data());
+      Info("CheckDataType", "%s", msg.Data());
       return;
    }
    useTags = slfn.Contains(".tag");
@@ -1399,7 +1441,7 @@ void AliAnalysisAlien::CheckDataType(const char *lfn, Bool_t &isCollection, Bool
    else                        msg += " type: unknown file;";
    if (useTags) msg += " using_tags: Yes";
    else          msg += " using_tags: No";
-   Info("CheckDataType", msg.Data());
+   Info("CheckDataType", "%s", msg.Data());
 }
 
 //______________________________________________________________________________
@@ -1423,6 +1465,50 @@ void AliAnalysisAlien::EnablePackage(const char *package)
    fPackages->Add(new TObjString(pkg));
 }      
 
+//______________________________________________________________________________
+TChain *AliAnalysisAlien::GetChainForTestMode(const char *treeName) const
+{
+// Make a tree from files having the location specified in fFileForTestMode. 
+// Inspired from JF's CreateESDChain.
+   if (fFileForTestMode.IsNull()) {
+      Error("GetChainForTestMode", "For proof test mode please use SetFileForTestMode() pointing to a file that contains data file locations.");
+      return NULL;
+   }
+   if (gSystem->AccessPathName(fFileForTestMode)) {
+      Error("GetChainForTestMode", "File not found: %s", fFileForTestMode.Data());
+      return NULL;
+   }   
+   // Open the file
+   ifstream in;
+   in.open(fFileForTestMode);
+   Int_t count = 0;
+    // Read the input list of files and add them to the chain
+    TString line;
+    TChain *chain = new TChain(treeName);
+    while (in.good())
+    {
+      in >> line;
+      if (line.IsNull()) continue;
+      if (count++ == fNtestFiles) break;
+      TString esdFile(line);
+      TFile *file = TFile::Open(esdFile);
+      if (file) {
+         if (!file->IsZombie()) chain->Add(esdFile);
+         file->Close();
+      } else {
+         Error("GetChainforTestMode", "Skipping un-openable file: %s", esdFile.Data());
+      }   
+    }
+    in.close();
+    if (!chain->GetListOfFiles()->GetEntries()) {
+       Error("GetChainForTestMode", "No file from %s could be opened", fFileForTestMode.Data());
+       delete chain;
+       return NULL;
+    }
+//    chain->ls();
+    return chain;
+}    
+
 //______________________________________________________________________________
 const char *AliAnalysisAlien::GetJobStatus(Int_t jobidstart, Int_t lastid, Int_t &nrunning, Int_t &nwaiting, Int_t &nerror, Int_t &ndone)
 {
@@ -1636,8 +1722,10 @@ Bool_t AliAnalysisAlien::CheckMergedFiles(const char *filename, const char *alie
       ::Error("GetNregisteredFiles", "You need to be connected to AliEn.");
       return kFALSE;
    }
-   printf("Checking directory <%s> for merged files <%s*> ...\n", aliendir, filename);
-   TString command = Form("find %s/ *%s", saliendir.Data(), filename);
+   sfilename = filename;
+   sfilename.ReplaceAll(".root", "*.root");
+   printf("Checking directory <%s> for merged files <%s> ...\n", aliendir, sfilename.Data());
+   TString command = Form("find %s/ *%s", saliendir.Data(), sfilename.Data());
    TGridResult *res = gGrid->Command(command);
    if (!res) {
       ::Error("GetNregisteredFiles","Error: No result for the find command\n");
@@ -1655,6 +1743,7 @@ Bool_t AliAnalysisAlien::CheckMergedFiles(const char *filename, const char *alie
       turl.ReplaceAll("alien://", "");
       turl.ReplaceAll(saliendir, "");
       sfilename = gSystem->BaseName(turl);
+      turl = turl.Strip(TString::kLeading, '/');
       // Now check to what the file corresponds to: 
       //    original output           - aliendir/%03d/filename
       //    merged file (which stage) - aliendir/filename-Stage%02d_%04d
@@ -1684,13 +1773,20 @@ Bool_t AliAnalysisAlien::CheckMergedFiles(const char *filename, const char *alie
       }
       if (doneFinal) {
          delete res;
+         printf("=> Removing files from previous stages...\n");
+         gGrid->Rm(Form("%s/*Stage*.root", aliendir));
+         for (i=1; i<stage; i++)
+            gGrid->Rm(Form("%s/*Stage%d*.zip", aliendir, i));
          return kTRUE;
       }               
    }
    delete res;
    // Compute number of jobs that were submitted for the current stage
    Int_t ntotstage = countOrig;
-   for (i=1; i<=stage; i++) ntotstage = (ntotstage/nperchunk)+1;
+   for (i=1; i<=stage; i++) {
+      if (ntotstage%nperchunk) ntotstage = (ntotstage/nperchunk)+1;
+      else                     ntotstage = (ntotstage/nperchunk);
+   }   
    // Now compare with the number of set bits in the chunksDone array
    Int_t nmissing = (stage>0)?(ntotstage - countStage):0;
    // Print the info
@@ -1698,7 +1794,6 @@ Bool_t AliAnalysisAlien::CheckMergedFiles(const char *filename, const char *alie
    if (stage==0) printf("*** No merging completed so far.\n");
    else          printf("*** Found %d out of %d files merged for stage %d\n", countStage, ntotstage, stage);
    if (nmissing) printf("*** Number of merged files missing for this stage: %d -> check merging job completion\n", nmissing);
-
    if (!submit) return doneFinal;
    // Sumbit merging jobs for all missing chunks for the current stage.
    TString query = Form("submit %s %s", jdl, aliendir);
@@ -1713,7 +1808,8 @@ Bool_t AliAnalysisAlien::CheckMergedFiles(const char *filename, const char *alie
    }
    // Submit next stage of merging
    if (stage==0) countStage = countOrig;
-   Int_t nchunks = (countStage/nperchunk)+1;
+   Int_t nchunks = (countStage/nperchunk);
+   if (countStage%nperchunk) nchunks += 1;
    for (i=0; i<nchunks; i++) {
       Int_t jobId = SubmitSingleJob(Form("%s %d %d", query.Data(), stage+1, i));
       if (!jobId) return kFALSE;
@@ -2062,14 +2158,194 @@ void AliAnalysisAlien::SetPreferedSE(const char */*se*/)
 Bool_t AliAnalysisAlien::StartAnalysis(Long64_t /*nentries*/, Long64_t /*firstEntry*/)
 {
 // Start remote grid analysis.
+   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+   Bool_t testMode = TestBit(AliAnalysisGrid::kTest);
+   if (!mgr || !mgr->IsInitialized()) {
+      Error("StartAnalysis", "You need an initialized analysis manager for this");
+      return kFALSE;
+   }
+   // Are we in PROOF mode ?
+   if (mgr->IsProofMode()) {
+      Info("StartAnalysis", "##### Starting PROOF analysis on via the plugin #####");
+      if (fProofCluster.IsNull()) {
+         Error("StartAnalysis", "You need to specify the proof cluster name via SetProofCluster");
+         return kFALSE;
+      }   
+      if (fProofDataSet.IsNull() && !testMode) {
+         Error("StartAnalysis", "You need to specify a dataset using SetProofDataSet()");
+         return kFALSE;
+      }   
+      // Set the needed environment
+      gEnv->SetValue("XSec.GSI.DelegProxy","2");
+      // Do we need to reset PROOF ? The success of the Reset operation cannot be checked
+      if (fProofReset && !testMode) {
+         if (fProofReset==1) {
+            Info("StartAnalysis", "Sending soft reset signal to proof cluster %s", fProofCluster.Data());
+            gROOT->ProcessLine(Form("TProof::Reset(\"%s\", kFALSE);", fProofCluster.Data()));
+         } else {         
+            Info("StartAnalysis", "Sending hard reset signal to proof cluster %s", fProofCluster.Data());
+            gROOT->ProcessLine(Form("TProof::Reset(\"%s\", kTRUE);", fProofCluster.Data()));
+         }
+         Info("StartAnalysis", "Stopping the analysis. Please use SetProofReset(0) to resume.");
+         return kFALSE;
+      }
+      // Do we need to change the ROOT version ? The success of this cannot be checked.
+      if (!fRootVersionForProof.IsNull() && !testMode) {
+         gROOT->ProcessLine(Form("TProof::Mgr(\"%s\")->SetROOTVersion(\"%s\");", 
+                            fProofCluster.Data(), fRootVersionForProof.Data()));
+      }
+      // Connect to PROOF and check the status
+      Long_t proof = 0;
+      TString sworkers;
+      if (fNproofWorkersPerSlave) sworkers = Form("workers=%dx", fNproofWorkersPerSlave);
+      else if (fNproofWorkers) sworkers = Form("workers=%d", fNproofWorkers);
+      if (!testMode) {
+         if (!sworkers.IsNull()) 
+            proof = gROOT->ProcessLine(Form("TProof::Open(\"%s\", \"%s\");", fProofCluster.Data(), sworkers.Data()));
+         else   
+            proof = gROOT->ProcessLine(Form("TProof::Open(\"%s\");", fProofCluster.Data()));
+      } else {
+         proof = gROOT->ProcessLine("TProof::Open(\"\");");
+         if (!proof) {
+            Error("StartAnalysis", "Could not start PROOF in test mode");
+            return kFALSE;
+         }   
+      }
+      if (!proof) {
+         Error("StartAnalysis", "Could not connect to PROOF cluster <%s>", fProofCluster.Data());
+         return kFALSE;
+      }   
+      if (fNproofWorkersPerSlave*fNproofWorkers > 0)
+         gROOT->ProcessLine(Form("gProof->SetParallel(%d);", fNproofWorkers));
+      // Is dataset existing ?
+      if (!testMode) {
+         TString dataset = fProofDataSet;
+         Int_t index = dataset.Index("#");
+         if (index>=0) dataset.Remove(index);
+//         if (!gROOT->ProcessLine(Form("gProof->ExistsDataSet(\"%s\");",fProofDataSet.Data()))) {
+//            Error("StartAnalysis", "Dataset %s not existing", fProofDataSet.Data());
+//            return kFALSE;
+//         }
+//         Info("StartAnalysis", "Dataset %s found", dataset.Data());
+      }
+      // Is ClearPackages() needed ?
+      if (TestSpecialBit(kClearPackages)) {
+         Info("StartAnalysis", "ClearPackages signal sent to PROOF. Use SetClearPackages(kFALSE) to reset this.");
+         gROOT->ProcessLine("gProof->ClearPackages();");
+      }
+      // Is a given aliroot mode requested ?
+      TList optionsList;
+      if (!fAliRootMode.IsNull()) {
+         TString alirootMode = fAliRootMode;
+         if (alirootMode == "default") alirootMode = "";
+         Info("StartAnalysis", "You are requesting AliRoot mode: %s", fAliRootMode.Data());
+         optionsList.SetOwner();
+         optionsList.Add(new TNamed("ALIROOT_MODE", alirootMode.Data()));
+         // Check the additional libs to be loaded
+         TString extraLibs;
+         if (!alirootMode.IsNull()) extraLibs = "ANALYSIS:ANALYSISalice";
+         // Parse the extra libs for .so
+         if (fAdditionalLibs.Length()) {
+            TObjArray *list = fAdditionalLibs.Tokenize(" ");
+            TIter next(list);
+            TObjString *str;
+            while((str=(TObjString*)next()) && str->GetString().Contains(".so")) {
+               TString stmp = str->GetName();
+               if (stmp.BeginsWith("lib")) stmp.Remove(0,3);
+               stmp.ReplaceAll(".so","");
+               if (!extraLibs.IsNull()) extraLibs += ":";
+               extraLibs += stmp;
+            }
+            if (list) delete list;            
+         }
+         if (!extraLibs.IsNull()) optionsList.Add(new TNamed("ALIROOT_EXTRA_LIBS",extraLibs.Data()));
+         // Check extra includes
+         if (!fIncludePath.IsNull()) {
+            TString includePath = fIncludePath;
+            includePath.ReplaceAll(" ",":");
+            includePath.Strip(TString::kTrailing, ':');
+            Info("StartAnalysis", "Adding extra includes: %s",includePath.Data()); 
+            optionsList.Add(new TNamed("ALIROOT_EXTRA_INCLUDES",includePath.Data()));
+         }
+         // Check if connection to grid is requested
+         if (TestSpecialBit(kProofConnectGrid)) 
+            optionsList.Add(new TNamed("ALIROOT_ENABLE_ALIEN", "1"));
+         // Enable AliRoot par
+         if (testMode) {
+         // Enable proof lite package
+            TString alirootLite = gSystem->ExpandPathName("$ALICE_ROOT/ANALYSIS/macros/AliRootProofLite.par");
+            for (Int_t i=0; i<optionsList.GetSize(); i++) {
+               TNamed *obj = (TNamed*)optionsList.At(i);
+               printf("%s  %s\n", obj->GetName(), obj->GetTitle());
+            }   
+            if (!gROOT->ProcessLine(Form("gProof->UploadPackage(\"%s\");",alirootLite.Data()))
+              && !gROOT->ProcessLine(Form("gProof->EnablePackage(\"%s\", (TList*)0x%lx);",alirootLite.Data(),(ULong_t)&optionsList))) {
+                  Info("StartAnalysis", "AliRootProofLite enabled");
+            } else {                      
+               Error("StartAnalysis", "There was an error trying to enable package AliRootProofLite.par");
+               return kFALSE;
+            }   
+         } else {
+            if (gROOT->ProcessLine(Form("gProof->EnablePackage(\"VO_ALICE@AliRoot::%s\", (TList*)0x%lx);", 
+                                   fAliROOTVersion.Data(), (ULong_t)&optionsList))) {
+               Error("StartAnalysis", "There was an error trying to enable package VO_ALICE@AliRoot::%s", fAliROOTVersion.Data());
+               return kFALSE;
+            }         
+         }
+      } else {
+         if (fAdditionalLibs.Contains(".so") && !testMode) {
+            Error("StartAnalysis", "You request additional libs to be loaded but did not enabled any AliRoot mode. Please refer to: \
+                   \n http://aaf.cern.ch/node/83 and use a parameter for SetAliRootMode()");
+            return kFALSE;       
+         }
+      }
+      // Enable par files if requested
+      if (fPackages && fPackages->GetEntries()) {
+         TIter next(fPackages);
+         TObject *package;
+         while ((package=next())) {
+            if (gROOT->ProcessLine(Form("gProof->UploadPackage(\"%s\");", package->GetName()))) {
+               if (gROOT->ProcessLine(Form("gProof->EnablePackage(\"%s\",kTRUE);", package->GetName()))) {
+                  Error("StartAnalysis", "There was an error trying to enable package %s", package->GetName());
+                  return kFALSE;
+               }
+            } else {
+               Error("StartAnalysis", "There was an error trying to upload package %s", package->GetName());
+               return kFALSE;
+            }
+         }
+      }
+      // Do we need to load analysis source files ?
+      // NOTE: don't load on client since this is anyway done by the user to attach his task.
+      if (fAnalysisSource.Length()) {
+         TObjArray *list = fAnalysisSource.Tokenize(" ");
+         TIter next(list);
+         TObjString *str;
+         while((str=(TObjString*)next())) {
+            gROOT->ProcessLine(Form("gProof->Load(\"%s+g\", kTRUE);", str->GetName()));
+         }   
+         if (list) delete list;
+      }
+      if (testMode) {
+      // Register dataset to proof lite.
+         if (fFileForTestMode.IsNull()) {
+            Error("GetChainForTestMode", "For proof test mode please use SetFileForTestMode() pointing to a file that contains data file locations.");
+            return kFALSE;
+         }
+         if (gSystem->AccessPathName(fFileForTestMode)) {
+            Error("GetChainForTestMode", "File not found: %s", fFileForTestMode.Data());
+            return kFALSE;
+         }   
+         TFileCollection *coll = new TFileCollection();
+         coll->AddFromFile(fFileForTestMode);
+         gROOT->ProcessLine(Form("gProof->RegisterDataSet(\"test_collection\", (TFileCollection*)0x%lx, \"OV\");", (ULong_t)coll));
+         gROOT->ProcessLine("gProof->ShowDataSets()");
+      }
+      return kTRUE;
+   }
    
    // Check if output files have to be taken from the analysis manager
    if (TestBit(AliAnalysisGrid::kDefaultOutputs)) {
-      AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
-      if (!mgr || !mgr->IsInitialized()) {
-         Error("StartAnalysis", "You need an initialized analysis manager for this");
-         return kFALSE;
-      }
       fOutputFiles = "";
       TIter next(mgr->GetOutputs());
       AliAnalysisDataContainer *output;
@@ -2145,7 +2421,7 @@ Bool_t AliAnalysisAlien::StartAnalysis(Long64_t /*nentries*/, Long64_t /*firstEn
    }   
    if (!CreateJDL()) return kFALSE;
    if (TestBit(AliAnalysisGrid::kOffline)) return kFALSE;
-   if (TestBit(AliAnalysisGrid::kTest)) {
+   if (testMode) {
       // Locally testing the analysis
       Info("StartAnalysis", "\n_______________________________________________________________________ \
       \n   Running analysis script in a daughter shell as on a worker node \
@@ -2254,9 +2530,17 @@ Bool_t AliAnalysisAlien::SubmitMerging()
          runOutDir = Form("%s/%03d",fGridOutputDir.Data(), i);
       }
       // Check now the number of merging stages.
-      TString outputFile = fOutputFiles;
-      Int_t index = outputFile.Index(",");
-      if (index>0) outputFile.Remove(index);
+      TObjArray *list = fOutputFiles.Tokenize(",");
+      TIter next(list);
+      TObjString *str;
+      TString outputFile;
+      while((str=(TObjString*)next())) {
+         outputFile = str->GetString();
+         Int_t index = outputFile.Index("@");
+         if (index > 0) outputFile.Remove(index);
+         if (!fMergeExcludes.Contains(outputFile)) break;
+      }
+      delete list;
       Bool_t done = CheckMergedFiles(outputFile, runOutDir, fMaxMergeFiles, kTRUE, mergeJDLName);
       if (!done) return kFALSE;
    }
@@ -2271,7 +2555,10 @@ Bool_t AliAnalysisAlien::SubmitMerging()
 //______________________________________________________________________________
 Bool_t AliAnalysisAlien::SubmitNext()
 {
-// Submit next bunch of master jobs if the queue is free.
+// Submit next bunch of master jobs if the queue is free. The first master job is
+// submitted right away, while the next will not be unless the previous was split.
+// The plugin will not submit new master jobs if there are more that 500 jobs in
+// waiting phase.
    static Bool_t iscalled = kFALSE;
    static Int_t firstmaster = 0;
    static Int_t lastmaster = 0;
@@ -2282,21 +2569,28 @@ Bool_t AliAnalysisAlien::SubmitNext()
    Int_t ntosubmit = 0;
    TGridResult *res;
    TString jobID = "";
-   if (!fNsubmitted) ntosubmit = 1;
-   else {
+   Int_t nmasterjobs = fInputFiles->GetEntries();
+   if (!fNsubmitted) {
+      ntosubmit = 1;
+      if (!IsUseSubmitPolicy()) {
+         if (ntosubmit>5)
+            Info("SubmitNext","### Warning submit policy not used ! Submitting too many jobs at a time may be prohibitted. \
+                \n### You can use SetUseSubmitPolicy() to enable if you have problems.");
+         ntosubmit = nmasterjobs;
+      }   
+   } else {
       TString status = GetJobStatus(firstmaster, lastmaster, nrunning, nwaiting, nerror, ndone);
       printf("=== master %d: %s\n", lastmaster, status.Data());
       // If last master not split, just return
       if (status != "SPLIT") {iscalled = kFALSE; return kTRUE;}
       // No more than 100 waiting jobs
-      if (nwaiting>100) {iscalled = kFALSE; return kTRUE;}
+      if (nwaiting>500) {iscalled = kFALSE; return kTRUE;}
       npermaster = (nrunning+nwaiting+nerror+ndone)/fNsubmitted;      
-      if (npermaster) ntosubmit = (100-nwaiting)/npermaster;
+      if (npermaster) ntosubmit = (500-nwaiting)/npermaster;
       if (!ntosubmit) ntosubmit = 1;
       printf("=== WAITING(%d) RUNNING(%d) DONE(%d) OTHER(%d) NperMaster=%d => to submit %d jobs\n", 
              nwaiting, nrunning, ndone, nerror, npermaster, ntosubmit);
    }
-   Int_t nmasterjobs = fInputFiles->GetEntries();
    for (Int_t i=0; i<ntosubmit; i++) {
       // Submit for a range of enumeration of runs.
       if (fNsubmitted>=nmasterjobs) {iscalled = kFALSE; return kTRUE;}
@@ -2424,6 +2718,9 @@ void AliAnalysisAlien::WriteAnalysisMacro()
       out << "// Automatically generated analysis steering macro executed in grid subjobs" << endl << endl;
       out << "   TStopwatch timer;" << endl;
       out << "   timer.Start();" << endl << endl;
+      // Change temp directory to current one
+      out << "// Set temporary merging directory to current one" << endl;
+      out << "   gSystem->Setenv(\"TMPDIR\", gSystem->pwd());" << endl << endl;   
       out << "// load base root libraries" << endl;
       out << "   gSystem->Load(\"libTree\");" << endl;
       out << "   gSystem->Load(\"libGeom\");" << endl;
@@ -2446,6 +2743,7 @@ void AliAnalysisAlien::WriteAnalysisMacro()
       if (fIncludePath.Length()) out << "   gSystem->AddIncludePath(\"" << fIncludePath.Data() << "\");" << endl;
       out << "   gSystem->AddIncludePath(\"-I$ALICE_ROOT/include\");" << endl << endl;
       out << "// Load analysis framework libraries" << endl;
+      TString setupPar = "AliAnalysisAlien::SetupPar";
       if (!fPackages) {
          out << "   gSystem->Load(\"libSTEERBase\");" << endl;
          out << "   gSystem->Load(\"libESD\");" << endl;
@@ -2457,7 +2755,6 @@ void AliAnalysisAlien::WriteAnalysisMacro()
          TIter next(fPackages);
          TObject *obj;
          TString pkgname;
-         TString setupPar = "AliAnalysisAlien::SetupPar";
          while ((obj=next())) {
             pkgname = obj->GetName();
             if (pkgname == "STEERBase" ||
@@ -2513,6 +2810,8 @@ void AliAnalysisAlien::WriteAnalysisMacro()
          while((str=(TObjString*)next())) {
             if (str->GetString().Contains(".so"))
                out << "   gSystem->Load(\"" << str->GetString().Data() << "\");" << endl;
+            if (str->GetString().Contains(".par"))
+               out << "   if (!" << setupPar << "(\"" << str->GetString() << "\")) return;" << endl;
          }
          if (list) delete list;
       }
@@ -2539,9 +2838,6 @@ void AliAnalysisAlien::WriteAnalysisMacro()
          out << "   gEnv->SetValue(\"XNet.ReconnectTimeout\",10);" << endl;
          out << "   gEnv->SetValue(\"XNet.FirstConnectMaxCnt\",1);" << endl << endl;
       }   
-      // Change temp directory to current one
-      out << "// Set temporary merging directory to current one" << endl;
-      out << "   gSystem->Setenv(\"TMPDIR\", gSystem->pwd());" << endl << endl;   
       out << "// connect to AliEn and make the chain" << endl;
       out << "   if (!TGrid::Connect(\"alien://\")) return;" << endl;
       if (IsUsingTags()) {
@@ -2615,7 +2911,7 @@ void AliAnalysisAlien::WriteAnalysisMacro()
             msg += "                      AliLHCTagCuts *lhcCuts,\n";
             msg += "                      AliDetectorTagCuts *detCuts,\n";
             msg += "                      AliEventTagCuts *evCuts)";
-            Info("WriteAnalysisMacro", msg.Data());
+            Info("WriteAnalysisMacro", "%s", msg.Data());
          }
       } 
       if (!IsUsingTags() || fFriendChainName!="") {
@@ -2874,7 +3170,7 @@ void AliAnalysisAlien::WriteMergingMacro()
       out << "   gSystem->Setenv(\"TMPDIR\", gSystem->pwd());" << endl << endl;   
       out << "// Connect to AliEn" << endl;
       out << "   if (!TGrid::Connect(\"alien://\")) return;" << endl;
-      out << "   Bool_t laststage = kTRUE;" << endl;
+      out << "   Bool_t laststage = kFALSE;" << endl;
       out << "   TString outputDir = dir;" << endl;  
       out << "   TString outputFiles = \"" << fOutputFiles << "\";" << endl;
       out << "   TString mergeExcludes = \"" << fMergeExcludes << "\";" << endl;
@@ -3075,7 +3371,7 @@ void AliAnalysisAlien::WriteExecutable()
       CdWork();
       TString workdir = gGrid->GetHomeDirectory();
       TString bindir = Form("%s/bin", workdir.Data());
-      if (!DirectoryExists(bindir)) gGrid->Mkdir(bindir);
+      if (!DirectoryExists(bindir)) gGrid->Mkdir(bindir,"-p");
       workdir += fGridWorkingDir;
       TString executable = Form("%s/bin/%s", gGrid->GetHomeDirectory(), fExecutable.Data());
       if (FileExists(executable)) gGrid->Rm(executable);
@@ -3133,7 +3429,7 @@ void AliAnalysisAlien::WriteMergeExecutable()
       CdWork();
       TString workdir = gGrid->GetHomeDirectory();
       TString bindir = Form("%s/bin", workdir.Data());
-      if (!DirectoryExists(bindir)) gGrid->Mkdir(bindir);
+      if (!DirectoryExists(bindir)) gGrid->Mkdir(bindir,"-p");
       workdir += fGridWorkingDir;
       TString executable = Form("%s/bin/%s", gGrid->GetHomeDirectory(), mergeExec.Data());
       if (FileExists(executable)) gGrid->Rm(executable);