Analyzer/HTRACKING/h_dump_cal.f - view

Return to h_dump_cal.f CVS log

Up to [HallC] / Analyzer / HTRACKING

File: [HallC] / Analyzer / HTRACKING / h_dump_cal.f (download)
Revision: 1.4, Thu Jun 10 15:48:04 1999 UTC (25 years, 3 months ago) by csa
Branch: MAIN
CVS Tags: spring03, sep-26-2002, sep-25-2002, sep-24-2002, sep-09-2002, sane, pionct, online07, online04, online03, oct1199, mduality, mar-24-2003, gep_online, gep3, fpi2, emc, e01004, bigcal, baryon, aug-12-2003, apr-02-2003, HEAD
Changes since 1.3: +19 -5 lines

(JRA) Added ycal, emeas calculations

      SUBROUTINE H_DUMP_CAL(ABORT,errmsg)
*--------------------------------------------------------
*-
*-   Purpose and Methods : Analyze scintillator information for each track 
*-
*-      Required Input BANKS     HMS_CALORIMETER
*-                               GEN_DATA_STRUCTURES
*-
*-   Output: ABORT           - success or failure
*-         : err             - reason for failure, if any
*- 
* author: John Arrington
* created: 9/7/95
*
* h_dump_cal writes out the raw calorimeterinformation for the final tracks.
* This data is analyzed by independent routines to fit the gains for each
* block.
*
* $Log: h_dump_cal.f,v $
* Revision 1.4  1999/06/10 16:48:04  csa
* (JRA) Added ycal, emeas calculations
*
* Revision 1.3  1998/12/17 22:02:38  saw
* Support extra set of tubes on HMS shower counter
*
* Revision 1.2  1996/01/17 18:17:47  cdaq
* (SAW) Remove extra () pair around implied do loop in write statement
*
* Revision 1.1  1995/10/09 20:16:45  cdaq
* Initial revision
*
*--------------------------------------------------------
      IMPLICIT NONE
*
      character*10 here
      parameter (here= 'H_DUMP_CAL')
*
      logical ABORT
      character*(*) errmsg
*
      INCLUDE 'hms_data_structures.cmn'
      include 'hms_calorimeter.cmn'

      integer*4 blk
      real*4 emeas,ycal

      save

*
*  Write out cal fitting data.
*
*  What should this do for new tubes?
*
      if (abs(hsdelta).le.10 .and. hcer_npe_sum.gt.2) then

! this correction comes from h_corect_cal_pos(neg).f

	ycal=hsy_fp + hcal_1pr_zpos*hsyp_fp
        ycal=min(35.,ycal)
        ycal=max(-35.,ycal)
	emeas=hsp*exp(ycal/165.4)/(1+ycal**2/50000.)

        write(35,'(1x,52(1x,f6.1),1x,e11.4)') 
     &      (hcal_realadc_pos(blk),blk=1,hmax_cal_blocks),emeas

!        if(hcal_num_neg_columns.gt.0) then
!        write(35,'(1x,52(1x,f6.1),1x,e11.4)') 
!     &       (hcal_realadc_neg(blk),blk=1,hmax_cal_blocks),hsp
!        endif
      endif

      RETURN
      END

Analyzer/Replay: Mark Jones, Documents: Stephen Wood