close 1
close 2
op nstat
op nbox
set ymgu 1
set ymgl 3
set xmgr 1.5
set xmgl 2.5
set xwid 8
set ywid 8
set lwid 8
set hwid 8
set bwid 8
set pwid 8
set fwid 8
SET GFON  22
SET VFON  22
SET LFON  132
SET TFON  22
SET CFON  22
SET VSIZ 0.7
SET ASIZ 0.7
set tsiz 0.7
set xsiz 19
set ksiz 0.36
set pmci 3
set plci 3
set hcol 6

set ndvx 505
*few usefull often used text symbol
* mostly page 122 and 123 paw user guide
** page
Percent = "Y#
Bracket = "P#
** greek letter 
Phi0     = [F]?0!
omega    = [w]

skip  = no
wwait = no
version = new
*format = postcript/encapsulated/gif
format  = gif
*mode    = mupair
*mode = singlemu
*mode = singlek
*mode = singlepi
mode = generic
dir = ~/html/fastsim/fastparam/TheDoc/ModelPlot
*** moose directory
if ([mode].eq.mupair.or.[mode].eq.singlemu) then
lumoose   = //lun2/4.
lufast    = //lun1/2.
moosehisto  = //lun2
fasthisto  = //lun1
elseif ([mode].eq.singlek) then
lumoose   = //lun2/6.
lufast    = //lun1/4.
moosehisto  = //lun2
fasthisto  = //lun1
else
lumoose   = //lun2/5.
lufast    = //lun1/3.
*lufast    = //lun1/5.
moosehisto  = //lun2
fasthisto  = //lun1
endif
*h/file 1 ~/sim_13.5.0/workdir/fast_generic_rescaled.hbook
h/file 1 ~/sim_13.5.0/workdir/test.hbook
h/file 2 ~/newmoose/workdir/generic_moose.hbook
*h/file 3 ~/sim_13.5.0//workdir/fast_generic_rescaled.hbook
*single pions
*h/file 2 ~/moose/workdir/sing_pi+_50-500Mev.hbook
*h/file 1 ~/sim/workdir/single_pi-_0.05-0.5.hbook
*h/file 1 ~/sim/workdir/sp500_new.hbook
*h/file 1 ~/sim/workdir/singlepi_new.hbook
*h/file 1 ~/sim/workdir/singlepi.hbook
**h/file 1 ~/sim/workdir/singlepi_test1.hbook
*h/file 2 ~/moose/workdir/sing_pi+_500Mev.hbook
**h/file 2 ~/moose/workdir/singlePI_notrigger.hbook
*mu pair
*h/file 2 ~/moose/workdir/mupair_notrigger.hbook
*h/file 1 ~/sim/workdir/mup_new.hbook
*h/file 1 ~/sim/workdir/mumu.hbook
*single mu
*h/file 1 ~/sim/workdir/singlek_test1.hbook
*h/file 2 ~/moose/workdir/singlek_notrigger.hbook
*h/file 2 ~/moose/workdir/single_mu_notrigger.hbook
*h/file 2 ~/moose/workdir/single_mu_dch.hbook
*h/file 1 ~/sim/workdir/singlemu_test1.hbook
*h/file 1 ~/sim/workdir/fast-param_test.hbook
*h/file 1 ~/sim/workdir/fast-param-low-mu.hbook
*h/file 2 ~/moose/workdir/single_mu_dch_low.hbook
*h/file 1 ~/sim/workdir/fast-param-low-mu-cor.hbook
*h/file 1 ~/sim/workdir/fast-param-low-mu-nocor.hbook
*apply different cuts to the two ntuple
*Ntuple/Cuts $1 chisq/nActive<20
*Ntuple/Cuts $1  abs(pt)>0.0
Ntuple/Cuts $1  abs(pt)>0.0.and.abs(sin(atan(tanDip)))<2.
*Ntuple/Cuts $1  abs(omegares)<0.2.and.abs(sin(atan(tanDip)))<2.
*.and.nActive>10
Ntuple/Cuts $2  abs(pt)>0.0
Ntuple/Cuts $10 abs(d0)<1.
Ntuple/Cuts $31 abs(omega)<0.01
Ntuple/Cuts $32 abs(curv0)<0.01
Ntuple/Cuts $42 abs(curv0/0.00453)<20.


