Plots Class Reference

set of generic plots More...

#include <computils.h>

Inheritance diagram for Plots:

Collaboration diagram for Plots:

Public Member Functions
	Plots (const std::string &s="", bool errors=false)
double	realmin (double lo=0, double hi=0)
double	realmax (double lo=0, double hi=0)
void	MaxScale (double scale=1.1, double lo=0, double hi=0)
void	MinScale (double scale, double lo=0, double hi=0)
void	Min (double scale)
void	Max (double scale)
std::vector< double >	findxrange (bool symmetric=false)
void	sortx (const AxisInfo &xinfo)
double	lo () const
double	hi () const
void	xrange (bool symmetric=false)
void	SetRangeUser (double lo, double hi)
void	limits ()
void	Draw (Legend &leg, bool means=false)
void	Draw_i (Legend &leg, bool means=false)
void	SetLogx (bool b=true)
void	SetLogy (bool b=true)
std::string	GetXaxisTitle ()
std::string	GetYaxisTitle ()
void	SetXaxisTitle (const std::string &s)
void	SetYaxisTitle (const std::string &s)
void	push_back (const Plotter &val)
	this is so that we can update the stats as we go along, but no updating isn't done at the moment

Static Public Member Functions
static void	setwatermark (bool b)

Public Attributes
T	elements
	STL member.

Private Attributes
std::string	m_name
bool	m_logx
	canvas log setting
bool	m_logy
bool	m_maxset
	yaxis range setting
double	m_max
bool	m_minset
double	m_min
bool	m_rangeset
	xaxis range setting
double	m_lo
double	m_hi
bool	m_trim_errors

Static Private Attributes
static bool	s_watermark = true

Detailed Description

set of generic plots

Definition at line 1210 of file computils.h.

Constructor & Destructor Documentation

Plots()

Plots::Plots ( const std::string & s = "", bool errors = false )

inline

Definition at line 1214 of file computils.h.

                                               : 
    m_name(s), 
    m_logx(false), m_logy(false), 
    m_maxset(false), m_max(0),
    m_minset(false), m_min(0),
    m_rangeset(false), 
    m_lo(0.0), m_hi(0.0),
    m_trim_errors(errors)
  { }

Member Function Documentation

Draw()

void Plots::Draw ( Legend & leg, bool means = false )

inline

Definition at line 1452 of file computils.h.

                                           {  
    Draw_i( leg, means );
    if ( LINES ) { 
      LINES=false;
      Draw_i( leg, means );
      LINES=true;
    }
  }

Draw_i()

void Plots::Draw_i ( Legend & leg, bool means = false )

inline

increase the number of log labels if only a few decades

still can't get this working correctly - for the time being leave these alternative loop orders in place but commented out until we can find a better solution ...

Definition at line 1461 of file computils.h.

                                             {  
 
    bool first = true;
 
#if 1
    if ( m_logy ) {      
      for ( unsigned i=0 ; i<size() ; i++, first=false ) { 
    double ymax = at(i).htest()->GetMaximum();
    double ymin = at(i).htest()->GetMinimum();
    at(i).htest()->GetYaxis()->SetMoreLogLabels(true);
    if ( ymax/ymin>1e6 ) at(i).htest()->GetYaxis()->SetMoreLogLabels(false);
    break;
      }
    }
#endif
 
#if 1
    for ( unsigned i=0 ; i<size() ; i++ ) at(i).trim_errors( m_trim_errors );

 
#endif
    
    // for ( unsigned i=0 ; i<size() ; i++,  first=false ) at(i).Draw( i, &leg, means, first, (i==size()-1) );
    //    for ( unsigned i=size() ; i-- ;  first=false ) at(i).Draw( i, &leg, means, first, i==0 );
 
    for ( unsigned i=size() ; i-- ;  first=false ) at(i).Draw( i, &leg, means, first, i==0 );
 
    // dev versions ... 
    //  for ( unsigned i=0 ; i<size() ; i++,  first=false ) at(i).DrawLegend( i, leg, means, first, (i==size()-1) );
    // for ( unsigned i=0 ; i<size() ; i++,  first=false ) at(i).Draw( i, &leg, means, first, i==size() );
 
 
    
    if ( s_watermark ) DrawLabel(0.1, 0.02, "built "+stime()+release, kBlack, 0.03 );
 
    gPad->SetLogy(m_logy);
    gPad->SetLogx(m_logx);
 
  }

