Analyzes the signal of a nanobeacon in a PMT. More...

#include <NBPulse.hh>

Public Member Functions
	NBPulse ()
	Default constructor. More...

	NBPulse (TH2D *h)
	Constructor. More...

	~NBPulse ()
	Destructor. More...

double	Kolmogorov_Smirnov_uniform (TH1D *h)
	Performs a Kolmogorov-Smirnov compatibility test between the hit time distribution and a uniform distribution. More...

void	initialize ()
	Leaves the pulse ready for analysis by projecting the time vs ToT distribution to the time axis and moving it out of the Root file directory. More...

void	analyzeFast ()
	Performs a fast analysis of the hit time distribution. More...

void	computeBaseline (TH1 *h)
	Computes baseline of the hit time distribution measured by the PMT. More...

void	fit ()
	Fits the hit nanobeacon signal to a model composed by a Landau and a Gauss functions. More...

bool	IsGood ()
	Check if the peak is good. More...

bool	IsFitted ()
	Check if the peak is fitted. More...

bool	IsSaturatedHit ()
	Check if the hit time distribution is saturated. More...

bool	IsWeak ()
	Check if the hit time distribution is weak. More...

TH1D *	getHtime_full ()
	Get the hit time distribution. More...

TH1D *	gettime_peak ()
	Get the hit time distribution. More...

TH1D *	gettot_peak ()
	Get the hit time distribution. More...

TH2D *	gettime_vs_tot ()
	Get the hit time vs ToT distribution. More...

RooWorkspace	getWorkspace ()
	Get all the information from the fit. More...

double	getPeakTime ()
	Get the arrival time of the nanobeacon pulse. More...

double	getMeanToT ()
	Get the mean of the ToT distribution of the hits produced by the nanobeacon. More...

double	getSigmaToT ()
	Get the standard deviation of the ToT distribution of the hits produced by the nanobeacon. More...

string	getName ()
	Get name of the hit time vs ToT histogram. More...

Private Attributes
TH2D *	time_vs_tot

TH1D *	htime_full

TH1D *	htime_peak

TH1D *	tot_peak

RooWorkspace	w

pair< double, double >	baseline

double	tlow

double	thigh

double	tpeak

double	nhits

double	meanToT

double	sigmaToT

double	max_ToT

int	n_bins_right

int	n_bins_left

int	n_bins_before_peak

int	n_bins_afterpulse

int	max_n_bins_baseline

int	min_n_bins_baseline

int	n_sigma_peak

int	n_sigma_before_peak

double	fit_region_min

double	fit_region_max

double	fit_range_min

double	fit_range_max

double	landau_m_min

double	landau_m_max

double	landau_s_min

double	landau_s_max

double	gaus_m_min

double	gaus_m_max

double	gaus_s_min

double	gaus_s_max

double	l_g_ratio_min

double	l_g_ratio_max

double	KS_threshold

bool	analyzed

bool	fitted

bool	is_good

bool	saturated_hit

bool	saturated_tot

bool	weak

Detailed Description

Analyzes the signal of a nanobeacon in a PMT.

Definition at line 33 of file NBPulse.hh.

Constructor & Destructor Documentation

◆ NBPulse() [1/2]

NBPulse::NBPulse ( )

inline

Default constructor.

Definition at line 66 of file NBPulse.hh.

66 {}

◆ NBPulse() [2/2]

NBPulse::NBPulse ( TH2D * h )

inline

Constructor.

This constructor initializes the nanobeacon pulse from a 2D histogram containing the time vs ToT hit distribution.

Parameters

h	time vs ToT histogram

