- Timestamp:
- 01/06/11 13:59:49 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosy/tpoint/plot.C
r9999 r10088 211 211 for (n=0; dates[n]>0; n++); 212 212 213 TH1F h0("Q 98", "", n-1, dates);214 TH1F h1("Q 95", "", n-1, dates);215 TH1F h2("Q 68", "", n-1, dates);216 TH1F h3("Q 50", "", 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); 217 217 218 218 Int_t first = 0; 219 219 220 cout << "n=" << n << endl; 221 220 cout << "Number of bins: " << n << endl; 221 cout << "Start Entries Median" << endl; 222 223 MTime date; 222 224 for (int i=0; dates[i+1]>0; i++) 223 225 { … … 225 227 Int_t bin1 = GetBin(g0, i+1); 226 228 227 cout << i << " " << bin0 << " " << bin1 << endl;229 //cout << i << " " << bin0 << " " << bin1 << endl; 228 230 if (bin0<0) 229 231 continue; … … 233 235 234 236 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; 235 240 236 241 h0.SetBinContent(i+1, quant[0]); … … 264 269 h2.DrawCopy(); 265 270 266 // .DrawCopy("same");271 //h3.DrawCopy("same"); 267 272 h1.DrawCopy("same"); 268 273 h0.DrawCopy("same"); … … 297 302 298 303 DrawGrid(h2); 299 300 301 } 304 }
Note:
See TracChangeset
for help on using the changeset viewer.