#include <CosmicEventParser.h>
Public Member Functions | |
CosmicEventParser () | |
CosmicEventParser () | |
const HepLorentzVector & | Momentum (void) |
const HepLorentzVector & | Momentum (void) |
int | pdgID (void) |
int | pdgID (void) |
const HepLorentzVector & | Vertex (void) |
const HepLorentzVector & | Vertex (void) |
Private Member Functions | |
std::istream & | read (std::istream &is) |
std::istream & | read (std::istream &is) |
std::ostream & | write (std::ostream &os) const |
std::ostream & | write (std::ostream &os) const |
Private Attributes | |
int | m_eventNumber |
HepLorentzVector | m_momentum |
int | m_pdgId |
HepLorentzVector | m_vertex |
Friends | |
std::ostream & | operator<< (std::ostream &os, const CosmicEventParser &ev) |
std::ostream & | operator<< (std::ostream &os, const CosmicEventParser &ev) |
std::istream & | operator>> (std::istream &is, CosmicEventParser &ev) |
std::istream & | operator>> (std::istream &is, CosmicEventParser &ev) |
|
00010 {};
|
|
00010 {};
|
|
00012 {return m_momentum;}
|
|
00012 {return m_momentum;}
|
|
00013 {return m_pdgId;}
|
|
00013 {return m_pdgId;}
|
|
|
|
00036 { 00037 int dummy; 00038 int id; 00039 Hep3Vector vert,mom; 00040 double v_x,v_y,v_z; 00041 is >> m_eventNumber >> dummy >> id >> v_x >> v_y >> v_z >> mom; 00042 00043 // 00044 // rotate over pi in x-z plane 00045 // 00046 vert.setX(-v_x); 00047 vert.setY( v_y); 00048 vert.setZ(-v_z); 00049 00050 // 00051 // convert to MeV's and mm units 00052 // 00053 mom = 1000*mom; 00054 vert = 10*vert; 00055 00056 00057 m_vertex.setVect(vert); 00058 m_vertex.setE(0.); 00059 00060 00061 m_momentum.setVect(mom); 00062 double energy = sqrt(pow(105.66,2)+mom.mag2()); 00063 m_momentum.setE(energy); 00064 00065 if(id == 5) m_pdgId = 13; 00066 else m_pdgId = -13; 00067 00068 return is; 00069 }
|
|
00011 {return m_vertex;}
|
|
00011 {return m_vertex;}
|
|
|
|
00073 { 00074 int dummy(1); 00075 int id(5); 00076 if(m_pdgId == -13) id = 6; 00077 os << m_eventNumber << " " << dummy << " " << id << " " 00078 << m_vertex.x() << " " << m_vertex.y() << " " << m_vertex.z() << " " 00079 << m_momentum.x() << " " << m_momentum.y() << " " << m_momentum.z(); 00080 00081 return os; 00082 }
|
|
00030 { 00031 return ev.write(os); 00032 }
|
|
00030 { 00031 return ev.write(os); 00032 }
|
|
00026 { 00027 return ev.read(is); 00028 }
|
|
00026 { 00027 return ev.read(is); 00028 }
|
|
|
|
|
|
|
|
|