00001 #ifndef bmptevent_cxx
00002 #define bmptevent_cxx
00003 #include "MessageService/MsgService.h"
00004 #include "MCReweight/BMPTEvent.h"
00005 #include <math.h>
00006
00007 CVSID("$Id: BMPTEvent.cxx,v 1.2 2005/03/07 11:21:57 tagg Exp $");
00008
00009
00010 BMPTEvent::BMPTEvent() :
00011 x(0),y(0),z(0),px(0),py(0),pz(0),pid(0),gen(0)
00012 {
00013 MSG("BMPT",Msg::kDebug) << "In BMPTEvent Constructor" << endl;
00014 }
00015
00016
00017 BMPTEvent::BMPTEvent(double X,double Y,double Z,
00018 double PX,double PY,double PZ,
00019 int PID,int GEN)
00020 {
00021 x=X; y=Y; z=Z;
00022 px=PX; py=PY; pz=PZ;
00023 pid=PID; gen=GEN;
00024 }
00025
00026
00027 BMPTEvent::~BMPTEvent()
00028 {
00029 MSG("BMPT",Msg::kDebug) << "In BMPTEvent Destructor" << endl;
00030 }
00031
00032
00033 double BMPTEvent::Momentum()
00034 {
00035 return sqrt(px*px + py*py + pz*pz);
00036 }
00037
00038
00039 double BMPTEvent::Mass()
00040 {
00041 if(pid==14) return 1.1;
00042 else if(pid==15) return 1.0;
00043 return 0;
00044 }
00045 #endif