1 cdaq  1.1.2.1       REAL FUNCTION coin_hbook()
  2                     IMPLICIT NONE
  3                     INCLUDE ?
  4               
  5                     LOGICAL FIRST
  6                     DATA FIRST /.TRUE./
  7                     SAVE FIRST
  8                     real haero_tot
  9                     real wlocut,whicut,q2min,thetacm,csum
 10                     logical good_electron,good_pion
 11                     logical good_pion1
 12                     logical coint_e_pi,coint_e_pi_delta

 13 cdaq  1.1.2.2       logical mmcut_relax,good_ran

 14 cdaq  1.1.2.1 c
 15               
 16                     IF(FIRST) THEN
 17                        FIRST=.FALSE.
 18                        PRINT *,  'BOOKING HISTOS...'
 19                        CALL HBOOK1(3,'scal E1 track t=scleantrack',
 20                    ^100,0.0,1.5,0)
 21                        CALL HBOOK1(4,'scal Etot track t=scleantrack',
 22                    ^100,0.0,2.0,0)
 23                        CALL HBOOK1(5,'scer npe t=scleantrack',100,0.,15.,0)  
 24                        CALL HBOOK1(6,'sbeta t=scleantrack',100,0.5,1.2,0)
 25                        CALL HBOOK2(7,'scal x vs Etot t=scleantrack',
 26                    ^100,0.0,2.0,100,-60.0,60.0,0)        
 27                        CALL HBOOK2(8,'scal y vs Etot t=scleantrack',
 28                    ^100,0.0,2.0,100,-30.0,30.0,0)        
 29               
 30               
 31               
 32                        CALL HBOOK1(13,'hbeta t=hclean',100,-1.,1.5,0)
 33                        CALL HBOOK1(14,'HMS s1xdedx t=hclean',100,0.0,800.,0)
 34                        CALL HBOOK1(15,'HMS TOFmass2 t=hclean',100,-1.,4.0,0)
 35 cdaq  1.1.2.1 
 36                        CALL HBOOK1(22,'haero-pos-npe-sum t=hclean',100,0.0,40.,0)
 37                        CALL HBOOK1(23,'haero-neg-npe-sum t=hclean',100,0.0,40.,0)
 38                        CALL HBOOK1(25,'haero-tot-npe-sum t=hclean',100,0.0,40.,0)
 39               
 40                        CALL HBOOK1(32,'hcer adc bottom t=hcelan',100,0.0,400.,0)
 41                        CALL HBOOK1(33,'hcer adc top t=hclean',100,0.0,400.,0)
 42                        CALL HBOOK1(34,'hcer npe sum t=hclean',100,0.0,40.,0)
 43                        CALL HBOOK2(35,'xcer vs hcer_npe t=hclean',100,0.0,20.,
 44                    ^100,-70.,70.,0)
 45               
 46                       
 47                        CALL HBOOK2(50,'HMS beta vs Cointime t=e pi',100,-10.,
 48                    ^10.,100,0.8,1.2,0)
 49                        CALL HBOOK2(51,'SOS beta vs Cointime t=e pi',100,-10.,
 50                    ^10.,100,0.8,1.2,0)
 51                        CALL HBOOK1(53,'SOS-HMS path corr(ns) t=e pi',100,-2.0,2.0,0)
 52                        CALL HBOOK2(54,'HMS ztgt vs SOS ztgt t=e pi',100,-8.,8.
 53                    ^,100,-8.,8.0,0)
 54                      
 55                      CALL HBOOK1(62,'Mm(pi prod) t=coint e pi delta',100,0.0,1.5,0)
 56 cdaq  1.1.2.1        CALL HBOOK1(63,'Em t=coint e pi delta',100,-0.4,2.4,0)
 57                      CALL HBOOK1(64,'Eexc t=coint e pi delta',100,-0.3,0.3,0)
 58               
 59                        CALL HBOOK1(65,'Q2 t=mmcut relax',100,1.0,4.,0)
 60                        CALL HBOOK1(66,'W t=mmcut relax',100,1.5,3.,0)
 61                        CALL HBOOK2(67,'W vs Q2 t=mmcut relax',100,1.4,
 62                    ^3.5,100,1.9,2.6,0)
 63                        CALL HBOOK2(68,'-t vs phipi t-mmcut relax',100,
 64                    ^-0.5,0.5,100,-0.5,0.5,0)
 65                         
 66               
 67                       CALL HBOOK1(101,'t, t=mmcut relax',100,0,0.8,0)
 68                       CALL HBOOK1(102,'phipi, t=mmcut relax',100,0.,6.28,0) 
 69               
 70                        CALL HBOOK1(2001,'hsdelta t=hclean',100,-20.0,20.0,0)
 71                        CALL HBOOK1(2002,'ssdelta,t=sclean',100,-20.0,20.0,0)
 72               

 73 cdaq  1.1.2.2          CALL HBOOK1(3000,'Cointime t=e pi',100,-10.,10.,0)

 74 cdaq  1.1.2.1         ENDIF
 75                 
 76                      haero_tot=haero_po+haero_ne
 77                 
 78                        good_electron = .false.

 79 cdaq  1.1.2.2          if(scer_npe.gt.1.0.and.ssshtrk.gt.0.7)

 80 cdaq  1.1.2.1      ^good_electron = .true.
 81                        good_pion1 = .false.

 82 cdaq  1.1.2.2          if(hcer_npe.lt.0.2.and.haero_tot.gt.4.)

 83 cdaq  1.1.2.1      ^good_pion1 = .true.
 84                        good_pion = .false.
 85                        if(good_pion1.and.hsmass2.lt.1.0) good_pion= .true.
 86                
 87                        coint_e_pi = .false.
 88                     if(abs(cointime).lt.3.5.and.good_pion.and.good_electron)
 89                    ^coint_e_pi=.true.
 90               
 91                       coint_e_pi_delta = .false.
 92                     if(coint_e_pi.and.abs(hsdelta).lt.8.5.and.abs(ssdelta).lt.15.)          
 93                    ^coint_e_pi_delta= .true.  
 94                         mmcut_relax = .false.
 95                     if(coint_e_pi_delta.and.abs(missmass-0.939).lt.0.05)
 96                    ^mmcut_relax=.true.    
 97               

 98 cdaq  1.1.2.2       good_ran=.false.
 99                     if(good_pion.and.good_electron.and.abs(hsdelta).lt.8.5.and.
100                    ^abs(ssdelta).lt.15.0.and.abs(missmass-0.939).lt.0.05) then
101               
102                        good_ran=.true.
103                     endif
104                     

