/usr/include/root/TMVA/ModulekNN.h is in libroot-tmva-dev 5.34.30-0ubuntu8.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 | // @(#)root/tmva $Id$
// Author: Rustem Ospanov
/**********************************************************************************
* Project: TMVA - a Root-integrated toolkit for multivariate data analysis *
* Package: TMVA *
* Class : ModulekNN *
* Web : http://tmva.sourceforge.net *
* *
* Description: *
* Module for k-nearest neighbor algorithm *
* *
* Author: *
* Rustem Ospanov <rustem@fnal.gov> - U. of Texas at Austin, USA *
* *
* Copyright (c) 2007: *
* CERN, Switzerland *
* MPI-K Heidelberg, Germany *
* U. of Texas at Austin, USA *
* *
* Redistribution and use in source and binary forms, with or without *
* modification, are permitted according to the terms listed in LICENSE *
* (http://tmva.sourceforge.net/LICENSE) *
**********************************************************************************/
#ifndef ROOT_TMVA_ModulekNN
#define ROOT_TMVA_ModulekNN
//______________________________________________________________________
/*
kNN::Event describes point in input variable vector-space, with
additional functionality like distance between points
*/
//______________________________________________________________________
// C++
#include <cassert>
#include <iosfwd>
#include <map>
#include <string>
#include <vector>
// ROOT
#ifndef ROOT_Rtypes
#include "Rtypes.h"
#endif
#ifndef ROOT_TRandom
#include "TRandom3.h"
#endif
#ifndef ROOT_TMVA_NodekNN
#include "TMVA/NodekNN.h"
#endif
namespace TMVA {
class MsgLogger;
namespace kNN {
typedef Float_t VarType;
typedef std::vector<VarType> VarVec;
class Event {
public:
Event();
Event(const VarVec &vec, Double_t weight, Short_t type);
Event(const VarVec &vec, Double_t weight, Short_t type, const VarVec &tvec);
~Event();
Double_t GetWeight() const;
VarType GetVar(UInt_t i) const;
VarType GetTgt(UInt_t i) const;
UInt_t GetNVar() const;
UInt_t GetNTgt() const;
Short_t GetType() const;
// keep these two function separate
VarType GetDist(VarType var, UInt_t ivar) const;
VarType GetDist(const Event &other) const;
void SetTargets(const VarVec &tvec);
const VarVec& GetTargets() const;
const VarVec& GetVars() const;
void Print() const;
void Print(std::ostream& os) const;
private:
VarVec fVar; // coordinates (variables) for knn search
VarVec fTgt; // targets for regression analysis
Double_t fWeight; // event weight
Short_t fType; // event type ==0 or == 1, expand it to arbitrary class types?
};
typedef std::vector<TMVA::kNN::Event> EventVec;
typedef std::pair<const Node<Event> *, VarType> Elem;
typedef std::list<Elem> List;
std::ostream& operator<<(std::ostream& os, const Event& event);
class ModulekNN
{
public:
typedef std::map<int, std::vector<Double_t> > VarMap;
public:
ModulekNN();
~ModulekNN();
void Clear();
void Add(const Event &event);
Bool_t Fill(const UShort_t odepth, UInt_t ifrac, const std::string &option = "");
Bool_t Find(Event event, UInt_t nfind = 100, const std::string &option = "count") const;
Bool_t Find(UInt_t nfind, const std::string &option) const;
const EventVec& GetEventVec() const;
const List& GetkNNList() const;
const Event& GetkNNEvent() const;
const VarMap& GetVarMap() const;
const std::map<Int_t, Double_t>& GetMetric() const;
void Print() const;
void Print(std::ostream &os) const;
private:
Node<Event>* Optimize(UInt_t optimize_depth);
void ComputeMetric(UInt_t ifrac);
const Event Scale(const Event &event) const;
private:
static TRandom3 fgRndm;
UInt_t fDimn;
Node<Event> *fTree;
std::map<Int_t, Double_t> fVarScale;
mutable List fkNNList; // latest result from kNN search
mutable Event fkNNEvent; // latest event used for kNN search
std::map<Short_t, UInt_t> fCount; // count number of events of each type
EventVec fEvent; // vector of all events used to build tree and analysis
VarMap fVar; // sorted map of variables in each dimension for all event types
mutable MsgLogger* fLogger; // message logger
MsgLogger& Log() const { return *fLogger; }
};
//
// inlined functions for Event class
//
inline VarType Event::GetDist(const VarType var1, const UInt_t ivar) const
{
const VarType var2 = GetVar(ivar);
return (var1 - var2) * (var1 - var2);
}
inline Double_t Event::GetWeight() const
{
return fWeight;
}
inline VarType Event::GetVar(const UInt_t i) const
{
return fVar[i];
}
inline VarType Event::GetTgt(const UInt_t i) const
{
return fTgt[i];
}
inline UInt_t Event::GetNVar() const
{
return fVar.size();
}
inline UInt_t Event::GetNTgt() const
{
return fTgt.size();
}
inline Short_t Event::GetType() const
{
return fType;
}
//
// inline functions for ModulekNN class
//
inline const List& ModulekNN::GetkNNList() const
{
return fkNNList;
}
inline const Event& ModulekNN::GetkNNEvent() const
{
return fkNNEvent;
}
inline const EventVec& ModulekNN::GetEventVec() const
{
return fEvent;
}
inline const ModulekNN::VarMap& ModulekNN::GetVarMap() const
{
return fVar;
}
inline const std::map<Int_t, Double_t>& ModulekNN::GetMetric() const
{
return fVarScale;
}
} // end of kNN namespace
} // end of TMVA namespace
#endif
|