KascadeGrandeLDF.cc

Go to the documentation of this file.

#include "KascadeGrandeLDF.h"
//
#include <TMath.h>

namespace MdLDFFinderAG {

  KascadeGrandeLDF::KascadeGrandeLDF(const double r) 
  {
    fName = "KascadeGrandeLDF";
    fReferenceDist = r;
  }
  
  double
  KascadeGrandeLDF::operator()(const double r, const double *par )
    const
  {
    //const double alpha = 1;
    //const double gamma = 1.85;
    //const double r0 = 150 ;
    const double alpha = 0.3;
    const double gamma = 4.6;
    const double r0 = 280;

    double rho  = par[0];
    double beta = par[1];

    double x, xx;

    //normalization factor such that f(r=fReferenceDist) = 1
    x = fReferenceDist/r0;
    xx = x*0.1*x*0.1;//(r/10r0)^2
    double k = TMath::Power( x, alpha) * TMath::Power( 1+x, beta ) * TMath::Power( 1+xx, gamma );
    
    x = r/r0;
    xx = x*0.1*x*0.1;//(r/10r0)^2

    return  rho * k * TMath::Power( x, -alpha) * TMath::Power( 1+x, -beta ) * TMath::Power( 1+xx, -gamma );

  }

}