void run(TString inFileName, TString plotName, TString sfile="bands", double mH = -1, int debug = 0){
_debug = debug;
double m2s = extractLimitAtQuantile(inFileName, plotName+"_-2sigma", 0.0275 );
double m2s_err = limitErr;
double m1s = extractLimitAtQuantile(inFileName, plotName+"_-1sigma", 0.16 );
double m1s_err = limitErr;
double med = extractLimitAtQuantile(inFileName, plotName+"_median", 0.5 );
double med_err = limitErr;
double p1s = extractLimitAtQuantile(inFileName, plotName+"_1sigma", 0.84 );
double p1s_err = limitErr;
double p2s = extractLimitAtQuantile(inFileName, plotName+"_2sigma", 0.975 );
double p2s_err = limitErr;
double dat = extractLimitAtQuantile(inFileName, plotName+"_observed", -1 );
double dat_err = limitErr;
cout<<"EXPECTED LIMIT BANDS from(obs, -2s,-1s,median,1s,2s) mass= "<<mH<<": ";
cout<<dat<<"+/-"<<dat_err<<", ";
cout<<m2s<<"+/-"<<m2s_err<<", "<<m1s<<"+/-"<<m1s_err<<", ";
cout<<med<<"+/-"<<med_err<<", ";
cout<<p1s<<"+/-"<<p1s_err<<", "<<p2s<<"+/-"<<p2s_err<<endl;
//cout<<"Observed data limit: "<<dat<<endl;
cout<<"Observed data limit: "<<dat<<" +/- "<<dat_err<<endl;
cout<<"expected median limit: "<<med<<" +/- "<<med_err<<endl;
SaveResults(sfile, mH, dat, dat_err, 0, 0, m2s, m1s, med, med_err, 0, p1s, p2s);
_m2s = m2s; _m1s=m1s;
}
BOOL APIENTRY ResultsDlgProc(
HWND hwnd,
UINT msg,
UINT wParam,
LONG lParam)
{
int i;
char szT[80];
BOOL fResults;
int aiT[1];
switch (msg) {
case WM_INITDIALOG:
aiT[0] = 190;
SendDlgItemMessage(hwnd, IDC_RESULTSLIST, LB_SETTABSTOPS, 1,
(LONG)(LPSTR)aiT);
fResults = FALSE;
for (i = 0; i < gcbmi; i++) {
if (gabmi[i].usPerIteration == 0)
continue;
wsprintf(szT, "%s\t%ld", (LPSTR)gabmi[i].pszBenchTitle,
gabmi[i].usPerIteration);
SendDlgItemMessage(hwnd, IDC_RESULTSLIST, LB_ADDSTRING, 0,
(LONG)(LPSTR)szT);
fResults = TRUE;
}
if (!fResults)
SendDlgItemMessage(hwnd, IDC_RESULTSLIST, LB_ADDSTRING, 0,
(LONG)(LPSTR)"No results have been generated yet!");
break;
case WM_COMMAND:
switch (LOWORD(wParam)) {
case IDOK:
EndDialog(hwnd, 0);
break;
case IDM_SAVERESULTS:
SaveResults(hwnd);
break;
default:
return FALSE;
}
break;
default:
return FALSE;
}
return TRUE;
}
int PrecForecast(stationData *stat, double *pobs, int *iteracao)
{
double P;
/* Calculos de previsao */
P = Georgakakos(*stat, *pobs, iteracao);
//pobs: regressao e etc
/* Gera arquivo de saida para o SWMM */
SaveResults(*stat, P);
(*iteracao)++;
return 0;
}
void CGUIControlProfiler::EndFrame(void)
{
m_iFrameCount++;
if (m_iFrameCount >= m_iMaxFrameCount)
{
const unsigned int dwSize = m_ItemHead.m_vecChildren.size();
for (unsigned int i=0; i<dwSize; ++i)
{
CGUIControlProfilerItem *p = m_ItemHead.m_vecChildren[i];
m_ItemHead.m_visTime += p->m_visTime;
m_ItemHead.m_renderTime += p->m_renderTime;
}
m_bIsRunning = false;
if (SaveResults())
m_ItemHead.Reset(this);
}
}
LONG APIENTRY FrameWndProc(
HWND hwnd,
UINT msg,
UINT wParam,
LONG lParam)
{
CLIENTCREATESTRUCT ccs;
char szT[80];
HMENU hmenu;
int i;
switch (msg) {
case WM_CREATE:
/*
* Add all the benchmark menu items to the benchmark menu.
*/
hmenu = GetSubMenu(GetMenu(hwnd), 2);
for (i = 0; i < gcbmi; i++) {
wsprintf(szT, "&%c. %s", 'A' + (i%26), (LPSTR)gabmi[i].pszBenchTitle);
if ((i > 0) && ((i%26) == 0))
AppendMenu(hmenu, MF_MENUBREAK, 200, NULL);
AppendMenu(hmenu, MF_STRING | MF_ENABLED, IDM_BENCHFIRST + i, szT);
gabmi[i].usPerIteration = 0;
}
/*
* Create the MDI client filling the client area
*/
ccs.hWindowMenu = NULL;
ccs.idFirstChild = 0;
ghwndMDIClient = CreateWindow("mdiclient",
NULL, WS_CHILD | WS_CLIPCHILDREN | WS_VSCROLL | WS_HSCROLL,
0, 0, 0, 0, hwnd, (HMENU)0xCAC, ghinst, (LPSTR)&ccs);
ShowWindow(ghwndMDIClient, SW_SHOW);
break;
case WM_COMMAND:
switch (LOWORD(wParam)) {
case IDM_ABOUT:
DialogBox(ghinst, (LPSTR)IDD_ABOUT, hwnd, (DLGPROC)AboutDlgProc);
break;
case IDM_VIEWRESULTS:
DialogBox(ghinst, (LPSTR)IDD_RESULTS, hwnd, (DLGPROC)ResultsDlgProc);
break;
case IDM_SAVERESULTS:
SaveResults(hwnd);
break;
case IDM_ALL:
gfAll = TRUE;
for (gibmi = 0; gibmi < gcbmi; gibmi++) {
gabmi[gibmi].pfnBenchmark(gabmi[gibmi].cIterations);
}
gfAll = FALSE;
DialogBox(ghinst, (LPSTR)IDD_RESULTS, hwnd, (DLGPROC)ResultsDlgProc);
break;
case IDM_PROFILEALL:
/*
* We don't want a million iterations of these guys when we're
* profiling because it takes too long. Each benchmark will
* iterate 5 times.
*/
gfAll = TRUE;
for (gibmi = 0; gibmi < gcbmi; gibmi++) {
gabmi[gibmi].pfnBenchmark(5);
}
gfAll = FALSE;
DialogBox(ghinst, (LPSTR)IDD_RESULTS, hwnd, (DLGPROC)ResultsDlgProc);
break;
default:
if (LOWORD(wParam) >= IDM_BENCHFIRST &&
LOWORD(wParam) < IDM_BENCHFIRST + gcbmi) {
gibmi = LOWORD(wParam) - IDM_BENCHFIRST;
gabmi[gibmi].pfnBenchmark(gabmi[gibmi].cIterations);
break;
}
return DefFrameProc(hwnd, ghwndMDIClient, WM_COMMAND, wParam,
lParam);
}
break;
case WM_CLOSE:
DestroyWindow(hwnd);
break;
case WM_DESTROY:
PostQuitMessage(0);
break;
default:
/*
* Use DefFrameProc() instead of DefWindowProc() since there are
* things that have to be handled differently because of MDI.
*/
return DefFrameProc(hwnd, ghwndMDIClient, msg, wParam, lParam);
}
return 0;
}