findxrange()

std::vector< double > Plots::findxrange ( bool symmetric = false )

inline

don't use empty histograms to find the bin limits unless they are all empty

Definition at line 1316 of file computils.h.

                                                     { 

    
    std::vector<double> v(2,0);
    
    TH1F* hf = at(0).htest();
 
    double vlo  =  1e21;
    double vhi  = -1e21;
    
    if ( hf->GetBinLowEdge(1)<vlo )                 vlo = hf->GetBinLowEdge(1);
    if ( hf->GetBinLowEdge(hf->GetNbinsX()+1)>vhi ) vhi = hf->GetBinLowEdge( hf->GetNbinsX()+1 );
      
    if ( size()>0 ) v = ::findxrangeuser( hf, symmetric );
 
    bool first = false;
 
    for ( unsigned i=1 ; i<size() ; i++ ) { 
  
      hf = at(i).htest();
 
      if ( ::empty( hf ) ) continue;
 
      if ( hf->GetBinLowEdge(1)<vlo )                 vlo = hf->GetBinLowEdge(1);
      if ( hf->GetBinLowEdge(hf->GetNbinsX()+1)>vhi ) vhi = hf->GetBinLowEdge( hf->GetNbinsX()+1 );
 
      std::vector<double> limits = ::findxrangeuser( hf, symmetric );
 
      double lo = limits[0];
      double hi = limits[1];
 
      if ( first ) { 
    v[0] = lo;
    v[1] = hi;
      }
      else { 
    if ( v[0]>lo ) v[0] = lo;
    if ( v[1]<hi ) v[1] = hi;
      }
 
      first = false;
    }
 
    double upper = ( v[1]-v[0] )*1.1 + v[0];
    double lower = v[0] - ( v[1]-v[0] )*0.1; 
 
    if ( m_logx ) {
      double dx = std::log10(v[1])-std::log10(v[0]);
      upper = std::pow(10,dx*1.1 + std::log10(v[0]));
      lower = std::pow(10,std::log10(v[0]) - dx*0.1); 
    }
 
    if ( lower<vlo ) lower = vlo;
    if ( upper>vhi ) upper = vhi;
 
    v[0] = lower;
    v[1] = upper;
    
    return v;
  }

GetXaxisTitle()

std::string Plots::GetXaxisTitle ( )

inline

Definition at line 1506 of file computils.h.

                            { 
    if ( size()>0 ) return at(0).htest()->GetXaxis()->GetTitle();
    return "";
  }

GetYaxisTitle()

std::string Plots::GetYaxisTitle ( )

inline

Definition at line 1511 of file computils.h.

                            { 
    if ( size()>0 ) return at(0).htest()->GetYaxis()->GetTitle();
    return "";
  }

hi()

double Plots::hi ( ) const

inline

Definition at line 1419 of file computils.h.

{ return m_hi; }

limits()

void Plots::limits ( )

inline

Definition at line 1442 of file computils.h.

                { 
    double rmax = realmax();
    double rmin = realmin();    
    if ( rmin<0 ) { 
      std::cout << "\tlimits \t" << m_name << "\tmin " << rmin << "\tmax " << rmax << std::endl; 
    }
  }

lo()

double Plots::lo ( ) const

inline

Definition at line 1418 of file computils.h.

{ return m_lo; }

Max()

void Plots::Max ( double scale )

inline

Definition at line 1308 of file computils.h.

                           { 
    m_maxset = true;
    for ( unsigned i=0 ; i<size() ; i++ ) {
      if ( at(i).href() ) at(i).href()->SetMaximum(scale);
      at(i).htest()->SetMaximum(scale);
    }
  }

