AMSTracker.h

Go to the documentation of this file.

#pragma once
 
#include <limits>
#include <vector>
 
#include <TH1F.h>
#include <TH2D.h>
 
#include "Line.h"
#include "Point.h"
 
const double TrkL1_ZPos = 158.920;          
const double TrkL1_ZPosCutoff = 2;          
const double closest_distance_cutoff = 10;  
 
 
namespace TrackerL1
{
    //___________________________________________________________________________
    bool HasL1XYHit(const TH2D * hTrackerL1_XY, const double& x, const double& y, const bool& debug = false)
    {
        int binX = hTrackerL1_XY->GetXaxis()->FindBin(x);
        int binY = hTrackerL1_XY->GetYaxis()->FindBin(y);
        
        if (debug == true)
            printf("\n\nFor (%f, %f):  (binX, binY) = (%d, %d). Bin content = %f\n\n", x, y, binX, binY, hTrackerL1_XY->GetBinContent(binX, binY));
        
        return (hTrackerL1_XY->GetBinContent(binX, binY) > 0);
    }
 
 
    //___________________________________________________________________________
    std::pair<Point::Point2DwID, bool> FindClosestTrkL1ClusterFrom3DLine(const TH2D * hTrackerL1_XY,
                                                                         const Line::Line3DParametricForm& line, 
                                                                         const std::vector<Point::Point2DwID>& gMCClustersXY,
                                                                         const bool& debug = false)
    {
        if (!hTrackerL1_XY)
        {
            std::cerr<<"\n\nWARNING! Bad pointer for hTrackerL1_XY\n\n";
            return std::make_pair(Point::Point2DwID(), false);
        }
        if (line.direction[2] == 0)
            return std::make_pair(Point::Point2DwID(), false);
      
        double t = (TrkL1_ZPos - line.point[2]) / line.direction[2];
 
        double x_extrapolated = line.point[0] + line.direction[0] * t;
        double y_extrapolated = line.point[1] + line.direction[1] * t;
        
        if (debug == true)
            std::cout<<"\nThe extrapolated X and Y at L1 for the 3D line are:  "<<x_extrapolated<<", "<<y_extrapolated;
 
        bool has_L1_hit = HasL1XYHit(hTrackerL1_XY, x_extrapolated, y_extrapolated);
        if (!has_L1_hit)
            return std::make_pair(Point::Point2DwID(), false);
 
        Point::Point2DwID closestPoint;
        double minDistance = std::numeric_limits<double>::max(); 
        
        for (const auto& MCCluster : gMCClustersXY)
        {
            double distance = Point::GetDistance2D(x_extrapolated, y_extrapolated, MCCluster.x, MCCluster.y); 
            
            if (debug == true)
                printf("\nDistance from (%7.2f, %7.2f) = %7.2f cm.", MCCluster.x, MCCluster.y, distance);
            
            if (distance < minDistance)
            {
                minDistance = distance;  
                closestPoint = MCCluster;
            }
        }
        
        
        
        if (minDistance <= closest_distance_cutoff)
        {
            closestPoint.closedis = minDistance;
            return std::make_pair(closestPoint, true);
        }
        
        else
            return std::make_pair(Point::Point2DwID(), false);
    }
}