sbndcode/sbndcode/CosmicId/Algs/GeometryCosmicIdAlg.cc

Go to the documentation of this file.

#include "GeometryCosmicIdAlg.h"

namespace sbnd{

GeometryCosmicIdAlg::GeometryCosmicIdAlg(const Config& config){

  this->reconfigure(config);

}


GeometryCosmicIdAlg::GeometryCosmicIdAlg(){

}


GeometryCosmicIdAlg::~GeometryCosmicIdAlg(){

}


void GeometryCosmicIdAlg::reconfigure(const Config& config){

  return;
}

// Remove any tracks in different TPC to beam activity
bool GeometryCosmicIdAlg::GeometryCosmicId(recob::Track track, std::vector<art::Ptr<recob::Hit>> hits, bool tpc0Flash, bool tpc1Flash){

  // Remove any tracks that are detected in one TPC and reconstructed in another
  int tpc = fTpcGeo.DetectedInTPC(hits);
  double startX = track.Start().X();
  double endX = track.End().X();

  // Check if track is stitched
  if(tpc == -1) return false;

  // Check the start/end points are in same TPC (track shifted into other TPC because time outside of beam)
  if(tpc == 0 && (startX>0 || endX>0)) return true;
  else if(tpc == 1 && (startX<0 || endX<0)) return true;

  // Check if track was detected in TPC where there was no beam activity
  if(tpc0Flash && !tpc1Flash && tpc == 1) return true;
  if(tpc1Flash && !tpc0Flash && tpc == 0) return true;
  
  return false;
  
}


}