50 printf(
"Generating histogram with %d peaks\n",npeaks);
65 for (
p=0;
p<npeaks;
p++) {
77 if (!
c1)
c1 =
new TCanvas(
"c1",
"c1",10,10,1000,700);
87 for (
p=0;
p<npeaks;
p++) {
88 for (pf=0;pf<nfound;pf++) {
91 if (diffx < 2*dx && diffy < 2*dy) ngood++;
94 if (ngood > nfound) ngood = nfound;
97 for (pf=0;pf<nfound;pf++) {
99 for (
p=0;
p<npeaks;
p++) {
102 if (diffx < 2*dx && diffy < 2*dy) nf++;
104 if (nf == 0) nghost++;
109 printf(
"Gener=%d, Found=%d, Good=%d, Ghost=%d\n",npeaks,nfound,ngood,nghost);
110 if (!
gROOT->IsBatch()) {
111 printf(
"\nDouble click in the bottom right corner of the pad to continue\n");
115void peaks2(
Int_t maxpeaks=50) {
117 for (
int i=0; i<10; ++i) {
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
R__EXTERN TRandom * gRandom
virtual void SetNpx(Int_t npx=100)
Set the number of points used to draw the function.
virtual void SetParameters(const Double_t *params)
A 2-Dim function with parameters.
virtual void SetNpy(Int_t npy=100)
Set the number of points used to draw the function.
virtual void SetStats(Bool_t stats=kTRUE)
Set statistics option on/off.
2-D histogram with a float per channel (see TH1 documentation)}
void FillRandom(const char *fname, Int_t ntimes=5000, TRandom *rng=nullptr) override
Fill histogram following distribution in function fname.
TObject * FindObject(const char *name) const override
Search if object named name is inside this pad or in pads inside this pad.
virtual Double_t Uniform(Double_t x1=1)
Returns a uniform deviate on the interval (0, x1).
Advanced 2-dimensional spectra processing.
Double_t * GetPositionY() const
Double_t * GetPositionX() const
void Print(Option_t *option="") const override
Print the array of positions.
virtual Int_t Search(const TH1 *hist, Double_t sigma=2, Option_t *option="", Double_t threshold=0.05)
This function searches for peaks in source spectrum in hin The number of found peaks and their positi...
Double_t Gaus(Double_t x, Double_t mean=0, Double_t sigma=1, Bool_t norm=kFALSE)
Calculates a gaussian function with mean and sigma.
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.