Background function?

[Wile_E_Coyote]

Simple “erfc” based background examples (the peak is a symmetric gaussian):

How to find the area error using peakfind.C

Using the "60Co_calibrated.root" file (attached in another post here), try this: { const Bool_t fit_with_erfc_background = kTRUE; // kFALSE or kTRUE TFile *f = TFile::Open("60Co_calibrated.root"); TCanvas *C2; f->GetObject("C2", C2); C2->Draw(); TH1 *h = (TH1*)C2->FindObject("h_energy"); h->SetDirectory(gROOT); h->GetListOfFunctions()->Delete(); // "remove" old stuff h->SetLineColor(kViolet); // h->Rebin(4); // try out the "average X bin width" != 1 case // Double_…

How to find the area error using peakfind.C

I’m afraid the way you create histograms in the “call_spectra_1.C” routine may create problems with getting proper “area” values later. Using the 05cm_60Co.Spe.txt file (extracted and renamed from the "root.zip" file, attached in another post here), try to play with the macro given below: Spe_60Co.cxx