Index: trunk/MagicSoft/Mars/mhist/MHRanForest.cc =================================================================== --- trunk/MagicSoft/Mars/mhist/MHRanForest.cc (revision 2060) +++ (revision ) @@ -1,186 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Hengstebeck 3/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// MHRanForest -// -// This histogram shows the evolution of the standard deviation -// of est. hadronness as the number of trees increases. -// It helps you to find out how many trees are really needed in g/h-sep. -// -//////////////////////////////////////////////////////////////////////////// -#include "MHRanForest.h" - -#include -#include -#include -#include -#include - -#include "MParList.h" -#include "MBinning.h" -#include "MRanForest.h" - -#include "MLog.h" -#include "MLogManip.h" - -#include "MMcEvt.hxx" - -ClassImp(MHRanForest); - -// -------------------------------------------------------------------------- -// -// Setup histograms, nbins is the number of bins used for the evaluation. -// The default is 100 bins. -// -MHRanForest::MHRanForest(Int_t nbins, const char *name, const char *title) -{ - // - // set the name and title of this object - // - fName = name ? name : "MHRanForest"; - fTitle = title ? title : "Histogram showing Standard deviation of estimated hadronness"; - - fGraphSigma = new TGraph; - fGraphSigma->SetTitle("Evolution of Standard deviation of estimated hadronness in tree combination"); - fGraphSigma->SetMaximum(1); - fGraphSigma->SetMarkerStyle(kFullDotSmall); -} - -// -------------------------------------------------------------------------- -// -// Delete the histograms. -// -MHRanForest::~MHRanForest() -{ - delete fGraphSigma; -} - -// -------------------------------------------------------------------------- -// -// Setup Filling of the histograms. It needs: -// MMcEvt and MRanForest -// -Bool_t MHRanForest::SetupFill(const MParList *plist) -{ - fMcEvt = (MMcEvt*)plist->FindObject("MMcEvt"); - if (!fMcEvt) - { - *fLog << err << dbginf << "MMcEvt not found... aborting." << endl; - return kFALSE; - } - - fRanForest = (MRanForest*)plist->FindObject("MRanForest"); - if (!fRanForest) - { - *fLog << err << dbginf << "MRanForest not found... aborting." << endl; - return kFALSE; - } - - fSigma.Set(fRanForest->GetNumTrees()); - fNumEvent=0; - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// -Bool_t MHRanForest::Fill(const MParContainer *par, const Stat_t w) -{ - fNumEvent++; - Double_t hest=0; - Double_t htrue=fMcEvt->GetPartId()==kGAMMA ? 0. : 1.; - - Int_t ntrees=fRanForest->GetNumTrees(); - - for (Int_t i=0;iGetTreeHad(j); - - hest/=i+1; - fSigma[i]+=(htrue-hest)*(htrue-hest); - } - - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// Finalize the histogram: -// calculate standard deviation and set histogram max and min -// -Bool_t MHRanForest::Finalize() -{ - Int_t n = fSigma.GetSize(); - - fGraphSigma->Set(n); - - Stat_t max=0.; - Stat_t min=0.; - for (Int_t i=0; iSetPoint(i,ip,ig); - } - fGraphSigma->SetMaximum(1.05*max); - fGraphSigma->SetMinimum(0.95*min); - - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// Draw histogram. (For the Meaning see class description) -// -void MHRanForest::Draw(Option_t *) -{ - if (fGraphSigma->GetN()==0) - return; - - TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this); - pad->SetBorderMode(0); - - AppendPad(""); - - fGraphSigma->Draw("ALP"); - pad->Modified(); - pad->Update(); - - TH1 *h=fGraphSigma->GetHistogram(); - if (!h) - return; - - h->GetXaxis()->SetRangeUser(0, 1); - h->SetXTitle("No.of Trees"); - h->SetYTitle("\\sigma of est.hadronness"); - - pad->Modified(); - pad->Update(); -} Index: trunk/MagicSoft/Mars/mhist/MHRanForest.h =================================================================== --- trunk/MagicSoft/Mars/mhist/MHRanForest.h (revision 2060) +++ (revision ) @@ -1,43 +1,0 @@ -#ifndef MARS_MHRanForest -#define MARS_MHRanForest - -#ifndef MARS_MH -#include "MH.h" -#endif - -#ifndef ROOT_TArrayF -#include -#endif - -class TH1D; -class TGraph; -class MParList; -class MMcEvt; -class MRanForest; - -class MHRanForest : public MH -{ -private: - const MMcEvt *fMcEvt; //! - const MRanForest *fRanForest; //! - - TArrayF fSigma; //! - Int_t fNumEvent; //! - TGraph *fGraphSigma; //-> - -public: - MHRanForest(Int_t nbins=100, const char *name=NULL, const char *title=NULL); - ~MHRanForest(); - - TGraph *GetGraphSigma() const { return fGraphSigma; } - - Bool_t SetupFill(const MParList *plist); - Bool_t Fill(const MParContainer *par, const Stat_t w=1); - Bool_t Finalize(); - - void Draw(Option_t *opt=""); - - ClassDef(MHRanForest, 1) // Histogram showing variance of estimated Hadronness -}; - -#endif Index: trunk/MagicSoft/Mars/mhist/MHRanForestGini.cc =================================================================== --- trunk/MagicSoft/Mars/mhist/MHRanForestGini.cc (revision 2060) +++ (revision ) @@ -1,169 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Hengstebeck 3/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// MHRanForest -// -// This histogram shows a measure of variable importance (mean decrease in -// Gini-index) -// -//////////////////////////////////////////////////////////////////////////// -#include "MHRanForestGini.h" - -#include -#include -#include -#include -#include - -#include "MParList.h" -#include "MBinning.h" -#include "MRanTree.h" -#include "MRanForest.h" - -#include "MLog.h" -#include "MLogManip.h" - -ClassImp(MHRanForestGini); - -// -------------------------------------------------------------------------- -// -// Setup histograms, nbins is the number of bins used for the evaluation. -// The default is 100 bins. -// -MHRanForestGini::MHRanForestGini(Int_t nbins, const char *name, const char *title) -{ - // - // set the name and title of this object - // - fName = name ? name : "MHRanForestGini"; - fTitle = title ? title : "Measure of importance of Random Forest-input parameters"; - - fGraphGini = new TGraph; - fGraphGini->SetTitle("Importance of RF-input parameters measured by mean Gini decrease"); - fGraphGini->SetMaximum(1); - fGraphGini->SetMarkerStyle(kFullDotSmall); -} - -// -------------------------------------------------------------------------- -// -// Delete the histograms. -// -MHRanForestGini::~MHRanForestGini() -{ - delete fGraphGini; -} - -// -------------------------------------------------------------------------- -// -// Setup Filling of the histograms. It needs: -// MMcEvt and MRanForest -// -Bool_t MHRanForestGini::SetupFill(const MParList *plist) -{ - fRanForest = (MRanForest*)plist->FindObject("MRanForest"); - if (!fRanForest) - { - *fLog << err << dbginf << "MRanForest not found... aborting." << endl; - return kFALSE; - } - - fGini.Set(fRanForest->GetNumDim()); - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// Fill the RanForest from a MRanForest container into the corresponding -// histogram dependant on the particle id. -// -// -Bool_t MHRanForestGini::Fill(const MParContainer *par, const Stat_t w) -{ - for (Int_t i=0;iGetNumDim();i++) - fGini[i]+=fRanForest->GetCurTree()->GetGiniDec(i); - - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// -Bool_t MHRanForestGini::Finalize() -{ - Int_t n = fGini.GetSize(); - - fGraphGini->Set(n); - - Stat_t max=0.; - Stat_t min=0.; - for (Int_t i=0; iGetNumTrees(); - fGini[i]/=fRanForest->GetNumData(); - - Stat_t ip = i+1.; - Stat_t ig = fGini[i]; - - max=TMath::Max(max,ig); - min=TMath::Min(min,ig); - - fGraphGini->SetPoint(i,ip,ig); - } - fGraphGini->SetMaximum(1.05*max); - fGraphGini->SetMinimum(0.95*min); - - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// Draw histogram. (For the Meaning see class description) -// -void MHRanForestGini::Draw(Option_t *) -{ - if (fGraphGini->GetN()==0) - return; - - TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this); - pad->SetBorderMode(0); - - AppendPad(""); - - fGraphGini->Draw("ALP"); - pad->Modified(); - pad->Update(); - - TH1 *h = fGraphGini->GetHistogram(); - if (!h) - return; - - h->GetXaxis()->SetRangeUser(0, 1); - h->SetXTitle("No.of RF-input parameter"); - h->SetYTitle("Mean decrease in Gini-index [a.u.]"); - - pad->Modified(); - pad->Update(); -} Index: trunk/MagicSoft/Mars/mhist/MHRanForestGini.h =================================================================== --- trunk/MagicSoft/Mars/mhist/MHRanForestGini.h (revision 2060) +++ (revision ) @@ -1,41 +1,0 @@ -#ifndef MARS_MHRanForestGini -#define MARS_MHRanForestGini - -#ifndef MARS_MH -#include "MH.h" -#endif - -#ifndef ROOT_TArrayF -#include -#endif - -class TH1D; -class TGraph; -class MParList; -class MRanForest; -class MRanTree; - -class MHRanForestGini : public MH -{ -private: - MRanForest *fRanForest; //! - - TArrayF fGini; //! - TGraph *fGraphGini; //-> - -public: - MHRanForestGini(Int_t nbins=100, const char *name=NULL, const char *title=NULL); - ~MHRanForestGini(); - - TGraph *GetGraphGini() const { return fGraphGini; } - - Bool_t SetupFill(const MParList *plist); - Bool_t Fill(const MParContainer *par, const Stat_t w=1); - Bool_t Finalize(); - - void Draw(Option_t *opt=""); - - ClassDef(MHRanForestGini, 1) -}; - -#endif