Definition at line 76 of file NBPulse.hh.

                    :
     max_ToT (26.4) ,
     n_bins_right (5) , 
     n_bins_left  (5) ,
     n_bins_before_peak(2) ,
     n_bins_afterpulse (10),
     max_n_bins_baseline (50),
     min_n_bins_baseline (10),
     n_sigma_peak (5) ,
     n_sigma_before_peak (4) ,
     fit_region_min (20.0), 
     fit_region_max (20.0) , 
     fit_range_min (6.0) , 
     fit_range_max (15.0) , 
     landau_m_min (2.0) , 
     landau_m_max (2.0) , 
     landau_s_min (1.0) , 
     landau_s_max (9.0) , 
     gaus_m_min (2.0) , 
     gaus_m_max (2.0) , 
     gaus_s_min (1.0) , 
     gaus_s_max (9.0) , 
     l_g_ratio_min (0.0) , 
     l_g_ratio_max (1.0) , 
     KS_threshold(1e-3) ,
     analyzed (false) ,
     fitted (false) ,
     is_good (false) ,
     saturated_hit (false) ,
     saturated_tot (false) , 
     weak (false) 
  
   {
  
     time_vs_tot = h ;
  
     initialize() ;
  
   }

◆ ~NBPulse()

NBPulse::~NBPulse ( )

inline

Destructor.

Definition at line 120 of file NBPulse.hh.

             {
  
     delete(htime_full) ;
  
     if (analyzed == true){
  
       delete(tot_peak) ; 
  
       delete(htime_peak) ;
  
     }
  
   }

Member Function Documentation

◆ Kolmogorov_Smirnov_uniform()

double NBPulse::Kolmogorov_Smirnov_uniform ( TH1D * h )

inline

Performs a Kolmogorov-Smirnov compatibility test between the hit time distribution and a uniform distribution.

Parameters

h	The hit time distribution

Returns: The probability that the hit time distribution is compatible with a uniform distribution.

Definition at line 141 of file NBPulse.hh.

                                             {
  
     TRandom3* r = new TRandom3() ;
  
     double low_x = 0 ;
  
     double high_x = h->GetBinCenter(h->GetNbinsX() - 1) ;
  
     int nentries = h->Integral(1 , h->GetNbinsX() - 1) ;
  
     if (h->GetBinCenter(1)>=0){
       
       low_x = h->GetBinCenter(1) ;
  
     }
  
     TH1D* u = (TH1D*)h->Clone("uniform") ;
  
     u->Reset() ;
  
     for (int i = 0 ; i<nentries ; i++){
  
       u->Fill(r->Uniform(low_x, high_x));
  
     }
     
     double proba = h->KolmogorovTest(u) ;
  
     delete(r) ;
  
     delete(u) ;
  
     return proba ;
  
   } ;

◆ initialize()

void NBPulse::initialize ( )

inline

Leaves the pulse ready for analysis by projecting the time vs ToT distribution to the time axis and moving it out of the Root file directory.

Definition at line 184 of file NBPulse.hh.

                    {
     
     htime_full = time_vs_tot->ProjectionX(time_vs_tot->GetName() , 1 , time_vs_tot->GetNbinsY());
   
     htime_full->SetDirectory(0) ;
  
     RooMsgService::instance().setGlobalKillBelow(RooFit::WARNING);
   
   }

◆ analyzeFast()

void NBPulse::analyzeFast ( )

inline

Performs a fast analysis of the hit time distribution.

