Return to phiplots_rear_normalized.kumac CVS log

Up to [HallC] / geant_gep

New kumacs for front and rear conetest plots

vec/create phi1e(180)
vec/create phi2e(180)
vec/read phi1e phi3_ejb.dat
vec/read phi2e phi4_ejb.dat
vec/create phi1b(180)
vec/create phi2b(180)
vec/read phi1b phi3_deblois.dat
vec/read phi2b phi4_deblois.dat

nbins=180
1d 5001 ' ' [nbins] -50 450
1d 5002 ' ' [nbins] -50 450
1d 6000 ' ' [nbins] -50 450
1d 7001 ' ' [nbins] -50 450
1d 7002 ' ' [nbins] -50 450
1d 8000 ' ' [nbins] -50 450

hist/put_vec/contents 5001 phi1e
hist/put_vec/contents 5002 phi2e
add 5001 5002 6000 1.0 1.0
normfac=$hinfo(5001,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
1d 5003 ' ' [nbins] -50 450
add 5001 5001 5003 [nfac] 0.0
hist/set/max 5003 1.02
hist/set/min 5003 0.98
1d 5004 ' ' [nbins] -50 450
normfac=$hinfo(5002,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
add 5002 5002 5004 [nfac] 0.0
hist/set/min 5004 0.75
hist/set/max 5004 1.25
normfac=$hinfo(6000,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
1d 6001 ' ' [nbins] -50 450
add 6000 6000 6001 [nfac] 0.0
hist/set/max 6001 1.02
hist/set/min 6001 0.98

hist/put_vec/contents 7001 phi1b
hist/put_vec/contents 7002 phi2b
add 7001 7002 8000 1.0 1.0
normfac=$hinfo(7001,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
1d 7003 ' ' [nbins] -50 450
add 7001 7001 7003 [nfac] 0.0
hist/set/max 7003 1.02
hist/set/min 7003 0.98
1d 7004 ' ' [nbins] -50 450
normfac=$hinfo(7002,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
add 7002 7002 7004 [nfac] 0.0
hist/set/min 7004 0.75
hist/set/max 7004 1.25
normfac=$hinfo(8000,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*360.0/500.0
1d 8001 ' ' [nbins] -50 450
add 8000 8000 8001 [nfac] 0.0
hist/set/max 8001 1.02
hist/set/min 8001 0.98

opt nstat
opt fit
fort/file 66 phiplots_rear_normalized.ps
meta 66

vec/create x(45)
vec/create y(45)
vec/create ex(45)
vec/create ey(45)
1d 5005 ' ' 45 -50 450
1d 5006 ' ' 45 -50 450
1d 6002 ' ' 45 -50 450
rebin 5003 x y ex ey 45 1 180
hist/put_vec/contents 5005 y
rebin 5004 x y ex ey 45 1 180
hist/put_vec/contents 5006 y
rebin 6001 x y ex ey 45 1 180
hist/put_vec/contents 6002 y
hist/set/min 5005 0.98
hist/set/max 5005 1.02
hist/set/min 5006 0.75
hist/set/max 5006 1.25
hist/set/min 6002 0.98
hist/set/max 6002 1.02


vec/create x(45)
vec/create y(45)
vec/create ex(45)
vec/create ey(45)
1d 7005 ' ' 45 -50 450
1d 7006 ' ' 45 -50 450
1d 8002 ' ' 45 -50 450
rebin 7003 x y ex ey 45 1 180
hist/put_vec/contents 7005 y
rebin 7004 x y ex ey 45 1 180
hist/put_vec/contents 7006 y
rebin 8001 x y ex ey 45 1 180
hist/put_vec/contents 8002 y
hist/set/min 7005 0.98
hist/set/max 7005 1.02
hist/set/min 7006 0.75
hist/set/max 7006 1.25
hist/set/min 8002 0.98
hist/set/max 8002 1.02

zone 2 3
v/create par1(5)
v/create errpar1(5)
hist/fit 5005(15.:345.) funq.for ' E' 5 par1 ! ! ! errpar1
v/create par4(5)
v/create errpar4(5)
hist/fit 7005(15.:345.) funq.for ' E' 5 par4 ! ! ! errpar4

v/create par2(5)
v/create errpar2(5)
hist/fit 5006(15.:345.) funq.for ' E' 5 par2 ! ! ! errpar2
v/create par5(5)
v/create errpar5(5)
hist/fit 7006(15.:345.) funq.for ' E' 5 par5 ! ! ! errpar5

v/create par3(5)
v/create errpar3(5)
hist/fit 6002(15.:345.) funq.for ' E' 5 par3 ! ! ! errpar3
v/create par6(5)
v/create errpar6(5)
hist/fit 8002(15.:345.) funq.for ' E' 5 par6 ! ! ! errpar6

close 66
shell 'ps2pdf phiplots_rear_normalized.ps'