---

ProcessedNormal.cxx

Go to the documentation of this file.

#include "FastShowerUtils/ProcessedNormal.h"
#include "FastShowerUtils/IDistProcessor.h"
#include "CLHEP/Random/RandGauss.h"
#include <assert.h>
#include <iostream>
namespace FastShower{
  double ProcessedNormal::sample() const{
    return RandGauss::shoot();
  }
  double ProcessedNormal::sample(const IDistProcessor* dp) const{
    return dp->process(RandGauss::shoot());
  }
  double ProcessedNormal::sample(double lowLim, 
                                 double upLim, 
                                 const IDistProcessor* dp) const{

    //cout<<"*** "<<text()<<": low/upLim = "<<lowLim<<" / "<<upLim<<endl;
    assert(lowLim<upLim);
    
    double value = dp->process( RandGauss::shoot() );

    //keep going until within truncation limits;
    int ncalls = 0;
    while(value < lowLim || value > upLim){
      value = dp->process( RandGauss::shoot() );
      ++ncalls;
      if( ncalls%2 == 1000){
        cout<<"Processed Normal called "<< ncalls<<" times"<<endl;
        cout<<"Processed Normal lowLim "
            << lowLim<<" upLim"<<upLim<<" value"<<value<<endl;
      }
    }
    
    return value;
  }
  IProcessedDist* ProcessedNormal::clone() const {
    IProcessedDist* iptr = new ProcessedNormal(*this);
    return iptr;
  } 

}

---