SndlhcMuonReco Namespace Reference

Classes
class	hough
class	MuonReco

Functions
	hit_finder (slope, intercept, box_centers, box_ds, tol=0.)
	numPlanesHit (systems, detector_ids)

Function Documentation

hit_finder()

SndlhcMuonReco.hit_finder	(		slope,
			intercept,
			box_centers,
			box_ds,
			tol = 0.
	)

 Finds hits intersected by Hough line 

Definition at line 9 of file SndlhcMuonReco.py.

def hit_finder(slope, intercept, box_centers, box_ds, tol = 0.) :
    """ Finds hits intersected by Hough line """
 
    # First check if track at center of box is within box limits
    d = np.abs(box_centers[0,:,1] - (box_centers[0,:,0]*slope + intercept))
    center_in_box = d < (box_ds[0,:,1]+tol)/2.
 
    # Now check if, assuming line is not in box at box center, the slope is large enough for line to clip the box at corner
    clips_corner = np.abs(slope) > np.abs((d - (box_ds[0,:,1]+tol)/2.)/(box_ds[0,:,0]+tol)/2.)
    
    # If either of these is true, line goes through hit:
    hit_mask = np.logical_or(center_in_box, clips_corner)
 
    # Return indices
    return np.where(hit_mask)[0]
 

numPlanesHit()

SndlhcMuonReco.numPlanesHit	(		systems,
			detector_ids
	)

Definition at line 206 of file SndlhcMuonReco.py.

def numPlanesHit(systems, detector_ids) :
    scifi_stations = []
    mufi_ds_planes = []
    mufi_us_planes = []
 
    scifi_stations.append( detector_ids[systems == 0]//1000000 )
    mufi_ds_planes.append( (detector_ids[systems == 3]%10000)//1000 )
    mufi_us_planes.append( (detector_ids[systems == 2]%10000)//1000 )
 
    return len(np.unique(scifi_stations)) + len(np.unique(mufi_ds_planes)) + len(np.unique(mufi_us_planes))