---

ProcessedFlatTester.cxx

Go to the documentation of this file.

// ================================================
// Test program for ProcessedFlat
// ================================================

#include "FastShowerUtils/ProcessedFlatTester.h"
#include "FastShowerUtils/ProcessedFlat.h"
#include "FastShowerUtils/IDistProcessor.h"
#include "FastShowerUtils/LinearProcessor.h"
#include "FastShowerUtils/TestStat.h"

namespace FastShower{
  //**********************************************************     
  TestStat ProcessedFlatTester::test(){
    m_ts+=this->test1();
    return m_ts;
  }
  //**********************************************************     
  TestStat ProcessedFlatTester::test1(){
    TestStat ts;
    const IDistProcessor* lp1 = new LinearProcessor(1., 0.);
    const IDistProcessor* lp2 = new LinearProcessor(2., 0.);
    const IDistProcessor* lp3 = new LinearProcessor(1., 1.);
    ProcessedFlat pf;

    bool pass = true;
    const int maxtry = 1000;

    int i = maxtry;
    while(i!=0 && pass){
      --i;
      double samp=pf.sample(lp1);
      //      cout<<"[0,1) "<<samp<<endl;
      if (samp <0. || samp>=1.) pass = false;
    }
    pass ? ts.incPass():ts.incFail();


    pass = true;
    i = maxtry;

    while(i!=0 && pass){
      --i;
      double samp=pf.sample(lp2);
      //      cout<<"[0,2) "<<samp<<endl;
      if (samp <0. || samp>=2.) pass = false;
    }
    pass ? ts.incPass():ts.incFail();


    pass = true;
    i = maxtry;

    while(i!=0 && pass){
      --i;
      double samp=pf.sample(lp3);
      //      cout<<"[1,2) "<<samp<<endl;
      if (samp <1. || samp>=2.) pass = false;
    }

    pass ? ts.incPass():ts.incFail();

    return ts;
  }
  //**********************************************************     
  void ProcessedFlatTester::results() const{
    cout<<"ProcessedFlatTester: "<<m_ts;
    if(!m_ts) cout<<" Failures! ";
    cout<<endl;
  }
    
  
}//namespace

---