47 #ifndef ROOT_TMVA_PDEFoamDiscriminantDensity
55 TMVA::PDEFoamDiscriminantDensity::PDEFoamDiscriminantDensity()
56 : PDEFoamDensityBase()
105 Log() <<
kFATAL <<
"<PDEFoamDiscriminantDensity::Density()> Binary tree not set!" <<
Endl;
108 std::vector<Double_t> lb(GetBox().size());
109 std::vector<Double_t> ub(GetBox().size());
112 const Double_t probevolume_inv = 1.0 / GetBoxVolume();
115 for (
UInt_t idim = 0; idim < GetBox().size(); ++idim) {
116 lb[idim] = xev[idim] - GetBox().at(idim) / 2.0;
117 ub[idim] = xev[idim] + GetBox().at(idim) / 2.0;
121 std::vector<const TMVA::BinarySearchTreeNode*> nodes;
124 const Double_t sumOfWeights = fBst->SearchVolume(&volume, &nodes);
127 event_density = nodes.size() * probevolume_inv;
131 for (std::vector<const TMVA::BinarySearchTreeNode*>::const_iterator it = nodes.begin();
132 it != nodes.end(); ++it) {
133 if ((*it)->GetClass() ==
fClass)
134 n_sig += (*it)->GetWeight();
138 return (n_sig / (sumOfWeights + 0.1)) * probevolume_inv;
MsgLogger & Endl(MsgLogger &ml)
void box(Int_t pat, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
TClass * fClass
pointer to the foreign object
PDEFoamDiscriminantDensity()
Abstract ClassifierFactory template that handles arbitrary types.
virtual Double_t Density(std::vector< Double_t > &Xarg, Double_t &event_density)
This function is needed during the foam buildup.
ClassImp(TMVA::PDEFoamDiscriminantDensity) TMVA