According to our discussion on last EGN software meeting
I made some changes in EC reconstruction code in RECSIS.
1. I wrote new subroutine for primary EC hit matching:
hit_match.F
2. Reformate ECHB BOS bank in order to include results of hit matching
as well as number of strips in the views for each hit. Now ECHB looks:
!-----------------------------------------------------------------------
! BANKname BANKtype ! Comments
TABLE ECHB ! create write display delete ! Forward calorimeter result bank
!
! ATTributes:
! -----------
!COL ATT-name FMT Min Max ! Comments
!
1 Sect I 0 0xFFFF ! Sector number & Layer number
2 E__hit F 0.0 6.0 ! energy found
3 dE_hit F 0.0 6.0 ! error on the energy found
4 t_hit F 0.0 9999.0 ! time found for
5 i_hit F 0.0 9999.0 ! sector rectangular coordinate
6 j_hit F 0.0 9999.0 ! sector rectangular coordinate
7 di_hit F 0.0 9999.0 ! sector rectangular coordinate error,
8 dj_hit F 0.0 9999.0 ! sector rectangular coordinate error,
9 x_hit F 0.0 9999.0 ! lab coordinate,
10 y_hit F 0.0 9999.0 ! lab coordinate,
11 z_hit F 0.0 9999.0 ! lab coordinate,
12 dx_hit F 0.0 9999.0 ! lab coordinate error,
13 dy_hit F 0.0 9999.0 ! lab coordinate error,
14 dz_hit F 0.0 9999.0 ! lab coordinate error,
15 u2_hit F 0.0 9999.0 ! second moment of u _hit pattern
16 v2_hit F 0.0 9999.0 ! second moment of v _hit pattern
17 w2_hit F 0.0 9999.0 ! second moment of w _hit pattern
18 u3_hit F -9999.0 9999.0 ! third moment of u _hit pattern
19 v3_hit F -9999.0 9999.0 ! third moment of v _hit pattern
20 w3_hit F -9999.0 9999.0 ! third moment of w _hit pattern
21 i2_hit F 0.0 9999.0 ! second moment of sector coordinates
22 j2_hit F 0.0 9999.0 ! second moment of sector coordinates
23 i3_hit F 0.0 9999.0 ! third moment of sector coordinates
24 j3_hit F 0.0 9999.0 ! third moment of sector coordinates
25 Nstrp_U I 0 36 ! Number of U strips in the hit
26 Nstrp_V I 0 36 ! Number of V strips in the hit
27 Nstrp_W I 0 36 ! Number of W strips in the hit
28 MatchID1 I 0 30 ! Id of matched hit in the layer1
29 CH21 F 0. 999. ! Quality measure of matching with layer1
30 MatchID2 I 0 30 ! Id of matched hit in the layer2
31 CH22 F 0. 999. ! Quality measure of matching with layer2
32 istat I 0 0xFFFF ! status word
!!
! For matching if current layer is WHOLE then layer1=INNER and layer2=OUTER
! if current layer is INNER then layer1=WHOLE and layer2=OUTER
! if current layer is OUTER then layer1=WHOLE and layer2=INNER
END TABLE
!
As ussual I crated a libraries on all platforms (SUN, HP & AIX).
It will be grate if users can run and test new EC as soon as possible.
If new bank will satisfy every one, than I think we can go to CVS with
one analyses flow and one result bank.
Below I made some studies on memory occupation of EC modules, EC alone
runs on both SUN's and JLABH1 with the same 4ms/event speed.
text data bss dec hex filename
368 132 20 520 208 ec_brun.o (ex ec.a)
1788 268 20 2076 81c ec_dalitz.o (ex ec.a)
220 132 36 388 184 ec_evnt.o (ex ec.a)
10852 452 68672 79976 13868 ec_fill_bos.o (ex ec.a)
1156 528 164 1848 738 ec_getcal_adc.o (ex ec.a)
1128 492 164 1784 6f8 ec_getcal_atten.o (ex ec.a)
3856 708 200 4764 129c ec_getcal_geom.o (ex ec.a)
1156 524 164 1844 734 ec_getcal_tdc.o (ex ec.a)
1136 132 40 1308 51c ec_ij_path.o (ex ec.a)
1388 268 84 1740 6cc ec_ijk_xyz.o (ex ec.a)
348 132 20 500 1f4 ec_ini_befor.o (ex ec.a)
3564 1304 104 4972 136c ec_init.o (ex ec.a)
4808 1296 280 6384 18f0 ec_read_ecca.o (ex ec.a)
4952 1308 272 6532 1984 ec_read_eccg.o (ex ec.a)
4892 1296 280 6468 1944 ec_read_eccl.o (ex ec.a)
4816 1296 280 6392 18f8 ec_read_ecct.o (ex ec.a)
976 132 8 1116 45c ec_sort.o (ex ec.a)
4692 272 6812 11776 2e00 ec_store_bos.o (ex ec.a)
5160 268 52 5480 1568 eccal_strip.o (ex ec.a)
384 268 20 672 2a0 ecfit_analyze.o (ex ec.a)
8636 268 29408 38312 95a8 ecfit_peaks_hits.o (ex ec.a)
4340 268 12928 17536 4480 ecfit_strips_peaks.o (ex ec.a)
10108 268 4060 14436 3864 ecfph.o (ex ec.a)
10020 268 3080 13368 3438 ecfsp.o (ex ec.a)
1460 276 4 1740 6cc ectcl_init.o (ex ec.a)
4664 404 200 5268 1494 hit_match.o (ex ec.a)
Note: ec_fill_bos is big since we filling two BOS banks.
ecfit_peaks_hits and ecfit_strips_peaks are for anal_opt=1
ecfph and ecfsp are for anal_opt=2
Regards Stepan.