ParticleID/RICHAnalysis/PidInfo_RICH_v1.cc

0001 #include "PidInfo_RICH_v1.h"
0002
0003 using namespace std;
0004
0005 PidInfo_RICH_v1::PidInfo_RICH_v1( const int id ):
0006   PidInfo( id ),
0007   _trackid( id ),
0008   _track_state( nullptr ),
0009   _pid_candidate_loglikelihood()
0010 {}
0011
0012 void
0013 PidInfo_RICH_v1::Reset()
0014 {
0015   _pid_candidate_loglikelihood.clear();
0016 }
0017
0018 float PidInfo_RICH_v1::get_likelihood(const PidInfo::PID_CANDIDATE candidate)
0019 {
0020   std::map<PidInfo::PID_CANDIDATE,float>::const_iterator citer = _pid_candidate_loglikelihood.find(candidate);
0021   if (citer == _pid_candidate_loglikelihood.end()) return NAN;
0022   return citer->second;
0023 }
0024
0025 /* //(for copy constructor:)
0026
0027 // copy over particle id probabilities
0028 std::vector<PID_CANDIDATE> pid_candidates;
0029 pid_candidates.push_back(SvtxTrack::ELECTRON);
0030 pid_candidates.push_back(SvtxTrack::CHARGEDPION);
0031 pid_candidates.push_back(SvtxTrack::CHARGEDKAON);
0032 pid_candidates.push_back(SvtxTrack::PROTON);
0033
0034 _pid_candidate_loglikelihood.clear();
0035
0036 for (unsigned int i=0; i<pid_candidates.size(); ++i) {
0037 if ( !isnan( track.get_pid_likelihood ( pid_candidates[i] ) ) )
0038 set_pid_likelihood( pid_candidates[i], track.get_pid_likelihood( pid_candidates[i] ) );
0039 }
0040 */