Up coming news, new revised version of EC on the tree!!!
Very Mary Christmas and Happy New Year!!!
The changes since 1st of December are:
==========================================================
New subroutines and parameters -
ec_getcal_geom.F
| Checks the tcl variable def_geom, if it is "T" then the
| ideal geometry is used, if it is "F" calls ec_read_eccg
| to read geometry parameters for given Run.
ec_read_eccg.F
| Reads calibration file containing the ECG BOS bank for
| geometry and fills the geometry common blocks. Geometry
| definitions are done according to Ralph Minehart's document.
| Concerning to geometry - 3 new tcl variables are implemented
| whole_surf, inner_surf and outer_surf, that defines surface of
| what layer to be taken as hit reconstruction plane for whole,
| inner and outer hits. Note that these parameters are not necessary
| to be full layer, it con be half a layer (7.5) or so.
| New common blocks for geometry:
ecgeom.inc | Main parameters
ecgeompar.inc | Intermediate parameters.
ec_ij_path.F
ec_ijk_xyz.F
| Some changes according to new geometry definitions.
| Also in the ec_ijk_xyz determination of errors on xyz
| coordinates are added.
ec_dalitz.F
| Determination of errors on ij are added.
| Also local coordinate i now determined by positions
| of UV and VW crossings, since i_UW==i_UV. According to
| simulations it is more accurate.
Ec_general.PAR
| New parameters, errors on defined hit coordinates are added.
ec_getcal_adc.F
ec_getcal_atten.F
ec_getcal_tdc.F
| Before reading the calibration files (calls to ec_read_****)
| check the tcl variables def_adc, def_tdc and def_atten,
| if def_adc is "T" then the - EcCal_Eo=0. and EcCal_Ech=0.0001
| if def_tdc is "T" then the - EcCal_To=0. and EcCal_Tch=0.05
| if def_atten is "T" then the - EcCal_atten=376.
| If any of them is "F" - calls ec_read_****
| to read calibration parameters for given Run.
ec_read_ecca.F
ec_read_eccl.F
ec_read_ecct.F
| Initializations (setting to the initial values) are taken
| to the ec_getcal_****.
| Also sequence of closing a Fpack file and freeing the reserved
| logical unit number are changed. Before it try to free logical
| unit number and then close the file, that did not work.
ec_ini_befor.F
| Some really unnecessary initialization loops are taken out.
| That makes code 6-7 times faster.
eccal_strip.F
| Some of the initialization of initial areas are moved here,
| more time effective
ecfph.F
| Determination of some new areas for more complete strip
| information analyzes are added.
| C- stripc is 1 if strip was in that hit and -1 otherwise
| real stripc(Maxstrips,0:Maxhits,Ec_U:Ec_W,whole:outer,Maxsectors)
| C- ppath is the path length from hit position to the edge
| real ppath(0:Maxhits,Ec_U:Ec_W,whole:outer,Maxsectors)
| C- nstripc is the number of strips in the peak on the edge at given hit
| integer nstripc(0:Maxhits,Ec_U:Ec_W,whole:outer,Maxsectors)
| C- ntstripc is the total number of strips in the given hit
| integer ntstripc(0:Maxhits,whole:outer,Maxsectors)
| Note that this information is not kept in the BOS bank.
| So if you wont to use them or dump into Ntuples you should include
| ec_strip.inc in your routine.
ec_strip.inc
| New areas for strip analyzes are included in the common block
ecfsp.F
| Before the allowed distance between strips in the peak was fixed - 2
| mistake! Now it is tcl controlled touch_id, as it should be.
ectcl_init.F
ec_init.F
ec_control.inc
| All new tcl variables are added, now tcl variables for EC are:
| set EcFitEdge_strip_threshold
| set EcFitEdge_peak_threshold
| set EcFitEdge_hit_threshold
| set def_atten
| set anal_opt
| set touch_id
| set def_geom
| set def_adc
| set def_tdc
| set whole_surf
| set inner_surf
| set outer_surf
ec_fill_bos.F
| Few changes are made:
| First - new BOS bank ECPI are filled in case of anal_opt=2.
| This bank will allow to transfer the hit information to CED
| for hit display.
| ECHB also will be filled only if anal_opt=2.
| Bank ECRB will be filled only if anal_opt=1.
| There are few reasons for this separation -
| first according to present cod in ECRB only the information
| of matching hits in both inner and outer planes are stored and
| scheme is so that if there are no any matching hits in both planes
| it stops the filling of any BOS banks;
| secondly in case of anal_opt=1 half of the information for ECHB
| are not determined, it also true for ECPI;
| and finely in case of anal_opt=2 more of obtained information
| from hit analyses are not stored in the ECRB.
| Of course one always can co RCS version and modify for his own.
================================================================
New version are checked in and libraries are built on HP and SUN
platforms. Corresponding links to xltest are were new libraries
are libec00010.a.
So good luck!!
Stepan.