TEvent.h

Go to the documentation of this file.

#ifndef RAD_TEvent
#define RAD_TEvent
 
#include "TConstants.h"
#include "TLorentzVector.h"
#include "TPhoton.h"
#include "TRandom.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