---

HadEcalHcalSharesEc.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Normalisers/HadEcalHcalSharesEc.h"

#include "FastShowerUtils/ParticleParameters.h"
//#include "FastShowerUtils/IDoubleShowererSelectorConfig.h"
#include "FastShowerUtils/ITripleShowererSelectorConfig.h"

#include <cmath>
#include <iostream>
namespace FastShower{
  //HadEcalHcalSharesEc::HadEcalHcalSharesEc(const IDoubleShowererSelectorConfig* c, 
  //                                         const std::string& s):
  HadEcalHcalSharesEc::HadEcalHcalSharesEc(const ITripleShowererSelectorConfig* c, 
                                           const std::string& s):
    IFnOfParticleParameters2(),HadEcalHcalSharesBase(c,s){}
  //
  double 
  HadEcalHcalSharesEc::mean(double energy) const{
    double a[2]={1.0, 0.820};
    double b[2]={-21.04, 0.420};
    double c[2]={-10.24, 0.048};

    return  (energy > 1.0)? a[0]+b[0]/std::pow((energy-c[0]),2) : 
                            a[1]*(1.-(1./(1.+std::exp((energy-b[1])/c[1])))) ;
  }
  //
  double 
  HadEcalHcalSharesEc::sigma(double energy) const{
    double alpha=0.49;
    double beta=0.08;

    return  std::sqrt( std::pow(alpha,2)/energy + std::pow(beta,2) );
  }
  //
  IFnOfParticleParameters2* HadEcalHcalSharesEc::clone() const {
    return new HadEcalHcalSharesEc(*this);
  }
}//namespace

---