Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
Background_width.C File Reference

Detailed Description

View in nbviewer Open in SWAN Example to illustrate the influence of the clipping window width on the estimated background.

void Background_width() {
Int_t i;
const Int_t nbins = 1024;
Double_t xmax = nbins;
Double_t source[nbins];
gROOT->ForceStyle();
TString dir = gROOT->GetTutorialDir();
TString file = dir+"/spectrum/TSpectrum.root";
TFile *f = new TFile(file.Data());
TH1F *back = (TH1F*) f->Get("back1");
TH1F *d1 = new TH1F("d1","",nbins,xmin,xmax);
TH1F *d2 = new TH1F("d2","",nbins,xmin,xmax);
TH1F *d3 = new TH1F("d3","",nbins,xmin,xmax);
back->GetXaxis()->SetRange(1,nbins);
back->SetTitle("Influence of clipping window width on the estimated background");
back->Draw("L");
TSpectrum *s = new TSpectrum();
for (i = 0; i < nbins; i++) source[i]=back->GetBinContent(i + 1);
for (i = 0; i < nbins; i++) d1->SetBinContent(i + 1,source[i]);
d1->Draw("SAME L");
for (i = 0; i < nbins; i++) source[i]=back->GetBinContent(i + 1);
for (i = 0; i < nbins; i++) d2->SetBinContent(i + 1,source[i]);
d2->Draw("SAME L");
for (i = 0; i < nbins; i++) source[i]=back->GetBinContent(i + 1);
for (i = 0; i < nbins; i++) d3->SetBinContent(i + 1,source[i]);
d3->Draw("SAME L");
}
#define f(i)
Definition RSha256.hxx:104
int Int_t
Definition RtypesCore.h:45
const Bool_t kFALSE
Definition RtypesCore.h:92
double Double_t
Definition RtypesCore.h:59
@ kRed
Definition Rtypes.h:66
@ kOrange
Definition Rtypes.h:67
@ kGreen
Definition Rtypes.h:66
float xmin
float xmax
#define gROOT
Definition TROOT.h:406
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition TAttLine.h:40
virtual void SetRange(Int_t first=0, Int_t last=0)
Set the viewing range for the axis using bin numbers.
Definition TAxis.cxx:920
A ROOT file is a suite of consecutive data records (TKey instances) with a well defined format.
Definition TFile.h:54
1-D histogram with a float per channel (see TH1 documentation)}
Definition TH1.h:575
virtual void SetTitle(const char *title)
See GetStatOverflows for more information.
Definition TH1.cxx:6678
TAxis * GetXaxis()
Get the behaviour adopted by the object about the statoverflows. See EStatOverflows for more informat...
Definition TH1.h:320
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:9062
virtual void Draw(Option_t *option="")
Draw this histogram with options.
Definition TH1.cxx:3073
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
Definition TH1.cxx:4993
Advanced Spectra Processing.
Definition TSpectrum.h:18
@ kBackOrder2
Definition TSpectrum.h:37
@ kBackDecreasingWindow
Definition TSpectrum.h:42
@ kBackSmoothing3
Definition TSpectrum.h:43
virtual TH1 * Background(const TH1 *hist, Int_t niter=20, Option_t *option="")
One-dimensional background estimation function.
Basic string class.
Definition TString.h:136
Definition file.py:1
Authors
Miroslav Morhac, Olivier Couet

Definition in file Background_width.C.