EdbPeakProb Class Reference

#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
TH2F *	eHD
TH2F *	eHpeak
Int_t	ePeakNmax
Int_t	ePeakNmin
Float_t	ePeakRMS
Float_t	ePeakRMSmax
Float_t	ePeakRMSmin
Float_t	ePeakX
Float_t	ePeakY
Float_t	eProb
Int_t	eVerbosity

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

eHD

TH2F* EdbPeakProb::eHD

eHpeak

TH2F* EdbPeakProb::eHpeak

ePeakNmax

Int_t EdbPeakProb::ePeakNmax

ePeakNmin

Int_t EdbPeakProb::ePeakNmin

ePeakRMS

Float_t EdbPeakProb::ePeakRMS

ePeakRMSmax

Float_t EdbPeakProb::ePeakRMSmax

ePeakRMSmin

Float_t EdbPeakProb::ePeakRMSmin

ePeakX

Float_t EdbPeakProb::ePeakX

ePeakY

Float_t EdbPeakProb::ePeakY

eProb

Float_t EdbPeakProb::eProb

eVerbosity

Int_t EdbPeakProb::eVerbosity

---

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

• src/libAlignment/EdbPeakProb.h
• src/libAlignment/EdbPeakProb.cxx