Resplot.cxx File Reference

#include <sstream>
#include <cmath>
#include "Directory.h"
#include "tagname.h"
#include "Resplot.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TMath.h"
#include "generate.h"
#include "rmsFrac.h"

Go to the source code of this file.

Functions
void	binwidth (TH1D *h)
void	ZeroErrors (TH1D *h)
void	unZeroErrors (TH1D *h)
double	getWeights (TH1D *h)
double	getWeights (TH2D *h)
void	GetStats (TH1D *h, std::vector< double > &stats)
std::string	number (const double &x)
double	NCounter (TH1D *h)
double	FindMean (TH1D *s, double frac=0.95)
int	findmax (TH1D *s)
Double_t	langaufun (Double_t *x_par, Double_t *par)

Variables
int	nexperiments_ = 0
int	nexperiments_max = 40
int	nexperiments_min = 20

Detailed Description

Author

M.Sutton

Date

Mon Jun 21 18:35:22 BST 2004

Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration

Definition in file Resplot.cxx.

Function Documentation

binwidth()

void binwidth

(

TH1D *

h

)

Definition at line 49 of file Resplot.cxx.

                       { 
  
  for (int i=1 ; i<=h->GetNbinsX() ; i++ ) { 
    double d =  h->GetBinLowEdge(i+1) - h->GetBinLowEdge(i);
    
    h->SetBinContent(i, h->GetBinContent(i)/d ); 
    h->SetBinError(i, h->GetBinError(i)/d ); 
  }
  
}

findmax()

int findmax

(

TH1D *

s

)

find maximum bin

Definition at line 1375 of file Resplot.cxx.

                     { 
  int imax = 1;
  for ( int i=2 ; i<=s->GetNbinsX() ; i++ ) {
    if ( s->GetBinContent(i)>s->GetBinContent(imax) ) imax = i;
  } 
  return imax;
}

FindMean()

double FindMean	(	TH1D *	s,
		double	frac = 0.95 )

Definition at line 1288 of file Resplot.cxx.

                                           {
  
  //  double entries = s->GetEntries() + s->GetBinContent(0) + s->GetBinContent(s->GetNbinsX()+1);
  double entries = generate::GetEntries(s,0,s->GetNbinsX()+1);
  
  //  std::cout << s->GetName() << "\tentries: " << entries << std::endl;
  
  if ( entries==0 ) return 0;  
  
  //  s->GetXaxis()->SetRangeUser(-100, 100);
  //  s->GetXaxis()->UnZoom();
  s->GetXaxis()->SetRange(1,s->GetNbinsX());
  
  double mean  = s->GetMean();
  double meane = s->GetMeanError();
  
  int upperbin = 0;
  int lowerbin = 0;
  
  
  for ( int it=0 ; it<20 ; it++ ) { 
    
    //    std::cout << it << "\tmean " << mean << " +- " << meane << std::endl; 
    
    int imax =  s->GetXaxis()->FindBin(mean);
    
    double sumn  = s->GetBinContent(imax); 
    
    //    double uppersum = 0; 
    //    double lowersum = 0; 
    
    upperbin = imax;
    lowerbin = imax;
    
    
    int i=1;
    while ( true ) { 
      
      const int upperbin_i = imax+i;
      const int lowerbin_i = imax-i;
      
      if ( upperbin_i>s->GetNbinsX() || lowerbin_i<1 ) break;
      
      double tsum = sumn + s->GetBinContent(upperbin_i) + s->GetBinContent(lowerbin_i);
      
      if ( tsum>entries*frac ) { 
        //  uppersum = tsum/entries; 
        break; 
      }
      
      //      lowersum = tsum/entries;
      
      sumn = tsum;
      
      upperbin = upperbin_i;
      lowerbin = lowerbin_i;
      
      i++;
    }
    
    s->GetXaxis()->SetRange(lowerbin, upperbin);
    
    double m  = s->GetMean();
    double me = s->GetMeanError(); 
    
    if ( it>0 ) { 
      if ( mean==m || 
          std::fabs(mean-m)<me*1e-5 || 
          std::fabs(mean-m)<meane*1e-5 ) { 
        mean = m;
        meane = me;
        break; 
      }
    }
    
    mean = m;
    meane = me;
    
  }
  
  //  std::cout << "mean " << mean << " +- " << meane << "\tentries: " << generate::GetEntries(s) << "\tlower " << lowerbin << " - " << upperbin << std::endl; 
  
  
  return mean;
}

GetStats()

void GetStats	(	TH1D *	h,
		std::vector< double > &	stats )

Definition at line 112 of file Resplot.cxx.

                                                   {
  
  stats.clear();
  stats.resize(4);
  
  TAxis& fXaxis = *h->GetXaxis();
  
  double s   = 0;
  double sx  = 0;
  double sx2 = 0;
  double sx4 = 0;
  
  
  //  if ( fXaxis.TestBit(TAxis::kAxisRange)  ) {
  //    for (bin=0;bin<4;bin++) stats[bin] = 0;
  
  Int_t firstBinX = fXaxis.GetFirst();
  Int_t lastBinX  = fXaxis.GetLast();
  // include underflow/overflow if TH1::StatOverflows(kTRUE) in case no range is set on the axis
  //      if (fgStatOverflows && !fXaxis.TestBit(TAxis::kAxisRange)) {
  //   if (firstBinX == 1) firstBinX = 0;
  //   if (lastBinX ==  fXaxis.GetNbins() ) lastBinX += 1;
  // }
  for (int binx = firstBinX; binx <= lastBinX; binx++) {
    double x   = fXaxis.GetBinCenter(binx);
    double w   = TMath::Abs(h->GetBinContent(binx));
    //    err = TMath::Abs(GetBinError(binx));
    s += w;
    // stats[1] += err*err;
    sx  += w*x;
    sx2 += w*x*x;
    sx4 += w*x*x*x*x;
  }
  
  if (s == 0){
    std::cerr<<"Denominator of zero in GetStats()\n";
    return;
  }
  double mu = sx/s;
  
  double v = sx2/s-mu*mu;
  
  double rms       = std::sqrt(v);
  double rmserror  = std::sqrt(0.25*(sx4/v-v*s))/s;
  //  double rmserror1 = std::sqrt(0.4*(sx4/s-v*v)/s);  << is this the correct one? it seems about 2* too large 
  
  stats[0] = mu;
  stats[1] = std::sqrt(v/s);
  
  stats[2] = rms;
  stats[3] = rmserror;
  
  //  double duff = std::sqrt(0.5*v/s);
  
#if 0
  std::cout << "GetStats() "  
  //   << "\tmean " << stats[0] << " +- " << stats[1]
  << "\trms "  << stats[2] << " +- " << stats[3] << "\t(" << rmserror1 << ")" << "\t root " << duff << std::endl;
#endif
  
} 

