TSpectrum.root file

ROOT
1

Hi everyone,

I’m beginner in ROOT. I have a gamma ray spectrum which needs to be deconvoluted, for this purpose I try to use TSpectrum class with its example .x Deconvolution.C (taken from this website root.cern.ch/root/html/TSpectrum.html). However, I receive an error message “Error in TFile::TFile: file spectra\TSpectrum.root does not exist”. My question is, doesn’t my root version 6.04.00 include TSpectrum class? What should I do?

Thanks in advance…

Nergis.

2

I solved my problem.

Nergis.

3

Hello!

Root version is 5.34/36.
My system: Scientific Linux release 7.2 (Nitrogen).
My gcc: gcc (GCC) 4.8.5 20150623 (Red Hat 4.8.5-4).

If I launch this code from root.cern.ch/doc/master/classTSpectrum.html

void Deconvolution() {
   Int_t i;
   const Int_t nbins = 256;
   Double_t xmin     = 0;
   Double_t xmax     = nbins;
   Double_t source[nbins];
   Double_t response[nbins];
   gROOT->ForceStyle();
   TString dir  = gROOT->GetTutorialsDir();
   TString file = dir+"/spectrum/TSpectrum.root";
   TFile *f     = new TFile(file.Data());
   TH1F *h = (TH1F*) f->Get("decon1");
   h->SetTitle("Deconvolution");
   TH1F *d = (TH1F*) f->Get("decon_response");
   for (i = 0; i < nbins; i++) source[i]=h->GetBinContent(i + 1);
   for (i = 0; i < nbins; i++) response[i]=d->GetBinContent(i + 1);
   h->SetMaximum(30000);
   h->Draw("L");
   TSpectrum *s = new TSpectrum();
   s->Deconvolution(source,response,256,1000,1,1);
   for (i = 0; i < nbins; i++) d->SetBinContent(i + 1,source[i]);
   d->SetLineColor(kRed);
   d->Draw("SAME L");
}

I have this error:

Error in <TFile::TFile>: file /home/darkside/Vlad_Programs/root/tutorials/spectrum/TSpectrum.root does not exist
Error: illegal pointer to class object h 0x0 163  unfolding.cpp:14:
*** Interpreter error recovered ***

So, where to find TSpectrum.root?

I do not have this file

pwd
/home/darkside/Vlad_Programs/root
find -name TSpectrum*
./hist/spectrum/inc/TSpectrum.h
./hist/spectrum/inc/TSpectrum2.h
./hist/spectrum/inc/TSpectrum2Fit.h
./hist/spectrum/inc/TSpectrum2Transform.h
./hist/spectrum/inc/TSpectrum3.h
./hist/spectrum/inc/TSpectrumFit.h
./hist/spectrum/inc/TSpectrumTransform.h
./hist/spectrum/src/TSpectrum.cxx
./hist/spectrum/src/TSpectrum2.cxx
./hist/spectrum/src/TSpectrum2Fit.cxx
./hist/spectrum/src/TSpectrum2Transform.cxx
./hist/spectrum/src/TSpectrum3.cxx
./hist/spectrum/src/TSpectrumFit.cxx
./hist/spectrum/src/TSpectrumTransform.cxx
./hist/spectrum/src/TSpectrum2.d
./hist/spectrum/src/TSpectrum2.o
./hist/spectrum/src/TSpectrum.d
./hist/spectrum/src/TSpectrum.o
./hist/spectrum/src/TSpectrum2Transform.d
./hist/spectrum/src/TSpectrum2Transform.o
./hist/spectrum/src/TSpectrumTransform.d
./hist/spectrum/src/TSpectrumTransform.o
./hist/spectrum/src/TSpectrum3.d
./hist/spectrum/src/TSpectrum3.o
./hist/spectrum/src/TSpectrumFit.d
./hist/spectrum/src/TSpectrumFit.o
./hist/spectrum/src/TSpectrum2Fit.d
./hist/spectrum/src/TSpectrum2Fit.o
./hist/spectrumpainter/inc/TSpectrum2Painter.h
./hist/spectrumpainter/src/TSpectrum2Painter.cxx
./hist/spectrumpainter/src/TSpectrum2Painter.d
./hist/spectrumpainter/src/TSpectrum2Painter.o
./include/TSpectrum2.h
./include/TSpectrum.h
./include/TSpectrum2Transform.h
./include/TSpectrumTransform.h
./include/TSpectrum3.h
./include/TSpectrumFit.h
./include/TSpectrum2Fit.h
./include/TSpectrum2Painter.h

Thank you in advance.
Best regards, Vladislav.

4

Here it is
TSpectrum.root (31.9 KB)

5

Perfetto!
Many thanks!

6

github.com/root-mirror/root/tre … s/spectrum

7

Note that these tutorials are for ROOT 6.x and you cannot use them with ROOT 5.34.

8

See also:
root.cern/doc/master/group__tut … ctrum.html