TPiKFun Class Reference

#include <TPiKFun.h>

Inheritance diagram for TPiKFun:

List of all members.

Public Member Functions
	TPiKFun ()
double	Eval_odd (const double &, const double &)
double	Eval_odd_ultra (const double &)
double	Eval_even (const double &, const double &)
double	Eval_even_odd (const double &, const double &)
void	Init (std::string, std::string)
double	Eval_int_odd (const double &, const double &)
double	Eval (const double &s, const double &c)
virtual double	Eval_even_ultra (const double &)
virtual double	Eval_Coloumb_Factor (const double &)
Static Public Member Functions
static double	Eval_a ()
Protected Attributes
float *	fm

---

Detailed Description

Definition at line 5 of file TPiKFun.h.

---

Constructor & Destructor Documentation

TPiKFun::TPiKFun

(

)

[inline]

Definition at line 7 of file TPiKFun.h.

:TMuKFun(){};

---

Member Function Documentation

double TMuKFun::Eval	(	const double &	s,
		const double &	c
	)			[inline, inherited]

Definition at line 15 of file TMuKFun.h.

References TMuKFun::Eval_a(), and TMuKFun::Eval_even_odd().

Referenced by TPiCrossPart::GetValue(), and TMuCrossPart::GetValue().

                                                      {
    return Eval_even_odd(s,c) + Eval_a();
  }

static double TMuKFun::Eval_a

(

)

[inline, static, inherited]

Definition at line 13 of file TMuKFun.h.

References gConst, and TConstants::Pi2().

Referenced by TMuKFun::Eval(), TPiCrossPart::GetValue(), and TMuCrossPart::GetValue().

{return gConst->Pi2()/6. - 1./4.;}

double TMuKFun::Eval_Coloumb_Factor

(

const double &

)

[virtual, inherited]

Definition at line 30 of file TMuKFun.C.

References TConstants::Alpha(), exp(), gConst, TRadGlobal::Get_MF2(), TRadGlobal::Get_s(), gGlobal, and TConstants::Pi().

Referenced by TMuCrossPart::GetValue(), and TMuCrossPart::TMuCrossPart().

                                                  {
  // Coloumb final state interaction, the first term of
  // expansion of Sakharov-Sommerfeld factor is taken into
  // account in formulae above, so it should be subtracted from K_Even

  double beta  = sqrt(1 - gGlobal->Get_MF2()*gGlobal->Get_s()/s);
  // relative velocity between final particles (inversed)
  double iv = (1+beta*beta)/(2*beta);
  double z = 2*gConst->Pi()*gConst->Alpha()*iv;
  // Sakharov-Sommerfeld factor
  return z/(1. - exp(-z));
}

double TPiKFun::Eval_even	(	const double &	,
		const double &
	)			[virtual]

Reimplemented from TMuKFun.

Definition at line 27 of file TPiKFun.C.

References EvtCyclic3::A, TLi2::Eval(), gConst, TRadGlobal::Get_MF2(), TRadGlobal::Get_s(), gGlobal, and TConstants::Pi2().

                                                         {
  // b-function from errata of JHEP10(1997) 006
  double r      = gGlobal->Get_MF2()*0.25*gGlobal->Get_s()/s;
  double rho    =-log(r);
  double beta   = sqrt(1-4*r);
  double beta1p = 1 + beta;
  double beta1m = 1 - beta;
  double rbeta  = beta1m/beta1p;
  double ibeta  = 1/beta;
  double A      = rho + gConst->Pi2()/6 - log(rbeta)*log(0.5*ibeta*ibeta*beta1p) + 2*TLi2::Eval(rbeta);
  double k_even = -1 + 0.5*beta1m*ibeta*rho + (2 + beta*beta)*ibeta*log(0.5*beta1p) + 0.5*(1+beta*beta)*ibeta*A;
  return k_even;
}

double TPiKFun::Eval_even_odd	(	const double &	,
		const double &
	)			[virtual]

Reimplemented from TMuKFun.

Definition at line 100 of file TPiKFun.C.

