PUBLIC float alifold(const char **strings, char *structure){
int length, energy, s, n_seq;
circ = 0;
length = (int) strlen(strings[0]);
if (length>init_length) init_alifold(length);
if ((P==NULL)||(fabs(P->temperature - temperature)>1e-6)) {
update_fold_params(); P = scale_parameters();
}
for (s=0; strings[s]!=NULL; s++);
n_seq = s;
alloc_sequence_arrays(strings, &S, &S5, &S3, &a2s, &Ss, circ);
make_pscores((const short **) S, strings, n_seq, structure);
energy = fill_arrays((const char **)strings);
backtrack((const char **)strings, 0);
parenthesis_structure(structure, length);
free_sequence_arrays(n_seq, &S, &S5, &S3, &a2s, &Ss);
if (backtrack_type=='C')
return (float) c[indx[length]+1]/(n_seq*100.);
else if (backtrack_type=='M')
return (float) fML[indx[length]+1]/(n_seq*100.);
else
return (float) f5[length]/(n_seq*100.);
}
float fold(const char *string, char *structure) {
int i, length, energy, bonus=0, bonus_cnt=0;
circ = 0;
length = (int) strlen(string);
if (length>init_length) initialize_fold(length);
if (fabs(P->temperature - temperature)>1e-6) update_fold_params();
encode_seq(string);
BP = (int *)space(sizeof(int)*(length+2));
make_ptypes(S, structure);
energy = fill_arrays(string);
backtrack(string, 0);
#ifdef PAREN
parenthesis_structure(structure, length);
#else
letter_structure(structure, length);
#endif
/* check constraints */
for(i=1;i<=length;i++) {
if((BP[i]<0)&&(BP[i]>-4)) {
bonus_cnt++;
if((BP[i]==-3)&&(structure[i-1]==')')) bonus++;
if((BP[i]==-2)&&(structure[i-1]=='(')) bonus++;
if((BP[i]==-1)&&(structure[i-1]!='.')) bonus++;
}
if(BP[i]>i) {
int l;
bonus_cnt++;
for(l=1; l<=base_pair[0].i; l++)
if((i==base_pair[l].i)&&(BP[i]==base_pair[l].j)) bonus++;
}
}
if (bonus_cnt>bonus) fprintf(stderr,"\ncould not enforce all constraints\n");
bonus*=BONUS;
free(S); free(S1); free(BP);
energy += bonus; /*remove bonus energies from result */
if (backtrack_type=='C')
return (float) c[indx[length]+1]/100.;
else if (backtrack_type=='M')
return (float) fML[indx[length]+1]/100.;
else
return (float) energy/100.;
}
PUBLIC float alifold(const char **strings, char *structure){
int length, energy, s, n_seq;
circular = 0;
length = (int) strlen(strings[0]);
#ifdef _OPENMP
/* always init everything since all global static variables are uninitialized when entering a thread */
init_alifold(length);
#else
if (length>init_length) init_alifold(length);
#endif
if (fabs(P->temperature - temperature)>1e-6) update_alifold_params();
for (s=0; strings[s]!=NULL; s++);
n_seq = s;
alloc_sequence_arrays(strings, &S, &S5, &S3, &a2s, &Ss, circular);
make_pscores((const short **) S, strings, n_seq, structure);
energy = fill_arrays((const char **)strings);
backtrack((const char **)strings, 0);
#ifdef PAREN
parenthesis_structure(structure, base_pair2, length);
#else
letter_structure(structure, base_pair2, length);
#endif
/*
* Backward compatibility:
* This block may be removed if deprecated functions
* relying on the global variable "base_pair" vanishs from within the package!
*/
base_pair = base_pair2;
/*
{
if(base_pair) free(base_pair);
base_pair = (bondT *)space(sizeof(bondT) * (1+length/2));
memcpy(base_pair, base_pair2, sizeof(bondT) * (1+length/2));
}
*/
free_sequence_arrays(n_seq, &S, &S5, &S3, &a2s, &Ss);
if (backtrack_type=='C')
return (float) c[indx[length]+1]/(n_seq*100.);
else if (backtrack_type=='M')
return (float) fML[indx[length]+1]/(n_seq*100.);
else
return (float) f5[length]/(n_seq*100.);
}
void generate_table1( const char* scan_results_file = "scan-results-19fb.txt" ) {
fill_arrays( scan_results_file ) ;
FILE* outfp ;
if ( (outfp=fopen( "tables.tex", "w" ))==NULL ) {
printf("\n\n *** can't open output file.\n\n") ; return ;
}
fprintf( outfp, "\\documentclass[11pt]{article}\n" ) ;
fprintf( outfp, "\\def\\vsmtvs {\\rule[-0.3cm]{0cm}{0.75cm}}\n" ) ;
fprintf( outfp, "\\begin{document}\n" ) ;
fprintf( outfp, "\n\n\n" ) ;
fprintf( outfp, "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n") ;
fprintf( outfp, " \\begin{tabular}{|l||c|c|c||c|}\n") ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " & HT2 & HT3 & HT4 & HT2-4 \\vsmtvs \\\\ \n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " Observed events & %.0f & %.0f & %.0f & %.0f \\vsmtvs \\\\ \n",
data_2b_obs[0], data_2b_obs[1], data_2b_obs[2], data_2b_obs[3] ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " Unbiased SM background predictions & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n",
ub_2b_val[0], ub_2b_p1s[0], ub_2b_m1s[0],
ub_2b_val[1], ub_2b_p1s[1], ub_2b_m1s[1],
ub_2b_val[2], ub_2b_p1s[2], ub_2b_m1s[2],
ub_2b_val[3], ub_2b_p1s[3], ub_2b_m1s[3]
) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " SM background from full fit & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n",
ff_2b_val[0], ff_2b_p1s[0], ff_2b_m1s[0],
ff_2b_val[1], ff_2b_p1s[1], ff_2b_m1s[1],
ff_2b_val[2], ff_2b_p1s[2], ff_2b_m1s[2],
ff_2b_val[3], ff_2b_p1s[3], ff_2b_m1s[3]
) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " MC value & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.0f \\pm %.0f$ \\vsmtvs \\\\ \n",
mc_2b_val[0], mc_2b_err[0],
mc_2b_val[1], mc_2b_err[1],
mc_2b_val[2], mc_2b_err[2],
mc_2b_val[3], mc_2b_err[3]
) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\end{tabular}\n" ) ;
fprintf( outfp, " \\label{tab:pred-and-observed-nb2}\n") ;
fprintf( outfp, "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n") ;
fprintf( outfp, "\n\n\n" ) ;
fprintf( outfp, " \\begin{tabular}{|c||c|c|c|c||c|}\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\multicolumn{6}{|c|}{ Observed number of events} \\\\ \n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " & HT1 & HT2 & HT3 & HT4 & HT1-4 \\\\ \n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2 & %.0f & %.0f & %.0f & %.0f & %.0f \\vsmtvs \\\\ \n",
data_3b_obs[0][0],
data_3b_obs[0][1],
data_3b_obs[0][2],
data_3b_obs[0][3],
data_3b_obs[0][4] ) ;
fprintf( outfp, " $N_b\\ge3$, MET3 & %.0f & %.0f & %.0f & %.0f & %.0f \\vsmtvs \\\\ \n",
data_3b_obs[1][0],
data_3b_obs[1][1],
data_3b_obs[1][2],
data_3b_obs[1][3],
data_3b_obs[1][4] ) ;
fprintf( outfp, " $N_b\\ge3$, MET4 & --- & %.0f & %.0f & %.0f & %.0f \\vsmtvs \\\\ \n",
data_3b_obs[2][1],
data_3b_obs[2][2],
data_3b_obs[2][3],
data_3b_obs[2][4] ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2-4 & %.0f & %.0f & %.0f & %.0f & %.0f \\vsmtvs \\\\ \n",
data_3b_obs[3][0],
data_3b_obs[3][1],
data_3b_obs[3][2],
data_3b_obs[3][3],
data_3b_obs[3][4] ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\multicolumn{6}{|c|}{ Unbiased background predictions, SM-only} \\\\ \n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " & HT1 & HT2 & HT3 & HT4 & HT1-4 \\\\ \n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2 & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n"
,ub_3b_val[0][0], ub_3b_p1s[0][0], ub_3b_m1s[0][0]
,ub_3b_val[0][1], ub_3b_p1s[0][1], ub_3b_m1s[0][1]
,ub_3b_val[0][2], ub_3b_p1s[0][2], ub_3b_m1s[0][2]
,ub_3b_val[0][3], ub_3b_p1s[0][3], ub_3b_m1s[0][3]
,ub_3b_val[0][4], ub_3b_p1s[0][4], ub_3b_m1s[0][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET3 & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ \\vsmtvs \\\\ \n"
,ub_3b_val[1][0], ub_3b_p1s[1][0], ub_3b_m1s[1][0]
,ub_3b_val[1][1], ub_3b_p1s[1][1], ub_3b_m1s[1][1]
,ub_3b_val[1][2], ub_3b_p1s[1][2], ub_3b_m1s[1][2]
,ub_3b_val[1][3], ub_3b_p1s[1][3], ub_3b_m1s[1][3]
,ub_3b_val[1][4], ub_3b_p1s[1][4], ub_3b_m1s[1][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET4 & --- & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ \\vsmtvs \\\\ \n"
,ub_3b_val[2][1], ub_3b_p1s[2][1], ub_3b_m1s[2][1]
,ub_3b_val[2][2], ub_3b_p1s[2][2], ub_3b_m1s[2][2]
,ub_3b_val[2][3], ub_3b_p1s[2][3], ub_3b_m1s[2][3]
,ub_3b_val[2][4], ub_3b_p1s[2][4], ub_3b_m1s[2][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2-4 & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n"
,ub_3b_val[3][0], ub_3b_p1s[3][0], ub_3b_m1s[3][0]
,ub_3b_val[3][1], ub_3b_p1s[3][1], ub_3b_m1s[3][1]
,ub_3b_val[3][2], ub_3b_p1s[3][2], ub_3b_m1s[3][2]
,ub_3b_val[3][3], ub_3b_p1s[3][3], ub_3b_m1s[3][3]
,ub_3b_val[3][4], ub_3b_p1s[3][4], ub_3b_m1s[3][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\multicolumn{6}{|c|}{ Full fit results, SM-only} \\\\ \n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " & HT1 & HT2 & HT3 & HT4 & HT1-4 \\\\ \n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2 & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n"
,ff_3b_val[0][0], ff_3b_p1s[0][0], ff_3b_m1s[0][0]
,ff_3b_val[0][1], ff_3b_p1s[0][1], ff_3b_m1s[0][1]
,ff_3b_val[0][2], ff_3b_p1s[0][2], ff_3b_m1s[0][2]
,ff_3b_val[0][3], ff_3b_p1s[0][3], ff_3b_m1s[0][3]
,ff_3b_val[0][4], ff_3b_p1s[0][4], ff_3b_m1s[0][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET3 & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ \\vsmtvs \\\\ \n"
,ff_3b_val[1][0], ff_3b_p1s[1][0], ff_3b_m1s[1][0]
,ff_3b_val[1][1], ff_3b_p1s[1][1], ff_3b_m1s[1][1]
,ff_3b_val[1][2], ff_3b_p1s[1][2], ff_3b_m1s[1][2]
,ff_3b_val[1][3], ff_3b_p1s[1][3], ff_3b_m1s[1][3]
,ff_3b_val[1][4], ff_3b_p1s[1][4], ff_3b_m1s[1][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET4 & --- & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ \\vsmtvs \\\\ \n"
,ff_3b_val[2][1], ff_3b_p1s[2][1], ff_3b_m1s[2][1]
,ff_3b_val[2][2], ff_3b_p1s[2][2], ff_3b_m1s[2][2]
,ff_3b_val[2][3], ff_3b_p1s[2][3], ff_3b_m1s[2][3]
,ff_3b_val[2][4], ff_3b_p1s[2][4], ff_3b_m1s[2][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2-4 & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.1f\\ ^{+%.1f}_{-%.1f}$ & $%.0f\\ ^{+%.0f}_{-%.0f}$ \\vsmtvs \\\\ \n"
,ff_3b_val[3][0], ff_3b_p1s[3][0], ff_3b_m1s[3][0]
,ff_3b_val[3][1], ff_3b_p1s[3][1], ff_3b_m1s[3][1]
,ff_3b_val[3][2], ff_3b_p1s[3][2], ff_3b_m1s[3][2]
,ff_3b_val[3][3], ff_3b_p1s[3][3], ff_3b_m1s[3][3]
,ff_3b_val[3][4], ff_3b_p1s[3][4], ff_3b_m1s[3][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\multicolumn{6}{|c|}{ Monte Carlo, SM-only} \\\\ \n" ) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " & HT1 & HT2 & HT3 & HT4 & HT1-4 \\\\ \n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2 & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ \\vsmtvs \\\\ \n"
,mc_3b_val[0][0], mc_3b_err[0][0]
,mc_3b_val[0][1], mc_3b_err[0][1]
,mc_3b_val[0][2], mc_3b_err[0][2]
,mc_3b_val[0][3], mc_3b_err[0][3]
,mc_3b_val[0][4], mc_3b_err[0][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET3 & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ \\vsmtvs \\\\ \n"
,mc_3b_val[1][0], mc_3b_err[1][0]
,mc_3b_val[1][1], mc_3b_err[1][1]
,mc_3b_val[1][2], mc_3b_err[1][2]
,mc_3b_val[1][3], mc_3b_err[1][3]
,mc_3b_val[1][4], mc_3b_err[1][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET4 & --- & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ & $%.1f \\pm %.1f$ \\vsmtvs \\\\ \n"
,mc_3b_val[2][1], mc_3b_err[2][1]
,mc_3b_val[2][2], mc_3b_err[2][2]
,mc_3b_val[2][3], mc_3b_err[2][3]
,mc_3b_val[2][4], mc_3b_err[2][4]
) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " \\hline\n" ) ;
fprintf( outfp, " $N_b\\ge3$, MET2-4 & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ & $%.0f \\pm %.0f$ \\vsmtvs \\\\ \n"
,mc_3b_val[3][0], mc_3b_err[3][0]
,mc_3b_val[3][1], mc_3b_err[3][1]
,mc_3b_val[3][2], mc_3b_err[3][2]
,mc_3b_val[3][3], mc_3b_err[3][3]
,mc_3b_val[3][4], mc_3b_err[3][4]
) ;
fprintf( outfp, " \\hline\n") ;
fprintf( outfp, " \\end{tabular}\n" ) ;
fprintf( outfp, " \\label{tab:pred-and-observed-nb3}\n") ;
fprintf( outfp, "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n") ;
fprintf( outfp, "\n\n\n" ) ;
fprintf( outfp, "\\end{document}\n" ) ;
fclose( outfp ) ;
}
void full_fit_hsbins_plots_components( bool doLogy=false ) {
TLatex* nblabel=new TLatex();
float header_x=0.37; float header_y=0.96;
if (usePublicStyle_) {
initCmsStyle();
gROOT->SetStyle("CMS");
gROOT->ForceStyle();
gStyle->SetPadTopMargin(0.06) ; //was 0.02
nblabel->SetTextSize(0.095);
nblabel->SetTextFont(42);
nblabel->SetNDC();
}
fill_arrays() ;
fill_frac_arrays("fitresults-ws-data-unblind-susyFixed0.0.txt") ;
/*
int fillStyle1sig = 1001 ;
int fillStyle2sig = 1001 ;
int fillColor1sig = kGreen-7 ;
int fillColor2sig = kGreen-10 ;
int lineColor = kGreen+3 ;
//appease sijin
fillColor2sig = kGreen-9; //or -7
fillColor1sig = kGreen+1;//-3;
*/
//int fillStyle1sig = 3344 ;
//int fillStyle2sig = 3644 ;
int fillStyle1sig = 3004 ;
int fillStyle2sig = 3005 ;
//int fillStyle1sig = 3001 ;
//int fillStyle2sig = 3003 ;
int fillColor1sig = kBlack ;
int fillColor2sig = kBlack ;
int lineColor = kBlack ;
gStyle->SetOptStat(0) ;
gStyle->SetLabelSize(usePublicStyle_? 0.17 : 0.12,"x") ;
gStyle->SetLabelSize(usePublicStyle_? 0.12 : 0.08,"y") ;
gStyle->SetTitleW(0.8) ;
gStyle->SetTitleH(0.08) ;
gStyle->SetTitleAlign(13) ;
gStyle->SetPadLeftMargin(usePublicStyle_ ? 0.35 : 0.14) ;
gStyle->SetPadRightMargin(0.02) ;
gStyle->SetLabelOffset(0.02,"y") ;
if (usePublicStyle_) gStyle->SetTitleSize(0.12,"y");
if (usePublicStyle_) gStyle->SetTitleOffset(1.5,"y");
gStyle->SetPadGridY(1) ;
if ( doLogy ) {
gStyle->SetOptLogy(1) ;
} else {
gStyle->SetOptLogy(0) ;
}
gDirectory->Delete("h*") ;
double xvals[4] = { 1.05, 2.05, 3.05, 4.05 } ;
double xel[4] = {0.5, 0.5, 0.5, 0.5} ;
double xeh[4] = {0.5, 0.5, 0.5, 0.5} ;
if (usePublicStyle_) {
for (int k=0;k<4;k++) xvals[k]-=0.05; //cancel out owen's fancy offset
}
TCanvas* cmip = (TCanvas*) gDirectory->FindObject("cmip") ;
if ( cmip == 0x0 ) {
cmip = new TCanvas("cmip", "MIP", 800, 600) ;
//cmip = new TCanvas("cmip", "MIP", 950, 600) ;
}
cmip->Clear() ;
//make pads by hand
/* original coordinates using divide command
pad 0.01 0.01 0.24 0.99
pad 0.26 0.01 0.49 0.99
pad 0.51 0.01 0.74 0.99
pad 0.76 0.01 0.99 0.99
*/
float ymax=usePublicStyle_ ? 0.85 : 0.99;
TPad* pad1 = new TPad("pad1","pad1", 0.01,0.01,0.24,ymax);
TPad* pad2 = new TPad("pad2","pad2", 0.26,0.01,0.49,ymax);
TPad* pad3 = new TPad("pad3","pad3", 0.51,0.01,0.74,ymax);
TPad* pad4 = new TPad("pad4","pad4", 0.76,0.01,0.99,ymax);
TPad* padh = usePublicStyle_ ? new TPad("padh","padh",0.01,ymax+0.02,0.99,0.99) : 0;
pad1->Draw(); pad2->Draw(); pad3->Draw(); pad4->Draw();
if (usePublicStyle_) padh->Draw();
TLatex* plotheader=0;
if (usePublicStyle_) { //copied from drawReducedTrees.h
padh->cd();
TString astring;
astring.Form("CMS, L = %.1f fb^{-1}, #sqrt{s} = 8 TeV",19.39);
plotheader = new TLatex(3.570061,23.08044,astring);
plotheader->SetNDC();
// plotheader->SetTextAlign(13);
plotheader->SetX(0.05);
// plotheader->SetY(0.4);
plotheader->SetY(0.25); //jmt now that _int is gone from L, realign this in Y
plotheader->SetTextFont(42);
plotheader->SetTextSizePixels(24);
plotheader->SetTextSize(0.35);
plotheader->Draw();
}
//--------
pad1->cd();
TH1F* h_data_nb2_m4 = new TH1F( "h_data_nb2_m4", "Data, nB=2, MET4", 4, 0.5, 4.5 ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb2_m4 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb2_m4 -> SetBinContent( 1, 0 ) ;
h_data_nb2_m4 -> SetBinContent( 2, data_2b_obs[0] ) ;
h_data_nb2_m4 -> SetBinContent( 3, data_2b_obs[1] ) ;
h_data_nb2_m4 -> SetBinContent( 4, data_2b_obs[2] ) ;
h_data_nb2_m4 -> SetMarkerStyle(20) ;
h_data_nb2_m4 -> SetMarkerSize(1.2) ;
h_data_nb2_m4 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb2_m4 -> SetMinimum( 10 ) ;
h_data_nb2_m4 -> SetMaximum( 90 ) ;
if(only1sigma_) h_data_nb2_m4 -> SetMaximum( 80 ) ;
if (usePublicStyle_) {
h_data_nb2_m4->SetYTitle("Events");
h_data_nb2_m4->GetXaxis()->LabelsOption("v");
}
h_data_nb2_m4 -> Draw("e") ;
double nb2_m4_yvals[4] = { -99., ff_2b_val[0], ff_2b_val[1], ff_2b_val[2] } ;
double nb2_m4_yerrh[4] = { 0.0, ff_2b_p1s[0], ff_2b_p1s[1], ff_2b_p1s[2] } ;
double nb2_m4_yerrl[4] = { 0.0, ff_2b_m1s[0], ff_2b_m1s[1], ff_2b_m1s[2] } ;
double nb2_m4_yerrh_2s[4] = { 0.0, ff_2b_p2s[0], ff_2b_p2s[1], ff_2b_p2s[2] } ;
double nb2_m4_yerrl_2s[4] = { 0.0, ff_2b_m2s[0], ff_2b_m2s[1], ff_2b_m2s[2] } ;
TGraphAsymmErrors* gr_nb2_m4 = new TGraphAsymmErrors( 4, xvals, nb2_m4_yvals, xel, xeh,
nb2_m4_yerrl, nb2_m4_yerrh ) ;
TGraphAsymmErrors* gr_nb2_m4_2s = new TGraphAsymmErrors( 4, xvals, nb2_m4_yvals, xel, xeh,
nb2_m4_yerrl_2s, nb2_m4_yerrh_2s ) ;
gr_nb2_m4 -> SetFillStyle(fillStyle1sig) ;
gr_nb2_m4 -> SetFillColor(fillColor1sig) ;
gr_nb2_m4 -> SetLineColor(lineColor) ;
gr_nb2_m4 -> SetLineWidth(1) ;
gr_nb2_m4 -> SetMarkerColor(4) ;
//gr_nb2_m4 -> SetMarkerStyle(24) ;
gr_nb2_m4_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb2_m4_2s -> SetFillColor(fillColor2sig) ;
gr_nb2_m4_2s -> SetLineColor(lineColor) ;
gr_nb2_m4_2s -> SetMarkerColor(4) ;
//gr_nb2_m4_2s -> SetMarkerStyle(24) ;
//Component Stack
int reorder[4] = {3,2,1,0}; //0 ttwj, 1 qcd, 2 znn, 3 vv
THStack* st_nb2_m4 = new THStack("h_stack_nb2_m4", "Stack, nB=2, MET4") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent = new TH1F("hcomponent_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent->SetFillColor(getFillColor(compName));
st_nb2_m4->Add(hcomponent);
hcomponent->SetBinContent(1,0);
hcomponent->SetBinContent(2,frac_2b_val[ci][0]*ff_2b_val[0]);
hcomponent->SetBinContent(3,frac_2b_val[ci][1]*ff_2b_val[1]);
hcomponent->SetBinContent(4,frac_2b_val[ci][2]*ff_2b_val[2]);
}
st_nb2_m4->Draw("SAME");
TString drawopt = usePublicStyle_? "2" : "p2";
if(!only1sigma_) gr_nb2_m4_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb2_m4_2s -> Draw("p") ;
gr_nb2_m4 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb2_m4_noyerr=(TGraphAsymmErrors*) gr_nb2_m4 ->Clone("gr_nb2_m4_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb2_m4_noyerr);
gr_nb2_m4_noyerr->Draw("pz");
}
else gr_nb2_m4 -> Draw("p") ;
h_data_nb2_m4 -> Draw("same e") ;
//--------
pad2->cd();
TH1F* h_data_nb3_m2 = new TH1F( "h_data_nb3_m2", "Data, nB=3, MET2", 4, 0.5, 4.5 ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m2 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m2 -> SetBinContent( 1, data_3b_obs[0][0] ) ;
h_data_nb3_m2 -> SetBinContent( 2, data_3b_obs[0][1] ) ;
h_data_nb3_m2 -> SetBinContent( 3, data_3b_obs[0][2] ) ;
h_data_nb3_m2 -> SetBinContent( 4, data_3b_obs[0][3] ) ;
h_data_nb3_m2 -> SetMarkerStyle(20) ;
h_data_nb3_m2 -> SetMarkerSize(1.2) ;
h_data_nb3_m2 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m2 -> SetMinimum( 4 ) ;
h_data_nb3_m2 -> SetMaximum( 220 ) ;
if(only1sigma_) h_data_nb3_m2 -> SetMaximum( 200 ) ;
if (usePublicStyle_) {
h_data_nb3_m2->SetYTitle("Events");
h_data_nb3_m2->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m2 -> Draw("e") ;
double nb3_m2_yvals[4] = { ff_3b_val[0][0], ff_3b_val[0][1], ff_3b_val[0][2], ff_3b_val[0][3] } ;
double nb3_m2_yerrh[4] = { ff_3b_p1s[0][0], ff_3b_p1s[0][1], ff_3b_p1s[0][2], ff_3b_p1s[0][3] } ;
double nb3_m2_yerrl[4] = { ff_3b_m1s[0][0], ff_3b_m1s[0][1], ff_3b_m1s[0][2], ff_3b_m1s[0][3] } ;
double nb3_m2_yerrh_2s[4] = { ff_3b_p2s[0][0], ff_3b_p2s[0][1], ff_3b_p2s[0][2], ff_3b_p2s[0][3] } ;
double nb3_m2_yerrl_2s[4] = { ff_3b_m2s[0][0], ff_3b_m2s[0][1], ff_3b_m2s[0][2], ff_3b_m2s[0][3] } ;
TGraphAsymmErrors* gr_nb3_m2 = new TGraphAsymmErrors( 4, xvals, nb3_m2_yvals, xel, xeh,
nb3_m2_yerrl, nb3_m2_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m2_2s = new TGraphAsymmErrors( 4, xvals, nb3_m2_yvals, xel, xeh,
nb3_m2_yerrl_2s, nb3_m2_yerrh_2s ) ;
gr_nb3_m2 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m2 -> SetFillColor(fillColor1sig) ;
gr_nb3_m2 -> SetLineColor(lineColor) ;
gr_nb3_m2 -> SetLineWidth(1) ;
gr_nb3_m2 -> SetMarkerColor(4) ;
//gr_nb3_m2 -> SetMarkerStyle(24) ;
gr_nb3_m2_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m2_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m2_2s -> SetLineColor(lineColor) ;
gr_nb3_m2_2s -> SetMarkerColor(4) ;
//gr_nb3_m2_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m2 = new THStack("h_stack_nb3_m2", "Stack, nB=3, MET2") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m2 = new TH1F("hcomponent_nb3_m2_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m2->SetFillColor(getFillColor(compName));
st_nb3_m2->Add(hcomponent_nb3_m2);
hcomponent_nb3_m2->SetBinContent(1,frac_3b_val[ci][0][0]*ff_3b_val[0][0]);
hcomponent_nb3_m2->SetBinContent(2,frac_3b_val[ci][0][1]*ff_3b_val[0][1]);
hcomponent_nb3_m2->SetBinContent(3,frac_3b_val[ci][0][2]*ff_3b_val[0][2]);
hcomponent_nb3_m2->SetBinContent(4,frac_3b_val[ci][0][3]*ff_3b_val[0][3]);
}
st_nb3_m2->Draw("SAME");
if(!only1sigma_) gr_nb3_m2_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m2_2s -> Draw("p") ;
gr_nb3_m2 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m2_noyerr=(TGraphAsymmErrors*) gr_nb3_m2 ->Clone("gr_nb3_m2_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m2_noyerr);
gr_nb3_m2_noyerr->Draw("pz");
}
else gr_nb3_m2 -> Draw("p") ;
h_data_nb3_m2 -> Draw("same e") ;
//--------
pad3->cd();
TH1F* h_data_nb3_m3 = new TH1F( "h_data_nb3_m3", "Data, nB=3, MET3", 4, 0.5, 4.5 ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m3 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m3 -> SetBinContent( 1, data_3b_obs[1][0] ) ;
h_data_nb3_m3 -> SetBinContent( 2, data_3b_obs[1][1] ) ;
h_data_nb3_m3 -> SetBinContent( 3, data_3b_obs[1][2] ) ;
h_data_nb3_m3 -> SetBinContent( 4, data_3b_obs[1][3] ) ;
h_data_nb3_m3 -> SetMarkerStyle(20) ;
h_data_nb3_m3 -> SetMarkerSize(1.2) ;
h_data_nb3_m3 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m3 -> SetMinimum( 0.3 ) ;
h_data_nb3_m3 -> SetMaximum( 45 ) ;
if (usePublicStyle_) {
h_data_nb3_m3->SetYTitle("Events");
h_data_nb3_m3->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m3 -> Draw("e") ;
double nb3_m3_yvals[4] = { ff_3b_val[1][0], ff_3b_val[1][1], ff_3b_val[1][2], ff_3b_val[1][3] } ;
double nb3_m3_yerrh[4] = { ff_3b_p1s[1][0], ff_3b_p1s[1][1], ff_3b_p1s[1][2], ff_3b_p1s[1][3] } ;
double nb3_m3_yerrl[4] = { ff_3b_m1s[1][0], ff_3b_m1s[1][1], ff_3b_m1s[1][2], ff_3b_m1s[1][3] } ;
double nb3_m3_yerrh_2s[4] = { ff_3b_p2s[1][0], ff_3b_p2s[1][1], ff_3b_p2s[1][2], ff_3b_p2s[1][3] } ;
double nb3_m3_yerrl_2s[4] = { ff_3b_m2s[1][0], ff_3b_m2s[1][1], ff_3b_m2s[1][2], ff_3b_m2s[1][3] } ;
TGraphAsymmErrors* gr_nb3_m3 = new TGraphAsymmErrors( 4, xvals, nb3_m3_yvals, xel, xeh,
nb3_m3_yerrl, nb3_m3_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m3_2s = new TGraphAsymmErrors( 4, xvals, nb3_m3_yvals, xel, xeh,
nb3_m3_yerrl_2s, nb3_m3_yerrh_2s ) ;
gr_nb3_m3 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m3 -> SetFillColor(fillColor1sig) ;
gr_nb3_m3 -> SetLineColor(lineColor) ;
gr_nb3_m3 -> SetLineWidth(1) ;
gr_nb3_m3 -> SetMarkerColor(4) ;
//gr_nb3_m3 -> SetMarkerStyle(24) ;
gr_nb3_m3_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m3_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m3_2s -> SetLineColor(lineColor) ;
gr_nb3_m3_2s -> SetMarkerColor(4) ;
//gr_nb3_m3_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m3 = new THStack("h_stack_nb3_m3", "Stack, nB=3, MET3") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m3 = new TH1F("hcomponent_nb3_m3_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m3->SetFillColor(getFillColor(compName));
st_nb3_m3->Add(hcomponent_nb3_m3);
hcomponent_nb3_m3->SetBinContent(1,frac_3b_val[ci][1][0]*ff_3b_val[1][0]);
hcomponent_nb3_m3->SetBinContent(2,frac_3b_val[ci][1][1]*ff_3b_val[1][1]);
hcomponent_nb3_m3->SetBinContent(3,frac_3b_val[ci][1][2]*ff_3b_val[1][2]);
hcomponent_nb3_m3->SetBinContent(4,frac_3b_val[ci][1][3]*ff_3b_val[1][3]);
}
st_nb3_m3->Draw("SAME");
if(!only1sigma_) gr_nb3_m3_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m3_2s -> Draw("p") ;
gr_nb3_m3 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m3_noyerr=(TGraphAsymmErrors*) gr_nb3_m3 ->Clone("gr_nb3_m3_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m3_noyerr);
gr_nb3_m3_noyerr->Draw("pz");
}
else gr_nb3_m3 -> Draw("p") ;
h_data_nb3_m3 -> Draw("same e") ;
//--------
pad4->cd();
TH1F* h_data_nb3_m4 = new TH1F( "h_data_nb3_m4", "Data, nB=3, MET4", 4, 0.5, 4.5 ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 1, "HT1" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 2, "HT2" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 3, "HT3" ) ;
h_data_nb3_m4 -> GetXaxis() -> SetBinLabel( 4, "HT4" ) ;
h_data_nb3_m4 -> SetBinContent( 1, 0 ) ;
h_data_nb3_m4 -> SetBinContent( 2, data_3b_obs[2][1] ) ;
h_data_nb3_m4 -> SetBinContent( 3, data_3b_obs[2][2] ) ;
h_data_nb3_m4 -> SetBinContent( 4, data_3b_obs[2][3] ) ;
h_data_nb3_m4 -> SetMarkerStyle(20) ;
h_data_nb3_m4 -> SetMarkerSize(1.2) ;
h_data_nb3_m4 -> SetLineWidth(2) ;
if ( doLogy) h_data_nb3_m4 -> SetMinimum( 0.1 ) ;
h_data_nb3_m4 -> SetMaximum( 14 ) ;
if(only1sigma_) h_data_nb3_m4 -> SetMaximum( 12 ) ;
if (usePublicStyle_) {
h_data_nb3_m4->SetYTitle("Events");
h_data_nb3_m4->GetXaxis()->LabelsOption("v");
}
h_data_nb3_m4 -> Draw("e") ;
double nb3_m4_yvals[4] = { -99., ff_3b_val[2][1], ff_3b_val[2][2], ff_3b_val[2][3] } ;
double nb3_m4_yerrh[4] = { 0.0, ff_3b_p1s[2][1], ff_3b_p1s[2][2], ff_3b_p1s[2][3] } ;
double nb3_m4_yerrl[4] = { 0.0, ff_3b_m1s[2][1], ff_3b_m1s[2][2], ff_3b_m1s[2][3] } ;
double nb3_m4_yerrh_2s[4] = { 0.0, ff_3b_p2s[2][1], ff_3b_p2s[2][2], ff_3b_p2s[2][3] } ;
double nb3_m4_yerrl_2s[4] = { 0.0, ff_3b_m2s[2][1], ff_3b_m2s[2][2], ff_3b_m2s[2][3] } ;
TGraphAsymmErrors* gr_nb3_m4 = new TGraphAsymmErrors( 4, xvals, nb3_m4_yvals, xel, xeh,
nb3_m4_yerrl, nb3_m4_yerrh ) ;
TGraphAsymmErrors* gr_nb3_m4_2s = new TGraphAsymmErrors( 4, xvals, nb3_m4_yvals, xel, xeh,
nb3_m4_yerrl_2s, nb3_m4_yerrh_2s ) ;
gr_nb3_m4 -> SetFillStyle(fillStyle1sig) ;
gr_nb3_m4 -> SetFillColor(fillColor1sig) ;
gr_nb3_m4 -> SetLineColor(lineColor) ;
gr_nb3_m4 -> SetLineWidth(1) ;
gr_nb3_m4 -> SetMarkerColor(4) ;
//gr_nb3_m4 -> SetMarkerStyle(24) ;
gr_nb3_m4_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_m4_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_m4_2s -> SetLineColor(lineColor) ;
gr_nb3_m4_2s -> SetMarkerColor(4) ;
//gr_nb3_m4_2s -> SetMarkerStyle(24) ;
//Component Stack
THStack* st_nb3_m4 = new THStack("h_stack_nb3_m4", "Stack, nB=3, MET4") ;
for ( int ici=0; ici<ncomps; ici++ ) {
int ci = reorder[ici];
TString compName = comp_fracfile_name[ci];
TH1F* hcomponent_nb3_m4 = new TH1F("hcomponent_nb3_m4_"+compName, "component "+compName, 4, 0.5, 4.5) ;
hcomponent_nb3_m4->SetFillColor(getFillColor(compName));
st_nb3_m4->Add(hcomponent_nb3_m4);
hcomponent_nb3_m4->SetBinContent(1,0);
hcomponent_nb3_m4->SetBinContent(2,frac_3b_val[ci][2][1]*ff_3b_val[2][1]);
hcomponent_nb3_m4->SetBinContent(3,frac_3b_val[ci][2][2]*ff_3b_val[2][2]);
hcomponent_nb3_m4->SetBinContent(4,frac_3b_val[ci][2][3]*ff_3b_val[2][3]);
}
st_nb3_m4->Draw("SAME");
if(!only1sigma_) gr_nb3_m4_2s -> Draw(drawopt) ;
if (!usePublicStyle_) gr_nb3_m4_2s -> Draw("p") ;
gr_nb3_m4 -> Draw(drawopt) ;
TGraphAsymmErrors* gr_nb3_m4_noyerr=(TGraphAsymmErrors*) gr_nb3_m4 ->Clone("gr_nb3_m4_noyerr");
if (usePublicStyle_ ) {
formatClone(gr_nb3_m4_noyerr);
gr_nb3_m4_noyerr->Draw("pz");
}
else gr_nb3_m4 -> Draw("p") ;
h_data_nb3_m4 -> Draw("same e") ;
TLegend * leg1=0;
//TLegend * leg2=0;
//dummy histograms for legend -- hardcoded :(
TH1F* hdummy_ttwj = new TH1F("hdummy_ttwj", "dummy ttwj", 4, 0.5, 4.5) ;
hdummy_ttwj->SetFillColor(getFillColor("ttwj"));
TH1F* hdummy_qcd = new TH1F("hdummy_qcd", "dummy qcd", 4, 0.5, 4.5) ;
hdummy_qcd->SetFillColor(getFillColor("qcd"));
TH1F* hdummy_znn = new TH1F("hdummy_znn", "dummy znn", 4, 0.5, 4.5) ;
hdummy_znn->SetFillColor(getFillColor("znn"));
TH1F* hdummy_vv = new TH1F("hdummy_vv", "dummy vv", 4, 0.5, 4.5) ;
hdummy_vv->SetFillColor(getFillColor("vv"));
TH1F* hdummy_blank = new TH1F("hdummy_blank", "", 4, 0.5, 4.5) ;
hdummy_blank->SetFillColor(kWhite);
hdummy_blank->SetLineColor(kWhite);
if (usePublicStyle_) {
pad1->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} = 2, MET4");
pad2->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET2");
pad3->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET3");
pad4->cd();
nblabel->DrawLatex(header_x,header_y,"N_{b-jet} #geq 3, MET4");
pad1->GetFrame()->Draw();
pad1->RedrawAxis();
pad2->GetFrame()->Draw();
pad2->RedrawAxis();
pad3->GetFrame()->Draw();
pad3->RedrawAxis();
pad4->GetFrame()->Draw();
pad4->RedrawAxis();
padh->cd();
float leg_x1=0.5,leg_y1=0.01,leg_x2=0.98,leg_y2=0.98;
leg1 = new TLegend(leg_x1,leg_y1,leg_x2,leg_y2);
leg1->SetBorderSize(0);
leg1->SetLineStyle(0);
leg1->SetTextFont(42);
leg1->SetFillStyle(0);
leg1->SetNColumns(2);
leg1->AddEntry(h_data_nb2_m4,"Data");
leg1->AddEntry(gr_nb2_m4_noyerr,"Full fit","lf");
//leg1->AddEntry(hdummy_znn,"Z#nu#nu","f");
leg1->AddEntry(hdummy_ttwj,"Top quark & W+jets","f");
leg1->AddEntry(hdummy_qcd,"QCD","f");
leg1->AddEntry(hdummy_znn,"Z+jets","f");
leg1->AddEntry(hdummy_vv, "Diboson & Drell-Yan", "f");
//leg1->AddEntry(hdummy_vv,"Diboson + Drell-Yan","f");
leg1->Draw();
/*
float leg2_x1=0.668,leg2_y1=0.01,leg2_x2=0.9,leg2_y2=0.333;
leg2 = new TLegend(leg2_x1,leg2_y1,leg2_x2,leg2_y2);
leg2->SetBorderSize(0);
leg2->SetLineStyle(0);
leg2->SetTextFont(42);
leg2->SetFillStyle(0);
//leg2->SetNColumns(2);
leg2->AddEntry(hdummy_ttwj,"top quark and W+jets","f");
leg2->Draw();
*/
}
cmip->Update() ;
cmip->Draw() ;
TString nsig = "2sig";
if(only1sigma_) nsig = "1sig";
if ( doLogy ) {
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb-logy_components.pdf" ) ;
} else {
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb_components.pdf" ) ;
cmip->SaveAs("fitresult-fullfit-hsbins-"+nsig+"-19fb_components.png" ) ;
}
}
void binsum_hsbins_plots( bool doLogy=false ) {
fill_arrays() ;
int fillStyle1sig = 1001 ;
int fillStyle2sig = 1001 ;
int fillColor1sig = kBlue-9 ;
int fillColor2sig = kBlue-10 ;
int fillColor1sig_ff = kGreen-7 ;
int fillColor2sig_ff = kGreen-10 ;
gStyle->SetOptStat(0) ;
gStyle->SetLabelSize(0.12,"x") ;
gStyle->SetLabelSize(0.08,"y") ;
gStyle->SetTitleW(0.8) ;
gStyle->SetTitleH(0.08) ;
gStyle->SetTitleAlign(13) ;
gStyle->SetPadLeftMargin(0.14) ;
gStyle->SetPadRightMargin(0.02) ;
gStyle->SetPadBottomMargin(0.30) ;
gStyle->SetLabelOffset(0.02,"y") ;
gStyle->SetPadGridY(1) ;
if ( doLogy ) {
gStyle->SetOptLogy(1) ;
} else {
gStyle->SetOptLogy(0) ;
}
double xvals[4] = { 1.05, 2.05, 3.05, 4.05 } ;
double xel[4] = {0.5, 0.5, 0.5, 0.5} ;
double xeh[4] = {0.5, 0.5, 0.5, 0.5} ;
TCanvas* cmipbs = (TCanvas*) gDirectory->FindObject("cmipbs") ;
if ( cmipbs == 0x0 ) {
cmipbs = new TCanvas("cmipbs", "MIP", 800, 600) ;
}
cmipbs->Clear() ;
cmipbs->Divide(4,1) ;
//--------
cmipbs->cd(1) ;
TH1F* h_data_nb3_hs = new TH1F( "h_data_nb3_hs", "Data, nB=3, HT bins", 4, 0.5, 4.5 ) ;
h_data_nb3_hs -> GetXaxis() -> SetBinLabel( 1, "HT1, sum of MET2,3" ) ;
h_data_nb3_hs -> GetXaxis() -> SetBinLabel( 2, "HT2, sum of MET2,3,4" ) ;
h_data_nb3_hs -> GetXaxis() -> SetBinLabel( 3, "HT3, sum of MET2,3,4" ) ;
h_data_nb3_hs -> GetXaxis() -> SetBinLabel( 4, "HT4, sum of MET2,3,4" ) ;
h_data_nb3_hs -> GetXaxis() -> LabelsOption("v") ;
h_data_nb3_hs -> SetBinContent( 1, data_3b_obs[3][0] ) ;
h_data_nb3_hs -> SetBinContent( 2, data_3b_obs[3][1] ) ;
h_data_nb3_hs -> SetBinContent( 3, data_3b_obs[3][2] ) ;
h_data_nb3_hs -> SetBinContent( 4, data_3b_obs[3][3] ) ;
h_data_nb3_hs -> SetMarkerStyle(20) ;
h_data_nb3_hs -> SetMarkerSize(1.2) ;
h_data_nb3_hs -> SetLineWidth(2) ;
if ( doLogy ) h_data_nb3_hs -> SetMinimum( 4 ) ;
h_data_nb3_hs -> SetMaximum( 440 ) ;
h_data_nb3_hs -> Draw("e") ;
double nb3_hs_yvals[4] = { ub_3b_val[3][0], ub_3b_val[3][1], ub_3b_val[3][2], ub_3b_val[3][3] } ;
double nb3_hs_yerrh[4] = { ub_3b_p1s[3][0], ub_3b_p1s[3][1], ub_3b_p1s[3][2], ub_3b_p1s[3][3] } ;
double nb3_hs_yerrl[4] = { ub_3b_m1s[3][0], ub_3b_m1s[3][1], ub_3b_m1s[3][2], ub_3b_m1s[3][3] } ;
double nb3_hs_yerrh_2s[4] = { ub_3b_p2s[3][0], ub_3b_p2s[3][1], ub_3b_p2s[3][2], ub_3b_p2s[3][3] } ;
double nb3_hs_yerrl_2s[4] = { ub_3b_m2s[3][0], ub_3b_m2s[3][1], ub_3b_m2s[3][2], ub_3b_m2s[3][3] } ;
TGraphAsymmErrors* gr_nb3_hs = new TGraphAsymmErrors( 4, xvals, nb3_hs_yvals, xel, xeh,
nb3_hs_yerrl, nb3_hs_yerrh ) ;
TGraphAsymmErrors* gr_nb3_hs_2s = new TGraphAsymmErrors( 4, xvals, nb3_hs_yvals, xel, xeh,
nb3_hs_yerrl_2s, nb3_hs_yerrh_2s ) ;
gr_nb3_hs -> SetFillStyle(fillStyle1sig) ;
gr_nb3_hs -> SetFillColor(fillColor1sig) ;
gr_nb3_hs -> SetLineColor(4) ;
gr_nb3_hs -> SetLineWidth(1) ;
gr_nb3_hs -> SetMarkerColor(4) ;
gr_nb3_hs_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_hs_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_hs_2s -> SetLineColor(4) ;
gr_nb3_hs_2s -> SetMarkerColor(4) ;
gr_nb3_hs_2s -> Draw("p2") ;
gr_nb3_hs_2s -> Draw("p") ;
gr_nb3_hs -> Draw("p2") ;
gr_nb3_hs -> Draw("p") ;
h_data_nb3_hs -> Draw("same e") ;
//--------
cmipbs->cd(2) ;
TH1F* h_data_nb3_ms = new TH1F( "h_data_nb3_ms", "Data, nB=3, MET bins", 3, 0.5, 3.5 ) ;
h_data_nb3_ms -> GetXaxis() -> SetBinLabel( 1, "MET2, sum of HT1,2,3,4" ) ;
h_data_nb3_ms -> GetXaxis() -> SetBinLabel( 2, "MET3, sum of HT1,2,3,4" ) ;
h_data_nb3_ms -> GetXaxis() -> SetBinLabel( 3, "MET4, sum of HT2,3,4" ) ;
h_data_nb3_ms -> GetXaxis() -> LabelsOption("v") ;
h_data_nb3_ms -> SetBinContent( 1, data_3b_obs[0][4] ) ;
h_data_nb3_ms -> SetBinContent( 2, data_3b_obs[1][4] ) ;
h_data_nb3_ms -> SetBinContent( 3, data_3b_obs[2][4] ) ;
h_data_nb3_ms -> SetMarkerStyle(20) ;
h_data_nb3_ms -> SetMarkerSize(1.2) ;
h_data_nb3_ms -> SetLineWidth(2) ;
if ( doLogy ) h_data_nb3_ms -> SetMinimum( 4 ) ;
h_data_nb3_ms -> SetMaximum( 580 ) ;
h_data_nb3_ms -> Draw("e") ;
double nb3_ms_yvals[3] = { ub_3b_val[0][4], ub_3b_val[1][4], ub_3b_val[2][4] } ;
double nb3_ms_yerrh[3] = { ub_3b_p1s[0][4], ub_3b_p1s[1][4], ub_3b_p1s[2][4] } ;
double nb3_ms_yerrl[3] = { ub_3b_m1s[0][4], ub_3b_m1s[1][4], ub_3b_m1s[2][4] } ;
double nb3_ms_yerrh_2s[3] = { ub_3b_p2s[0][4], ub_3b_p2s[1][4], ub_3b_p2s[2][4] } ;
double nb3_ms_yerrl_2s[3] = { ub_3b_m2s[0][4], ub_3b_m2s[1][4], ub_3b_m2s[2][4] } ;
TGraphAsymmErrors* gr_nb3_ms = new TGraphAsymmErrors( 4, xvals, nb3_ms_yvals, xel, xeh,
nb3_ms_yerrl, nb3_ms_yerrh ) ;
TGraphAsymmErrors* gr_nb3_ms_2s = new TGraphAsymmErrors( 4, xvals, nb3_ms_yvals, xel, xeh,
nb3_ms_yerrl_2s, nb3_ms_yerrh_2s ) ;
gr_nb3_ms -> SetFillStyle(fillStyle1sig) ;
gr_nb3_ms -> SetFillColor(fillColor1sig) ;
gr_nb3_ms -> SetLineColor(4) ;
gr_nb3_ms -> SetLineWidth(1) ;
gr_nb3_ms -> SetMarkerColor(4) ;
gr_nb3_ms_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_ms_2s -> SetFillColor(fillColor2sig) ;
gr_nb3_ms_2s -> SetLineColor(4) ;
gr_nb3_ms_2s -> SetMarkerColor(4) ;
gr_nb3_ms_2s -> Draw("p2") ;
gr_nb3_ms_2s -> Draw("p") ;
gr_nb3_ms -> Draw("p2") ;
gr_nb3_ms -> Draw("p") ;
h_data_nb3_ms -> Draw("same e") ;
//++++++++++++++++
cmipbs->cd(3) ;
h_data_nb3_hs -> Draw("e") ;
double nb3_hs_ff_yvals[4] = { ff_3b_val[3][0], ff_3b_val[3][1], ff_3b_val[3][2], ff_3b_val[3][3] } ;
double nb3_hs_ff_yerrh[4] = { ff_3b_p1s[3][0], ff_3b_p1s[3][1], ff_3b_p1s[3][2], ff_3b_p1s[3][3] } ;
double nb3_hs_ff_yerrl[4] = { ff_3b_m1s[3][0], ff_3b_m1s[3][1], ff_3b_m1s[3][2], ff_3b_m1s[3][3] } ;
double nb3_hs_ff_yerrh_2s[4] = { ff_3b_p2s[3][0], ff_3b_p2s[3][1], ff_3b_p2s[3][2], ff_3b_p2s[3][3] } ;
double nb3_hs_ff_yerrl_2s[4] = { ff_3b_m2s[3][0], ff_3b_m2s[3][1], ff_3b_m2s[3][2], ff_3b_m2s[3][3] } ;
TGraphAsymmErrors* gr_nb3_hs_ff = new TGraphAsymmErrors( 4, xvals, nb3_hs_ff_yvals, xel, xeh,
nb3_hs_ff_yerrl, nb3_hs_ff_yerrh ) ;
TGraphAsymmErrors* gr_nb3_hs_ff_2s = new TGraphAsymmErrors( 4, xvals, nb3_hs_ff_yvals, xel, xeh,
nb3_hs_ff_yerrl_2s, nb3_hs_ff_yerrh_2s ) ;
gr_nb3_hs_ff -> SetFillStyle(fillStyle1sig) ;
gr_nb3_hs_ff -> SetFillColor(fillColor1sig_ff) ;
gr_nb3_hs_ff -> SetLineColor(kGreen+2) ;
gr_nb3_hs_ff -> SetLineWidth(1) ;
gr_nb3_hs_ff -> SetMarkerColor(4) ;
gr_nb3_hs_ff_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_hs_ff_2s -> SetFillColor(fillColor2sig_ff) ;
gr_nb3_hs_ff_2s -> SetLineColor(kGreen+2) ;
gr_nb3_hs_ff_2s -> SetMarkerColor(4) ;
gr_nb3_hs_ff_2s -> Draw("p2") ;
gr_nb3_hs_ff_2s -> Draw("p") ;
gr_nb3_hs_ff -> Draw("p2") ;
gr_nb3_hs_ff -> Draw("p") ;
h_data_nb3_hs -> Draw("same e") ;
//--------
cmipbs->cd(4) ;
h_data_nb3_ms -> Draw("e") ;
double nb3_ms_ff_yvals[3] = { ff_3b_val[0][4], ff_3b_val[1][4], ff_3b_val[2][4] } ;
double nb3_ms_ff_yerrh[3] = { ff_3b_p1s[0][4], ff_3b_p1s[1][4], ff_3b_p1s[2][4] } ;
double nb3_ms_ff_yerrl[3] = { ff_3b_m1s[0][4], ff_3b_m1s[1][4], ff_3b_m1s[2][4] } ;
double nb3_ms_ff_yerrh_2s[3] = { ff_3b_p2s[0][4], ff_3b_p2s[1][4], ff_3b_p2s[2][4] } ;
double nb3_ms_ff_yerrl_2s[3] = { ff_3b_m2s[0][4], ff_3b_m2s[1][4], ff_3b_m2s[2][4] } ;
TGraphAsymmErrors* gr_nb3_ms_ff = new TGraphAsymmErrors( 4, xvals, nb3_ms_ff_yvals, xel, xeh,
nb3_ms_ff_yerrl, nb3_ms_ff_yerrh ) ;
TGraphAsymmErrors* gr_nb3_ms_ff_2s = new TGraphAsymmErrors( 4, xvals, nb3_ms_ff_yvals, xel, xeh,
nb3_ms_ff_yerrl_2s, nb3_ms_ff_yerrh_2s ) ;
gr_nb3_ms_ff -> SetFillStyle(fillStyle1sig) ;
gr_nb3_ms_ff -> SetFillColor(fillColor1sig_ff) ;
gr_nb3_ms_ff -> SetLineColor(kGreen+2) ;
gr_nb3_ms_ff -> SetLineWidth(1) ;
gr_nb3_ms_ff -> SetMarkerColor(4) ;
gr_nb3_ms_ff_2s -> SetFillStyle(fillStyle2sig) ;
gr_nb3_ms_ff_2s -> SetFillColor(fillColor2sig_ff) ;
gr_nb3_ms_ff_2s -> SetLineColor(kGreen+2) ;
gr_nb3_ms_ff_2s -> SetMarkerColor(4) ;
gr_nb3_ms_ff_2s -> Draw("p2") ;
gr_nb3_ms_ff_2s -> Draw("p") ;
gr_nb3_ms_ff -> Draw("p2") ;
gr_nb3_ms_ff -> Draw("p") ;
h_data_nb3_ms -> Draw("same e") ;
if ( doLogy ) {
cmipbs->SaveAs("fitresult-binsums-hsbins-2sig-19fb-logy.pdf") ;
} else {
cmipbs->SaveAs("fitresult-binsums-hsbins-2sig-19fb.pdf") ;
}
}
