Changeset 10088


Ignore:
Timestamp:
Jan 6, 2011, 1:59:49 PM (10 years ago)
Author:
tbretz
Message:
Added output of filename for outliers.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Cosy/tpoint/plot.C

    r9999 r10088  
    211211    for (n=0; dates[n]>0; n++);
    212212
    213     TH1F h0("Q98", "", n-1, dates);
    214     TH1F h1("Q95", "", n-1, dates);
    215     TH1F h2("Q68", "", n-1, dates);
    216     TH1F h3("Q50", "", n-1, dates);
     213    TH1F h0("Q50", "", n-1, dates);
     214    TH1F h1("Q68", "", n-1, dates);
     215    TH1F h2("Q95", "", n-1, dates);
     216    TH1F h3("Q98", "", n-1, dates);
    217217
    218218    Int_t first = 0;
    219219
    220     cout << "n=" << n << endl;
    221 
     220    cout << "Number of bins: " << n << endl;
     221    cout << "Start      Entries Median" << endl;
     222
     223    MTime date;
    222224    for (int i=0; dates[i+1]>0; i++)
    223225    {
     
    225227        Int_t bin1 = GetBin(g0, i+1);
    226228
    227         cout << i << " " << bin0 << " " << bin1 << endl;
     229        //cout << i << " " << bin0 << " " << bin1 << endl;
    228230        if (bin0<0)
    229231            continue;
     
    233235
    234236        TMath::Quantiles(bin1-bin0, 4, g2.GetY()+bin0, quant, prob, kFALSE, NULL, 7);
     237
     238        date.SetAxisTime(dates[i]);
     239        cout << date.GetStringFmt("%x", "de_DE") << "  " << Form("%4d", bin1-bin0) << "   " << Form("%4.1f", quant[0]*60) << endl;
    235240
    236241        h0.SetBinContent(i+1, quant[0]);
     
    264269    h2.DrawCopy();
    265270
    266     //.DrawCopy("same");
     271    //h3.DrawCopy("same");
    267272    h1.DrawCopy("same");
    268273    h0.DrawCopy("same");
     
    297302
    298303    DrawGrid(h2);
    299 
    300 
    301 }
     304}
Note: See TracChangeset for help on using the changeset viewer.