References EvtCyclic3::A, TLi2::Eval(), TMuKFun::Eval_int_odd(), gConst, TRadGlobal::Get_MF2(), TRadGlobal::Get_s(), gGlobal, and TConstants::Pi2().

                                                             {
  // b-function from errata of JHEP10(1997) 006
  double r      = gGlobal->Get_MF2()*0.25*gGlobal->Get_s()/s;
  double rho    =-log(r);
  double beta2  = 1 - 4*r;
  double beta   = sqrt(beta2);
  double beta1p = 1 + beta;
  double beta1m = 1 - beta;
  double rbeta  = beta1m/beta1p;
  double ibeta  = 1/beta;
  double lgbeta = log(0.5*beta1p);
  double A      = rho + gConst->Pi2()/6.0 - log(rbeta)*log(0.5*ibeta*ibeta*beta1p) + 2*TLi2::Eval(rbeta);
  double k_even = -1 + 0.5*beta1m*ibeta*rho + (2 + beta*beta)*ibeta*lgbeta + 0.5*(1+beta*beta)*ibeta*A;

  double betas2 = beta2*(1-c*c);
  double i2betas2 = 1/(2*betas2);

  double d      =-0.25*ibeta*beta1m*beta1m*(rho*rho + gConst->Pi2()/3)
    + ibeta*(1 + beta2)*(-rho*lgbeta - TLi2::Eval(-rbeta) + 2*TLi2::Eval(0.5*beta1m));

  double betac  = beta*c;
  double betc1m = 1 - betac;
  double l_m    = log(0.5*betc1m);
  double l_m2   = 0.5*l_m*l_m;
  double t1     = (1-beta2)/(2*betc1m);
  double a      = l_m2 - log(1 - t1)*(rho + l_m ) + TLi2::Eval(t1);
  double li2t1  = TLi2::Eval(1 - t1);
  double li2t2  = TLi2::Eval(betas2/(1 + beta2 - 2*betac));

  double res_plus = l_m2 - li2t1 + li2t2 + i2betas2*(a*(1 - beta2) + betc1m*(-2*a + d));

  //changing sign of theta cosine
  betac  = -betac;
  betc1m = 1 - betac;
  l_m    = log(0.5*betc1m);
  l_m2   = 0.5*l_m*l_m;
  t1     = (1-beta2)/(2*betc1m);
  a      = l_m2 - log(1 - t1)*(rho + l_m ) + TLi2::Eval(t1);
  li2t1  = TLi2::Eval(1 - t1);
  li2t2  = TLi2::Eval(betas2/(1 + beta2 - 2*betac));
  double res_minus= l_m2 - li2t1 + li2t2 + i2betas2*(a*(1 - beta2) + betc1m*(-2*a + d));

  return k_even + res_plus - res_minus + Eval_int_odd(beta,c);
}

double TMuKFun::Eval_even_ultra

(

const double &

)

[virtual, inherited]

Definition at line 43 of file TMuKFun.C.

References gConst, TRadGlobal::Get_MF2(), TRadGlobal::Get_s(), gGlobal, and TConstants::Pi2().

                                              {
  return 3./4.*log(s/(gGlobal->Get_MF2()*0.25*gGlobal->Get_s())) + gConst->Pi2()/6. - 1.;
}

double TMuKFun::Eval_int_odd	(	const double &	,
		const double &
	)			[inherited]

Definition at line 214 of file TMuKFun.C.

References genRecEmupikp::i, ganga-rec::j, NX, and T.

Referenced by Eval_even_odd(), TMuKFun::Eval_even_odd(), Eval_odd(), and TMuKFun::Eval_odd().

                                                               {
  int i = int(floor((beta - 0.01)*50));
  int j = int(floor((c + 1)*50));
  double wy = 50*c - (j - 50);
  double res = 0;
  if(i >= 0 && i < NX-1){
    double wx = 50*beta - i - 0.5;
    res = 
      (T(i,j  )*(1-wx) + T(i+1,j  )*wx)*(1-wy) + 
      (T(i,j+1)*(1-wx) + T(i+1,j+1)*wx)*wy;
  } else if (i == -1){
    double wx = 100*beta;
     res = 
      (T(0,j  )*(1-wy) + T(0,  j+1)*wy)*wx;
  } else if (i == NX - 1){
    double wx = 1 - 100*(beta-0.99);
    res = 
      (T(NX-1,j)*(1-wy)+ T(NX-1,j+1)*wy)*wx;
  }
  return -res;
}

double TPiKFun::Eval_odd	(	const double &	,
		const double &
	)			[virtual]

Reimplemented from TMuKFun.

Definition at line 41 of file TPiKFun.C.

