---

ProcessedNormalTester.cxx

Go to the documentation of this file.

// ================================================
// Test program for ProcessedNormal
// ================================================

#include "FastShowerUtils/ProcessedNormalTester.h"
#include "FastShowerUtils/ProcessedNormal.h"
#include "FastShowerUtils/IDistProcessor.h"
#include "FastShowerUtils/LinearProcessor.h"
#include "FastShowerUtils/TestStat.h"
#include <pair.h>
#include <cmath>

namespace FastShower{
  using std::pair;
  //**********************************************************     
  TestStat ProcessedNormalTester::test(){
    m_ts+=this->test1();
    return m_ts;
  }
  //**********************************************************     
  TestStat ProcessedNormalTester::test1(){
    TestStat ts;
    std::vector<pair<double, double> > pars;
    pars.push_back( pair<double, double>(1., 0.) );
    pars.push_back( pair<double, double>(2., 0.) );
    pars.push_back( pair<double, double>(1., 1.) );
    
    std::vector<pair<double, double> >::const_iterator pptr = pars.begin();
    std::vector<pair<double, double> >::const_iterator pend = pars.end();
    
    const int maxtry = 100000;
    const double dmaxtry = static_cast<double>(maxtry);

    for(; pptr!=pend; ++pptr){
      double a = (*pptr).first;
      double b = (*pptr).second;
    
      const IDistProcessor* lp = new LinearProcessor(a, b);
      ProcessedNormal pn;

      int i        = 0;
      double sum   = 0.;
      double sumsq = 0.;
      double samp;

      for(; i!=maxtry; ++i){
        samp=pn.sample(lp);
        //      cout<< "a, b, samp "<<a<<" "<<b<<" "<<samp<<endl;
        sum+=samp;
        sumsq+=samp*samp;
      }

      double mean = sum/dmaxtry;
      double var= (sumsq/dmaxtry)-( mean*mean );
      
      //      cout<<"a,b "<<a<<" "<<b<<endl;
      //      cout<<"mean, var: "<<mean<<" "<<var<<endl;
      
      (abs(mean-b)<0.01 && abs(var-(a*a))<0.1) ? ts.incPass():ts.incFail();
    }
      return ts;
  }
  //**********************************************************     
  void ProcessedNormalTester::results() const{
    cout<<"ProcessedNormalTester: "<<m_ts;
    if(!m_ts) cout<<" Failures! ";
    cout<<endl;
  }
  
}//namespace

---