*======================================================================= subroutine s_init_cal(abort,errmsg) *======================================================================= *- *- Purpose: SOS Calorimeter Initialization *- *- Created: 20 Mar 1994 Tsolak A. Amatuni * $Log: s_init_cal.f,v $ * Revision 1.3 1995/05/22 19:45:40 cdaq * (SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts" * * Revision 1.2 1994/06/14 04:30:27 cdaq * (DFG) Remove hardwired parameters * * Revision 1.1 1994/04/13 18:18:40 cdaq * Initial revision * *----------------------------------------------------------------------- * * implicit none save * logical abort character*(*) errmsg character*10 here parameter (here='S_INIT_CAL') * integer*4 block !Block number integer*4 row !Row number integer*4 column !Column number * include 'sos_data_structures.cmn' include 'sos_calorimeter.cmn' * *-----Initialize the positions do column=1,smax_cal_columns do row=1,smax_cal_rows block=row+smax_cal_rows*(column-1) * if(column.eq.1) then scal_block_xc(block)=scal_1pr_top(row)+0.5*scal_1pr_thick scal_block_yc(block)=0.5*(scal_1pr_left+scal_1pr_right) scal_block_zc(block)=scal_1pr_zpos+0.5*scal_1pr_thick else if(column.eq.2) then scal_block_xc(block)=scal_2ta_top(row)+0.5*scal_2ta_thick scal_block_yc(block)=0.5*(scal_2ta_left+scal_2ta_right) scal_block_zc(block)=scal_2ta_zpos+0.5*scal_2ta_thick else if(column.eq.3) then scal_block_xc(block)=scal_3ta_top(row)+0.5*scal_3ta_thick scal_block_yc(block)=0.5*(scal_3ta_left+scal_3ta_right) scal_block_zc(block)=scal_3ta_zpos+0.5*scal_3ta_thick else scal_block_xc(block)=scal_4ta_top(row)+0.5*scal_4ta_thick scal_block_yc(block)=0.5*(scal_4ta_left+scal_4ta_right) scal_block_zc(block)=scal_4ta_zpos+0.5*scal_4ta_thick endif enddo !End loop over rows enddo !End loop over columns * * scal_block_xsize= scal_4ta_top(2) - scal_4ta_top(1) scal_block_ysize= scal_4ta_left - scal_4ta_right scal_block_zsize= scal_4ta_thick scal_xmax= scal_4ta_top(scal_4ta_nr) + scal_block_xsize scal_xmin= scal_4ta_top(1) scal_ymax= scal_4ta_left scal_ymin= scal_4ta_right scal_zmin= scal_1pr_zpos scal_zmax= scal_4ta_zpos scal_fv_xmin=scal_xmin+5. scal_fv_xmax=scal_xmax-5. scal_fv_ymin=scal_ymin+5. scal_fv_ymax=scal_ymax-5. scal_fv_zmin=scal_zmin scal_fv_zmax=scal_zmax * return end