References TLi2::Eval(), TMuKFun::Eval_int_odd(), gConst, TRadGlobal::Get_MF2(), TRadGlobal::Get_s(), gGlobal, and TConstants::Pi2().

                                                        {
  double r      = gGlobal->Get_MF2()*0.25*gGlobal->Get_s()/s;
  double rho    =-log(r);
  double beta2  = 1 - 4*r;
  double betas2 = beta2*(1-c*c);
  double i2betas2 = 1/(2*betas2);
  double beta   = sqrt(beta2);
  double ibeta  = 1/beta;

  double d      =-0.25*ibeta*(1-beta)*(1-beta)*(rho*rho + gConst->Pi2()/3)
    + ibeta*(1 + beta2)*(-rho*log(0.5*(1 + beta)) - 
                         TLi2::Eval(-(1-beta)/(1+beta)) + 2*TLi2::Eval(0.5*(1-beta)));

  double betac  = beta*c;
  double betc1m = 1 - betac;
  double l_m    = log(0.5*betc1m);
  double l_m2   = 0.5*l_m*l_m;
  double t1     = (1-beta2)/(2*betc1m);
  double a      = l_m2 - log(1 - t1)*(rho + l_m ) + TLi2::Eval(t1);
  double li2t1  = TLi2::Eval(1 - t1);
  double li2t2  = TLi2::Eval(betas2/(1 + beta2 - 2*betac));

  double res_plus = l_m2 - li2t1 + li2t2 + i2betas2*(a*(1 - beta2) + betc1m*(-2*a + d));

  //changing sign of theta cosine
  betac  = -betac;
  betc1m = 1 - betac;
  l_m    = log(0.5*betc1m);
  l_m2   = 0.5*l_m*l_m;
  t1     = (1-beta2)/(2*betc1m);
  a      = l_m2 - log(1 - t1)*(rho + l_m ) + TLi2::Eval(t1);
  li2t1  = TLi2::Eval(1 - t1);
  li2t2  = TLi2::Eval(betas2/(1 + beta2 - 2*betac));
  double res_minus= l_m2 - li2t1 + li2t2 + i2betas2*(a*(1 - beta2) + betc1m*(-2*a + d));

  return res_plus - res_minus + Eval_int_odd(beta,c);
}

double TPiKFun::Eval_odd_ultra

(

const double &

)

[virtual]

Reimplemented from TMuKFun.

Definition at line 79 of file TPiKFun.C.

References TLi2::Eval().

                                             {
  // Ultra-relativistic limit of odd part of K-function 
  double st_h2 = 0.5*(1-c);
  double st_h  = sqrt(st_h2);
  double ct_h2 = 0.5*(1+c);
  double ct_h  = sqrt(ct_h2);
  double lst_h = log(st_h);
  double lct_h = log(ct_h);
  /*
  std::cout<<c<<" "
           <<st_h2<<" "
           <<ct_h2<<" "
           <<st_h<<" "
           <<ct_h<<" "
           <<lst_h<<" "
           <<lct_h<<" "
           <<std::endl;
  */
  return (2-1/ct_h2)*lst_h*lst_h - (2-1/st_h2)*lct_h*lct_h - TLi2::Eval(st_h2) + TLi2::Eval(ct_h2);
}

void TMuKFun::Init	(	std::string	,
		std::string
	)			[inherited]

Definition at line 166 of file TMuKFun.C.

References false, TMuKFun::fm, fname, genRecEmupikp::i, ganga-rec::j, NX, NY, deljobs::string, and T.

Referenced by TMuCrossPart::TMuCrossPart().

                                                {
  if(!fm){
    fm = new float[NX*NY];
/*
    char *intd = getenv("INTEG_DIR");
    char *intf = getenv("INTEG_FNAME");
    std::string fname;
    if(intf == NULL){
      if(intd == NULL){
        fname = "integ.dat";
      }else{
        std::ostringstream temp;
        temp<<intd<<"/"<<"integ.dat";
        fname = temp.str();
      }
    } else {
      fname = intf;
    }
*/
    std::string fname;
    if(! intf.empty()) fname = intf;
    else fname = "integ.dat";

    if(! intd.empty()) fname = intd + "/" + fname;
    
    std::ifstream IN(fname.c_str());
    if(IN.is_open() == false ){
      std::cout<<"Can't read \""<<fname<<"\". "
               <<"Angular distribution will be partially incorrect!"
               <<std::endl;
      std::cout<<"Probably you should tune INTEG_DIR or INTEG_FNAME variables."<<std::endl;
      for(int i=0;i<NX;i++){
        for(int j=0;j<NY;j++)
          T(i,j) = 0;
      }
      return;
    }
    float t1,t2;
    for(int i=0;i<NX;i++){
      for(int j=0;j<NY-1;j++)
        IN>>t1>>t2>>T(i,j);
      T(i,0) = 0;
      T(i,NY-1) = 0;
    }
    IN.close();
  }
}

---

Member Data Documentation

float* TMuKFun::fm [protected, inherited]

Definition at line 8 of file TMuKFun.h.

Referenced by TMuKFun::Init(), TMuKFun::TMuKFun(), and TMuKFun::~TMuKFun().

---