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 1204 of file computils.h.

Constructor & Destructor Documentation

Plots()

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

inline

Definition at line 1208 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 1446 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 1455 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 1310 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 1500 of file computils.h.

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

GetYaxisTitle()

std::string Plots::GetYaxisTitle ( )

inline

Definition at line 1505 of file computils.h.

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

hi()

double Plots::hi ( ) const

inline

Definition at line 1413 of file computils.h.

{ return m_hi; }

limits()

void Plots::limits ( )

inline

Definition at line 1436 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 1412 of file computils.h.

{ return m_lo; }

Max()

void Plots::Max ( double scale )

inline

Definition at line 1302 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 1243 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 1289 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 1262 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 1521 of file computils.h.

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

realmax()

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

inline

Definition at line 1229 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 1218 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 1497 of file computils.h.

{ m_logx=b; } 

SetLogy()

void Plots::SetLogy ( bool b = true )

inline

Definition at line 1498 of file computils.h.

{ m_logy=b; } 

SetRangeUser()

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

inline

Definition at line 1424 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 1527 of file computils.h.

{ s_watermark = b; }

SetXaxisTitle()

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

inline

Definition at line 1510 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 1514 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 1375 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 1416 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 1547 of file computils.h.

m_lo

double Plots::m_lo

private

Definition at line 1546 of file computils.h.

m_logx

bool Plots::m_logx

private

canvas log setting

Definition at line 1534 of file computils.h.

m_logy

bool Plots::m_logy

private

Definition at line 1535 of file computils.h.

m_max

double Plots::m_max

private

Definition at line 1539 of file computils.h.

m_maxset

bool Plots::m_maxset

private

yaxis range setting

Definition at line 1538 of file computils.h.

m_min

double Plots::m_min

private

Definition at line 1542 of file computils.h.

m_minset

bool Plots::m_minset

private

Definition at line 1541 of file computils.h.

m_name

std::string Plots::m_name

private

Definition at line 1531 of file computils.h.

m_rangeset

bool Plots::m_rangeset

private

xaxis range setting

Definition at line 1545 of file computils.h.

m_trim_errors

bool Plots::m_trim_errors

private

Definition at line 1549 of file computils.h.

s_watermark

bool Plots::s_watermark = true

staticprivate

Definition at line 1553 of file computils.h.

---

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

• computils.h
• computils.cxx