Return to
h_reconstruction.f
CVS log
Up to [HallC] / Analyzer / HTRACKING
|
Return to
h_reconstruction.f
CVS log
|
Up to [HallC] / Analyzer / HTRACKING
|
|
File: [HallC] / Analyzer / HTRACKING / h_reconstruction.f
(download)
Revision: 1.13.24.9, Mon Apr 27 20:11:34 2009 UTC (15 years, 4 months ago) by puckett Branch: gep_online Changes since 1.13.24.8: +12 -4 lines latest updates |
SUBROUTINE H_reconstruction(ABORT,err)
*--------------------------------------------------------
*- Prototype C analysis routine
*-
*-
*- Purpose and Methods : reconstruction of HMS quantities
*-
*- Output: ABORT - success or failure
*- : err - reason for failure, if any
*-
*- Created 8-Nov-1993 Kevin B. Beard, HU
*- Modified 20-Nov-1993 KBB for new errors
* $Log: h_reconstruction.f,v $
* Revision 1.13.24.9 2009/04/27 21:11:34 puckett
* latest updates
*
* Revision 1.13.24.8 2008/10/09 15:53:29 puckett
* added best track selection method for FPP1 and FPP2
*
* Revision 1.13.24.7 2008/09/29 16:01:20 puckett
* added checking for existence of histograms, prevents analyzer crash when histograms are turned off to save memory
*
* Revision 1.13.24.6 2007/11/02 19:52:52 cdaq
* Implementation of simple "constant drift velocity" type model. (ejb)
*
* Revision 1.13.24.5 2007/11/01 19:14:51 cdaq
* added wire--HMS track distance to FPP Ntuple
*
* Revision 1.13.24.4 2007/10/22 18:14:53 cdaq
* commented out print statement
*
* Revision 1.13.24.3 2007/10/17 19:38:50 cdaq
* FPP fixes
*
* Revision 1.13.24.2 2007/09/12 14:40:03 brash
* *** empty log message ***
*
* Revision 1.13.24.1 2007/08/22 19:09:30 frw
* added FPP
*
*
* Revision 1.20 2004/04/26 19:53:33 frw
* inserted calls for FPP
*
* Revision 1.13 2002/12/20 21:53:33 jones
* Modified by Hamlet for new HMS aerogel
*
* Revision 1.13 2002/09/30
* (Hamlet) Add call HMS Aerogel (Took cp from Vardan)
*
* Revision 1.12 1996/08/30 20:33:42 saw
* (DVW?) Add hbypass_track_eff
*
* Revision 1.11 1995/10/10 17:33:45 cdaq
* (JRA) Don't make an error just because no track is found
*
* Revision 1.10 1995/08/31 14:46:06 cdaq
* (JRA) Add call to h_trans_cer
*
* Revision 1.9 1995/05/22 19:39:25 cdaq
* (SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts"
*
* Revision 1.8 1995/05/11 19:05:48 cdaq
* (JRA) Add call to h_trans_misc
*
* Revision 1.7 1995/02/02 13:06:13 cdaq
* (SAW) Add call to h_select_best_track
*
* Revision 1.6 1994/06/06 16:49:49 cdaq
* (DFG) add h_recon_num and bypass switches
*
* Revision 1.5 1994/05/12 21:18:13 cdaq
* (DFG) Put h_prt_track_tests here. Remove from h_tof
*
* Revision 1.4 1994/04/13 16:06:00 cdaq
* (DFG) Add consolidated call to h_raw_dump_all
* Commented out returns after ABORT's
*
* Revision 1.3 1994/02/22 15:51:33 cdaq
* (DFG) Replace with real version
* (SAW) Move to TRACKING directory
*
* Revision 1.2 1994/02/04 20:49:31 cdaq
* Print out some raw hit data
*
* Revision 1.1 1994/02/04 20:47:59 cdaq
* Initial revision
*
*-
*- All standards are from "Proposal for Hall C Analysis Software
*- Vade Mecum, Draft 1.0" by D.F.Geesamn and S.Wood, 7 May 1993
*-
*--------------------------------------------------------
IMPLICIT NONE
SAVE
*
character*16 here
parameter (here= 'H_reconstruction')
*
logical ABORT
character*(*) err
*
INCLUDE 'hms_data_structures.cmn'
INCLUDE 'gen_constants.par'
INCLUDE 'gen_units.par'
include 'hms_scin_parms.cmn'
include 'hms_bypass_switches.cmn'
include 'hms_statistics.cmn'
include 'gep_data_structures.cmn'
*
* Local variables
integer*4 istat
*--------------------------------------------------------
*
ccc ABORT= .TRUE.
ccc err= ':no events analyzed!'
*
* increment reconstructed number
c h_recon_num= h_recon_num + 1
*
* dump
call h_raw_dump_all(ABORT,err)
if(ABORT) then
call g_add_path(here,err)
return
endif ! end test on h_raw_dump_all ABORT
*
* TRANSLATE SCINTILATORS AND CALCULATE START TIME
* HMS_RAW_SCIN ====> HMS_DECODED_SCIN
*
If(hbypass_trans_scin.eq.0) then
call H_TRANS_SCIN(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on SCIN ABORT
endif ! end test on hbypass_trans_scin
*
* TRANSLATE HMISC TDC HITS.
* H_RAW_MISC ====> HMS_DECODED_MISC
*
If(hbypass_trans_scin.eq.0) then
call H_TRANS_MISC(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on SCIN ABORT
endif ! end test on hbypass_trans_scin
*
* TRANSLATE CERENKOV
* HMS_RAW_CER ====> HMS_DECODED_CER
*
If(hbypass_trans_cer.eq.0) then
call H_TRANS_CER(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on CER ABORT
endif ! end test on hbypass_trans_cer
*
* TRANSLATE FPP
* HMS_RAW_FPP ====> HMS_FPP_event
*
c write(*,*) 'ENTERING H_TRANS_FPP'
If(hbypass_trans_fpp.eq.0.or.hbypass_trans_fpp.eq.4) then
call h_trans_fpp(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on FPP ABORT
endif
c write(*,*) 'h_trans_fpp successful'
*
* Next Aerogel Cerenkov information
* HMS_DECODED_AERO====> HMS_TRACK_TESTS
*
if(hbypass_haero.eq.0) then
call H_AERO(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test of H_AERO ABORT
endif ! end test on hbypass_aero
*
* TRANSLATE CALORIMETER
* HMS_RAW_CAL ====> HMS_DECODED_CAL
*
if(hbypass_trans_cal.eq.0) then
call H_TRANS_CAL(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on CAL ABORT
endif ! end test on hbypass_trans_cal
*
* TRANLATE DRIFT CHAMBERS
* HMS_RAW_DC + HMS_DECODED_SCIN ====> HMS_DECODED_DC
if(hbypass_trans_dc.eq.0) then
call H_TRANS_DC(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
return
endif ! end test on H_TRANS_DC ABORT
endif ! end test on hbypass_trans_dc
c write(*,*) 'h_trans_dc successful'
if(hbypass_track.eq.0) then
call H_TRACK(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
return
endif ! end test on H_TRACK ABORT
c write(*,*) 'h_track successful'
if(hbypass_track_eff.eq.0) then
call h_track_tests
endif ! end test on hbypass_trackeff
c write(*,*) 'h_track_tests successful'
endif ! end test on hbypass_track
* only proceed if the number of tracks is greater than one
*
c write(*,*) 'h_track successful'
if(HNTRACKS_FP .lt. 1) then
c don't want error message every time a track is not found.
c ABORT=.FALSE.
c err=":no tracks found!"
else
c print *,' HMS track found!'
* Proceed if one or more track has been found
*
* Project tracks back to target
* HMS_FOCAL_PLANE ====> HMS_TARGET
*
if(hbypass_targ_trans.eq. 0) then
call H_TARG_TRANS(ABORT,err,istat)
if(ABORT) then
call G_add_path(here,err)
return
endif ! end test on H_TARG_TRANS ABORT
endif ! end test on hbypass_target_trans
c write(*,*) 'h_targ_trans successful'
*
* Now begin to process particle identification information
* First scintillator and time of flight
* HMS_RAW_SCIN ====> HMS_TRACK_TESTS
*
if(hbypass_tof.eq.0) then
call H_TOF(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
return
endif ! end test of H_TOF ABORT
endif ! end test on hbypass_tof
* Next Calorimeter information
* HMS_DECODED_CAL ====> HMS_TRACK_TESTS
*
c write(*,*) 'h_tof successful'
if(hbypass_cal.eq.0) then
call H_CAL(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test of H_CAL ABORT
endif ! end test on hbypass_cal
* Next Cerenkov information
* HMS_DECODED_CER ====> HMS_TRACK_TESTS
*
if(hbypass_cer.eq.0) then
call H_CER(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test of H_CER ABORT
endif ! end test on hbypass_cer
*
* Dump HMS_TRACK_TESTS if hdebugprinttracktests is set
if( hdebugprinttracktests .ne. 0 ) then
call h_prt_track_tests
endif
c write(*,*) 'h_prt_track_tests successful'
*
* Combine results in HMS physics analysis
* HMS_TARGET + HMS_TRACK_TESTS ====> HMS_PHYSICS
*
if(hbypass_track.eq.0) then
call h_select_best_track(abort,err)
if(ABORT) then
call G_add_path(here,err)
return
endif
endif
c write(*,*) 'h_select_best_track successful'
*
if(hbypass_physics.eq.0) then
call h_physics(abort,err)
if(ABORT) then
call G_add_path(here,err)
return
endif ! end test of H_PHYSICS ABORT
endif ! end test on hbypass_physics
*
* process FPP information
*
c write(*,*) 'All HMS reconstruction successful'
if(hbypass_trans_fpp.eq.2) then !option to use HMS tracks to fill FPP variables
call h_trans_fpp_hms(ABORT,err)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test on FPP ABORT
endif ! end test on hbypass_trans_fpp
c write(*,*) 'h_trans_fpp_hms successful'
if(hbypass_fpp.eq.0) then
call h_fpp(ABORT,err) ! now includes "best track" selection for each FPP (AJP)
if(ABORT) then
call G_add_path(here,err)
* return
endif ! end test of h_fpp ABORT
endif ! end test on hbypass_fpp
c write(*,*) 'h_fpp successful'
*
endif ! end test no tracks found
* * fill FPP histogramms even if no HMS track
if(hbypass_fpp.eq.0) then
call h_fill_fpp(ABORT,err)
if (ABORT) then
call g_add_path(here,err)
return
endif
endif ! end test on hbypass_fpp
c write(*,*) 'h_fill_fpp successful'
*
* Successful return
ABORT=.FALSE.
RETURN
END
| Analyzer/Replay: Mark Jones, Documents: Stephen Wood |
Powered by ViewCVS 0.9.2-cvsgraph-1.4.0 |