Definition at line 198 of file NBPulse.hh.

                     {
  
     if (analyzed == true) return ;
  
     computeBaseline(htime_full) ;
     
     int peakbin = htime_full->GetMaximumBin() ;
  
     double height = htime_full->GetBinContent(peakbin) ;
  
     int leftbin = peakbin - n_bins_left ;
  
     int rightbin = peakbin + n_bins_right ;
  
     tpeak = htime_full->GetBinCenter(peakbin) ;
  
     tlow = htime_full->GetBinCenter(leftbin) ;
  
     thigh = htime_full->GetBinCenter(rightbin) ;
  
     htime_peak = (TH1D*)htime_full->Clone();
  
     for (int i=0 ; i<htime_peak->GetNbinsX() ; i++){
  
       if (i<leftbin || i>rightbin){
  
         htime_peak->SetBinContent(i,0);
  
       }
  
     }
  
     htime_peak->SetName((std::string("htime_")+htime_peak->GetName()).c_str());
  
     htime_peak->SetDirectory(0) ;
   
     tot_peak = (TH1D*)time_vs_tot->ProjectionY((std::string("ToT_") + time_vs_tot->GetName()).c_str() , leftbin , rightbin) ; 
  
     tot_peak->SetDirectory(0);
  
     meanToT = tot_peak->GetMean() ;
   
     if (meanToT > max_ToT){
     
       saturated_tot = true ;
     
     }
  
     sigmaToT = tot_peak->GetStdDev() ;
  
     double p_KS = Kolmogorov_Smirnov_uniform(htime_full) ;
  
     if (baseline.first==0 && baseline.second==0) weak = true ;
  
     if (p_KS > KS_threshold) weak = true ;
  
     if (height < (baseline.first + n_sigma_peak * baseline.second)) weak = true ;
  
     for (int i = 1 ; i<=n_bins_before_peak ; i++){
  
       if (htime_full->GetBinContent(peakbin - n_bins_before_peak) < (baseline.first + n_sigma_before_peak * baseline.second)) weak = true ;
  
     }
  
     if (weak == false && htime_full->GetBinContent(peakbin + n_bins_afterpulse) == 0) saturated_hit = true ;
  
     if (weak == false && saturated_hit == false) is_good = true ;
  
     nhits = htime_peak->Integral() ;
  
     analyzed = true ;
  
   }

◆ computeBaseline()

void NBPulse::computeBaseline ( TH1 * h )

inline

Computes baseline of the hit time distribution measured by the PMT.

For this the first bins of the distribution (where the nanobeacon pulse is absent) are used.

Parameters

h	Histogram with the hit time distribution

Definition at line 279 of file NBPulse.hh.

                               {
  
     int N = 0;
  
     double sum = 0 ;
  
     vector<double> bins ;
  
     for (int i=1 ; i<max_n_bins_baseline + 1 ; i++){
  
       if(h->GetBinContent(i) > 0 ){
  
         sum += h->GetBinContent(i);
  
         bins.push_back(h->GetBinContent(i));
  
         N += 1 ;
  
       }
  
     }
  
     if(N < min_n_bins_baseline){
  
       bins.clear();
  
       sum = 0 ;
  
       N = 0 ;
  
       for (int i=1 ; i<max_n_bins_baseline + 1 ; i++){
  
         if(h->GetBinContent(h->GetNbinsX() - i)>0){
  
           sum += h->GetBinContent(h->GetNbinsX() - i);
  
           bins.push_back(h->GetBinContent(h->GetNbinsX() - i));
  
           N += 1 ;
  
         }
  
       }
  
     }
  
     if(N < min_n_bins_baseline){
  
       baseline.first = 0 ;
  
       baseline.second = 0 ;
  
       return ;
  
     }
  
     double bg_mean = sum/N ;
  
     baseline.first = bg_mean ;
  
     double sdev = 0;
  
     for (int i=0 ; i<N ; i++){
     
       sdev += pow((bins[i] - bg_mean) , 2) ;
     
     }
  
     baseline.second = sqrt(sdev/N);
  
   }

◆ fit()

void NBPulse::fit ( )

inline

Fits the hit nanobeacon signal to a model composed by a Landau and a Gauss functions.