105 cdaq  1.1.2.1 
106                     CALL HFILL(3,ssprtrk,0.,1.)
107                     CALL HFILL(4,ssshtrk,0.,1.)
108                     CALL HFILL(5,scer_npe,0.,1.)
109                     CALL HFILL(6,ssbeta,0.,1.)
110                     CALL HFILL(7,ssshtrk,scal_x,1.0)
111                     CALL HFILL(8,ssshtrk,scal_y,1.0)
112               
113                     CALL HFILL(13,hsbeta,0.,1.)
114                     CALL HFILL(14,hsdedx1,0.,1.)
115                     CALL HFILL(15,hsmass2,0.,1.)
116                   
117                     CALL HFILL(22,haero_po,0.,1.)
118                     CALL HFILL(23,haero_ne,0.,1.)
119                     CALL HFILL(25,haero_tot,0.,1.)
120               
121                     CALL HFILL(32,hceradc1,0.,1.)
122                     CALL HFILL(33,hceradc2,0.,1.)
123                     CALL HFILL(34,hcer_npe,0.,1.)
124                     CALL HFILL(35,hcer_npe,hcersx,1.)
125                    
126 cdaq  1.1.2.1       CALL HFILL(2001,hsdelta,0.,1.)
127                     CALL HFILL(2002,ssdelta,0.,1.)
128                   
129                     if(good_pion.and.good_electron)then
130                     CALL HFILL(50, cointime,hsbeta,1.)
131                     CALL HFILL(51, cointime,ssbeta,1.)
132                     CALL HFILL(53, coinpath,0.,1.)
133                     CALL HFILL(54,sszbeam,hszbeam,1.)
134                     endif     
135               
136                     if(coint_e_pi_delta)then
137                     CALL HFILL(62, missmass,0.,1.)
138                     CALL HFILL(63, Em,0.,1.)
139                     CALL HFILL(64, Eexc,0.,1.)
140               

141 cdaq  1.1.2.2       endif
142                     
143                     if(good_ran) then
144                        CALL HFILL(3000,cointime,0.,1.)
145                     endif

146 cdaq  1.1.2.1       if(mmcut_relax)then
147                     CALL HFILL(65, Q2,0.,1.)
148                     CALL HFILL(66, W,0.,1.)
149                     CALL HFILL(67, Q2,W,1.)
150                     CALL HFILL(68,ctphix,ctphiy,1.)
151                     CALL HFILL(101,t,0.,1.)
152               c      print*, 'phi_pq=',phi_pq 
153                     CALL HFILL(102,phi_pq,0.,1.)
154               
155                     endif
156               
157                     coin_hbook = 1.
158                     return
159                     end
160               
161               
162               
163               
164               
165               
166               
167 cdaq  1.1.2.1 
168