MaxScale()

void Plots::MaxScale ( double scale = 1.1, double lo = 0, double hi = 0 )

inline

Definition at line 1249 of file computils.h.

                                                      {  
    
    if ( size()<1 ) return;
 
    double tmax = realmax(lo,hi);
    double tmin = realmin(lo,hi);
   
    m_max = scale*tmax;
    
    if ( m_logy ) m_min = tmin;
 
    for ( unsigned i=0 ; i<size() ; i++ ) {
      if ( at(i).href() ) at(i).href()->SetMaximum(scale*tmax);
      at(i).htest()->SetMaximum(scale*tmax);
    } 
 
  }

Min()

void Plots::Min ( double scale )

inline

Definition at line 1295 of file computils.h.

                           { 
    m_minset = true;
    for ( unsigned i=0 ; i<size() ; i++ ) {
      if ( at(i).href() ) at(i).href()->SetMinimum(scale);
      at(i).htest()->SetMinimum(scale);
      if ( m_logy ) { 
    if ( at(i).href() ) if ( at(i).href()->GetMinimum()<=0 )  at(i).href()->GetMinimum(1e-4);
    if ( at(i).htest()->GetMinimum()<=0 ) at(i).htest()->GetMinimum(1e-4);
      }
    }
  }

MinScale()

void Plots::MinScale ( double scale, double lo = 0, double hi = 0 )

inline

Definition at line 1268 of file computils.h.

                                                      { 
 
    if ( size()<1 ) return;
 
    if ( scale==0 ) { 
      for ( unsigned i=0 ; i<size() ; i++ ) {
    if ( at(i).href() ) at(i).href()->SetMinimum(0);
    at(i).htest()->SetMinimum(0);
      }  
      return;
    }
 
    
    double tmin = 0;
 
    if ( lo!=hi ) tmin = realmin( lo, hi );
    else          tmin = realmin();
 
    m_min = scale*tmin;
 
    for ( unsigned i=0 ; i<size() ; i++ ) {
      if ( at(i).href() ) at(i).href()->SetMinimum(scale*tmin);
      at(i).htest()->SetMinimum(scale*tmin);
    }
  }

push_back()

void Plots::push_back ( const Plotter & val )

inline

this is so that we can update the stats as we go along, but no updating isn't done at the moment

Definition at line 1527 of file computils.h.

                                     {
    std::vector<Plotter>::push_back( val );
  }

realmax()

double Plots::realmax ( double lo = 0, double hi = 0 )

inline

Definition at line 1235 of file computils.h.

                                       {
    bool first = true;
    double max = 0;
    for ( unsigned i=0 ; i<size() ; i++ ) {
      //      double rmref  = realmin( at(i).href(), false );
      double rmtest = ::realmax( at(i).htest(), false, lo, hi );
      if ( rmtest!=0 && ( first || max<rmtest ) ) max = rmtest;
      if ( rmtest!=0 ) first = false;
    }
    return max;
  }

realmin()

double Plots::realmin ( double lo = 0, double hi = 0 )

inline

Definition at line 1224 of file computils.h.

                                         {
    bool first = true;
    double min = 1000;
    for ( unsigned i=0 ; i<size() ; i++ ) {
      double rmtest = ::realmin( at(i).htest(), false, lo, hi );
      if ( rmtest!=0 && ( first || min>rmtest ) ) min = rmtest;
      if ( rmtest!=0 ) first = false;
    }
    return min;
  }

SetLogx()

void Plots::SetLogx ( bool b = true )

inline

Definition at line 1503 of file computils.h.

{ m_logx=b; } 

SetLogy()

void Plots::SetLogy ( bool b = true )

inline

Definition at line 1504 of file computils.h.

{ m_logy=b; } 

SetRangeUser()

void Plots::SetRangeUser ( double lo, double hi )

inline

