doc/master/Background__width2_8C.html

void Background_width2()

{

   Int_t i;

   const Int_t nbins = 4096;

   Double_t xmin = 0;

   Double_t xmax = 4096;

   Double_t source[nbins];

   gROOT->ForceStyle();


   TH1F *d1 = new TH1F("d1", "", nbins, xmin, xmax);

   TH1F *d2 = new TH1F("d2", "", nbins, xmin, xmax);

   TH1F *d3 = new TH1F("d3", "", nbins, xmin, xmax);

   TH1F *d4 = new TH1F("d4", "", nbins, xmin, xmax);


   TString dir = gROOT->GetTutorialDir();

   TString file = dir + "/legacy/spectrum/TSpectrum.root";

   TFile *f = new TFile(file.Data());

   TH1F *back = (TH1F *)f->Get("back2");


   back->SetTitle("Influence of clipping window width on the estimated background");

   back->SetAxisRange(0, 1000);

   back->SetMaximum(7000);

   back->Draw("L");


   TSpectrum *s = new TSpectrum();


   for (i = 0; i < nbins; i++)

      source[i] = back->GetBinContent(i + 1);

   s->Background(source, nbins, 10, TSpectrum::kBackDecreasingWindow, TSpectrum::kBackOrder2, kFALSE,

                 TSpectrum::kBackSmoothing3, kFALSE);

   for (i = 0; i < nbins; i++)

      d1->SetBinContent(i + 1, source[i]);

   d1->SetLineColor(kRed);

   d1->Draw("SAME L");


   for (i = 0; i < nbins; i++)

      source[i] = back->GetBinContent(i + 1);

   s->Background(source, nbins, 20, TSpectrum::kBackDecreasingWindow, TSpectrum::kBackOrder2, kFALSE,

                 TSpectrum::kBackSmoothing3, kFALSE);

   for (i = 0; i < nbins; i++)

      d2->SetBinContent(i + 1, source[i]);

   d2->SetLineColor(kBlue);

   d2->Draw("SAME L");


   for (i = 0; i < nbins; i++)

      source[i] = back->GetBinContent(i + 1);

   s->Background(source, nbins, 30, TSpectrum::kBackDecreasingWindow, TSpectrum::kBackOrder2, kFALSE,

                 TSpectrum::kBackSmoothing3, kFALSE);

   for (i = 0; i < nbins; i++)

      d3->SetBinContent(i + 1, source[i]);

   d3->SetLineColor(kGreen);

   d3->Draw("SAME L");


   for (i = 0; i < nbins; i++)

      source[i] = back->GetBinContent(i + 1);

   s->Background(source, nbins, 40, TSpectrum::kBackDecreasingWindow, TSpectrum::kBackOrder2, kFALSE,

                 TSpectrum::kBackSmoothing3, kFALSE);

   for (i = 0; i < nbins; i++)

      d4->SetBinContent(i + 1, source[i]);

   d4->SetLineColor(kMagenta);

   d4->Draw("SAME L");

}

f
#define f(i)
Definition RSha256.hxx:104

Int_t
int Int_t
Definition RtypesCore.h:45

kFALSE
constexpr Bool_t kFALSE
Definition RtypesCore.h:94

Double_t
double Double_t
Definition RtypesCore.h:59

kRed
@ kRed
Definition Rtypes.h:66

kGreen
@ kGreen
Definition Rtypes.h:66

kMagenta
@ kMagenta
Definition Rtypes.h:66

kBlue
@ kBlue
Definition Rtypes.h:66

xmin
float xmin
Definition THbookFile.cxx:95

xmax
float xmax
Definition THbookFile.cxx:95

gROOT
#define gROOT
Definition TROOT.h:406

TAttLine::SetLineColor
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition TAttLine.h:40

TFile
A ROOT file is an on-disk file, usually with extension .root, that stores objects in a file-system-li...
Definition TFile.h:53

TH1F
1-D histogram with a float per channel (see TH1 documentation)
Definition TH1.h:623

TH1::SetTitle
void SetTitle(const char *title) override
Change/set the title.
Definition TH1.cxx:6747

TH1::SetMaximum
virtual void SetMaximum(Double_t maximum=-1111)
Definition TH1.h:405

TH1::Draw
void Draw(Option_t *option="") override
Draw this histogram with options.
Definition TH1.cxx:3068

TH1::SetBinContent
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:9255

TH1::GetBinContent
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
Definition TH1.cxx:5090

TH1::SetAxisRange
virtual void SetAxisRange(Double_t xmin, Double_t xmax, Option_t *axis="X")
Set the "axis" range.
Definition Haxis.cxx:201

TSpectrum
Advanced Spectra Processing.
Definition TSpectrum.h:18

TSpectrum::Background
virtual TH1 * Background(const TH1 *hist, Int_t niter=20, Option_t *option="")
One-dimensional background estimation function.
Definition TSpectrum.cxx:145

TSpectrum::kBackOrder2
@ kBackOrder2
Definition TSpectrum.h:37

TSpectrum::kBackDecreasingWindow
@ kBackDecreasingWindow
Definition TSpectrum.h:42

TSpectrum::kBackSmoothing3
@ kBackSmoothing3
Definition TSpectrum.h:43

TString
Basic string class.
Definition TString.h:139

TString::Data
const char * Data() const
Definition TString.h:376