Logo ROOT  
Reference Guide
No Matches
Src3.C File Reference

Detailed Description

View in nbviewer Open in SWAN Example to illustrate high resolution peak searching function (class TSpectrum2).

Found 13 candidate peaks
posx= 46, posy= 44, value=1962
posx= 39, posy= 4, value=1829
posx= 46, posy= 19, value=1762
posx= 50, posy= 32, value=1745
posx= 23, posy= 56, value=1710
posx= 35, posy= 48, value=1546
posx= 19, posy= 3, value=1524
posx= 39, posy= 58, value=1499
posx= 15, posy= 22, value=1376
posx= 38, posy= 16, value=1170
posx= 4, posy= 21, value=1166
posx= 10, posy= 16, value=1013
posx= 8, posy= 35, value=708
#include <TSpectrum2.h>
void Src3() {
Int_t i, j, nfound;
const Int_t nbinsx = 64;
const Int_t nbinsy = 64;
Double_t xmax = (Double_t)nbinsx;
Double_t ymax = (Double_t)nbinsy;
Double_t** source = new Double_t*[nbinsx];
for (i=0;i<nbinsx;i++)
source[i]=new Double_t[nbinsy];
Double_t** dest = new Double_t*[nbinsx];
for (i=0;i<nbinsx;i++)
dest[i]=new Double_t[nbinsy];
TString dir = gROOT->GetTutorialDir();
TString file = dir+"/spectrum/TSpectrum2.root";
TFile *f = new TFile(file.Data());
auto search = (TH2F*) f->Get("search1");
auto s = new TSpectrum2();
for (i = 0; i < nbinsx; i++){
for (j = 0; j < nbinsy; j++){
source[i][j] = search->GetBinContent(i + 1,j + 1);
nfound = s->SearchHighRes(source, dest, nbinsx, nbinsy, 2, 2, kFALSE, 3, kFALSE, 1);//3, 10, 100
printf("Found %d candidate peaks\n",nfound);
Double_t *PositionX = s->GetPositionX();
Double_t *PositionY = s->GetPositionY();
auto m = new TMarker();
for (i=0;i<nfound;i++) {
printf("posx= %d, posy= %d, value=%d\n",(Int_t)(PositionX[i]+0.5), (Int_t)(PositionY[i]+0.5),
#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
float xmin
float ymin
float xmax
float ymax
#define gROOT
Definition TROOT.h:406
R__EXTERN TStyle * gStyle
Definition TStyle.h:412
A ROOT file is a suite of consecutive data records (TKey instances) with a well defined format.
Definition TFile.h:54
2-D histogram with a float per channel (see TH1 documentation)}
Definition TH2.h:251
Manages Markers.
Definition TMarker.h:22
Advanced 2-dimensional spectra processing.
Definition TSpectrum2.h:18
Basic string class.
Definition TString.h:136
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
Definition TStyle.cxx:1589
Definition file.py:1
auto * m
Definition textangle.C:8
#define dest(otri, vertexptr)
Definition triangle.c:1040
Miroslav Morhac, Olivier Couet

Definition in file Src3.C.