Skip to content

Correct assignment of energy to tracks in the HMS. #14

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
vardant wants to merge 2 commits into
base: master
Choose a base branch
from
Open

Correct assignment of energy to tracks in the HMS. #14

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a767491
Move hit clusters with zero energy deposition in the HMS calorimeter
vardant Oct 15, 2014
a5893f6
Correct assignment of hit clusters in the HMS calorimeter to the
vardant Oct 15, 2014
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 5 additions & 2 deletions HTRACKING/h_clusters_cal.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ subroutine h_clusters_cal(abort,errmsg)
*- Modified 25 Mar 1994 DFG
*- Change name of print routine
* 10 Apr 1994 DFG Protect for Et=0 division
* $Log$
* $Log: h_clusters_cal.f,v $
* Revision 1.5 1999/02/03 21:13:23 saw
* Code for new Shower counter tubes
*
Expand Down Expand Up @@ -204,12 +204,15 @@ subroutine h_clusters_cal(abort,errmsg)
if(hcluster_et(nc).gt.0.) then
hcluster_xc(nc)=hcluster_xc(nc)/hcluster_et(nc)
else
hcluster_xc(nc)= -1.0 ! Set fraction negative for bad et
hcluster_xc(nc)= -999999. ! Set out of the calorimeter volume
endif
enddo
*
100 continue
if(hdbg_clusters_cal.gt.0) call h_prt_cal_clusters
c
write(88,*) hnclusters_cal
c
*
return
end
18 changes: 7 additions & 11 deletions HTRACKING/h_tracks_cal.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -67,8 +67,8 @@ subroutine h_tracks_cal(abort,errmsg)
real*4 delta_x !Distance between track & cluster in X projection
logical*4 track_in_fv

integer*4 t_nt,t_nc
real*4 t_minx,temp_x
integer*4 t_nc
real*4 t_minx

include 'hms_data_structures.cmn'
include 'hms_calorimeter.cmn'
Expand Down Expand Up @@ -113,24 +113,20 @@ subroutine h_tracks_cal(abort,errmsg)
if( (hcal_fv_test.ne.0.and.track_in_fv) .or. hcal_fv_test.eq.0) then

if(hnclusters_cal.gt.0) then
t_minx=9
t_nt = 1
t_nc = 1
t_minx=999999.
t_nc = -1
do nc=1,hnclusters_cal
delta_x=abs(xf-hcluster_xc(nc))
if(delta_x.le.(0.5*hcal_block_xsize + hcal_slop)) then
temp_x = delta_x
if(temp_x.lt.t_minx) then
t_minx = temp_x
t_nt = nt
if(delta_x.lt.t_minx) then
t_minx = delta_x
t_nc = nc
endif
! hcluster_track(nt)=nc !Track matches cluster #nc
hntracks_cal =hntracks_cal+1
endif !End ... if matched
enddo !End loop over clusters

hcluster_track(t_nt)=t_nc !Track matches cluster #nc with min deviation
hcluster_track(nt)=t_nc !Track matches cluster #nc with min deviation
endif !End ... if number of clusters > 0
endif
enddo !End loop over detector tracks
Expand Down