]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDCalibraExbAltFit.cxx
clean up the log summary creation
[u/mrichter/AliRoot.git] / TRD / AliTRDCalibraExbAltFit.cxx
index 250fc4f719bc30b19380b61840c9674db8cbd623..430f43b8196530a2bab70b255393918056454d1d 100644 (file)
@@ -417,6 +417,7 @@ void AliTRDCalibraExbAltFit::FillPEArray2()
       
       Int_t nEntries=0;
       TGraphErrors *gg=DrawMS(fitterhisto,nEntries);
+      if (!gg) continue;
       // printf("N: %i\n",gg->GetN());
       // printf("entries: %i\n",nEntries);
 
@@ -471,6 +472,11 @@ void AliTRDCalibraExbAltFit::FillPEArray2()
 //_________Helper function__________________________________________________
 TGraphErrors* AliTRDCalibraExbAltFit::DrawMS(const TH2 *const h2, Int_t &nEntries)
 {
+  //
+  // Debug function
+  //
+
+
   TF1 fg("fg", "gaus", -10., 30.);
   TGraphErrors *gp = new TGraphErrors();
 
@@ -492,7 +498,7 @@ TGraphErrors* AliTRDCalibraExbAltFit::DrawMS(const TH2 *const h2, Int_t &nEntrie
       //Warning("drawMS()", Form("reject x[%d]=%f on NDF=%d", jpt, x, fg.GetNDF()));
       continue;
     }
-    if(fg.GetParameter(1)+fg.GetParameter(2)/2>ay->GetXmax() || fg.GetParameter(1)-fg.GetParameter(2)/2<ay->GetXmin() || !fg.GetParameter(0)) continue;
+    if(((fg.GetParameter(1)+fg.GetParameter(2)/2)>ay->GetXmax()) || ((fg.GetParameter(1)-fg.GetParameter(2)/2)<ay->GetXmin()) || (TMath::Abs(fg.GetParameter(0))< 0.00001)) continue;
     gp->SetPoint(ig, ax->GetBinCenter(jpt), fg.GetParameter(1));
     gp->SetPointError(ig, 0, TMath::Sqrt(pow(fg.GetParError(1),2) + (1/pow(fg.GetParameter(0),2))));
     ig++;