TripleGaussCollFit.cxx File Reference

#include "dqm_algorithms/TripleGaussCollFit.h"
#include <dqm_core/AlgorithmManager.h>
#include <dqm_core/exceptions.h>
#include <dqm_algorithms/tools/AlgorithmHelper.h>
#include "ers/ers.h"
#include "TH2.h"
#include "TH1.h"
#include "TF1.h"
#include "TClass.h"
#include <iostream>
#include <cmath>

Go to the source code of this file.

Namespaces
namespace	dqm_algorithms

Macros
#define	PARAM_IS_AREA
#define	ADJUST_VAL   startVal = startVal * hist1D->GetBinWidth(1) / m_gaus3_fn->GetParameter(2) / 2.5066

Typedefs
typedef std::map< std::string, double >	param_map
typedef std::map< std::string, double >	err_map

Functions
Double_t	getStartValue (TH1 *thehist, Double_t minx, Double_t maxx)
void	multLB_writeTags (std::vector< int > *goodLB, dqm_core::Result *result, const TH1 *thehist)

Variables
const char *const	fmt = "%0*d"

Macro Definition Documentation

ADJUST_VAL

#define ADJUST_VAL   startVal = startVal * hist1D->GetBinWidth(1) / m_gaus3_fn->GetParameter(2) / 2.5066

Definition at line 380 of file TripleGaussCollFit.cxx.

PARAM_IS_AREA

#define PARAM_IS_AREA

Definition at line 165 of file TripleGaussCollFit.cxx.

Typedef Documentation

err_map

typedef std::map<std::string, double> err_map

Definition at line 163 of file TripleGaussCollFit.cxx.

param_map

typedef std::map<std::string, double> param_map

Definition at line 162 of file TripleGaussCollFit.cxx.

Function Documentation

getStartValue()

Double_t getStartValue	(	TH1 *	thehist,
		Double_t	minx,
		Double_t	maxx )

Definition at line 70 of file TripleGaussCollFit.cxx.

                                                                   {
 
  Double_t result = 0.0;
  int minidx = thehist->FindBin(minx), maxidx = thehist->FindBin(maxx);
  for (int idx = minidx; idx <= maxidx; idx++) {
    Double_t val = thehist->GetBinContent(idx);
    if (val > result)
      result = val;
  }
 
  return result;
}

multLB_writeTags()

void multLB_writeTags	(	std::vector< int > *	goodLB,
		dqm_core::Result *	result,
		const TH1 *	thehist )

Definition at line 86 of file TripleGaussCollFit.cxx.

                                                                                          {
 
  /* This method assumes that the entries in goodLB are in ascending order!! */
 
  char buf[10]; 
  std::string lb_str = "LB_";
  int firstLB = (*goodLB)[0], offset = 1; 
  unsigned int c = 1, sz = goodLB->size(), psz = static_cast<unsigned int>( std::log10(sz))+1;
  Double_t lo, hi;
    
  while (c < sz) {
    if ((*goodLB)[c] == firstLB + offset) 
      offset++;
      
    else {
      std::string tagname = lb_str;
      //find the low LB
      
      lo = thehist->GetXaxis()->GetBinLowEdge(firstLB);
 
      if (lo == std::floor(lo))
    sprintf(buf, fmt, psz, static_cast<Int_t>(lo));
      else 
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(lo))+1);
 
      tagname +=  buf;
      tagname += "_";
 
      hi = thehist->GetXaxis()->GetBinUpEdge(firstLB+offset-1);
 
      if (hi == std::floor(hi))
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(hi))-1);
      else
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(hi)));
 
      tagname += buf;
      result->tags_ [tagname] = 1.0;
      firstLB = (*goodLB)[c];
      offset = 1;
    }
 
    c++;  //I've always wanted to write this
  }
 
  lo = thehist->GetXaxis()->GetBinLowEdge(firstLB);
 
  if (lo == std::floor(lo))
    sprintf(buf, fmt, psz, static_cast<Int_t>(lo));
  else 
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(lo))+1);
 
  lb_str += buf;
  lb_str += "_";
  hi = thehist->GetXaxis()->GetBinUpEdge(firstLB+offset-1);
  if (hi == std::floor(hi))
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(hi))-1);
  else
    sprintf(buf, fmt, psz, static_cast<Int_t>(std::floor(hi)));
  lb_str += buf;
  result->tags_ [lb_str] = 1.0;
 
}

Variable Documentation

fmt

const char* const fmt = "%0*d"

Definition at line 84 of file TripleGaussCollFit.cxx.