BOSS_7.0.2: /home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/Mcgpj/Mcgpj-00-01-04/src/code/src/TPhoton.C Source File

00001 #include <iostream>
00002 #include "TPhoton.h"
00003 #define DB(x) #x<<"="<<x<<", "
00004 using namespace std;
00005 void TPhotonI::Init(){
00006   // Initializing photons in Initial state
00007   // by distribution 1/x/(1-x)
00008   fP0 = gGlobal->Get_XMin()/(1-gGlobal->Get_XMin());
00009   fP1 = gGlobal->Get_XMax()/(1-gGlobal->Get_XMax())/fP0;
00010   fNorm = log(fP1);
00011   fP1 = fNorm; // to increase speed of exponent in GetEnergy
00012   fBetaI = gGlobal->Get_BetaI();
00013   fLnD   = gGlobal->Get_LnD();
00014   //  cout<<DB(fP0)<<DB(fP1)<<DB(fNorm)<<DB(fBetaI)<<DB(fLnD)<<endl<<flush;
00015 }
00016 
00017 void TPhotonF::Init(){
00018   // Initializing photons in final state 
00019   // by distribution 1/x
00020   fP0 = gGlobal->Get_XMin();
00021   fP1 = gGlobal->Get_XMax()/fP0;
00022   fNorm = log(fP1);
00023   fP1 = fNorm; // to increase speed of exponent in GetEnergy
00024   fBetaI = gGlobal->Get_BetaI();
00025   fLnD   = gGlobal->Get_LnD();
00026 
00027   fiBetaI = 1/fBetaI;
00028   double c1 = gGlobal->Get_CosTheta0();
00029   double c2 = -1;// + 0.5*(1-c1);
00030   fBt    = (1 + fBetaI*c1)/(1 - fBetaI*c1);
00031   fLnD_n = log(1./fBt*(1. + fBetaI*c2)/(1. - fBetaI*c2));
00032 
00033   fBt2    = 1 - fBetaI*c1;
00034   fLnD_n2 = log((1. - fBetaI*c2)/(1 - fBetaI*c1));
00035   //  cout<<DB(fBt)<<DB(fLnD_n)<<DB(fLnD)<<endl<<flush;
00036   //  exit(0);
00037 }
00038 
00039 void TPhotonS::Init(){
00040   // Initializing photons in final state 
00041   // by distribution 1/x
00042   Init(gGlobal->Get_XMin(), gGlobal->Get_Beta());
00043   fBetaI = gGlobal->Get_BetaI();
00044   fLnD   = gGlobal->Get_LnD();
00045 }
00046 
00047 void TPhotonS::Init(const double & delta, const double & beta){
00048   // Initializing photons in final state 
00049   // by distribution 1/x
00050   fP0 = delta;
00051   fP1 = 2/beta;
00052   fP2 = pow(2/beta*pow(delta,-beta/2),4);
00053   fP3 = 1.-beta/2;
00054   fNorm = pow(pow(delta,beta/2),4);
00055   fBetaI = gGlobal->Get_BetaI();
00056   fLnD   = gGlobal->Get_LnD();
00057 }
00058 
00059 void TPhotonD::Init(){
00060   // Initializing photons in final state 
00061   // by distribution 1/x
00062   Init(gGlobal->Get_XMin(), gGlobal->Get_Beta());
00063   fBetaI = gGlobal->Get_BetaI();
00064   fLnD   = gGlobal->Get_LnD();
00065 }
00066 
00067 void TPhotonD::Init(const double & delta, const double & beta){
00068   // Initializing photons in final state 
00069   // by distribution 1/x
00070   fP0 = delta;
00071   fP1 = 2/beta;
00072   fNorm = pow(delta,beta/2);
00073   fBetaI = gGlobal->Get_BetaI();
00074   fLnD   = gGlobal->Get_LnD();
00075 }