* Makes the 4x4 plot for the MSSM paper that has the most important optimized
* cuts for hZ (Y34 and Min(Mij)) and hA (Btag and Chi_mass).
* Signals are normalized to 50 x efficiency
 
 exec style_2by2
 set xlab 1.75
 set ylab -0.4   | Yes, bizzare but true
 set ndvx 10508
 set ndvy 3
 set hmax 0.1
 set yhti 0.9
 set tsiz 0.58
 
 close 1
 close 2
 close 3
 close 4 

* histogram/file 1 ~/temp/hists/black_events_inc_391-391_70-70.hist
 histogram/file 1 ~/temp/hist/maximal_2_060-095_391_391.hist
 
 v/crea err(40) r
 
 h/pl //lun1/hasg/50310
 h/pl //lun1/hasg/50311
 h/pl //lun1/hasg/50312
 
 h/pl //lun1/hasg/50330
 h/pl //lun1/hasg/50331
 h/pl //lun1/hasg/50332
 
*$$$ h/pl //lun1/hzsg/51420
*$$$ h/pl //lun1/hzsg/51421
*$$$ h/pl //lun1/hzsg/51422
*$$$ 
*$$$ h/pl //lun1/hzsg/51440
*$$$ h/pl //lun1/hzsg/51441
*$$$ h/pl //lun1/hzsg/51442
 
 cd //pawc
 get/con 50310 err
 sigma err=sqrt(err)
 put/err 50310 err
 get/con 50330 err
 sigma err=sqrt(err)
 put/err 50330 err

 zone 2 2
 opt utit
 opt logy
* opt liny

* AZHA Analysis on hA->bbbb
* tit 'HA Analysis on hA"5#bb&^-!bb&^-!' u
* tit 'hA?bbbb! Analysis on hA"5#bb&^_!bb&^_!' u
 tit 'F(m?A!=80 GeV, m?h!=80 GeV)' u
* nul 0 5 0.001 500
 set hmax 0.001
 set hcol 1303
 set htyp 235
 norm=$eval($hinfo(50312,'SUM')*50.0)
 norm 50312 [norm]
 h/pl 50312 hist
 set hcol 1
 set htyp 0
 h/pl 50311 shist
 tit '' u
 set hmax 0.1
 h/pl 50310 se0
* atit 'hA?bbbb! Analysis' 'Events/0.125' ! 220 T
* atit '' 'Events/0.125' ! 220
* exec stathist#tristat 10 11 12 0.25 0.35 0.75
 set chhe 0.8
 opt liny
 null 0 1 0 1 sab
 itx 0.80 0.85 'a)' 
 exec data   0.2 0.95 0.05 20 0.6
 exec signal 0.2 0.85 0.6 555 'BG'
 exec signal 0.2 0.75 0.6 235 'hA"5#bb&^[-]!bb&^[-]!' 3 
 opt logy
 
* HA Analysis on hA->bbbb
* tit 'hA"\242#?bbbb! Analysis on hA"5#bb&^_!bb&^_!' u
 set hmax 0.01
 tit 'F(m?A!=80 GeV, m?h!=80 GeV)' u
* nul 0 5 0.001 60
 set hcol 1303
 set htyp 235
 norm=$eval($hinfo(50332,'SUM')*50.0)
 norm 50332 [norm]
 h/pl 50332 hist
 tit '' u
 set hcol 1
 set htyp 0
 h/pl 50331 shist
 h/pl 50330 se0
* atit 'hA"\242#?bbbb! Analysis' 'Events/0.125' ! 220 T
* atit '' 'Events/0.125' ! 220 
* exec stathist#tristat 20 21 22 0.25 0.35 0.75
 set chhe 0.8
 opt liny
 null 0 1 0 1 sab
 itx 0.8 0.85 'b)'
 opt logy
 set hmax 0.1
 
 return 
*$$$* hZ analyses
*$$$ 
*$$$ get/con 51420 err
*$$$ sigma err=sqrt(err)
*$$$ put/err 51420 err
*$$$ get/con 51440 err
*$$$ sigma err=sqrt(err)
*$$$ put/err 51440 err
*$$$
*$$$* AZHZ Analysis on hZ->bbqq
*$$$ tit 'F(m?Z!, m?h!=97 GeV)' u
*$$$* nul 0 5 0.001 500
*$$$ set hmax 0.01
*$$$ set hcol 1303
*$$$ set htyp 235
*$$$ norm=$eval($hinfo(51422,'SUM')*50.0)
*$$$ norm 51422 [norm]
*$$$ h/pl 51422 hist
*$$$ set hcol 1
*$$$ set htyp 0
*$$$ h/pl 51421 shist
*$$$ tit '' u
*$$$ set hmax 0.1
*$$$ h/pl 51420 se0
*$$$* atit 'hA?bbbb! Analysis' 'Events/0.125' ! 220 T
*$$$* atit '' 'Events/0.125' ! 220
*$$$* exec stathist#tristat 10 11 12 0.25 0.35 0.75
*$$$ set chhe 0.8
*$$$ opt liny
*$$$ null 0 1 0 1 sab
*$$$ itx 0.80 0.85 'c)' 
*$$$ exec data   0.2 0.95 0.05 20 0.6
*$$$ exec signal 0.2 0.85 0.6 555 'BG'
*$$$ exec signal 0.2 0.75 0.6 235 'hZ"5#bb&^[-]!qq&^_!' 3 
*$$$ opt logy
*$$$ 
*$$$* HZ Analysis on hZ->bbqq
*$$$ set hmax 0.01
*$$$ tit 'F(m?Z!, m?h!=97 GeV)' u
*$$$* nul 0 5 0.001 60
*$$$ set hcol 1303
*$$$ set htyp 235
*$$$ norm=$eval($hinfo(51442,'SUM')*50.0)
*$$$ norm 51442 [norm]
*$$$ h/pl 51442 hist
*$$$ tit '' u
*$$$ set hcol 1
*$$$ set htyp 0
*$$$ h/pl 51441 shist
*$$$ h/pl 51440 se0
*$$$* atit 'hA"\242#?bbbb! Analysis' 'Events/0.125' ! 220 T
*$$$* atit '' 'Events/0.125' ! 220 
*$$$* exec stathist#tristat 20 21 22 0.25 0.35 0.75
*$$$ set chhe 0.8
*$$$ opt liny
*$$$ null 0 1 0 1 sab
*$$$ itx 0.8 0.85 'd)'
*$$$ opt logy
*$$$ set hmax 0.1
*$$$ 
*$$$ return