---

HadEcalHcalSharesBar.cxx

Go to the documentation of this file.

#include "FastShowerUtils/Normalisers/HadEcalHcalSharesBar.h"

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

#include <cmath>
#include <iostream>
namespace FastShower{
  //HadEcalHcalSharesBar::HadEcalHcalSharesBar(const IDoubleShowererSelectorConfig* c, 
  //                                           const std::string& s):
  HadEcalHcalSharesBar::HadEcalHcalSharesBar(const ITripleShowererSelectorConfig* c, 
                                             const std::string& s):
    IFnOfParticleParameters2(),HadEcalHcalSharesBase(c,s){}
  //
  double 
  HadEcalHcalSharesBar::mean(double energy) const{
    double a[2]={1.0, 0.820};
    double b[2]={-16.779, 0.420};
    double c[2]={-8.374, 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 
  HadEcalHcalSharesBar::sigma(double energy) const{
    double alpha=0.46;
    double beta=0.08;

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

---