Changeset 10088 for trunk/Cosy/tpoint

Timestamp:

01/06/11 13:59:49 (14 years ago)

Author:

tbretz

Message:

Added output of filename for outliers.

File:

• trunk/Cosy/tpoint/plot.C (modified) (5 diffs)

trunk/Cosy/tpoint/plot.C

@@ -211 +211 @@
     for (n=0; dates[n]>0; n++);
 
-    TH1F h0("Q98", "", n-1, dates);
-    TH1F h1("Q95", "", n-1, dates);
-    TH1F h2("Q68", "", n-1, dates);
-    TH1F h3("Q50", "", n-1, dates);
+    TH1F h0("Q50", "", n-1, dates);
+    TH1F h1("Q68", "", n-1, dates);
+    TH1F h2("Q95", "", n-1, dates);
+    TH1F h3("Q98", "", n-1, dates);
 
     Int_t first = 0;
 
-    cout << "n=" << n << endl;
-
+    cout << "Number of bins: " << n << endl;
+    cout << "Start      Entries Median" << endl;
+
+    MTime date;
     for (int i=0; dates[i+1]>0; i++)
     {
@@ -225 +227 @@
         Int_t bin1 = GetBin(g0, i+1);
 
-        cout << i << " " << bin0 << " " << bin1 << endl;
+        //cout << i << " " << bin0 << " " << bin1 << endl;
         if (bin0<0)
             continue;
@@ -233 +235 @@
 
         TMath::Quantiles(bin1-bin0, 4, g2.GetY()+bin0, quant, prob, kFALSE, NULL, 7);
+
+        date.SetAxisTime(dates[i]);
+        cout << date.GetStringFmt("%x", "de_DE") << "  " << Form("%4d", bin1-bin0) << "   " << Form("%4.1f", quant[0]*60) << endl;
 
         h0.SetBinContent(i+1, quant[0]);
@@ -264 +269 @@
     h2.DrawCopy();
 
-    //.DrawCopy("same");
+    //h3.DrawCopy("same");
     h1.DrawCopy("same");
     h0.DrawCopy("same");
@@ -297 +302 @@
 
     DrawGrid(h2);
-
-
-}
+}