FPGATrackSimHitFilteringTool.cxx File Reference

#include "FPGATrackSimHitFilteringTool.h"
#include "PathResolver/PathResolver.h"
#include "geometry_constants.h"
#include "FPGATrackSimObjects/FPGATrackSimMultiTruth.h"
#include "FPGATrackSimObjects/FPGATrackSimCluster.h"
#include "FPGATrackSimClusteringTool.h"
#include <AsgMessaging/MessageCheck.h>
#include <vector>
#include <fstream>
#include <cmath>
#include <iostream>
#include <sstream>
#include <boost/algorithm/string.hpp>

Go to the source code of this file.

Functions
void	fill_nearby_hits (const FPGATrackSimHit &, std::vector< std::pair< FPGATrackSimHit, unsigned >> &, std::vector< std::pair< FPGATrackSimHit, unsigned >> &)
int	countHitsFromHeader (FPGATrackSimLogicalEventInputHeader &)
void	printHitsFromHeader (FPGATrackSimLogicalEventInputHeader &)
void	fill_nearby_hits (const FPGATrackSimHit &innerHit, std::vector< std::pair< FPGATrackSimHit, unsigned > > &outerHits, std::vector< std::pair< FPGATrackSimHit, unsigned > > &nearby_outer_hits)

Function Documentation

countHitsFromHeader()

int countHitsFromHeader

(

FPGATrackSimLogicalEventInputHeader &

header

)

Definition at line 501 of file FPGATrackSimHitFilteringTool.cxx.

{
    int nHits = 0;
    for (int i = 0; i < header.nTowers(); ++i) {
        FPGATrackSimTowerInputHeader &tower = *header.getTower(i);
        std::vector<FPGATrackSimHit>  hits = tower.hits();
        nHits += hits.size();
    }
    return nHits;
}

fill_nearby_hits() [1/2]

void fill_nearby_hits	(	const FPGATrackSimHit &	,
		std::vector< std::pair< FPGATrackSimHit, unsigned >> &	,
		std::vector< std::pair< FPGATrackSimHit, unsigned >> &
	)

fill_nearby_hits() [2/2]

void fill_nearby_hits	(	const FPGATrackSimHit &	innerHit,
		std::vector< std::pair< FPGATrackSimHit, unsigned > > &	outerHits,
		std::vector< std::pair< FPGATrackSimHit, unsigned > > &	nearby_outer_hits
	)

Definition at line 467 of file FPGATrackSimHitFilteringTool.cxx.

{
  for (const auto& outerHit : outerHits) {
    // layer needs to be inner hit physLayer + 1
    if(outerHit.first.getPhysLayer() != innerHit.getPhysLayer() + 1)
      continue;
 
    // want to do just +-1 in eta/phi as appropriate
    // but need to interpret module IDs first
 
    // for now, just do by z in barrel and rho in EC
    // strips are at most 50mm long in barrel, 60mm long in EC
    // adjacent one is at most a bit over 60mm away
 
    if(innerHit.isBarrel()) {
      if(fabs(innerHit.getZ() - outerHit.first.getZ()) > 80)
        continue;
    }
    else {
      if(fabs(innerHit.getR() - outerHit.first.getR()) > 80)
        continue;
    }
 
    // now distance in x-y plane
    if(fabs(innerHit.getX() - outerHit.first.getX()) > 80)
      continue;
    if( fabs(innerHit.getY() - outerHit.first.getY()) > 80)
      continue;
 
    nearby_outer_hits.push_back(outerHit);
  }
}

printHitsFromHeader()

void printHitsFromHeader

(

FPGATrackSimLogicalEventInputHeader &

header

)

Definition at line 513 of file FPGATrackSimHitFilteringTool.cxx.

{
    for (int i = 0; i < header.nTowers(); ++i) {
        FPGATrackSimTowerInputHeader &tower = *header.getTower(i);
        std::vector<FPGATrackSimHit>  hits = tower.hits();
        ANA_MSG_DEBUG("tower " << i << " with " << hits.size() << " hits:");
        for(const auto& hit : hits) {
          ANA_MSG_DEBUG("  isStrip, side, x, y, z, rho, phi, physLayer, layer, section, phiM, etaM = " << hit.isStrip() << ", " << hit.getPhysLayer()%2 << ", " << hit.getX() << ", " << hit.getY() << ", " << hit.getZ() << ", " << hit.getR() << ", " << hit.getGPhi() << ", " << hit.getPhysLayer() << ", " << hit.getLayer() << ", " << hit.getSection() << ", " << hit.getPhiModule() << ", " << hit.getEtaModule());
        }
    }
}