prob the peak probability (2-dim only for the moment) More...

#include <EdbPeakProb.h>

Inheritance diagram for EdbPeakProb:

Collaboration diagram for EdbPeakProb:

Public Member Functions
	EdbPeakProb ()

Float_t	Prob ()

Float_t	Probability2D (TH2F *hd, float &xpeak, float &ypeak)

virtual	~EdbPeakProb ()

Public Attributes
TH1F *	eHbin
	service spectrum histo More...

TH2F *	eHD
	input histogram to be analysed More...

TH2F *	eHpeak
	service peak histo More...

Int_t	ePeakNmax
	apriori limits in the number of entries to accept the found peak More...

Int_t	ePeakNmin
	apriori limits in the number of entries to accept the found peak More...

Float_t	ePeakRMS
	found peak parameters More...

Float_t	ePeakRMSmax
	apriori limits in RMS to accept the found peak More...

Float_t	ePeakRMSmin
	apriori limits in RMS to accept the found peak More...

Float_t	ePeakX

Float_t	ePeakY

Float_t	eProb
	the probability the found peak is not the BG fluctuation More...

Int_t	eVerbosity
	0-no any message, 1-print, 2-plot; default=2 More...

Detailed Description

prob the peak probability (2-dim only for the moment)

Constructor & Destructor Documentation

◆ EdbPeakProb()

EdbPeakProb::EdbPeakProb ( )

◆ ~EdbPeakProb()

EdbPeakProb::~EdbPeakProb ( )

virtual

{
  SafeDelete(eHbin);
  SafeDelete(eHpeak);
}

Member Function Documentation

◆ Prob()

Float_t EdbPeakProb::Prob ( )

{
  float prob = Probability2D( eHD, ePeakX, ePeakY);
  return prob;
}

◆ Probability2D()

float EdbPeakProb::Probability2D	(	TH2F *	hd,
		float &	xpeak,
		float &	ypeak
	)

estimate the probability that the highest bin of the histogram do not belong to the background
Method: estimate the deviation of the tail of the coincidence distribution from exponential shape

{
 
  float prob=0;
 
  float maxbin = hd->GetMaximum();
  eHbin = new TH1F("hbin","hbin", (int)maxbin+1, 0, maxbin+1);
  int nx = hd->GetNbinsX();
  int ny = hd->GetNbinsY();
  for(int ix=1; ix<=nx; ix++)
    for(int iy=1; iy<=ny; iy++)
      eHbin->Fill( hd->GetBinContent(ix,iy) );
 
  float tailbin =  eHbin->GetMean()+eHbin->GetRMS();
  eHbin->Fit("expo","QR0","",tailbin,maxbin);
  float a = eHbin->GetFunction("expo")->GetParameter(0);
  float b = eHbin->GetFunction("expo")->GetParameter(1);
 
  float lowlim = (int)((-a)/b+0.99999);
  if(lowlim<1) lowlim=1;
 
  float real = eHbin->GetBinContent((int)maxbin+1);
  if( lowlim>maxbin )    prob=0;
  else                   prob = 1.- eHbin->GetFunction("expo")->Integral(maxbin,2*(maxbin+1)+10) / real;
 
  // calculate mean and rms for bins>lowlim to check the peak integrity
 
  eHpeak = (TH2F*)hd->Clone("hdpeak");
  eHpeak->Reset();
 
  float bin;
  for(int ix=1; ix<=nx; ix++)
    for(int iy=1; iy<=ny; iy++) {
      bin  = hd->GetBinContent(hd->GetBin(ix, iy));
      if(bin<lowlim) continue;
      eHpeak->Fill( hd->GetXaxis()->GetBinCenter( ix ),
            hd->GetYaxis()->GetBinCenter( iy ),
            bin );
    }
  Log(3, "AlignmentProbability", "peak: ( %f %f )   rms: ( %f %f )",  
     eHpeak->GetMean(1), eHpeak->GetMean(2), eHpeak->GetRMS(1),eHpeak->GetRMS(2));
 
  if(eVerbosity>1) {
    TCanvas *c = new TCanvas("aa","aa");
    c->Divide(2,2);
    c->cd(1); hd->Draw("colz");
    c->GetPad(2)->SetLogy();
    c->cd(2); eHbin->Draw();    eHbin->GetFunction("expo")->Draw("same");
    //c->cd(3); hd->SetMinimum(lowlim); hd->Draw();
    c->cd(3);
    eHpeak->Draw("colz");
  } else {
    SafeDelete(eHbin);
    SafeDelete(eHpeak);
  }
 
  Log(2,"Probability2D","%f   BGlim=%6.1f, maxbin=%6.1f",prob,lowlim,maxbin);
  return prob;
}

Member Data Documentation

◆ eHbin

TH1F* EdbPeakProb::eHbin

service spectrum histo

◆ eHD

TH2F* EdbPeakProb::eHD

input histogram to be analysed

◆ eHpeak

TH2F* EdbPeakProb::eHpeak

service peak histo

◆ ePeakNmax

Int_t EdbPeakProb::ePeakNmax

apriori limits in the number of entries to accept the found peak

◆ ePeakNmin

Int_t EdbPeakProb::ePeakNmin

apriori limits in the number of entries to accept the found peak

◆ ePeakRMS

Float_t EdbPeakProb::ePeakRMS

found peak parameters

◆ ePeakRMSmax

Float_t EdbPeakProb::ePeakRMSmax

apriori limits in RMS to accept the found peak

◆ ePeakRMSmin

Float_t EdbPeakProb::ePeakRMSmin

apriori limits in RMS to accept the found peak

◆ ePeakX

Float_t EdbPeakProb::ePeakX

◆ ePeakY

Float_t EdbPeakProb::ePeakY

◆ eProb

Float_t EdbPeakProb::eProb

the probability the found peak is not the BG fluctuation

◆ eVerbosity

Int_t EdbPeakProb::eVerbosity

0-no any message, 1-print, 2-plot; default=2

The documentation for this class was generated from the following files:

/home/antonio/fedra_doxygen/src/libAlignment/EdbPeakProb.h
/home/antonio/fedra_doxygen/src/libAlignment/EdbPeakProb.cxx

Public Member Functions

Public Attributes