TMom Class Reference

#include <TMom.h>

Inheritance diagram for TMom:

Public Member Functions
void	addEntry (double val)
void	addEntry (double val, const std::vector< double > &valcut)
double	getM3 ()
double	getMax ()
double	getMean ()
double	getMean2 ()
double	getMean3 ()
double	getMin ()
int	getNevt ()
std::vector< double >	getPeak ()
double	getRMS ()
void	setCut (const std::vector< double > &, const std::vector< double > &)
void	setCut (double, double)
	TMom ()
	TMom (const std::vector< double > &, const std::vector< double > &)
	TMom (double, double)
	~TMom () override

Public Attributes
int	_dimCut

Private Member Functions
void	init (const std::vector< double > &, const std::vector< double > &)
void	init (double, double)

Private Attributes
std::vector< double >	_ampl
std::vector< double >	_cutHigh
std::vector< double >	_cutLow
int	bing [101]
double	M3
double	max
double	mean
double	mean2
double	mean3
double	min
int	nevt
double	peak
double	rms
double	sum
double	sum2
double	sum3

Detailed Description

Definition at line 7 of file TMom.h.

Constructor & Destructor Documentation

TMom() [1/3]

TMom::TMom

(

)

Definition at line 18 of file TMom.cc.

{ init(0.0, 10.0e9); }

References init.

TMom() [2/3]

TMom::TMom	(	double	cutlow,
		double	cuthigh
	)

Definition at line 21 of file TMom.cc.

{ init(cutlow, cuthigh); }

References init.

TMom() [3/3]

TMom::TMom	(	const std::vector< double > &	cutlow,
		const std::vector< double > &	cuthigh
	)

Definition at line 22 of file TMom.cc.

{ init(cutlow, cuthigh); }

References init.

~TMom()

TMom::~TMom ( )

override

Definition at line 25 of file TMom.cc.

{}

Member Function Documentation

addEntry() [1/2]

void TMom::addEntry

(

double

val

)

Definition at line 88 of file TMom.cc.

                              {
  std::vector<double> dumb;
  dumb.push_back(val);
  addEntry(val, dumb);
}

References heppy_batch::val.

Referenced by EcalABAnalyzer::analyze(), EcalLaserAnalyzer2::analyze(), EcalLaserAnalyzer::analyze(), and EcalTestPulseAnalyzer::endJob().

addEntry() [2/2]

void TMom::addEntry	(	double	val,
		const std::vector< double > &	valcut
	)

Definition at line 94 of file TMom.cc.

                                                               {
  int passingAllCuts = 1;
 
  for (int iCut = 0; iCut < _dimCut; iCut++) {
    int passing;
    if (valcut.at(iCut) > _cutLow.at(iCut) && valcut.at(iCut) <= _cutHigh.at(iCut)) {
      passing = 1;
    } else
      passing = 0;
    passingAllCuts *= passing;
  }
 
  if (passingAllCuts == 1) {
    nevt += 1;
    sum += val;
    sum2 += val * val;
    sum3 += val * val * val;
    if (val > max)
      max = val;
    if (val < min)
      min = val;
 
    // for peak stuff
    _ampl.push_back(val);
  }
}

References SiStripPI::max, min(), nevt, TagProbeFitTreeAnalyzer_cfi::passing, combinedConstraintHelpers::sum2(), and heppy_batch::val.

getM3()

double TMom::getM3

(

)

Definition at line 156 of file TMom.cc.

                   {
  double m = getMean();
  double m2 = getMean2();
  double m3 = getMean3();
  double sig = getRMS();
 
  if (nevt != 0 && sig != 0)
    M3 = (m3 - 3.0 * m * (m2 - m * m) - m * m * m) / (sig * sig * sig);
  else
    M3 = 0.0;
  return M3;
}

References visualization-live-secondInstance_cfg::m, callgraph::m2, and nevt.

Referenced by EcalTestPulseAnalyzer::endJob().

getMax()

double TMom::getMax

(

)

Definition at line 170 of file TMom.cc.

{ return max; }

References SiStripPI::max.

Referenced by EcalTestPulseAnalyzer::endJob().

getMean()

double TMom::getMean

(

)

Definition at line 121 of file TMom.cc.

                     {
  if (nevt != 0)
    mean = sum / nevt;
  else
    mean = 0.0;
  return mean;
}

References SiStripPI::mean, and nevt.

Referenced by EcalABAnalyzer::endJob(), EcalTestPulseAnalyzer::endJob(), EcalLaserAnalyzer2::endJob(), and EcalLaserAnalyzer::endJob().

getMean2()

double TMom::getMean2

(

)

Definition at line 129 of file TMom.cc.

                      {
  if (nevt != 0)
    mean2 = sum2 / nevt;
  else
    mean2 = 0.0;
  return mean2;
}

References nevt, and combinedConstraintHelpers::sum2().

getMean3()

double TMom::getMean3

(

)

Definition at line 136 of file TMom.cc.

                      {
  if (nevt != 0)
    mean3 = sum3 / nevt;
  else
    mean3 = 0.0;
  return mean3;
}

References nevt.

getMin()

double TMom::getMin

(

)

Definition at line 169 of file TMom.cc.

{ return min; }

References min().

Referenced by EcalTestPulseAnalyzer::endJob().

getNevt()

int TMom::getNevt

(

)

Definition at line 144 of file TMom.cc.

{ return nevt; }

References nevt.

Referenced by EcalTestPulseAnalyzer::endJob().

getPeak()

std::vector< double > TMom::getPeak

(

)