Definition at line 355 of file NBPulse.hh.

             {
  
     if (fitted==true)   return ;
  
     if (analyzed==false) analyzeFast() ;
  
     w.SetName((std::string("W_") + htime_full->GetName()).c_str());
  
     double tpeak = htime_full->GetBinCenter(htime_full->GetMaximumBin()) ;
  
     RooRealVar x("time" , "hit time" , tpeak - fit_region_min , tpeak + fit_region_max) ;
  
     RooRealVar ml("mean landau" , "mean landau" , tpeak - landau_m_min , tpeak + landau_m_max);
  
     RooRealVar sl("sigma landau" , "sigma landau" , landau_s_min , landau_s_max);
   
     RooLandau L1("L1" , "Main_Landau" , x , ml , sl) ;
     
     RooRealVar mg("mean Gauss" , "mean Gauss" , tpeak - gaus_m_min , tpeak + gaus_m_max);
  
     RooRealVar sg("sigma Gauss" , "sigma Gauss" , gaus_s_min , gaus_s_max);
  
     RooGaussian G1("G1" , "Main_Gaussian" , x , mg , sg) ;
   
     RooRealVar c1("c1" , "c1" , l_g_ratio_min , l_g_ratio_max);
  
     RooAddPdf M("Model" , "Model" , RooArgList(G1,L1) , RooArgList(c1));
  
     w.import(M) ;
  
     RooDataHist Data ("data" , "data" , *w.var("time") , Import(*htime_full));
   
     w.pdf("Model")->fitTo(Data , Range(tpeak - fit_range_min  , tpeak + fit_range_max) , RooFit::PrintLevel(-1000) , RooFit::PrintEvalErrors(-1000)) ;
  
     RooRealVar Tmin ("tmin" , "tmin" , tpeak - fit_range_min) ;
   
     RooRealVar Tmax ("tmax" , "tmax" , tpeak + fit_range_max) ;
   
     RooRealVar Tpeak ("tpeak" , "tpeak" , tpeak) ;
  
     RooRealVar Mpv("mpv" , "mpv" , Tmin.getVal()) ;
  
     RooRealVar Diff("peak-mpv" , "peak-mpv" , Tpeak.getVal() - Tmin.getVal()) ;
  
     double max_pdfval = 0 ;
  
     RooRealVar *var = (RooRealVar*)w.function("Model")->getObservables(RooArgSet(*w.var("time")))->first();
  
     for(double x = Tmin.getVal() ; x<Tmax.getVal() ; x+=0.1){
  
       var->setVal(x) ;
  
       double pdfval = w.function("Model")->getVal(RooArgSet(*var)) ;
  
       if (pdfval>max_pdfval){
  
         max_pdfval = pdfval ;
  
         Mpv.setVal(x) ;
  
       }
  
     }
  
     w.import(Mpv) ;
   
     w.import(Tmin) ;
  
     w.import(Tmax) ;
  
     w.import(Tpeak) ;
  
     w.import(Diff) ;
  
     fitted = true ;
  
   }

◆ IsGood()

bool NBPulse::IsGood ( )

inline

Check if the peak is good.

Returns: true if the peak is classified as good.

Definition at line 439 of file NBPulse.hh.

                 {
  
     return is_good ;
  
   }

◆ IsFitted()

bool NBPulse::IsFitted ( )

inline

Check if the peak is fitted.

Returns: true if the peak has been fitted.

Definition at line 451 of file NBPulse.hh.

                   {
  
     return fitted ;
  
   }

◆ IsSaturatedHit()

bool NBPulse::IsSaturatedHit ( )

inline

Check if the hit time distribution is saturated.

Returns: true if the hit time distributon is classified as saturated.

Definition at line 463 of file NBPulse.hh.

                        {
     
     return saturated_hit ;
   
   }

◆ IsWeak()

bool NBPulse::IsWeak ( )

inline

Check if the hit time distribution is weak.

Returns: true if the hit time distributon is classified as weak.

Definition at line 475 of file NBPulse.hh.

                {
  
     return weak ;
     
   }

◆ getHtime_full()

TH1D* NBPulse::getHtime_full ( )

inline

Get the hit time distribution.

Returns: histogram with hit time distribution

Definition at line 488 of file NBPulse.hh.

                        {
     
     return htime_full ;
     
   } 

◆ gettime_peak()

TH1D* NBPulse::gettime_peak ( )

inline

Get the hit time distribution.

Returns: histogram with hit time distribution

Definition at line 501 of file NBPulse.hh.

                       {
     
     return htime_peak ;
     
   }

◆ gettot_peak()

TH1D* NBPulse::gettot_peak ( )

inline

Get the hit time distribution.

Returns: histogram with hit time distribution

Definition at line 514 of file NBPulse.hh.

                      {
     
     return tot_peak ;
     
   }

◆ gettime_vs_tot()

TH2D* NBPulse::gettime_vs_tot ( )

inline

