*  eep.kumac 
*  analyse GEANT simulations
* K. Aniol, May 17, 2005
*
read ikin ' enter kinematics number 1-7'
 read eloss ' do energy loss correction?, yes = 1, no = 0 '
 read i1  ' i1 = first event number to analyse'
 read i2  ' i2 = last event number to analyse'
read ntid ' ntid = ntuple id to analyse '
* nevent is the total number of events to analyse
 nevent = [i2] - [i1] + 1

*
*create histograms from the ntuple
read hid0 ' hid0 = base histogram id number for ntuple projections '
hid1 = [hid0]
hi/cr/1d [hid1] 'ztgtex' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgtey' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgtpx' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgtpy' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgte' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgtp' 400 -4. 4.
hid1 = [hid1] + 1
hi/cr/1d [hid1] 'ztgtep' 400 -4. 4.
*hid1 = [hid1]+1
*hidp = [hid1]
*hi/cr/1d [hid1] 'deltap/p smearing function ' 100 -5.e-4 5.e-4
hid1=[hid1]+1
hi/cr/1d [hid1] 'e3 excitation energy for tgt1' 600 -2.0 10.
hid1=[hid1]+1
hi/cr/1d [hid1] 'e3 excitation energy for tgt2' 600 -2.0 10.
hid1=[hid1]+1
hi/cr/1d [hid1] 'e3 excitation energy for tgt3' 600 -2.0 10.
*
ip = [hid0] + 1000
ih = [hid0] + 1001
iv = [hid0] + 1002
ie = [hid0] + 1003
ipr =[hid0] + 1004
iep =[hid0] + 1005
i12 =[hid0] + 1006
hi/cr/1d [ip] 'momentum smearing' 100 -0.0004 0.0004
hi/cr/1d [ih] ' horizontal angle smearing' 100 -0.002 0.002
hi/cr/1d [iv] ' vertical angle smearing' 100 -0.006 0.006
hi/cr/1d [ie] ' electron momentum at sp07' 500 1.7 2.1
hi/cr/1d [ipr] ' proton momentum at collimator exit' 500 0.910 1.01
hi/cr/2d [iep] 'emiss vs pmiss' 600 -2. 10. 600 -600. 600.
hi/cr/2d [i12] 'emom vs pmom' 120 1900. 2120. 120 940. 1060.
*

nt/loop [ntid] eeppb.f([hid0],[eloss],[ikin]) [nevent] [i1]
*opt sta
*hi/pl [hidp]
*hi/fit [hidp] g

zone 2 2
opt sta
*hi/pl [ip]
hi/fit [ip] g
*hi/pl [ih]
hi/fit [ih] g
*hi/pl [iv]
hi/fit [iv] g
opt nsta
hi/pl [ie]
hi/pl [ipr]
opt sta
zone 1 1
fort/close 1
vec/cr a(600) r
hi/get_vect/cont 108 a
vec/dr a