A class for generating a r-t and t-r graphs by binning the 2D histograms in Calibrator::rtHists in r ot t bins. More...

#include <Calibrator.h>

Collaboration diagram for RtGraph:

Public Member Functions
	RtGraph (TH2F , int, const char , bool, TDirectory *)
	the constructor More...

	~RtGraph ()
	the destructor More...

	RtGraph (const RtGraph &other)=delete
	copy constructor More...

RtGraph &	operator= (const RtGraph &other)=delete
	assignment More...

Public Attributes
std::vector< double >	rval
	array of the histograms for all bins More...

std::vector< double >	tval
	the r values More...

TGraphErrors *	trgr
	the t(r) graph More...

TGraphErrors *	rtgr
	the r(t) graph More...

int	npoints
	the number of graph points More...

float	mintime
	the minimum t-value More...

Private Types
enum	bintype { LOW, GOOD, HIGH, EMPTY }

Private Attributes
std::string	m_chname
	the histogram name More...

std::string	m_chtit
	the histogram title More...

bintype *	m_btype

float *	m_rightsig

float *	m_leftsig

float *	m_maxval

float *	m_leftval

float *	m_rightval

float	m_mean

float	m_mindistance

int *	m_maxbin

int	m_ipoint

float	m_t

float	m_d

float	m_et

float	m_ed

float *	m_tv

float *	m_dv

float *	m_etv

float *	m_edv

Detailed Description

A class for generating a r-t and t-r graphs by binning the 2D histograms in Calibrator::rtHists in r ot t bins.

Definition at line 34 of file Calibrator.h.

Member Enumeration Documentation

◆ bintype

enum RtGraph::bintype

private

Enumerator
LOW
GOOD
HIGH
EMPTY

Definition at line 51 of file Calibrator.h.

51 {LOW, GOOD, HIGH, EMPTY};

Constructor & Destructor Documentation

◆ RtGraph() [1/2]

RtGraph::RtGraph	(	TH2F *	rtHist,
		int	binvar,
		const char *	binlabel,
		bool	pflag,
		TDirectory *	dir
	)

the constructor