Get the hit time vs ToT distribution.

Returns: a 2d histogram with hit time vs ToT distribution

Definition at line 528 of file NBPulse.hh.

                         {
  
     return  time_vs_tot ;
     
   }

◆ getWorkspace()

RooWorkspace NBPulse::getWorkspace ( )

inline

Get all the information from the fit.

Returns: Roofit Workspace with all the information and parameters from the fit

Definition at line 541 of file NBPulse.hh.

◆ getPeakTime()

double NBPulse::getPeakTime ( )

inline

Get the arrival time of the nanobeacon pulse.

Returns: pulse arrival time

Definition at line 554 of file NBPulse.hh.

                        {
  
     return tpeak ;
     
   }

◆ getMeanToT()

double NBPulse::getMeanToT ( )

inline

Get the mean of the ToT distribution of the hits produced by the nanobeacon.

Returns: mean ToT value

Definition at line 567 of file NBPulse.hh.

                       {
     
     return meanToT ;
     
   }

◆ getSigmaToT()

double NBPulse::getSigmaToT ( )

inline

Get the standard deviation of the ToT distribution of the hits produced by the nanobeacon.

Returns: standard deviation from ToT distribution

Definition at line 580 of file NBPulse.hh.

                        {
  
     return sigmaToT ;  
     
   }

◆ getName()

string NBPulse::getName ( )

inline

Get name of the hit time vs ToT histogram.

Returns: The name of the histogram.

Definition at line 593 of file NBPulse.hh.

                    {
     
     return time_vs_tot->GetName();
     
   }

Member Data Documentation

◆ time_vs_tot

TH2D* NBPulse::time_vs_tot

private

Definition at line 35 of file NBPulse.hh.

◆ htime_full

TH1D* NBPulse::htime_full

private

Definition at line 37 of file NBPulse.hh.

◆ htime_peak

TH1D * NBPulse::htime_peak

private

Definition at line 37 of file NBPulse.hh.

◆ tot_peak

TH1D * NBPulse::tot_peak

private

Definition at line 37 of file NBPulse.hh.

◆ w

RooWorkspace NBPulse::w

private

Definition at line 39 of file NBPulse.hh.

◆ baseline

pair<double , double> NBPulse::baseline

private

Definition at line 41 of file NBPulse.hh.

◆ tlow

double NBPulse::tlow

private

Definition at line 43 of file NBPulse.hh.

◆ thigh

double NBPulse::thigh

private

Definition at line 43 of file NBPulse.hh.

◆ tpeak

double NBPulse::tpeak

private

Definition at line 43 of file NBPulse.hh.

◆ nhits

double NBPulse::nhits

private

Definition at line 43 of file NBPulse.hh.

◆ meanToT

double NBPulse::meanToT

private

Definition at line 43 of file NBPulse.hh.

◆ sigmaToT

double NBPulse::sigmaToT

private

Definition at line 43 of file NBPulse.hh.

◆ max_ToT

double NBPulse::max_ToT

private

Definition at line 43 of file NBPulse.hh.

◆ n_bins_right

int NBPulse::n_bins_right

private

Definition at line 45 of file NBPulse.hh.

◆ n_bins_left

int NBPulse::n_bins_left

private

Definition at line 45 of file NBPulse.hh.

◆ n_bins_before_peak

int NBPulse::n_bins_before_peak

private

Definition at line 45 of file NBPulse.hh.

◆ n_bins_afterpulse

int NBPulse::n_bins_afterpulse

private

Definition at line 45 of file NBPulse.hh.

◆ max_n_bins_baseline

int NBPulse::max_n_bins_baseline

private

Definition at line 45 of file NBPulse.hh.

◆ min_n_bins_baseline

int NBPulse::min_n_bins_baseline

private

Definition at line 45 of file NBPulse.hh.

◆ n_sigma_peak

int NBPulse::n_sigma_peak

private

Definition at line 45 of file NBPulse.hh.

◆ n_sigma_before_peak

int NBPulse::n_sigma_before_peak

private

Definition at line 45 of file NBPulse.hh.