Example10.cxx

Go to the documentation of this file.

// $Id$
//-----------------------------------------------------------------------
//       The GSI Online Offline Object Oriented (Go4) Project
//         Experiment Data Processing at EE department, GSI
//-----------------------------------------------------------------------
// Copyright (C) 2000- GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
//                     Planckstr. 1, 64291 Darmstadt, Germany
// Contact:            http://go4.gsi.de
//-----------------------------------------------------------------------
// This software can be used under the license agreements as stated
// in Go4License.txt file which is part of the distribution.
//-----------------------------------------------------------------------

 
 
#ifndef __CINT__
 
#include "TH1.h"
#include "TFile.h"
#include "TCanvas.h"
#include "TApplication.h"
 
#include "TGo4FitMinuit.h"
#include "TGo4Fitter.h"
#include "TGo4FitDataHistogram.h"
#include "TGo4FitModelPolynom.h"
#include "TGo4FitModelGauss1.h"
#include "TGo4FitModelFromData.h"
 
void Example10();
 
int main(int argc, char **argv)
{
   TApplication theApp("Application", nullptr, nullptr);
 
   Example10();
 
   theApp.Run();
 
   return 0;
}
int main(int argc, char **argv) {…}
 
#endif
 
// routine to read histogram from examples file
TH1D* GetHistogram(const char *HistogramName)
{
   TFile *f = TFile::Open("histograms.root");
   if (!f) return nullptr;
   TH1D *histo = nullptr;
   f->GetObject(HistogramName, histo);
   if (histo) histo->SetDirectory(nullptr);
   delete f;
   return histo;
}
TH1D* GetHistogram(const char *HistogramName) {…}
 
void Example10()
{
// create fitter and select function to fit
   TGo4Fitter fitter("Fitter", TGo4Fitter::ff_ML_Poisson, kTRUE);
 
   TH1D* histo1 = GetHistogram("hDeg120_P_c");            // get histogram from file
   TH1D* histo2 = GetHistogram("hDeg120_CND");            // get histogram from file
 
// add histogram to fitter, which should be fitted
   fitter.AddH1("data2", histo2, kTRUE, 2200., 2900.);
 
// create polynom of first order
   fitter.AddPolynomX("data2", "Pol", 1);
 
// add first histogram as component of second histogram
   fitter.AddModel("data2", new TGo4FitModelFromData("Model1",histo1,kTRUE) );
 
// create two gaussians
   fitter.AddGauss1("data2", "Gauss3", 2597., 14.);
   fitter.AddGauss1("data2", "Gauss4", 2717., 14.);
 
// perform all actions
   fitter.DoActions();
 
   // call fitter Draw method with options #data1,Gauss1,Gauss2
   fitter.Draw("#data2,Model1,Gauss3,Gauss4");
   fitter.Print("Pars");
}
void Example10() {…}