Logo ROOT   6.12/07
Reference Guide
Background_incr.C
Go to the documentation of this file.
1 /// \file
2 /// \ingroup tutorial_spectrum
3 /// \notebook
4 /// Example to illustrate the background estimator (class TSpectrum).
5 ///
6 /// \macro_image
7 /// \macro_code
8 ///
9 /// \authors Miroslav Morhac, Olivier Couet
10 
11 void Background_incr() {
12  Int_t i;
13  const Int_t nbins = 1024;
14  Double_t xmin = 0;
15  Double_t xmax = nbins;
16  Double_t source[nbins];
17  gROOT->ForceStyle();
18 
19  TH1F *d = new TH1F("d","",nbins,xmin,xmax);
20 
21  TString dir = gROOT->GetTutorialDir();
22  TString file = dir+"/spectrum/TSpectrum.root";
23  TFile *f = new TFile(file.Data());
24  TH1F *back = (TH1F*) f->Get("back1");
25  back->GetXaxis()->SetRange(1,nbins);
26  back->SetTitle("Estimation of background with increasing window");
27  back->Draw("L");
28 
29  TSpectrum *s = new TSpectrum();
30 
31  for (i = 0; i < nbins; i++) source[i] = back->GetBinContent(i + 1);
32 
33  // Estimate the background
37 
38  // Draw the estimated background
39  for (i = 0; i < nbins; i++) d->SetBinContent(i + 1,source[i]);
40  d->SetLineColor(kRed);
41  d->Draw("SAME L");
42 }
float xmin
Definition: THbookFile.cxx:93
Definition: Rtypes.h:59
THist< 1, float, THistStatContent, THistStatUncertainty > TH1F
Definition: THist.hxx:285
A ROOT file is a suite of consecutive data records (TKey instances) with a well defined format...
Definition: TFile.h:46
virtual TObject * Get(const char *namecycle)
Return pointer to object identified by namecycle.
#define gROOT
Definition: TROOT.h:402
Basic string class.
Definition: TString.h:125
1-D histogram with a float per channel (see TH1 documentation)}
Definition: TH1.h:567
int Int_t
Definition: RtypesCore.h:41
virtual void Draw(Option_t *option="")
Default Draw method for all objects.
Definition: TObject.cxx:195
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition: TAttLine.h:40
virtual void Draw(Option_t *option="")
Draw this histogram with options.
Definition: TH1.cxx:2969
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...
Definition: TH1.cxx:8477
float xmax
Definition: THbookFile.cxx:93
const Bool_t kFALSE
Definition: RtypesCore.h:88
double Double_t
Definition: RtypesCore.h:55
virtual TH1 * Background(const TH1 *hist, Int_t niter=20, Option_t *option="")
One-dimensional background estimation function.
Definition: TSpectrum.cxx:152
static constexpr double s
Advanced Spectra Processing.
Definition: TSpectrum.h:18
Definition: file.py:1
const char * Data() const
Definition: TString.h:345