New classes for event plane calculation (Johanna)

[u/mrichter/AliRoot.git] / T0 / AliT0CalibSeasonTimeShift.cxx

diff --git a/T0/AliT0CalibSeasonTimeShift.cxx b/T0/AliT0CalibSeasonTimeShift.cxx
index 7c688c396d4a67ca8fcec2556426069307ece933..828c4da30816ce082c97bc10afe624e575e86718 100644 (file)
--- a/T0/AliT0CalibSeasonTimeShift.cxx
+++ b/T0/AliT0CalibSeasonTimeShift.cxx
@@ -118,8 +118,15 @@ Bool_t AliT0CalibSeasonTimeShift::SetT0Par(const char* filePhys)
     {
       gFile->ls();
     TDirectory *dr = (TDirectory*) gFile->Get("T0Calib");
-    // dr->Dump();
+      if (!dr) {
+      AliWarning(Form("no Tzero Directory in file collected "));
+      return ok;
+    }
     TObjArray * TzeroObj = (TObjArray*) dr->Get("fTzeroObject");
+    if(!TzeroObj) { 
+      AliWarning(Form("no Tzero Object in file collected "));
+      return ok;
+    }
     TString histname[4]={"fTzeroORAplusORC", "fTzeroORA", "fTzeroORC",  "fResolution"};
     for (Int_t i=0; i<4; i++)
        {
@@ -127,12 +134,14 @@ Bool_t AliT0CalibSeasonTimeShift::SetT0Par(const char* filePhys)
          gFile->Get(histname[i].Data());
          if(!cfd) AliWarning(Form("no histograms collected for %s", histname[i].Data()));
          if(cfd) {
+           //      printf(" T0s: %i %s mean %f  rms %f  \n", i, histname[i].Data(),cfd->GetMean , cfd ->GetRMS() ); 
            GetMeanAndSigma(cfd, mean, sigma);
            if (sigma == 0 || sigma > 500) ok = false;
            else
              { 
-               fMeanPar[i] = mean;
+               fMeanPar[i] =  mean;
                fSigmaPar[i] = sigma;
+               printf(" T0s: %i %s %f %f  \n", i, histname[i].Data(), fMeanPar[i],fSigmaPar[i] ); 
                ok=true;
              }
          }