#include <iostream>
#include <cmath>
#include "TEvent.h"
#include "TRadGlobal.h"
#include "TKinemCut.h"
Go to the source code of this file.
Defines | |
#define | INRANGE(x) if ( x < 0 ) x += 2*gConst->Pi(); else if( x > 2*gConst->Pi() ) x -= 2*gConst->Pi() |
#define | swap(x) {double t = x[Ne]; x[Ne] = x[Np]; x[Np] = t;} |
Functions | |
void | lorbck (double *pcm, double *pincm, double *pinlab) |
#define INRANGE | ( | x | ) | if ( x < 0 ) x += 2*gConst->Pi(); else if( x > 2*gConst->Pi() ) x -= 2*gConst->Pi() |
Definition at line 11 of file TEvent.C.
Referenced by TEvent::MakeEvent(), and TEvent::MakeEventgg().
Definition at line 593 of file TEvent.C.
Referenced by MdcHoughFinder::execute(), HoughValidUpdate::execute(), BbEmc::execute(), KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), KalFitTrack::order_hits(), KalFitTrack::order_wirhit(), MdcBbEmcEff::selectBbByEmcShower(), and TEvent::Swap().
void lorbck | ( | double * | pcm, | |
double * | pincm, | |||
double * | pinlab | |||
) |
Definition at line 33 of file TEvent.C.
References genRecEmupikp::i.
00033 { 00034 double beta[3]; 00035 for(int i=0;i<3;i++) 00036 beta[i] = pcm[i]/pcm[3]; 00037 double gamma = 1./sqrt(1 - beta[0]*beta[0] - beta[1]*beta[1] - beta[2]*beta[2]); 00038 double bdotp = beta[0]*pincm[0] + beta[1]*pincm[1] + beta[2]*pincm[2]; 00039 double fact = gamma*bdotp/(1 + gamma) + pincm[3]; 00040 00041 for(int i=0; i<3; i++) 00042 pinlab[i] = pincm[i] + beta[i]*fact*gamma; 00043 pinlab[3] = gamma*(pincm[3] + bdotp); 00044 }