* 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 
* TanB 50 hA dominates
 histogram/file 1 /hp/sd01d/data/ifisk/max_mssm_391-391_70-85.hist 8190 X
* TanB 1.5 hZ dominates
 histogram/file 2 /hp/sd01d/data/ifisk/max_mssm_041-041_326-370.hist 8190 X
 
 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 //lun2/hzsg/50420
 h/pl //lun2/hzsg/50421
 h/pl //lun2/hzsg/50422
 
 h/pl //lun2/hzsg/50340
 h/pl //lun2/hzsg/50341
 h/pl //lun2/hzsg/50342
 
 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
 
 
* hZ analyses
 
 get/con 50420 err
 sigma err=sqrt(err)
 put/err 50420 err
 get/con 50340 err
 sigma err=sqrt(err)
 put/err 50340 err

* AZHZ Analysis on hZ->bbqq
 tit 'F(m?Z!, m?h!=95 GeV)' u
* nul 0 5 0.001 500
 set hmax 0.01
 set hcol 1303
 set htyp 235
 norm=$eval($hinfo(50422,'SUM')*50.0)
 norm 50422 [norm]
 h/pl 50422 hist
 set hcol 1
 set htyp 0
 h/pl 50421 shist
 tit '' u
 set hmax 0.1
 h/pl 50420 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 'hA"5#bb&^[-]!bb&^[-]!' 3 
 opt logy
 
* HZ Analysis on hZ->bbqq
 set hmax 0.01
 tit 'F(m?Z!, m?h!=95 GeV)' u
* nul 0 5 0.001 60
 set hcol 1303
 set htyp 235
 norm=$eval($hinfo(50342,'SUM')*50.0)
 norm 50342 [norm]
 h/pl 50342 hist
 tit '' u
 set hcol 1
 set htyp 0
 h/pl 50341 shist
 h/pl 50340 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 
 
*$$$* HZ Analysis on hZ->bbqq
*$$$* tit 'hZ?bbqq! Analysis on hZ"5#bb&^_!qq&^_!' u
*$$$ tit 'F(m?Z!, m?h!=85 GeV)' u
*$$$ h/pl 31 hist
*$$$ tit '' u
*$$$ set hcol 1303
*$$$ set htyp 235
*$$$ norm=$eval($hinfo(32,'SUM')*50.0)
*$$$ norm 32 [norm]
*$$$ h/pl 32 shist
*$$$ set hcol 1
*$$$ set htyp 0
*$$$* atit 'hZ?bbqq! Analysis' 'Events/0.5' ! 220 T
*$$$* atit '' 'Events/0.5' ! 220 
*$$$* exec stathist#tristat 30 31 32 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
*$$$* tit 'hZ"\242#?bbqq! Analysis on hZ"5#bb&^_!qq&^_!' u
*$$$ tit 'F(m?Z!, m?h!=85 GeV)' u
*$$$ nul 0 5 0.005 150
*$$$ h/pl 41 shist
*$$$ tit '' u
*$$$ h/pl 40 se0
*$$$ set hcol 1303
*$$$ set htyp 235
*$$$ norm=$eval($hinfo(42,'SUM')*50.0)
*$$$ norm 42 [norm]
*$$$ h/pl 42 shist
*$$$ set hcol 1
*$$$ set htyp 0
*$$$ h/pl 40 se0
*$$$* atit 'hZ"\242#?bbqq! Analysis' 'Events/0.125' ! 220 T
*$$$* atit '' 'Events/0.125' ! 220
*$$$* exec stathist#tristat 40 41 42 0.25 0.35 0.75
*$$$ set chhe 0.8
*$$$ opt liny
*$$$ null 0 1 0 1 sab
*$$$ itx 0.80 0.85 'd)'
*$$$ opt logy
*$$$
*$$$ 
*$$$ 
*$$$ opt liny
*$$$ set hmax 0.8
*$$$ set tsiz 0.7