RtGraph Class Reference

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
	~RtGraph ()
	the destructor
	RtGraph (const RtGraph &other)=delete
	copy constructor
RtGraph &	operator= (const RtGraph &other)=delete
	assignment

Public Attributes
std::vector< double >	rval
	array of the histograms for all bins
std::vector< double >	tval
	the r values
TGraphErrors *	trgr
	the t(r) graph
TGraphErrors *	rtgr
	the r(t) graph
int	npoints
	the number of graph points
float	mintime
	the minimum t-value

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

Private Attributes
std::string	m_chname
	the histogram name
std::string	m_chtit
	the histogram title
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.

{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:

• Calibrator.h
• Calibrator.cxx