source: fact/tools/rootmacros/FPedestal.c@ 13637

#include <cstdio>
int FPedestal(char *title_base, fits &datafile, vector<int16_t> &data, vector<int16_t> &data_offset, vector<float> &drs_basemean, vector<float> &drs_gainmean, vector<float> &drs_triggeroffsetmean, UInt_t data_roi, size_t pixelnr, float &pedestal_mean, float &pedestal_rms)
{
//      TCanvas *canv_base = new TCanvas( "canv_base", "Determine the baseline", 100, 520, 700, 500 );
        TH1F *baseline = new TH1F("baseline", title_base, 1000,-10,40);
//      canv_base->SetLogy();
        baseline->GetXaxis()->SetTitle("Sample value (mV)");
        baseline->GetYaxis()->SetTitle("Entries");
        
//-------------------------------------------
//Find the baseline
//-------------------------------------------
        for (size_t i=0; i<datafile.GetNumRows(); i++)
        {
                datafile.GetRow(i);
                
                //Iterate over the slices
                for (UInt_t k=0; k<data_roi; k++)
                {
                        UInt_t drs_calib_offset = (k+data_offset[pixelnr])%data_roi;
                        float sample = (data[pixelnr*data_roi+k]*2000/4096.-drs_basemean[pixelnr*data_roi+drs_calib_offset]-drs_triggeroffsetmean[pixelnr*data_roi+k])/drs_gainmean[pixelnr*data_roi+drs_calib_offset]*1907.35;
                        baseline->Fill(sample);
                }
        }
        pedestal_mean = baseline->GetXaxis()->GetBinCenter(baseline->GetMaximumBin());
        pedestal_rms = baseline->GetRMS();
        std::cout << "Value of maximal probability: " << pedestal_mean << " +- " << pedestal_rms << std::endl;
        
//      canv_base->cd();
//      baseline->Draw();
//      canv_base->Modified();
//      canv_base->Update();
        
        return 0;
}