/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/TrackingTools/GsfTools/interface/GSUtilities.h

Go to the documentation of this file.

#ifndef GSUtilities_h
#define GSUtilities_h

class GSUtilities {
public:
  GSUtilities (const unsigned nComp, const float* weights,
               const float* parameters, const float* errors) :
    theNComp(nComp),
    theWeights(0),
    theParameters(0),
    theErrors(0)
  {
    if ( theNComp ) {
      theWeights = new float[theNComp];
      theParameters = new float[theNComp];
      theErrors = new float[theNComp];
    }
    const float* wPtr1(weights);
    const float* pPtr1(parameters);
    const float* ePtr1(errors);
    float* wPtr2(theWeights);
    float* pPtr2(theParameters);
    float* ePtr2(theErrors);
    for ( unsigned i=0; i<theNComp; i++ ) {
      *(wPtr2++) = weights ? *(wPtr1++) : 1.;
      *(pPtr2++) = *(pPtr1++);
      *(ePtr2++) = *(ePtr1++);
    }
  } 
  ~GSUtilities () 
  {
    delete [] theWeights;
    delete [] theParameters;
    delete [] theErrors;
  }
  float quantile (const float) const;
  float mode () const;
  double pdf (const double&) const;
  double cdf (const double&) const;
  double dpdf1 (const double&) const;
  double dpdf2 (const double&) const;

  double combinedMean() const;
  //  mean value of errors 
  double errorCombinedMean() const;
  //  error for the highest weight
  float errorHighestWeight() const;
  // max weight component - chiara
  float maxWeight() const;
  // mode error + some utilities functions
  float errorMode();
  float getMax(float);
  float getMin(float);

private:
  double gauss (const double&, const double&, const double&) const;
  double gaussInt (const double&, const double&, const double&) const;
  double findMode (const double) const;

private:
  unsigned theNComp;
  float* theWeights;
  float* theParameters;
  float* theErrors;
};
#endif