Definition at line 118 of file Calibrator.cxx.

                                                                                            {
  
   int fitresult;
   npoints = binvar==0 ? rtHist->GetNbinsX() : rtHist->GetNbinsY() ;
  
  
   m_mean = -20.0; 
   m_t = -20.0;
   m_d = -20.0;
   m_et = -20.0;
   m_ed = -20.0;
   //  m_ff =  nullptr ;
  
   
   TH1D** hslizes = new TH1D*[npoints];
   m_btype = new bintype[npoints];
   m_tv = new float[npoints];
   m_dv = new float[npoints];
   m_etv = new float[npoints];
   m_edv = new float[npoints];
   m_rightsig = new float[npoints];
   m_leftsig = new float[npoints];
   m_leftval = new float[npoints];
   m_rightval = new float[npoints];
   m_maxbin = new int[npoints];
   m_maxval = new float[npoints];
   
   m_ipoint=0;
   mintime=999;
   m_mindistance=0;
   
   if (pflag){
     TDirectory* binhistdir = dir->mkdir("binhist");
     binhistdir->cd();
     std::cout << " Calibrator::RtGraph: " << " Directory " << std::string(binhistdir->GetPath()) << " binlabel " << std::string(binlabel) << std::endl;
   }
  
   std::cout << " Calibrator::RtGraph: " << " number of points " << npoints << " binvar " << binvar << std::endl;
   // check and classify the bin-histograms
   for (int i=0;i<npoints;i++) {
     
     m_chname = std::string(Form("%s slize_%i",binlabel,i));
     m_chtit = binvar==0 ? std::string(Form("bin %i : %4.2f < %s < %4.2f  [ns]",i,rtHist->GetXaxis()->GetBinLowEdge(i+1),binlabel,rtHist->GetXaxis()->GetBinUpEdge(i+1))) :
                         std::string(Form("bin %i : %4.2f < %s < %4.2f  [mm]",i,rtHist->GetYaxis()->GetBinLowEdge(i+1),binlabel,rtHist->GetYaxis()->GetBinUpEdge(i+1))) ;
     hslizes[i] = binvar==0 ? rtHist->ProjectionY(m_chname.data(),i+1,i+1) :
                              rtHist->ProjectionX(m_chname.data(),i+1,i+1);
     hslizes[i]->SetTitle(m_chtit.data());
     if (!pflag){
       hslizes[i]->SetDirectory(nullptr);
     }
     
     m_maxbin[i]=1; m_maxval[i]=0;
     for (int j=1;j<=hslizes[i]->GetNbinsX();j++) {
       float val=hslizes[i]->GetBinContent(j);
       if (val>m_maxval[i]){
         m_maxval[i]=val;
         m_maxbin[i]=j;
         m_leftval[i] = hslizes[i]->GetBinContent(j-1);
         m_rightval[i] = hslizes[i]->GetBinContent(j+1);
       }
     }
     
     if(m_maxval[i]==0) m_maxval[i]=1.0e-4;
     m_rightsig[i] = m_rightval[i]==0 ? 100 : std::sqrt(m_maxval[i])/m_maxval[i] + std::sqrt(m_rightval[i])/m_rightval[i];
     m_leftsig[i]  = m_leftval[i]==0 ? 100 : std::sqrt(m_maxval[i])/m_maxval[i] + std::sqrt(m_leftval[i])/m_leftval[i];
     
     m_btype[i]=EMPTY;
     if(true){
       m_btype[i]=GOOD;
       m_btype[i-1]=GOOD;
     } 
     else{
       if (hslizes[i]->GetEntries()>0){
         if (i<(float)hslizes[i]->GetNbinsX()/2)  m_btype[i]=LOW;
         else m_btype[i]=HIGH;
       }
     }
     int rtEntries = rtHist->GetEntries();
     if(rtEntries==0) rtEntries=1;
     printf("%s ... %8f %8i %4i %8f %8f %4i\n", m_chtit.data(), (float)hslizes[i]->GetEntries()/(float)rtEntries, (int)hslizes[i]->GetEntries(), m_maxbin[i], m_leftsig[i], m_rightsig[i], m_btype[i]);
   }
   
   float frmin=0,frmax=0;
   for (int i=0;i<npoints;i++) {
     
     // prepare initial fit parameters   
     m_mean=hslizes[i]->GetBinCenter(m_maxbin[i]);
     if (m_btype[i]==LOW) {frmin = hslizes[i]->GetBinCenter(1);        frmax = hslizes[i]->GetBinCenter(m_maxbin[i]+4);}
     if (m_btype[i]==GOOD) {frmin = hslizes[i]->GetBinCenter(m_maxbin[i]-3); frmax = hslizes[i]->GetBinCenter(m_maxbin[i]+3);}
     if (m_btype[i]==HIGH) {frmin = hslizes[i]->GetBinCenter(m_maxbin[i]-3); frmax = hslizes[i]->GetBinCenter(m_maxbin[i]+4);}
     
     std::unique_ptr<TF1> ff(new  TF1("dtfitfunc","gaus"));
     
     ff->SetRange(frmin,frmax);
  
     m_t = binvar==0 ? rtHist->GetXaxis()->GetBinCenter(i+1) : rtHist->GetYaxis()->GetBinCenter(i+1);
     m_et = binvar==0 ? rtHist->GetXaxis()->GetBinWidth(1)/std::sqrt(12.0) : rtHist->GetYaxis()->GetBinWidth(1)/std::sqrt(12.0);
     m_d = 0;
     m_ed = 0;
     
     if (m_btype[i]==LOW){
       m_d =  0;
       m_ed = 0;
     }
     if (m_btype[i]==HIGH){
       m_d =  2;
       m_ed = 0;
     }
     if (m_btype[i]==EMPTY){
       m_d =  0;
       m_ed = 0;
     }
     if (m_btype[i]==GOOD && hslizes[i]->GetEntries()>3){
       if (m_btype[i-1]==GOOD){
         fitresult=hslizes[i]->Fit("dtfitfunc","QR");
         ff->SetRange(m_mean-1.0*ff->GetParameter(2), m_mean+1.0*ff->GetParameter(2));
         fitresult=hslizes[i]->Fit("dtfitfunc","QR");
         m_d =  ff->GetParameter(1);
         m_ed = ff->GetParError(1);
         std::cout << fitresult << " " << m_ed << std::endl;
       }
       else{
         m_d =  0;
         m_ed = 0;
       }
     }
  
  
     if (m_btype[i]==GOOD && m_ed!=0){
       m_tv[m_ipoint]=m_t;
       m_dv[m_ipoint]=m_d;
       m_etv[m_ipoint]=m_et;
       m_edv[m_ipoint]=m_ed;
       m_ipoint++;
       if (binvar==0) {
         rval.push_back(m_d);
         tval.push_back(m_t);
       }
       else {
         rval.push_back(m_t);
         tval.push_back(m_d);
       }
     }
  
     
     if(m_d>=m_mindistance && m_t<mintime) mintime = m_t; 
  
   }  
  
   rtgr = binvar==0 ? new TGraphErrors(m_ipoint,m_tv,m_dv,m_etv,m_edv) :  new TGraphErrors(m_ipoint,m_dv,m_tv,m_edv,m_etv);
   rtgr->SetTitle((std::string("binning in ") + binlabel).data());
   trgr = binvar==0 ? new TGraphErrors(m_ipoint,m_dv,m_tv,m_edv,m_etv) :  new TGraphErrors(m_ipoint,m_tv,m_dv,m_etv,m_edv);
   trgr->SetTitle((std::string("binning in ") + binlabel).data());
  
   rtgr->SetMarkerStyle(1) ;
   rtgr->SetMarkerColor(2) ;
   rtgr->SetLineColor(2) ;
   rtgr->SetName("rtgraph") ;
   trgr->SetMarkerStyle(1) ;
   trgr->SetMarkerColor(2) ;
   trgr->SetLineColor(2) ;
   trgr->SetName("trgraph") ;
  
   dir->cd();
   delete [] hslizes;
 }