if ([mode].eq.mupair.or.[mode].eq.singlemu.or.[mode].eq.singlek.or.[mode].eq.singlepi) then
Ntuple/Cuts $51 abs(pt)<6.
Ntuple/Cuts $41 abs(1/pt)<1.5
Ntuple/Cuts $70 abs(sigmad0)<0.015
Ntuple/Cuts $71 sigmaphi0<0.01
Ntuple/Cuts $72 sigmaomega<0.00005
Ntuple/Cuts $73 sigmaz0<0.06
Ntuple/Cuts $74 sigmatanDip<0.01
Ntuple/Cuts $75 abs(d0res)<5
Ntuple/Cuts $76 abs(phi0res)<5
Ntuple/Cuts $77 abs(omegares)<5.
Ntuple/Cuts $78 abs(z0res)<5.
Ntuple/Cuts $79 abs(tanDipres)<5

Ntuple/Cuts $80 abs(d0res*sigmad0)<0.05
Ntuple/Cuts $81 abs(phi0res*sigmaphi0)<0.02
Ntuple/Cuts $82 abs(omegares*sigmaomega)<0.00025
Ntuple/Cuts $83 abs(z0res*sigmaz0)<0.1
Ntuple/Cuts $84 abs(tanDipres*sigmatanDip)<0.05

else
Ntuple/Cuts $51 abs(pt)<3.
Ntuple/Cuts $41 abs(1/pt)<3.

Ntuple/Cuts $70 abs(sigmad0)<0.16
Ntuple/Cuts $71 sigmaphi0<0.1
Ntuple/Cuts $72 sigmaomega<0.0004
Ntuple/Cuts $73 sigmaz0<1.
Ntuple/Cuts $74 sigmatanDip<.1
Ntuple/Cuts $75 abs(d0res)<5
Ntuple/Cuts $76 abs(phi0res)<5
Ntuple/Cuts $77 abs(omegares)<5.
Ntuple/Cuts $78 abs(z0res)<5.
Ntuple/Cuts $79 abs(tanDipres)<5

Ntuple/Cuts $80 abs(d0res*sigmad0)<0.25
Ntuple/Cuts $81 abs(phi0res*sigmaphi0)<0.5
Ntuple/Cuts $82 abs(omegares*sigmaomega)<0.0025
Ntuple/Cuts $83 abs(z0res*sigmaz0)<0.5
Ntuple/Cuts $84 abs(tanDipres*sigmatanDip)<.5

endif

Ntuple/Cuts $52 abs(0.00453/curv0)<6.
Ntuple/Cuts $60 abs(tanDip)<10.
Ntuple/Cut  $20 z0<7.and.z0>-5
* abs(d0)
if ([format].eq.postscript) then
f/close 09
f/file 9 ~mangeol/html/fastsim/fastparam/TheDoc/fast-param_[mode]_[version].ps
meta 9 -111
endif

if ([skip].eq.no) then

if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/d0_bb[version].eps
meta 10 -113
endif
*zone 2 4
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose](d0) $1.and.$10  100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose](d0) $1.and.$10 100000   1 100001   'ES' 1000000
*itx 0 1500 [Percent]
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast](d0) $1.and.$2.and.$10 100000   1 100001   'E1S' 1000000
set hcol 6
if ([format].eq.gif) then 
pict/print [dir]/d0_bb[version].gif
endif

if ([wwait].eq.yes) then 
wait 
endif
* phi0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/phi0_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]phi0 $1 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]phi0 $1 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3 
Ntuple/Plot [lufast]phi0 $1.and.$2 100000   1 100001   'E1S' 1000000
set hcol 6
if ([format].eq.gif) then 
pict/print [dir]/phi0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* omega
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/omega_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]omega $1.and.$31 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]omega $1.and.$31 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]omega $1.and.$2.and.$31 100000   1 100001   'E1S' 1000000
set hcol 6
if ([format].eq.gif) then 
pict/print [dir]/omega_bb[version].gif
endif

if ([wwait].eq.yes) then 
wait 
endif
* z0
if ([format].eq.encapsulated) then

f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/z0_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]z0 $1.and.$20 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]z0 $1.and.$20 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]z0 $1.and.$2.and.$20 100000   1 100001   'E1S' 1000000
set hcol 6
if ([format].eq.gif) then 
pict/print [dir]/z0_bb[version].gif
endif