Definition at line 1430 of file computils.h.

                                           { 
    m_rangeset = true;
    m_lo = lo;
    m_hi = hi;
    for ( unsigned i=0 ; i<size() ; i++ ) { 
      if ( at(i).href() ) at(i).href()->GetXaxis()->SetRangeUser( m_lo, m_hi );
      at(i).htest()->GetXaxis()->SetRangeUser( m_lo, m_hi );
      at(i).SetRange(m_lo, m_hi);
    }
  }

setwatermark()

void Plots::setwatermark ( bool b )

inlinestatic

Definition at line 1533 of file computils.h.

{ s_watermark = b; }

SetXaxisTitle()

void Plots::SetXaxisTitle ( const std::string & s )

inline

Definition at line 1516 of file computils.h.

                                         { 
    if ( size()>0 ) at(0).htest()->GetXaxis()->SetTitle(s.c_str());
  }

SetYaxisTitle()

void Plots::SetYaxisTitle ( const std::string & s )

inline

Definition at line 1520 of file computils.h.

                                         { 
    if ( size()>0 ) at(0).htest()->GetYaxis()->SetTitle(s.c_str());
  }

sortx()

void Plots::sortx ( const AxisInfo & xinfo )

inline

Definition at line 1381 of file computils.h.

                                      {
    
    if ( xinfo.rangeset() ) {       
      std::cout << "xr1: " << xinfo << std::endl;
      m_lo = xinfo.lo();
      m_hi = xinfo.hi();
    }
    
    if ( xinfo.autoset() && size() > 0 ) {
      std::vector<double> limits = findxrange( xinfo.symmetric() );
      if ( xinfo.rangeset() ) { 
    std::cout << "xr2: " << xinfo << std::endl;
    if ( limits[0]<m_lo ) m_lo = limits[0];
    if ( limits[1]>m_hi ) m_hi = limits[1];
      }
      else { 
    std::cout << "xr3: " << xinfo << std::endl;
    m_lo = limits[0];
    m_hi = limits[1];
      }
    }
    else if (size() == 0)
    {
        std::cout<<"Warning in computils.h::sortx() size=0.  Setting m_lo/m_hi to 0/1.  You will probably have empty figures."<<std::endl;
        m_lo = 0;
        m_hi = 1;
    }
 
    if ( xinfo.rangeset() || xinfo.autoset() ) { 
      SetRangeUser( m_lo, m_hi );
      if ( xinfo.log() && m_lo>0 ) SetLogx(true);
      else                         SetLogx(false);
    }
 
  }

xrange()

void Plots::xrange ( bool symmetric = false )

inline

Definition at line 1422 of file computils.h.

                                  { 
    m_rangeset = false;
    for ( unsigned i=0 ; i<size() ; i++ ) { 
      if ( at(i).href() ) ::xrange( at(i).href(), symmetric );
      ::xrange( at(i).htest(), symmetric );
    }
  }

Member Data Documentation

elements

T std::vector< T >::elements

inherited

STL member.

m_hi

double Plots::m_hi

private

Definition at line 1553 of file computils.h.

m_lo

double Plots::m_lo

private

Definition at line 1552 of file computils.h.

m_logx

bool Plots::m_logx

private

canvas log setting

Definition at line 1540 of file computils.h.

m_logy

bool Plots::m_logy

private

Definition at line 1541 of file computils.h.

m_max

double Plots::m_max

private

Definition at line 1545 of file computils.h.

m_maxset

bool Plots::m_maxset

private

yaxis range setting

Definition at line 1544 of file computils.h.

m_min

double Plots::m_min

private

Definition at line 1548 of file computils.h.

m_minset

bool Plots::m_minset

private

Definition at line 1547 of file computils.h.

m_name

std::string Plots::m_name

private

Definition at line 1537 of file computils.h.

m_rangeset

bool Plots::m_rangeset

private

xaxis range setting

Definition at line 1551 of file computils.h.

m_trim_errors

bool Plots::m_trim_errors

private

Definition at line 1555 of file computils.h.

s_watermark

bool Plots::s_watermark = true

staticprivate

Definition at line 1559 of file computils.h.

---

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

• computils.h
• computils.cxx