◆ ~RtGraph()

RtGraph::~RtGraph ( )

the destructor

Definition at line 287 of file Calibrator.cxx.

                  {  
  
  delete [] m_btype     ;
  delete [] m_tv        ;
  delete [] m_dv        ;
  delete [] m_etv       ;
  delete [] m_edv       ;
  delete [] m_rightsig  ;
  delete [] m_leftsig   ;
  delete [] m_leftval   ;
  delete [] m_rightval  ;
  delete [] m_maxbin    ;
  delete [] m_maxval    ;
 }

◆ RtGraph() [2/2]

RtGraph::RtGraph ( const RtGraph & other )

delete

copy constructor

Member Function Documentation

◆ operator=()

RtGraph& RtGraph::operator= ( const RtGraph & other )

delete

assignment

Member Data Documentation

◆ m_btype

bintype* RtGraph::m_btype

private

Definition at line 52 of file Calibrator.h.

◆ m_chname

std::string RtGraph::m_chname

private

the histogram name

Definition at line 49 of file Calibrator.h.

◆ m_chtit

std::string RtGraph::m_chtit

private

the histogram title

Definition at line 50 of file Calibrator.h.

◆ m_d

float RtGraph::m_d

private

Definition at line 63 of file Calibrator.h.

◆ m_dv

float* RtGraph::m_dv

private

Definition at line 67 of file Calibrator.h.

◆ m_ed

float RtGraph::m_ed

private

Definition at line 65 of file Calibrator.h.

◆ m_edv

float* RtGraph::m_edv

private

Definition at line 69 of file Calibrator.h.

◆ m_et

float RtGraph::m_et

private

Definition at line 64 of file Calibrator.h.

◆ m_etv

float* RtGraph::m_etv

private

Definition at line 68 of file Calibrator.h.

◆ m_ipoint

int RtGraph::m_ipoint

private

Definition at line 61 of file Calibrator.h.

◆ m_leftsig

float* RtGraph::m_leftsig

private

Definition at line 54 of file Calibrator.h.

◆ m_leftval

float* RtGraph::m_leftval

private

Definition at line 56 of file Calibrator.h.

◆ m_maxbin

int* RtGraph::m_maxbin

private

Definition at line 60 of file Calibrator.h.

◆ m_maxval

float* RtGraph::m_maxval

private

Definition at line 55 of file Calibrator.h.

◆ m_mean

float RtGraph::m_mean

private

Definition at line 58 of file Calibrator.h.

◆ m_mindistance

float RtGraph::m_mindistance

private

Definition at line 59 of file Calibrator.h.

◆ m_rightsig

float* RtGraph::m_rightsig

private

Definition at line 53 of file Calibrator.h.

◆ m_rightval

float* RtGraph::m_rightval

private

Definition at line 57 of file Calibrator.h.

◆ m_t

float RtGraph::m_t

private

Definition at line 62 of file Calibrator.h.

◆ m_tv

float* RtGraph::m_tv

private

Definition at line 66 of file Calibrator.h.

◆ mintime

float RtGraph::mintime

the minimum t-value

Definition at line 47 of file Calibrator.h.

◆ npoints

int RtGraph::npoints

the number of graph points

Definition at line 46 of file Calibrator.h.

◆ rtgr

TGraphErrors* RtGraph::rtgr

the r(t) graph

Definition at line 45 of file Calibrator.h.

◆ rval

std::vector<double> RtGraph::rval

array of the histograms for all bins

the r values

Definition at line 42 of file Calibrator.h.

◆ trgr

TGraphErrors* RtGraph::trgr

the t(r) graph

Definition at line 44 of file Calibrator.h.

◆ tval

std::vector<double> RtGraph::tval

the r values

Definition at line 43 of file Calibrator.h.

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

Public Member Functions

Public Attributes

Private Types

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ bintype

Constructor & Destructor Documentation

◆ RtGraph() [1/2]

◆ ~RtGraph()

◆ RtGraph() [2/2]

Member Function Documentation

◆ operator=()

Member Data Documentation

◆ m_btype

◆ m_chname

◆ m_chtit

◆ m_d

◆ m_dv

◆ m_ed

◆ m_edv

◆ m_et

◆ m_etv

◆ m_ipoint

◆ m_leftsig

◆ m_leftval

◆ m_maxbin

◆ m_maxval

◆ m_mean

◆ m_mindistance

◆ m_rightsig

◆ m_rightval

◆ m_t

◆ m_tv

◆ mintime

◆ npoints

◆ rtgr

◆ rval

◆ trgr

◆ tval