if ([wwait].eq.yes) then 
wait 
endif
* tanDip
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/tanDip_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]tanDip $1.and.$60 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]tanDip $1.and.$60 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]tanDip $1.and.$2.and.$60 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/tanDip_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* tanDip
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/costheta_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]sin(atan(tanDip)) $1.and.$60 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]sin(atan(tanDip)) $1.and.$60 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]sin(atan(tanDip)) $1.and.$2.and.$60 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/costheta_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
*** pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/pt_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]abs(pt) $1.and.$51 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]abs(pt) $1.and.$51 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]abs(pt) $1.and.$2.and.$51 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/pt_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
*** 1/pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/invpt_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]1/pt $1.and.$41 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]1/pt $1.and.$41 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]1/pt $1.and.$2.and.$41 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/invpt_bb[version].gif
endif

******* hits *****************

if ([wwait].eq.yes) then 
wait 
endif
** active hits
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/active_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]nActive $1 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]nActive $1 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]nActive $1.and.$2 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/active_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
** svt hits
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/svt_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]nSvt $1 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]nSvt $1 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]nSvt $1.and.$2 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/svt_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
** dch hits
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/dch_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]nDch $1 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]nDch $1 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]nDch $1.and.$2 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/dch_bb[version].gif
endif



if ([wwait].eq.yes) then 
wait 
endif
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dd0_bb[version].eps
meta 10 -113
endif
* abs(d0)
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]abs(sigmad0) $1.and.$2.and.$70 100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]abs(sigmad0) $1.and.$70  100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose]abs(sigmad0) $1.and.$70 100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Sigmad0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* phi0
if ([format].eq.encapsulated) then

f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dphi0_bb[version].eps
meta 10 -113
endif
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]sigmaphi0 $1.and.$2.and.$71 100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]sigmaphi0 $1.and.$71 100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose]sigmaphi0 $1.and.$71  100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Sigmaphi0_bb[version].gif
endif



if ([wwait].eq.yes) then 
wait 
endif
* omega
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Domega_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]sigmaomega $1.and.$72 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]sigmaomega $1.and.$72 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3 
Ntuple/Plot [lufast]sigmaomega $1.and.$2.and.$72 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/Sigmaomega_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* z0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dz0_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]sigmaz0 $1.and.$73 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]sigmaz0 $1.and.$73 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]sigmaz0 $1.and.$2.and.$73 100000   1 100001   'E1S' 1000000
set hcol 6


if ([format].eq.gif) then 
pict/print [dir]/Sigmaz0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* tanDip
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/DtanDip_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]sigmatanDip $1.and.$74 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]sigmatanDip $1.and.$74 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]sigmatanDip $1.and.$2.and.$74 100000   1 100001   'E1S' 1000000
set hcol 6


if ([format].eq.gif) then 
pict/print [dir]/SigmatanDip_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dd0_bb[version].eps
meta 10 -113
endif
* abs(d0)
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]d0res $1.and.$2.and.$75 100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]d0res $1.and.$75  100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose]d0res $1.and.$75 100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Pulld0_bb[version].gif
endif



if ([wwait].eq.yes) then 
wait 
endif
* phi0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dphi0_bb[version].eps
meta 10 -113
endif
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]phi0res $1.and.$2.and.$76  100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]phi0res $1.and.$76 100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose]phi0res $1.and.$76  100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Pullphi0_bb[version].gif
endif



if ([wwait].eq.yes) then 
wait 
endif
* omega
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Domega_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]omegares $1.and.$77 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]omegares $1.and.$77 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]omegares $1.and.$2.and.$77 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/Pullomega_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* z0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dz0_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]z0res $1.and.$78 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]z0res $1.and.$78 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]z0res $1.and.$2.and.$78 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/Pullz0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* tanDip
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/DtanDip_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]tanDipres $1.and.$79 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]tanDipres $1.and.$79 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]tanDipres $1.and.$2.and.$79 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/PulltanDip_bb[version].gif
endif

if ([wwait].eq.yes) then 
wait 
endif
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dd0_bb[version].eps
meta 10 -113
endif
* abs(d0)
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast](d0res*sigmad0) $1.and.$2.and.$80 100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose](d0res*sigmad0) $1.and.$80  100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose](d0res*sigmad0) $1.and.$80 100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Dd0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* phi0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dphi0_bb[version].eps
meta 10 -113
endif
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]phi0res*sigmaphi0 $1.and.$2.and.$81  100000   1 100001   'E1' 1000000
set hcol 6
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]phi0res*sigmaphi0 $1.and.$81 100000   1 100001   'S' 1000000
set dmod 3
Ntuple/Plot [lumoose]phi0res*sigmaphi0 $1.and.$81  100000   1 100001   'ES' 1000000


if ([format].eq.gif) then 
pict/print [dir]/Dphi0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* omega
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Domega_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]omegares*sigmaomega $1.and.$82 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]omegares*sigmaomega $1.and.$82 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]omegares*sigmaomega $1.and.$2.and.$82 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/Domega_bb[version].gif
endif



