106 static const
Float_t kHigh= FLT_MAX;
137 fFillFoamWithOrigWeights(
kFALSE),
138 fDTSeparation(
kFoam),
141 fTimer(new
Timer(0, "PDEFoam", kTRUE)),
147 fVariableNames->SetOwner(kTRUE);
176 fFillFoamWithOrigWeights(
kFALSE),
177 fDTSeparation(
kFoam),
184 if(strlen(name) > 128)
197 delete fVariableNames;
199 if (fDistr)
delete fDistr;
200 if (fPseRan)
delete fPseRan;
201 if (fXmin)
delete [] fXmin; fXmin=0;
202 if (fXmax)
delete [] fXmax; fXmax=0;
206 for(
Int_t i=0; i<fNCells; i++)
delete fCells[i];
243 , fFillFoamWithOrigWeights(
kFALSE)
244 , fDTSeparation(
kFoam)
251 Log() <<
kFATAL <<
"COPY CONSTRUCTOR NOT IMPLEMENTED" <<
Endl;
264 Log() <<
kFATAL <<
"<SetDim>: Dimension is zero or negative!" <<
Endl;
267 if (fXmin)
delete [] fXmin;
268 if (fXmax)
delete [] fXmax;
278 if (idim<0 || idim>=GetTotDim())
279 Log() <<
kFATAL <<
"<SetXmin>: Dimension out of bounds!" <<
Endl;
289 if (idim<0 || idim>=GetTotDim())
290 Log() <<
kFATAL <<
"<SetXmax>: Dimension out of bounds!" <<
Endl;
308 if(fPseRan==0)
Log() <<
kFATAL <<
"Random number generator not set" <<
Endl;
309 if(fDistr==0)
Log() <<
kFATAL <<
"Distribution function not set" <<
Endl;
310 if(fDim==0)
Log() <<
kFATAL <<
"Zero dimension not allowed" <<
Endl;
317 if(fRvec==0)
Log() <<
kFATAL <<
"Cannot initialize buffer fRvec" <<
Endl;
321 if(fAlpha==0)
Log() <<
kFATAL <<
"Cannot initialize buffer fAlpha" <<
Endl;
326 fInhiDiv =
new Int_t[fDim];
327 for(
Int_t i=0; i<fDim; i++) fInhiDiv[i]=0;
331 fMaskDiv =
new Int_t[fDim];
332 for(
Int_t i=0; i<fDim; i++) fMaskDiv[i]=1;
336 for(
Int_t i=0; i<fDim; i++){
338 hname = fName+
TString(
"_HistEdge_");
340 htitle = TString(
"Edge Histogram No. ");
342 (*fHistEdg)[i] =
new TH1D(hname.
Data(),htitle.
Data(),fNBin,0.0, 1.0);
343 ((
TH1D*)(*fHistEdg)[i])->Sumw2();
371 for(
Int_t i=0; i<fNCells; i++)
delete fCells[i];
377 Log() <<
kFATAL <<
"not enough memory to create " << fNCells
380 for(
Int_t i=0; i<fNCells; i++){
382 fCells[i]->SetSerial(i);
391 for(
Long_t iCell=0; iCell<=fLastCe; iCell++){
392 Explore( fCells[iCell] );
403 if (fLastCe==fNCells){
408 cell = fCells[fLastCe];
410 cell->
Fill(status, parent, 0, 0);
463 cell->
GetHcub(cellPosi,cellSize);
473 for (
Int_t idim = 0; idim < fDim; ++idim)
474 vol_scale *= fXmax[idim] - fXmin[idim];
480 toteventsOld = GetCellElement(cell, 0);
491 for (i=0;i<fDim;i++) ((
TH1D *)(*fHistEdg)[i])->Reset();
495 for (iev=0;iev<fNSampl;iev++){
498 if (fDim>0)
for (j=0; j<fDim; j++) xRand[j]= cellPosi[j] +fAlpha[j]*(cellSize[j]);
500 wt = dx*
Eval(xRand, event_density);
501 totevents += event_density;
505 for (k=0; k<fDim; k++) {
507 ((
TH1D *)(*fHistEdg)[nProj])->Fill(xproj,wt);
515 if (ceSum[3]>wt) ceSum[3]=wt;
516 if (ceSum[4]<wt) ceSum[4]=wt;
518 if (ceSum[1]>0) nevEff = ceSum[0]*ceSum[0]/ceSum[1];
520 if ( nevEff >= fNBin*fEvPerBin)
break;
524 if (fNSampl>0) totevents /= fNSampl;
528 if (cell==fCells[0] && ceSum[0]<=0.0){
530 Log() <<
kFATAL <<
"No events were found during exploration of "
531 <<
"root cell. Please check PDEFoam parameters nSampl "
532 <<
"and VolFrac." <<
Endl;
534 Log() <<
kWARNING <<
"Negative number of events found during "
535 <<
"exploration of root cell" <<
Endl;
540 for (k=0; k<fDim;k++){
542 if ( fInhiDiv[k]==1) fMaskDiv[k] =0;
549 Double_t intTrue = ceSum[0]/(nevMC+0.000001);
552 if (kBest == -1) Varedu(ceSum,kBest,xBest,yBest);
553 intDriv =
sqrt(ceSum[1]/nevMC) -intTrue;
560 SetCellElement(cell, 0, totevents);
564 for (parent = cell->
GetPare(); parent!=0; parent = parent->
GetPare()){
567 parent->
SetIntg( parIntg +intTrue -intOld );
568 parent->
SetDriv( parDriv +intDriv -driOld );
569 SetCellElement( parent, 0, GetCellElement(parent, 0) + totevents - toteventsOld);
594 for(
Int_t kProj=0; kProj<fDim; kProj++) {
595 if( fMaskDiv[kProj]) {
602 for(
Int_t jLo=1; jLo<=fNBin; jLo++) {
604 for(
Int_t jUp=jLo; jUp<=fNBin;jUp++) {
605 aswIn += ((
TH1D *)(*fHistEdg)[kProj])->GetBinContent(jUp);
606 asswIn += Sqr(((
TH1D *)(*fHistEdg)[kProj])->GetBinError( jUp));
612 else sswIn =
sqrt(asswIn) /
sqrt(nent*(xUp-xLo)) *(xUp-xLo);
615 else sswOut=
sqrt(sswAll-asswIn)/
sqrt(nent*(1.0-xUp+xLo)) *(1.0-xUp+xLo);
616 if( (sswIn+sswOut) < sswtBest) {
617 sswtBest = sswIn+sswOut;
634 if(iLo == 0) xBest=yBest;
635 if(iUp == fNBin) yBest=xBest;
640 if( (kBest >= fDim) || (kBest<0) )
651 fPseRan->RndmArray(fDim,fRvec);
652 for(
Int_t k=0; k<fDim; k++) fAlpha[k] = fRvec[k];
672 for(i=0; i<=fLastCe; i++) {
673 if( fCells[i]->GetStat() == 1 ) {
675 driv = fCells[i]->GetDriv();
686 if (GetMaxDepth() > 0)
687 bCutMaxDepth = fCells[i]->GetDepth() < GetMaxDepth();
691 bCutNmin = GetCellElement(fCells[i], 0) > GetNmin();
694 if(driv > drivMax && bCutNmin && bCutMaxDepth) {
703 Log() <<
kVERBOSE <<
"Warning: No cell with more than "
704 << GetNmin() <<
" events found!" <<
Endl;
705 else if (!bCutMaxDepth)
706 Log() <<
kVERBOSE <<
"Warning: Maximum depth reached: "
707 << GetMaxDepth() <<
Endl;
709 Log() <<
kWARNING <<
"<PDEFoam::PeekMax>: no more candidate cells (drivMax>0) found for further splitting." <<
Endl;
730 if(fLastCe+1 >= fNCells)
Log() <<
kFATAL <<
"Buffer limit is reached, fLastCe=fnBuf" <<
Endl;
737 if( kBest<0 || kBest>=fDim )
Log() <<
kFATAL <<
"Wrong kBest" <<
Endl;
743 Int_t d1 = CellFill(1, cell);
744 Int_t d2 = CellFill(1, cell);
748 Explore( (fCells[d1]) );
749 Explore( (fCells[d2]) );
764 std::vector<Double_t> xvec;
765 xvec.reserve(GetTotDim());
766 for (
Int_t idim = 0; idim < GetTotDim(); ++idim)
767 xvec.push_back( VarTransformInvers(idim, xRand[idim]) );
769 return GetDistr()->Density(xvec, event_density);
780 fTimer->Init(fNCells);
785 while ( (fLastCe+2) < fNCells ) {
788 if ( (iCell<0) || (iCell>fLastCe) ) {
789 Log() <<
kVERBOSE <<
"Break: "<< fLastCe+1 <<
" cells created" <<
Endl;
791 for (
Long_t jCell=fLastCe+1; jCell<fNCells; jCell++)
792 delete fCells[jCell];
796 newCell = fCells[iCell];
800 if ( Divide( newCell )==0)
break;
805 Log() <<
kVERBOSE << GetNActiveCells() <<
" active cells created" <<
Endl;
817 fInhiDiv =
new Int_t[ fDim ];
818 for(
Int_t i=0; i<fDim; i++) fInhiDiv[i]=0;
821 if( ( 0<=iDim) && (iDim<fDim)) {
822 fInhiDiv[iDim] = inhiDiv;
839 errors = 0; warnings = 0;
840 if (level==1)
Log() <<
kVERBOSE <<
"Performing consistency checks for created foam" <<
Endl;
841 for(iCell=1; iCell<=fLastCe; iCell++) {
842 cell = fCells[iCell];
847 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d has only one daughter " << iCell <<
Endl;
851 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d has no daughter and is inactive " << iCell <<
Endl;
855 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d has two daughters and is active " << iCell <<
Endl;
859 if( (cell->
GetPare())!=fCells[0] ) {
862 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d parent not pointing to this cell " << iCell <<
Endl;
868 if(cell != (cell->
GetDau0())->GetPare()) {
870 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d daughter 0 not pointing to this cell " << iCell <<
Endl;
874 if(cell != (cell->
GetDau1())->GetPare()) {
876 if (level==1)
Log() <<
kFATAL <<
"ERROR: Cell's no %d daughter 1 not pointing to this cell " << iCell <<
Endl;
881 if(level==1)
Log() <<
kFATAL <<
"ERROR: Cell no. " << iCell <<
" has Volume of <1E-50" <<
Endl;
886 for(iCell=0; iCell<=fLastCe; iCell++) {
887 cell = fCells[iCell];
890 if(level==1)
Log() <<
kFATAL <<
"ERROR: Cell no. " << iCell <<
" is active but Volume is 0 " <<
Endl;
895 Log() <<
kVERBOSE <<
"Check has found " << errors <<
" errors and " << warnings <<
" warnings." <<
Endl;
898 Info(
"CheckAll",
"Check - found total %d errors \n",errors);
908 if (iCell < 0 || iCell > fLastCe) {
910 <<
")>: cell number " << iCell <<
" out of bounds!"
916 fCells[iCell]->GetHcub(cellPosi,cellSize);
917 Int_t kBest = fCells[iCell]->GetBest();
918 Double_t xBest = fCells[iCell]->GetXdiv();
920 Log() <<
"Cell[" << iCell <<
"]={ ";
921 Log() <<
" " << fCells[iCell] <<
" " <<
Endl;
922 Log() <<
" Xdiv[abs. coord.]="
923 << VarTransformInvers(kBest,cellPosi[kBest] + xBest*cellSize[kBest])
925 Log() <<
" Abs. coord. = (";
926 for (
Int_t idim=0; idim<fDim; idim++) {
927 Log() <<
"dim[" << idim <<
"]={"
928 << VarTransformInvers(idim,cellPosi[idim]) <<
","
929 << VarTransformInvers(idim,cellPosi[idim] + cellSize[idim])
935 fCells[iCell]->
Print(
"1");
937 Log() <<
"Elements: [";
941 if (i>0)
Log() <<
", ";
942 Log() << GetCellElement(fCells[iCell], i);
955 for(
Long_t iCell=0; iCell<=fLastCe; iCell++)
970 std::vector<Float_t> tvalues = VarTransform(values);
975 SetCellElement(cell, 0, GetCellElement(cell, 0) + wt);
976 SetCellElement(cell, 1, GetCellElement(cell, 1) + wt*wt);
987 for (
Long_t iCell = 0; iCell < fNCells; ++iCell) {
988 TObject* elements = fCells[iCell]->GetElement();
991 fCells[iCell]->SetElement(
NULL);
1029 std::vector<Float_t> txvec(VarTransform(xvec));
1031 return GetCellValue(FindCell(txvec), cv);
1033 return kernel->
Estimate(
this, txvec, cv);
1057 std::map<Int_t,Float_t> txvec;
1058 for (std::map<Int_t,Float_t>::const_iterator it=xvec.begin(); it!=xvec.end(); ++it)
1059 txvec.insert(std::pair<Int_t, Float_t>(it->first, VarTransform(it->first, it->second)));
1062 std::vector<PDEFoamCell*> cells = FindCells(txvec);
1065 std::vector<Float_t> cell_values;
1066 cell_values.reserve(cells.size());
1067 for (std::vector<PDEFoamCell*>::const_iterator cell_it=cells.begin();
1068 cell_it != cells.end(); ++cell_it)
1069 cell_values.push_back(GetCellValue(*cell_it, cv));
1093 PDEFoamVect cellPosi0(GetTotDim()), cellSize0(GetTotDim());
1101 cell0->
GetHcub(cellPosi0,cellSize0);
1103 if (xvec.at(idim)<=cellPosi0[idim]+cellSize0[idim])
1131 PDEFoamVect cellPosi0(GetTotDim()), cellSize0(GetTotDim());
1139 map<Int_t, Float_t>::const_iterator it = txvec.find(idim);
1141 if (it != txvec.end()){
1145 cell0->
GetHcub(cellPosi0,cellSize0);
1147 if (it->second <= cellPosi0[idim] + cellSize0[idim])
1153 FindCells(txvec, cell->
GetDau0(), cells);
1154 FindCells(txvec, cell->
GetDau1(), cells);
1158 cells.push_back(cell);
1179 std::map<Int_t, Float_t> txvec_map;
1180 for (
UInt_t i=0; i<txvec.size(); ++i)
1181 txvec_map.insert(std::pair<Int_t, Float_t>(i, txvec.at(i)));
1184 std::vector<PDEFoamCell*> cells(0);
1187 FindCells(txvec_map, fCells[0], cells);
1210 std::vector<PDEFoamCell*> cells(0);
1213 FindCells(txvec, fCells[0], cells);
1232 if ( GetTotDim()!=1 )
1233 Log() <<
kFATAL <<
"<Draw1Dim>: function can only be used for 1-dimensional foams!"
1239 h1=
new TH1D(hname,
"1-dimensional Foam", nbin, fXmin[0], fXmax[0]);
1241 if (!h1)
Log() <<
kFATAL <<
"ERROR: Can not create histo" << hname <<
Endl;
1246 std::vector<Float_t> txvec;
1247 txvec.push_back( VarTransform(0, h1->
GetBinCenter(ibinx)) );
1249 if (kernel !=
NULL) {
1251 val = kernel->
Estimate(
this, txvec, cell_value);
1253 val = GetCellValue(FindCell(txvec), cell_value);
1284 if ((idim1>=GetTotDim()) || (idim1<0) ||
1285 (idim2>=GetTotDim()) || (idim2<0) ||
1287 Log() <<
kFATAL <<
"<Project2>: wrong dimensions given: "
1288 << idim1 <<
", " << idim2 <<
Endl;
1294 Log() <<
kWARNING <<
"Warning: number of bins too big: " << nbin
1295 <<
" Using 1000 bins for each dimension instead." <<
Endl;
1297 }
else if (nbin<1) {
1299 <<
"; set nbin=50" <<
Endl;
1309 h1=
new TH2D(hname.
Data(),
Form(
"var%d vs var%d",idim1,idim2), nbin, fXmin[idim1], fXmax[idim1], nbin, fXmin[idim2], fXmax[idim2]);
1311 if (!h1)
Log() <<
kFATAL <<
"ERROR: Can not create histo" << hname <<
Endl;
1319 std::map<Int_t, Float_t> txvec;
1325 std::vector<TMVA::PDEFoamCell*> cells = FindCells(txvec);
1330 for (std::vector<TMVA::PDEFoamCell*>::const_iterator it = cells.begin();
1331 it != cells.end(); ++it) {
1333 PDEFoamVect cellPosi(GetTotDim()), cellSize(GetTotDim());
1334 (*it)->GetHcub(cellPosi,cellSize);
1337 std::vector<Float_t> tvec;
1338 for (
Int_t i=0; i<GetTotDim(); ++i) {
1339 if ( i != idim1 && i != idim2 )
1340 tvec.push_back(cellPosi[i] + 0.5*cellSize[i]);
1342 tvec.push_back(txvec[i]);
1344 if (kernel !=
NULL) {
1346 sum_cv += kernel->
Estimate(
this, tvec, cell_value);
1348 sum_cv += GetCellValue(FindCell(tvec), cell_value);
1373 return GetCellElement(cell, 0);
1376 return GetCellElement(cell, 1);
1382 return GetCellValue(cell,
kValue)/volume;
1386 Log() <<
kWARNING <<
"<GetCellDensity(cell)>: ERROR: cell volume"
1387 <<
" negative or zero!"
1388 <<
" ==> return cell density 0!"
1389 <<
" cell volume=" << volume
1390 <<
" cell entries=" << GetCellValue(cell,
kValue) <<
Endl;
1392 Log() <<
kWARNING <<
"<GetCellDensity(cell)>: WARNING: cell volume"
1393 <<
" close to zero!"
1394 <<
" cell volume: " << volume <<
Endl;
1416 Log() <<
kFATAL <<
"<GetCellValue>: unknown cell value" <<
Endl;
1458 Log() <<
kFATAL <<
"<SetCellElement> ERROR: cell element is not a TVectorD*" <<
Endl;
1474 Log() <<
kINFO <<
"Elapsed time: " + fTimer->GetElapsedTime()
1481 if (fNCells >= 100) modulo =
Int_t(fNCells/100);
1482 if (fLastCe%modulo == 0) fTimer->DrawProgressBar( fLastCe );
1525 if (GetTotDim() != 2)
1526 Log() <<
kFATAL <<
"RootPlot2dim() can only be used with "
1527 <<
"two-dimensional foams!" <<
Endl;
1535 }
else if (opt.
Contains(
"rms_ov_mean")){
1543 plotcellnumber =
kTRUE;
1548 outfile<<
"{" << std::endl;
1553 outfile <<
"TColor *graycolors[100];" << std::endl;
1554 outfile <<
"for (Int_t i=0.; i<100; i++)" << std::endl;
1555 outfile <<
" graycolors[i]=new TColor(1000+i, 1-(Float_t)i/100.,1-(Float_t)i/100.,1-(Float_t)i/100.);"<< std::endl;
1558 outfile <<
"cMap = new TCanvas(\"" << fName <<
"\",\"Cell Map for "
1559 << fName <<
"\",600,600);" << std::endl;
1561 outfile<<
"TBox*a=new TBox();"<<std::endl;
1562 outfile<<
"a->SetFillStyle(0);"<<std::endl;
1563 outfile<<
"a->SetLineWidth(4);"<<std::endl;
1564 outfile<<
"TBox *b1=new TBox();"<<std::endl;
1565 outfile<<
"TText*t=new TText();"<<std::endl;
1567 outfile << (colors ?
"gStyle->SetPalette(1, 0);" :
"gStyle->SetPalette(0);")
1569 outfile <<
"b1->SetFillStyle(1001);"<<std::endl;
1570 outfile<<
"TBox *b2=new TBox();"<<std::endl;
1571 outfile <<
"b2->SetFillStyle(0);"<<std::endl;
1574 outfile <<
"b1->SetFillStyle(0);"<<std::endl;
1586 for (
Long_t iCell=1; iCell<=fLastCe; iCell++) {
1587 if ( fCells[iCell]->GetStat() == 1) {
1588 Float_t value = GetCellValue(fCells[iCell], cell_value);
1595 outfile <<
"// observed minimum and maximum of distribution: " << std::endl;
1596 outfile <<
"// Float_t zmin = "<< zmin <<
";" << std::endl;
1597 outfile <<
"// Float_t zmax = "<< zmax <<
";" << std::endl;
1600 outfile <<
"// used minimum and maximum of distribution (taking into account log scale if applicable): " << std::endl;
1601 outfile <<
"Float_t zmin = "<< zmin <<
";" << std::endl;
1602 outfile <<
"Float_t zmax = "<< zmax <<
";" << std::endl;
1608 Float_t scale = (ncolors-1)/(zmax - zmin);
1609 PDEFoamVect cellPosi(GetTotDim()), cellSize(GetTotDim());
1613 outfile <<
"// =========== Rectangular cells ==========="<< std::endl;
1614 for (
Long_t iCell=1; iCell<=fLastCe; iCell++) {
1615 if ( fCells[iCell]->GetStat() == 1) {
1616 fCells[iCell]->GetHcub(cellPosi,cellSize);
1617 x1 = offs+lpag*(cellPosi[0]);
1618 y1 = offs+lpag*(cellPosi[1]);
1619 x2 = offs+lpag*(cellPosi[0]+cellSize[0]);
1620 y2 = offs+lpag*(cellPosi[1]+cellSize[1]);
1624 Float_t value = GetCellValue(fCells[iCell], cell_value);
1631 color = 1000+(
Int_t((value-zmin)*scale));
1634 outfile <<
"b1->SetFillColor(" << color <<
");" << std::endl;
1638 outfile<<
"b1->DrawBox("<<x1<<
","<<y1<<
","<<x2<<
","<<y2<<
");"<<std::endl;
1640 outfile<<
"b2->DrawBox("<<x1<<
","<<y1<<
","<<x2<<
","<<y2<<
");"<<std::endl;
1643 if (plotcellnumber) {
1644 outfile<<
"t->SetTextColor(4);"<<std::endl;
1646 outfile<<
"t->SetTextSize(0.025);"<<std::endl;
1647 else if(fLastCe<251)
1648 outfile<<
"t->SetTextSize(0.015);"<<std::endl;
1650 outfile<<
"t->SetTextSize(0.008);"<<std::endl;
1651 x = offs+lpag*(cellPosi[0]+0.5*cellSize[0]);
1652 y = offs+lpag*(cellPosi[1]+0.5*cellSize[1]);
1653 outfile<<
"t->DrawText("<<x<<
","<<y<<
","<<
"\""<<iCell<<
"\""<<
");"<<std::endl;
1657 outfile<<
"// ============== End Rectangles ==========="<< std::endl;
1659 outfile <<
"}" << std::endl;
1670 GetDistr()->FillBinarySearchTree(ev);
1679 if(fDistr)
delete fDistr;
virtual void FillFoamCells(const Event *ev, Float_t wt)
This function fills a weight 'wt' into the PDEFoam cell, which corresponds to the given event 'ev'...
TVectorT< Element > & ResizeTo(Int_t lwb, Int_t upb)
Resize the vector to [lwb:upb] .
Random number generator class based on M.
MsgLogger & Endl(MsgLogger &ml)
void SetElement(TObject *fobj)
Long_t PeekMax()
Internal subprogram used by Create.
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
TObjArray * fVariableNames
timer for graphical output
void OutputGrow(Bool_t finished=false)
message logger
R__EXTERN TStyle * gStyle
virtual void SetOwner(Bool_t enable=kTRUE)
Set whether this collection is the owner (enable==true) of its content.
void Fill(Int_t, PDEFoamCell *, PDEFoamCell *, PDEFoamCell *)
Fills in certain data into newly allocated cell.
std::vector< double > values
void Grow()
Internal subrogram used by Create.
static const char * filename()
static Bool_t AddDirectoryStatus()
static function: cannot be inlined on Windows/NT
void SetXmin(Int_t idim, Double_t wmin)
set lower foam bound in dimension idim
TObject * GetElement() const
virtual Int_t GetNbinsX() const
void SetCellElement(PDEFoamCell *cell, UInt_t i, Double_t value)
Set cell element i of cell to value.
Int_t Divide(PDEFoamCell *)
Internal subrogram used by Create.
TH1D * Draw1Dim(ECellValue cell_value, Int_t nbin, PDEFoamKernelBase *kernel=NULL)
Draws 1-dimensional foam (= histogram)
void CheckAll(Int_t)
User utility, miscellaneous and debug.
Double_t GetVolume() const
std::vector< TMVA::PDEFoamCell * > FindCells(const std::vector< Float_t > &) const
Find all cells, that contain txvec.
static void AddDirectory(Bool_t add=kTRUE)
Sets the flag controlling the automatic add of histograms in memory.
Int_t CellFill(Int_t, PDEFoamCell *)
Internal subprogram used by Create.
const char * Data() const
void MakeAlpha()
Internal subrogram used by Create.
virtual Bool_t CellValueIsUndefined(PDEFoamCell *)
Returns true, if the value of the given cell is undefined.
static const double x2[5]
virtual Float_t Estimate(PDEFoam *, std::vector< Float_t > &, ECellValue)=0
void RootPlot2dim(const TString &filename, TString opt, Bool_t createCanvas=kTRUE, Bool_t colors=kTRUE)
Debugging tool which plots the cells of a 2-dimensional PDEFoam as rectangles in C++ format readable ...
void SetDau1(PDEFoamCell *Daug)
void Info(const char *location, const char *msgfmt,...)
Int_t GetNumberOfColors() const
Return number of colors in the color palette.
void ResetCellElements()
Remove the cell elements from all cells.
virtual void Explore(PDEFoamCell *Cell)
Internal subprogram used by Create.
TVectorT< Double_t > TVectorD
void Print(Option_t *option) const
Printout of the cell geometry parameters for the debug purpose.
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
virtual void Print(Option_t *option="") const
This method must be overridden when a class wants to print itself.
void InitCells()
Internal subprogram used by Create.
void PrintCells()
Prints geometry of ALL cells of the FOAM.
virtual Double_t GetBinCenter(Int_t bin) const
return bin center for 1D historam Better to use h1.GetXaxis().GetBinCenter(bin)
void Varedu(Double_t[], Int_t &, Double_t &, Double_t &)
Internal subrogram used by Create.
void SetXdiv(Double_t Xdiv)
void SetXmax(Int_t idim, Double_t wmax)
set upper foam bound in dimension idim
Double_t GetCellElement(const PDEFoamCell *cell, UInt_t i) const
Returns cell element i of cell 'cell'.
PDEFoamCell * GetPare() const
void CalcVolume()
Calculates volume of the cell using size params which are calculated.
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content see convention for numbering bins in TH1::GetBin In case the bin number is greater th...
TVectorT< Element > & Zero()
Set vector elements to zero.
char * Form(const char *fmt,...)
void SetDriv(Double_t Driv)
void SetDim(Int_t kDim)
Sets dimension of cubical space.
1-D histogram with a double per channel (see TH1 documentation)}
void PrintCell(Long_t iCell=0)
Prints geometry of and elements of 'iCell', as well as relations to parent and daughter cells...
Int_t GetColorPalette(Int_t i) const
Return color number i in current palette.
Double_t Eval(Double_t *xRand, Double_t &event_density)
Internal subprogram.
static const double x1[5]
void Create()
Basic initialization of FOAM invoked by the user.
virtual TH2D * Project2(Int_t idim1, Int_t idim2, ECellValue cell_value=kValue, PDEFoamKernelBase *kernel=NULL, UInt_t nbin=50)
Project foam variable idim1 and variable idim2 to histogram.
ClassImp(TMVA::PDEFoam) static const Float_t kHigh
virtual Float_t GetCellValue(const std::vector< Float_t > &xvec, ECellValue cv, PDEFoamKernelBase *)
This function finds the cell, which corresponds to the given untransformed event vector 'xvec' and re...
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
Mother of all ROOT objects.
virtual Int_t GetNbinsY() const
void SetIntg(Double_t Intg)
void DeleteBinarySearchTree()
Delete the foam's density estimator, which contains the binary search tree.
std::vector< Float_t > & GetValues()
void FillBinarySearchTree(const Event *ev)
Insert event to internal foam's density estimator PDEFoamDensityBase.
PDEFoamCell * GetDau1() const
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
virtual ~PDEFoam()
Default destructor.
PDEFoamCell * FindCell(const std::vector< Float_t > &) const
Find cell that contains 'xvec' (in foam coordinates [0,1]).
std::vector< double > errors
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content.
void SetDau0(PDEFoamCell *Daug)
static const Float_t kVlow
void GetHcub(PDEFoamVect &, PDEFoamVect &) const
Provides size and position of the cell These parameter are calculated by analyzing information in all...
void SetPalette(Int_t ncolors=kBird, Int_t *colors=0, Float_t alpha=1.)
See TColor::SetPalette.
void SetInhiDiv(Int_t, Int_t)
This can be called before Create, after setting kDim It defines which variables are excluded in the p...
PDEFoam()
Default constructor for streamer, user should not use it.
PDEFoamCell * GetDau0() const
2-D histogram with a double per channel (see TH1 documentation)}