Definition at line 172 of file TMom.cc.

                                {
  std::vector<double> p;
  double wbin = (max - min) / 100.;
  int bung;
 
  for (unsigned int i = 0; i < _ampl.size(); i++) {
    if (wbin <= 0.0)
      bung = 1;
    else
      bung = (int)((_ampl.at(i) - min) / wbin) + 1;
    if (1 <= bung && bung <= 100)
      bing[bung]++;
  }
 
  TMarkov* peakM = new TMarkov();
 
  int popmax = 0;
 
  for (int k = 1; k < 101; k++) {
    if (bing[k] > popmax) {
      popmax = bing[k];
    }
  }
 
  peakM->peakFinder(&bing[0]);
  p.push_back(peakM->getPeakValue(0));
  p.push_back(peakM->getPeakValue(1));
 
  return p;
}

References TMarkov::getPeakValue(), mps_fire::i, createfilelist::int, dqmdumpme::k, SiStripPI::max, min(), AlCaHLTBitMon_ParallelJobs::p, and TMarkov::peakFinder().

getRMS()

double TMom::getRMS

(

)

Definition at line 146 of file TMom.cc.

                    {
  double m = getMean();
  double m2 = getMean2();
  if (nevt != 0)
    rms = TMath::Sqrt(m2 - m * m);
  else
    rms = 0.0;
  return rms;
}

References visualization-live-secondInstance_cfg::m, callgraph::m2, nevt, and SiStripPI::rms.

Referenced by EcalTestPulseAnalyzer::endJob().

init() [1/2]

void TMom::init ( const std::vector< double > &  cutlow, const std::vector< double > &  cuthigh  )

private

Definition at line 49 of file TMom.cc.

                                                                                 {
  nevt = 0;
  mean = 0;
  mean2 = 0;
  mean3 = 0;
  sum = 0;
  sum2 = 0;
  sum3 = 0;
  rms = 0;
  M3 = 0;
  peak = 0;
  min = 10.0e9;
  max = 0.;
  assert(cutlow.size() == cuthigh.size());
  _cutLow.clear();
  _cutHigh.clear();
  _dimCut = cutlow.size();
  _cutLow = cutlow;
  _cutHigh = cuthigh;
  for (int i = 0; i < 101; i++) {
    bing[i] = 0;
  }
}

References cms::cuda::assert(), mps_fire::i, SiStripPI::max, SiStripPI::mean, min(), nevt, SiStripPI::rms, and combinedConstraintHelpers::sum2().

init() [2/2]

void TMom::init ( double  cutlow, double  cuthigh  )

private

Definition at line 27 of file TMom.cc.

                                             {
  nevt = 0;
  mean = 0;
  mean2 = 0;
  mean3 = 0;
  sum = 0;
  sum2 = 0;
  sum3 = 0;
  rms = 0;
  M3 = 0;
  peak = 0;
  min = 10.0e9;
  max = 0.;
  _cutLow.clear();
  _cutHigh.clear();
  _dimCut = 1;
  _cutLow.push_back(cutlow);
  _cutHigh.push_back(cuthigh);
  for (int i = 0; i < 101; i++) {
    bing[i] = 0;
  }
}

References mps_fire::i, SiStripPI::max, SiStripPI::mean, min(), nevt, SiStripPI::rms, and combinedConstraintHelpers::sum2().

setCut() [1/2]

void TMom::setCut	(	const std::vector< double > &	cutlow,
		const std::vector< double > &	cuthigh
	)

Definition at line 79 of file TMom.cc.

                                                                                   {
  assert(cutlow.size() == cuthigh.size());
  _cutLow.clear();
  _cutHigh.clear();
  _dimCut = cutlow.size();
  _cutLow = cutlow;
  _cutHigh = cuthigh;
}

References cms::cuda::assert().

setCut() [2/2]

void TMom::setCut	(	double	cutlow,
		double	cuthigh
	)

Definition at line 72 of file TMom.cc.

                                               {
  _cutLow.clear();
  _cutHigh.clear();
  _dimCut = 1;
  _cutLow.push_back(cutlow);
  _cutHigh.push_back(cuthigh);
}

Member Data Documentation

_ampl

std::vector<double> TMom::_ampl

private

Definition at line 24 of file TMom.h.

_cutHigh

std::vector<double> TMom::_cutHigh

private

Definition at line 23 of file TMom.h.

_cutLow

std::vector<double> TMom::_cutLow

private

Definition at line 22 of file TMom.h.

_dimCut

int TMom::_dimCut

Definition at line 30 of file TMom.h.

bing

int TMom::bing[101]

private

Definition at line 21 of file TMom.h.

M3

double TMom::M3

private

Definition at line 17 of file TMom.h.

max

double TMom::max

private

Definition at line 20 of file TMom.h.

mean

double TMom::mean

private

Definition at line 10 of file TMom.h.

mean2

double TMom::mean2

private

Definition at line 11 of file TMom.h.

mean3

double TMom::mean3

private

Definition at line 12 of file TMom.h.

min

double TMom::min

private

Definition at line 19 of file TMom.h.

nevt

int TMom::nevt

private

Definition at line 9 of file TMom.h.

peak

double TMom::peak

private

Definition at line 18 of file TMom.h.

rms

double TMom::rms

private

Definition at line 16 of file TMom.h.

sum

double TMom::sum

private

Definition at line 13 of file TMom.h.

sum2

double TMom::sum2

private

Definition at line 14 of file TMom.h.

sum3

double TMom::sum3

private

Definition at line 15 of file TMom.h.