TEvent.h

Go to the documentation of this file.

#ifndef RAD_TEvent
#define RAD_TEvent
 
#include "TConstants.h"
#include "TPhoton.h"
#include "TRandom.h"
#include "TLorentzVector.h"
#include <iostream>
#define Ne 0
#define Np 1
#define Ng1 2
 
class TEvent{
protected:
  double fen[3];  // energy of final particle 
  double fp[3];   // momentum of final particle 
  double fY[3];   // energy of particles before interaction for cross
          // section relativistic shift
 
  double fth[3];  // theta of final particle 
  double fcth[3]; // cos(theta) of final particle 
  double fsth[3]; // sin(theta) of final particle 
 
  double fph[3];  // phi of final particle 
  double fcph[3]; // cos(phi) of final particle 
  double fsph[3]; // sin(phi) of final particle 
  bool fIsSmear;  // to smear or not to smear that is a question
  double fdPar[9]; // smearing parameters for delta theta and delta phi
  double fxg[4];
  double fCosPsi;
  unsigned int fNSelStat;
  unsigned int *fSelStat; //[fNSelStat]
public:
  TEvent();
  virtual ~TEvent();
  void Init();
  // two photons radiation on arbitrary angles
  bool MakeEvent(const double &, const double &, const double &, 
         const double &, const double &, const double &, 
         const double &, const double &);
  // two photons radiation in collinear region for initial particles with masses
  bool MakeEvent(const double &, const double &, const double &);
  bool MakeEventgg(const double &, const double &, const double &);
  // two photons radiation in collinear region for final particles in massless approximation 
  bool MakeEventN(const double &, const double &);
  // four photons radiation in collinear region for initial&final particles in massless approximation 
  bool MakeEvent(const double &, const double &, const double &,
         const double &, const double &);
  // Collinear event
  bool MakeEvent(const double &);
  // One hard photon on big angle
  bool MakeEvent(const double &, TPhoton*, const unsigned int);
  // Kinematics cut
  bool Select();
  bool Selectgg();
  inline void GetEvent(double &p0, double &t0, double &ph0, 
               double &p1, double &t1, double &ph1){
    p0 = fp[Ne]; t0 = fth[Ne]; ph0 = fph[Ne];
    p1 = fp[Np]; t1 = fth[Np]; ph1 = fph[Np];
  }
  inline double GetPhi(){return gRandom->Rndm()*2*gConst->Pi();}
  inline double GetY(unsigned int i){return fY[i];}
  void Print();
  void PrintSelectStat(){
    for(unsigned int i=0; i<fNSelStat; i++)std::cout<<"Cut #"<<i<<" "<<fSelStat[i]<<std::endl;
  }
  inline double GetCosPsi(){return fCosPsi;};
  inline void GetPPar(unsigned int i, double &e, double &p, double &c){
    e = fen[i]; p = fp[i]; c = fcth[i];
  }
  inline TLorentzVector Get4Vector(unsigned int i){
    return TLorentzVector(fp[i]*fsth[i]*fcph[i],fp[i]*fsth[i]*fsph[i],fp[i]*fcth[i],fen[i]);
  }
  inline double GetChiPrime(unsigned int i){
    return fen[i]*fen[Ng1] - fp[i]*fp[Ng1]*(fsth[i]*fsth[Ng1]*(fcph[i]*fcph[Ng1]+fsph[i]*fsph[Ng1])+fcth[i]*fcth[Ng1]);
  }
  inline double GetSumComponent(unsigned int i){
    switch(i){
    case 0:
      return fen[Ne] + fen[Np] + fen[Ng1];
      break;
    case 1:
      return fp[Ne]*fsth[Ne]*fcph[Ne] + fp[Np]*fsth[Np]*fcph[Np] + fp[Ng1]*fsth[Ng1]*fcph[Ng1];
      break;
    case 2:
      return fp[Ne]*fsth[Ne]*fsph[Ne] + fp[Np]*fsth[Np]*fsph[Np] + fp[Ng1]*fsth[Ng1]*fsph[Ng1];
      break;
    case 3:
      return fp[Ne]*fcth[Ne] + fp[Np]*fcth[Np] + fp[Ng1]*fcth[Ng1];
      break;
    }
    return 0;
  }
  void Swap();
  void CosPrint();
  void SetdPar(double *p){
    fIsSmear = true;
    for(int i=0; i<9; i++) fdPar[i] = p[i];
  }
 
  inline double dDeltaTheta(){
    if ( !fIsSmear ) return 0;
    double sigma = 0;
    if(gRandom->Rndm()<fdPar[0])
      sigma = fdPar[1];
    else
      sigma = fdPar[2];
    return gRandom->Gaus(0,sigma);
  }
  inline double dDeltaPhi(){
    if ( !fIsSmear ) return 0;
    double sigma = 0;
    if(gRandom->Rndm()<fdPar[3])
      sigma = fdPar[4];
    else
      sigma = fdPar[5];
    return gRandom->Gaus(0,sigma);
  }
  inline double dATheta(){
    if ( !fIsSmear ) return 0;
    double sigma = 0;
    if(gRandom->Rndm()<fdPar[6])
      sigma = fdPar[7];
    else
      sigma = fdPar[8];
    return gRandom->Gaus(0,sigma);
  }
  void GetEvent(double*, int&);
};
 
#endif //#ifndef RAD_TEvent