if ([wwait].eq.yes) then 
wait 
endif
* z0
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/Dz0_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]z0res*sigmaz0 $1.and.$83 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]z0res*sigmaz0 $1.and.$83 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]z0res*sigmaz0 $1.and.$2.and.$83 100000   1 100001   'E1S' 1000000
set hcol 6


if ([format].eq.gif) then 
pict/print [dir]/Dz0_bb[version].gif
endif


if ([wwait].eq.yes) then 
wait 
endif
* tanDip
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/DtanDip_bb[version].eps
meta 10 -113
endif
igset mtyp 0
set dmod 1
Ntuple/Plot [lumoose]tanDipres*sigmatanDip $1.and.$84 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose]tanDipres*sigmatanDip $1.and.$84 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast]tanDipres*sigmatanDip $1.and.$2.and.$84 100000   1 100001   'E1S' 1000000
set hcol 6
if ([format].eq.gif) then 
pict/print [dir]/DtanDip_bb[version].gif
endif
if ([wwait].eq.yes) then 
wait 
endif
endif

************************************* truc vs pt *********************************
*********************************************************************************
Det = DCH

******************* Sigma Doca vs Pt
TheName = SDocPt
ToPlot =  [Det]pulls/121

if ([format].eq.encapsulated) then

f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif

TheName = SZ0Pt
ToPlot =  [Det]pulls/122
******************* Sigma Z0 vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif
TheName = STDipPt
ToPlot =  [Det]pulls/123
******************* Sigma Tan(Dip) vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif
TheName = SOmPt
ToPlot =  [Det]pulls/124
******************* Sigma Omega vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************

if ([wwait].eq.yes) then 
wait 
endif
TheName = SPhiPt
ToPlot =  [Det]pulls/125
******************* Sigma Phi vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
*********************************************************************************
Det = SVT
if ([wwait].eq.yes) then 
wait 
endif
TheName = SDocPt
ToPlot =  [Det]pulls/121
******************* Sigma Doca vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif
TheName = SZ0Pt
ToPlot =  [Det]pulls/122
******************* Sigma Z0 vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif
TheName = STDipPt
ToPlot =  [Det]pulls/123
******************* Sigma Tan(Dip) vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************
if ([wwait].eq.yes) then 
wait 
endif
TheName = SOmPt
ToPlot =  [Det]pulls/124
******************* Sigma Omega vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************

if ([wwait].eq.yes) then 
wait 
endif
TheName = SPhiPt
ToPlot =  [Det]pulls/125
******************* Sigma Phi vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif
igset mtyp 21
set dmod 1
set hcol 6
set pmci 6
set plci 6
h/Plot [moosehisto]/trkqamerge/[ToPlot] L 
set dmod 1
igset mtyp 20
set hcol 3
set pmci 3
set plci 3
h/Plot [fasthisto]/trkqamerge/[ToPlot] PS

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_[Det]_bb[version].gif
endif
******************************************************

***************************************************************************************
***************************************************************************************
*****                                    Matrix elements                          *****
***************************************************************************************
***************************************************************************************
*   1   d0
*   2   phi0
*   3   omega
*   4   z0
*   5   tanDip 
**************************************************************************************
* will put M12, M13, M14, M15, M23, M24, M25, M34, M45
*
op logy
set xval 0.8
do i=1,4
do j=[i]+1,5
if ([wwait].eq.yes) then 
wait 
endif
TheName = M[i][j]
Ntuple/Cuts $3 abs(M[i][j])<0.002
******************* Sigma Phi vs Pt
if ([format].eq.encapsulated) then
f/close 10
f/file 10 ~mangeol/html/fastsim/fastparam/[TheName]_[Det]_bb[version].eps
meta 10 -113
endif


igset mtyp 0
set dmod 1
set hcol 6
Ntuple/Plot [lumoose][TheName] $1.and.$3 100000   1 100001   '' 1000000
set dmod 3
Ntuple/Plot [lumoose][TheName] $1.and.$3 100000   1 100001   'ES' 1000000
set dmod 1
igset mtyp 20
set hcol 3
Ntuple/Plot [lufast][TheName] $1.and.$2.and.$3 100000   1 100001   'E1S' 1000000
set hcol 6

if ([format].eq.gif) then 
pict/print [dir]/[TheName]_bb[version].gif
endif
enddo
enddo
op liny
set xval 0.4
f/close 10
f/close 09