getWeights() [1/2]

double getWeights

(

TH1D *

h

)

Definition at line 78 of file Resplot.cxx.

                           { 
  double sum   = 0;
  double sume2 = 0;
  for ( int i=0 ; i<h->GetNbinsX()+1 ; i++ ) { 
    sum   += h->GetBinContent(i);
    sume2 += h->GetBinError(i)*h->GetBinError(i);
  }
  if ( sume2!=0 ) return sum*sum/sume2;
  return 0;
}

getWeights() [2/2]

double getWeights

(

TH2D *

h

)

Definition at line 90 of file Resplot.cxx.

                           { 
  
  double sum   = 0;
  double sume2 = 0;
  
  for ( int i=1 ; i<=h->GetNbinsX() ; i++ ) { 
    for ( int j=1 ; j<=h->GetNbinsY() ; j++ ) { 
      sum   += h->GetBinContent(i,j);
      sume2 += h->GetBinError(i,j)*h->GetBinError(i,j);
    }
  }
  
  if ( sume2!=0 ) return sum*sum/sume2;
  
  return 0;
}

langaufun()

Double_t langaufun	(	Double_t *	x_par,
		Double_t *	par )

Definition at line 1997 of file Resplot.cxx.

                                                   {
  
  //Fit parameters:
  //par[0]=Total area (integral -inf to inf, normalization constant)
  //par[1]=Most Probable (MP, location) parameter of Landau density
  //par[2]=Width (scale) parameter of Landau density
  //par[3]=Width (sigma) of convoluted Gaussian function
  //
  //In the Landau distribution (represented by the CERNLIB approximation),
  //the maximum is located at x=-0.22278298 with the location parameter=0.
  //This shift is corrected within this function, so that the actual
  //maximum is identical to the MP parameter.
  
  // Numeric constants
  Double_t invsq2pi = 0.3989422804014;   // (2 pi)^(-1/2)
  Double_t mpshift  = -0.22278298;       // Landau maximum location
  
  // Control constants
  Double_t np = 500;   // number of convolution steps
  Double_t sc =   5;   // convolution extends to +-sc Gaussian sigmas
  
  double& x = x_par[0];
  
  // Variables
  Double_t xx;
  Double_t mpc;
  Double_t fland;
  Double_t sum = 0.0;
  Double_t xlow,xupp;
  Double_t step;
  Double_t i;
  
  //      if ( par[2]<0.1*par[3] ) return 0;
  
  // correct the most probable bin position as described above
  mpc = par[1] - mpshift*par[2];
  
  // Range of convolution integral
  
  //      if ( par[3]<par[2] ) np = np*par[2]/par[3];
  
  xlow = x - sc * par[3];
  xupp = x + sc * par[3];
  
  step = (xupp-xlow) / np;
  
  // Convolution integral of Landau and Gaussian by sum
  for(i=1.0; i<=np/2; i++) {
    xx = xlow + (i-0.5)*step;
    fland = TMath::Landau(xx,mpc,par[2]) / par[2];
    sum += fland * TMath::Gaus(x,xx,par[3]);
    
    xx = xupp - (i-0.5)*step;
    fland = TMath::Landau(xx,mpc,par[2]) / par[2];
    sum += fland * TMath::Gaus(x,xx,par[3]);
  }
  
  return (par[0] * step * sum * invsq2pi / par[3]);
}

NCounter()

double NCounter

(

TH1D *

h

)

Definition at line 381 of file Resplot.cxx.

                         { 
  double N=0;
  for ( int i=1 ; i<=h->GetNbinsX() ; i++ ) N += h->GetBinContent(i);
  return N; 
}

number()

std::string number

(

const double &

x

)

Definition at line 374 of file Resplot.cxx.

                                  { 
  std::ostringstream s;
  s << x;
  return s.str();
}

unZeroErrors()

void unZeroErrors

(

TH1D *

h

)

Definition at line 72 of file Resplot.cxx.

                           {
  for (int i=1 ; i<=h->GetNbinsX() ; i++ ) if (h->GetBinContent(i)==0) h->SetBinError(i,0);
}

ZeroErrors()

void ZeroErrors

(

TH1D *

h

)

Definition at line 66 of file Resplot.cxx.

                         {
  for (int i=1 ; i<=h->GetNbinsX() ; i++ ) if (h->GetBinContent(i)==0) h->SetBinError(i,1);
}

Variable Documentation

nexperiments_

int nexperiments_ = 0

Definition at line 1660 of file Resplot.cxx.

nexperiments_max

int nexperiments_max = 40

Definition at line 1661 of file Resplot.cxx.

nexperiments_min

int nexperiments_min = 20

Definition at line 1662 of file Resplot.cxx.