BPFPId.cxx

Go to the documentation of this file.

///////////////////////////////////////////////////////////
//
// \file    BPFPId.cxx
// \brief   PId object created by BPFPId_module.cc.
// \version
// \author  Michael Baird - mbaird42@fnal.gov & Jonathan Davies jpdavies@sussex.ac.uk
//
///////////////////////////////////////////////////////////

#include "BreakPointFitter/DataProducts/BPFPId.h"

#include "canvas/Persistency/Common/FindManyP.h"
#include "cetlib_except/exception.h"
#include <iostream>

namespace bpfit
{

///////////////////////////////////////////////////////////

  BPFPId::BPFPId() : rb::PID()
  {
    SetDefaults();
  }

///////////////////////////////////////////////////////////

  BPFPId::BPFPId(int pdg, double val) : rb::PID(pdg, val) 
  {
    SetDefaults();
  }

///////////////////////////////////////////////////////////

  void BPFPId::SetDefaults()
  {
    SetLength(-5.0);
    SetChi2T(1.0e9);
    SetdEdXLL(-1.0e9);
    SetHitRatio(-5.0);
  }

///////////////////////////////////////////////////////////

  unsigned int HighestPIDTrack(
                const std::vector< art::Ptr<rb::Track> >& sliceTracks,
                const std::string& bpfpidModuleLabel, 
                const int & pdg,
                const art::Event& e)
  {
    int    bestTrack        = -1;
    double highestPID       = -2.0;
    // double bestTrackLength  = 0.0;
    // double bestStartZ       = std::numeric_limits<double>::max();
    // double bestADC          = std::numeric_limits<double>::max();

    art::FindManyP<bpfit::BPFPId> tracksToBPFPId(sliceTracks, e, bpfpidModuleLabel);

    for(size_t c = 0; c < sliceTracks.size(); ++c){
      geo::View_t view = sliceTracks[c]->View();
      //Only use 3D tracks
      if (view != geo::kXorY){  
        continue;
      } //End of loop over 2d tracks

      if(tracksToBPFPId.isValid()){
        std::vector<art::Ptr<bpfit::BPFPId > > trackBPFPId = tracksToBPFPId.at(c);   
        // If the PID is lower, we can bail, nothing to update 
        if(trackBPFPId.size() == 0) continue;
        if(trackBPFPId.at(0)->Pdg() != pdg) continue;
        if(trackBPFPId.at(0)->Value() > highestPID){
          
          highestPID = trackBPFPId.at(0)->Value();
          // bestTrackLength = sliceTracks[c]->TotalLength();
          // bestStartZ = sliceTracks[c]->Start().Z(); 
          // bestADC    = sliceTracks[c]->TotalADC(); 
          bestTrack = c;
        }
        //else -- handle a tie between two BPFPId scores //FIXME

      
      }//End of valid BPFPId object
      
    }//End of for loop over tracks


    if (bestTrack < 0)
      return 999;

    return (unsigned int)bestTrack;
    
  }
///////////////////////////////////////////////////////////

}//namespace bpfit