Class Reference for E1039 Core & Analysis Software
AnaPileup.cc
Go to the documentation of this file.
1 
10 #include "AnaPileup.h"
11 
12 
13 #include <interface_main/SQHitVector.h>
14 #include <interface_main/SQEvent.h>
15 #include <interface_main/SQMCEvent.h>
16 #include <interface_main/SQTrackVector.h>
17 #include <interface_main/SQDimuonVector.h>
18 
19 #include <ktracker/SRecEvent.h>
20 #include <geom_svc/GeomSvc.h>
21 
22 #include <fun4all/Fun4AllReturnCodes.h>
23 #include <fun4all/PHTFileServer.h>
24 #include <phool/PHNodeIterator.h>
25 #include <phool/PHIODataNode.h>
26 #include <phool/getClass.h>
27 
28 
29 #include <TFile.h>
30 #include <TTree.h>
31 #include <TClonesArray.h>
32 
33 
34 #define LogError(exp) std::cout<<"ERROR: " <<__FILE__<<": "<<__LINE__<<": "<< exp << std::endl
35 #define LogWarning(exp) std::cout<<"WARNING: "<<__FILE__<<": "<<__LINE__<<": "<< exp << std::endl
36 
37 using namespace std;
38 
39 AnaPileup::AnaPileup(const std::string& name) :
40  SubsysReco(name),
41  _event(0),
42  _event_header(nullptr),
43  _hit_vector(nullptr),
44  _out_name("eval.root"),
45  legacyContainer(true)
46 {
47 
48 }
49 AnaPileup::~AnaPileup()
50 {
51 
52  delete pos1;
53  delete pos2;
54  delete pos3;
55  delete posvtx;
56  delete mom1;
57  delete mom2;
58  delete mom3;
59  delete momvtx;
60  delete rec_mom1;
61  delete rec_momvtx;
62  delete rec_posvtx;
63  delete rec_momtgt;
64  delete rec_postgt;
65 
66  delete rec_pmom;
67  delete rec_nmom;
68  delete rec_ppos;
69  delete rec_npos;
70  delete rec_vtx;
71 
72 }
73 
74 
75 int AnaPileup::Init(PHCompositeNode* topNode) {
76 
77  pos1 = new TVector3();
78  pos2 = new TVector3();
79  pos3 = new TVector3();
80  posvtx = new TVector3();
81  mom1 = new TVector3();
82  mom2 = new TVector3();
83  mom3 = new TVector3();
84  momvtx = new TVector3();
85  rec_mom1 = new TVector3();
86  rec_momvtx = new TVector3();
87  rec_posvtx = new TVector3();
88  rec_momtgt = new TVector3();
89  rec_postgt = new TVector3();
90 
91  rec_pmom = new TVector3();
92  rec_nmom = new TVector3();
93  rec_ppos = new TVector3();
94  rec_npos = new TVector3();
95  rec_vtx = new TVector3();
96 
97  return Fun4AllReturnCodes::EVENT_OK;
98 }
99 
100 int AnaPileup::InitRun(PHCompositeNode* topNode) {
101 
102  ResetEvalVars();
103  InitEvalTree();
104 
105  p_geomSvc = GeomSvc::instance();
106 
107  int ret = GetNodes(topNode);
108  if(ret != Fun4AllReturnCodes::EVENT_OK) return ret;
109 
110  return Fun4AllReturnCodes::EVENT_OK;
111 }
112 
113 int AnaPileup::process_event(PHCompositeNode* topNode) {
114 
115  fpga1 = _event_header->get_trigger(SQEvent::MATRIX1);
116  fpga2 = _event_header->get_trigger(SQEvent::MATRIX2);
117  fpga3 = _event_header->get_trigger(SQEvent::MATRIX3);
118  fpga4 = _event_header->get_trigger(SQEvent::MATRIX4);
119  fpga5 = _event_header->get_trigger(SQEvent::MATRIX5);
120 
121  nim1 = _event_header->get_trigger(SQEvent::NIM1);
122  nim2 = _event_header->get_trigger(SQEvent::NIM2);
123 
124  n_bucketsize = mi_evt_true->get_cross_section();
125  eventID = _event_header->get_event_id();
126 
127  Eval_sqhit(topNode);
128  Eval_trkdim(topNode);
129 
130  evt_tree->Fill();
131 
132  ++_event;
133 
134  return Fun4AllReturnCodes::EVENT_OK;
135 }
136 
137 
139 int AnaPileup::End(PHCompositeNode* topNode) {
140  if(Verbosity() >= Fun4AllBase::VERBOSITY_A_LOT)
141  std::cout << "AnaPileup::End" << std::endl;
142 
143  PHTFileServer::get().cd(_out_name.c_str());
144  _sqhits_tree->Write();
145  trkTree->Write();
146  dimTree->Write();
147  evt_tree->Write();
148 
149  return Fun4AllReturnCodes::EVENT_OK;
150 }
151 //Evaluation fuction==============================================
152 int AnaPileup::Eval_sqhit(PHCompositeNode* topNode)
153 {
154  ResetEvalVars();
155 
156  h1Lhit = 0;
157  h1Rhit = 0;
158 
159  h1Bhit = 0;
160  h1Thit = 0;
161 
162  GeomSvc* geomSvc = GeomSvc::instance();
163 
164  for(int ihit=0; ihit<_hit_vector->size(); ++ihit) {
165  SQHit *sqhit = _hit_vector->at(ihit);
166  int sq_detid = sqhit->get_detector_id();
167 
169  if (sq_detid == geomSvc->getDetectorID("H1B")){
170  new((*pos_H1B)[h1Bhit]) TVector3(sqhit->get_truth_x(), sqhit->get_truth_y(), sqhit->get_truth_z());
171  new((*mom_H1B)[h1Bhit]) TVector3(sqhit->get_truth_px(), sqhit->get_truth_py(), sqhit->get_truth_pz());
172  eID_H1B[h1Bhit] = sqhit->get_element_id();
173  ++h1Bhit;
174  }
175  if (sq_detid == geomSvc->getDetectorID("H1T")){
176  new((*pos_H1T)[h1Thit]) TVector3(sqhit->get_truth_x(), sqhit->get_truth_y(), sqhit->get_truth_z());
177  new((*mom_H1T)[h1Thit]) TVector3(sqhit->get_truth_px(), sqhit->get_truth_py(), sqhit->get_truth_pz());
178  eID_H1T[h1Thit] = sqhit->get_element_id();
179  ++h1Thit;
180  }
181  if (sq_detid == geomSvc->getDetectorID("H1L")){
182  eID_H1L[h1Lhit] = sqhit->get_element_id();
183  new((*pos_H1L)[h1Lhit]) TVector3(sqhit->get_truth_x(), sqhit->get_truth_y(), sqhit->get_truth_z());
184  new((*mom_H1L)[h1Lhit]) TVector3(sqhit->get_truth_px(), sqhit->get_truth_py(), sqhit->get_truth_pz());
185  ++h1Lhit;
186  }
187  if (sq_detid == geomSvc->getDetectorID("H1R")){
188  eID_H1R[h1Rhit] = sqhit->get_element_id();
189  new((*pos_H1R)[h1Rhit]) TVector3(sqhit->get_truth_x(), sqhit->get_truth_y(), sqhit->get_truth_z());
190  new((*mom_H1R)[h1Rhit]) TVector3(sqhit->get_truth_px(), sqhit->get_truth_py(), sqhit->get_truth_pz());
191  ++h1Rhit;
192  }
193 
194 
195  D0Xhit = 0;
197  if (sq_detid == geomSvc->getDetectorID("D0X") ) {
198  eID_D0X[D0Xhit] = sqhit->get_element_id();
199  new((*pos_D0X)[D0Xhit]) TVector3(sqhit->get_truth_x(), sqhit->get_truth_y(), sqhit->get_truth_z());
200  new((*mom_D0X)[D0Xhit]) TVector3(sqhit->get_truth_px(), sqhit->get_truth_py(), sqhit->get_truth_pz());
201  ++D0Xhit;
202  }
203 
204 
205 
206  }
207 
208  _sqhits_tree->Fill();
209 
210  pos_H1B->Delete();
211  pos_H1R->Delete();
212  pos_H1L->Delete();
213  pos_H1T->Delete();
214 
215  mom_H1B->Delete();
216  mom_H1L->Delete();
217  mom_H1R->Delete();
218  mom_H1T->Delete();
219 
220  pos_D0X->Delete();
221  mom_D0X->Delete();
222 
223  return Fun4AllReturnCodes::EVENT_OK;
224 }
225 
227 int AnaPileup::Eval_trkdim(PHCompositeNode* topNode)
228 {
229  int nTracks = trackVector->size();
230  int nRecTracks = legacyContainer ? _recEvent->getNTracks() : recTrackVector->size();
231  for(int i = 0; i < nTracks; ++i)
232  {
233  SQTrack* track = trackVector->at(i);
234  *pos1 = track->get_pos_st1();
235  *mom1 = track->get_mom_st1().Vect();
236  *pos3 = track->get_pos_st3();
237  *mom3 = track->get_mom_st3().Vect();
238  *posvtx = track->get_pos_vtx();
239  *momvtx = track->get_mom_vtx().Vect();
240 
241  int recid = track->get_rec_track_id();
242  if(recid >= 0 && recid < nRecTracks)
243  {
244  SRecTrack* recTrack = legacyContainer ? &(_recEvent->getTrack(recid)) : dynamic_cast<SRecTrack*>(recTrackVector->at(recid));
245  *rec_mom1 = recTrack->getMomentumVecSt1();
246  *rec_momvtx = recTrack->getVertexMom();
247  *rec_posvtx = recTrack->getVertexPos();
248  *rec_momtgt = recTrack->getTargetMom();
249  *rec_postgt = recTrack->getTargetPos();
250  }
251  else
252  {
253  rec_mom1->SetXYZ(-999., -999., -999.);
254  rec_momvtx->SetXYZ(-999., -999., -999.);
255  rec_posvtx->SetXYZ(-999., -999., -999.);
256  }
257 
258  trkTree->Fill();
259  }
260 
261 
263  int nDimuons = dimuonVector->size();
264  int nRecDimuons = legacyContainer ? _recEvent->getNDimuons() : (recDimuonVector ? recDimuonVector->size() : -1);
265  for(int i = 0; i < nRecDimuons; ++i)
266  {
267  SRecDimuon *recDimuon = &(_recEvent->getDimuon(i));
268  rec_mass = recDimuon->mass;
269  *rec_pmom = recDimuon->p_pos.Vect();
270  *rec_nmom = recDimuon->p_neg.Vect();
271  *rec_ppos = recDimuon->vtx_pos;
272  *rec_npos = recDimuon->vtx_neg;
273  *rec_vtx = recDimuon->vtx;
274 
275  dimTree->Fill();
276  }
277  return Fun4AllReturnCodes::EVENT_OK;
278 }
279 
280 
281 int AnaPileup::InitEvalTree() {
282 
283  PHTFileServer::get().open(_out_name.c_str(), "RECREATE");
284 
285  mom_H1T = new TClonesArray("TVector3");
286  pos_H1T = new TClonesArray("TVector3");
287  mom_H1B = new TClonesArray("TVector3");
288  pos_H1B = new TClonesArray("TVector3");
289 
290  mom_H1L = new TClonesArray("TVector3");
291  pos_H1L = new TClonesArray("TVector3");
292  mom_H1R = new TClonesArray("TVector3");
293  pos_H1R = new TClonesArray("TVector3");
294 
295  mom_D0X = new TClonesArray("TVector3");
296  pos_D0X = new TClonesArray("TVector3");
297 
298 
299  evt_tree = new TTree("pileup_evt", "Event level info from piling up bkg interactions");
300 
301  evt_tree->Branch("n_bucketsize", &n_bucketsize, "n_bucketsize/I");
302  evt_tree->Branch("eventID", &eventID, "eventID/I");
303 
304  evt_tree->Branch("fpga1", &fpga1, "fpga1/I");
305  evt_tree->Branch("fpga2", &fpga2, "fpga2/I");
306  evt_tree->Branch("fpga3", &fpga3, "fpga3/I");
307  evt_tree->Branch("fpga4", &fpga4, "fpga4/I");
308  evt_tree->Branch("fpga5", &fpga5, "fpga5/I");
309 
310  evt_tree->Branch("nim1", &nim1, "nim1/I");
311  evt_tree->Branch("nim2", &nim2, "nim2/I");
312 
314  _sqhits_tree = new TTree("pileup_sqhits", "Analsysis of sqhits from piling up bkg interactions");
315 
316  _sqhits_tree->Branch("n_bucketsize", &n_bucketsize, "n_bucketsize/I");
317 
319  _sqhits_tree->Branch("fpga1", &fpga1, "fpga1/I");
320  _sqhits_tree->Branch("fpga2", &fpga2, "fpga2/I");
321  _sqhits_tree->Branch("fpga3", &fpga3, "fpga3/I");
322  _sqhits_tree->Branch("fpga4", &fpga4, "fpga4/I");
323  _sqhits_tree->Branch("fpga5", &fpga5, "fpga5/I");
324 
325  _sqhits_tree->Branch("nim1", &nim1, "nim1/I");
326  _sqhits_tree->Branch("nim2", &nim2, "nim2/I");
327 
329  _sqhits_tree->Branch("h1Bhit", &h1Bhit, "h1Bhit/I");
330  _sqhits_tree->Branch("h1Thit", &h1Thit, "h1Thit/I");
331  _sqhits_tree->Branch("eID_H1T", eID_H1T, "eID_H1T[h1Thit]/I");
332  _sqhits_tree->Branch("eID_H1B", eID_H1B, "eID_H1B[h1Bhit]/I");
333 
334  _sqhits_tree->Branch("h1Lhit", &h1Lhit, "h1Lhit/I");
335  _sqhits_tree->Branch("h1Rhit", &h1Rhit, "h1Rhit/I");
336  _sqhits_tree->Branch("eID_H1L", eID_H1L, "eID_H1L[h1Lhit]/I");
337  _sqhits_tree->Branch("eID_H1R", eID_H1R, "eID_H1R[h1Rhit]/I");
338 
339  _sqhits_tree->Branch("mom_H1T", &mom_H1T, 256000, 99);
340  _sqhits_tree->Branch("pos_H1T", &pos_H1T, 256000, 99);
341  _sqhits_tree->Branch("mom_H1B", &mom_H1B, 256000, 99);
342  _sqhits_tree->Branch("pos_H1B", &pos_H1B, 256000, 99);
343 
344  _sqhits_tree->Branch("mom_H1L", &mom_H1L, 256000, 99);
345  _sqhits_tree->Branch("pos_H1L", &pos_H1L, 256000, 99);
346  _sqhits_tree->Branch("mom_H1R", &mom_H1R, 256000, 99);
347  _sqhits_tree->Branch("pos_H1R", &pos_H1R, 256000, 99);
348 
350  _sqhits_tree->Branch("D0Xhit", &D0Xhit, "D0Xhit/I");
351  _sqhits_tree->Branch("eID_D0X", eID_D0X, "eID_D0X[D0Xhit]/I");
352  _sqhits_tree->Branch("mom_D0X", &mom_D0X, 256000, 99);
353  _sqhits_tree->Branch("pos_D0X", &pos_D0X, 256000, 99);
354 
355  mom_H1T->BypassStreamer();
356  pos_H1T->BypassStreamer();
357  mom_H1B->BypassStreamer();
358  pos_H1B->BypassStreamer();
359 
360  mom_H1L->BypassStreamer();
361  pos_H1L->BypassStreamer();
362  mom_H1R->BypassStreamer();
363  pos_H1R->BypassStreamer();
364 
365  mom_D0X->BypassStreamer();
366  pos_D0X->BypassStreamer();
367 
368  trkTree = new TTree("pileup_trk", "Track Tree Created by AnaPileup Module");
369  trkTree->Branch("pos1", &pos1, 256000, 99);
370  trkTree->Branch("pos2", &pos2, 256000, 99);
371  trkTree->Branch("pos3", &pos3, 256000, 99);
372  trkTree->Branch("posvtx", &posvtx, 256000, 99);
373  trkTree->Branch("mom1", &mom1, 256000, 99);
374  trkTree->Branch("mom2", &mom2, 256000, 99);
375  trkTree->Branch("mom3", &mom3, 256000, 99);
376  trkTree->Branch("momvtx", &momvtx, 256000, 99);
377  trkTree->Branch("rec_mom1", &rec_mom1, 256000, 99);
378  trkTree->Branch("rec_momvtx", &rec_momvtx, 256000, 99);
379  trkTree->Branch("rec_posvtx", &rec_posvtx, 256000, 99);
380  trkTree->Branch("rec_momtgt", &rec_momtgt, 256000, 99);
381  trkTree->Branch("rec_postgt", &rec_postgt, 256000, 99);
382 
383  dimTree = new TTree("pileup_dim", "Dimuon Tree Created by AnaPileup Module");
384  dimTree->Branch("rec_mass", &rec_mass, "rec_mass/D");
385  dimTree->Branch("rec_pmom", &rec_pmom, 256000, 99);
386  dimTree->Branch("rec_nmom", &rec_nmom, 256000, 99);
387  dimTree->Branch("rec_ppos", &rec_ppos, 256000, 99);
388  dimTree->Branch("rec_npos", &rec_npos, 256000, 99);
389  dimTree->Branch("rec_vtx", &rec_vtx, 256000, 99);
390 
391  return 0;
392 }
393 
394 int AnaPileup::ResetEvalVars() {
395 
396  run_id = std::numeric_limits<int>::max();
397  spill_id = std::numeric_limits<int>::max();
398  event_id = std::numeric_limits<int>::max();
399  emu_trigger = 0;
400 
401 
402  return 0;
403 }
404 
405 int AnaPileup::GetNodes(PHCompositeNode* topNode) {
406 
407 
408  _event_header = findNode::getClass<SQEvent>(topNode, "SQEvent");
409  if (!_event_header) {
410  LogError("!_event_header");
411  //return Fun4AllReturnCodes::ABORTEVENT;
412  }
413 
414  mi_evt_true = findNode::getClass<SQMCEvent >(topNode, "SQMCEvent");
415 
416 
417  _hit_vector = findNode::getClass<SQHitVector>(topNode, "SQHitVector");
418  if (!_hit_vector) {
419  LogError("!_hit_vector");
420  return Fun4AllReturnCodes::ABORTEVENT;
421  }
422 
423 
424  trackVector = findNode::getClass<SQTrackVector>(topNode, "SQTruthTrackVector");
425  dimuonVector = findNode::getClass<SQDimuonVector>(topNode, "SQTruthDimuonVector");
426  if(!trackVector || !dimuonVector)
427  {
428  return Fun4AllReturnCodes::ABORTEVENT;
429  }
430 
431  if(legacyContainer)
432  {
433  _recEvent = findNode::getClass<SRecEvent>(topNode, "SRecEvent");
434  if (!_recEvent)
435  {
436  LogError("!_recEvent");
437  //return Fun4AllReturnCodes::ABORTEVENT;
438  }
439  }
440  else
441  {
442  recTrackVector = findNode::getClass<SQTrackVector>(topNode, "SQRecTrackVector");
443  recDimuonVector = findNode::getClass<SQDimuonVector>(topNode, "SQRecDimuonVector");
444  if(!recTrackVector)
445  {
446  return Fun4AllReturnCodes::ABORTEVENT;
447  }
448  }
449 
450 
451  return Fun4AllReturnCodes::EVENT_OK;
452 }
453 
LogError
#define LogError(exp)
Definition: AnaPileup.cc:34
PHTFileServer.h
TFile clean handling.
AnaPileup::AnaPileup
AnaPileup(const std::string &name="AnaPileup.root")
Definition: AnaPileup.cc:39
AnaPileup::InitRun
int InitRun(PHCompositeNode *topNode)
Definition: AnaPileup.cc:100
AnaPileup::End
int End(PHCompositeNode *topNode)
===========================
Definition: AnaPileup.cc:139
AnaPileup::~AnaPileup
virtual ~AnaPileup()
Definition: AnaPileup.cc:49
AnaPileup::process_event
int process_event(PHCompositeNode *topNode)
Definition: AnaPileup.cc:113
AnaPileup::Init
int Init(PHCompositeNode *topNode)
Definition: AnaPileup.cc:75
AnaPileup::ResetEvalVars
int ResetEvalVars()
Definition: AnaPileup.cc:394
AnaPileup::InitEvalTree
int InitEvalTree()
Definition: AnaPileup.cc:281
Fun4AllBase::VERBOSITY_A_LOT
@ VERBOSITY_A_LOT
Output a lot of messages.
Definition: Fun4AllBase.h:48
Fun4AllBase::Verbosity
virtual int Verbosity() const
Gets the verbosity of this module.
Definition: Fun4AllBase.h:64
GeomSvc
User interface class about the geometry of detector planes.
Definition: GeomSvc.h:164
GeomSvc::getDetectorID
int getDetectorID(const std::string &detectorName) const
Get the plane position.
Definition: GeomSvc.h:219
GeomSvc::instance
static GeomSvc * instance()
singlton instance
Definition: GeomSvc.cxx:212
PHTFileServer::get
static PHTFileServer & get(void)
return reference to class singleton
Definition: PHTFileServer.h:36
PHTFileServer::open
void open(const std::string &filename, const std::string &type="RECREATE")
open a SafeTFile. If filename is not found in the map, create a new TFile and append to the map; incr...
Definition: PHTFileServer.cxx:29
PHTFileServer::cd
bool cd(const std::string &filename)
change to directory of TFile matching filename
Definition: PHTFileServer.cxx:76
SQDimuonVector::size
virtual size_t size() const =0
SQEvent::MATRIX2
@ MATRIX2
Definition: SQEvent.h:28
SQEvent::MATRIX1
@ MATRIX1
Definition: SQEvent.h:27
SQEvent::NIM2
@ NIM2
Definition: SQEvent.h:23
SQEvent::NIM1
@ NIM1
Definition: SQEvent.h:22
SQEvent::MATRIX3
@ MATRIX3
Definition: SQEvent.h:29
SQEvent::MATRIX4
@ MATRIX4
Definition: SQEvent.h:30
SQEvent::MATRIX5
@ MATRIX5
Definition: SQEvent.h:31
SQEvent::get_trigger
virtual bool get_trigger(const SQEvent::TriggerMask i) const =0
Return the trigger bit (fired or not) of the selected trigger channel.
SQEvent::get_event_id
virtual int get_event_id() const =0
Return the event ID, which is unique per run.
SQHitVector::at
virtual const SQHit * at(const size_t idkey) const =0
SQHitVector::size
virtual size_t size() const =0
SQHit
An SQ interface class to hold one detector hit.
Definition: SQHit.h:20
SQHit::get_truth_z
virtual float get_truth_z() const
Return the true z-position of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:78
SQHit::get_truth_pz
virtual float get_truth_pz() const
Return the true z-momentum of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:87
SQHit::get_truth_y
virtual float get_truth_y() const
Return the true y-position of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:75
SQHit::get_truth_px
virtual float get_truth_px() const
Return the true x-momentum of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:81
SQHit::get_element_id
virtual short get_element_id() const
Return the element ID of this hit.
Definition: SQHit.h:45
SQHit::get_truth_py
virtual float get_truth_py() const
Return the true y-momentum of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:84
SQHit::get_truth_x
virtual float get_truth_x() const
Return the true x-position of this hit. Meaningful only if this hit is of MC.
Definition: SQHit.h:72
SQHit::get_detector_id
virtual short get_detector_id() const
Return the detector ID of this hit.
Definition: SQHit.h:42
SQMCEvent::get_cross_section
virtual double get_cross_section() const =0
Return the cross section.
SQTrackVector::at
virtual const SQTrack * at(const size_t id) const =0
SQTrackVector::size
virtual size_t size() const =0
SQTrack
An SQ interface class to hold one true or reconstructed track.
Definition: SQTrack.h:8
SQTrack::get_pos_vtx
virtual TVector3 get_pos_vtx() const =0
Return the track position at vertex.
SQTrack::get_rec_track_id
virtual int get_rec_track_id() const =0
Return the track ID of associated reconstructed track. Valid only if this object holds truth track in...
SQTrack::get_mom_vtx
virtual TLorentzVector get_mom_vtx() const =0
Return the track momentum at vertex.
SQTrack::get_mom_st1
virtual TLorentzVector get_mom_st1() const =0
Return the track momentum at Station 1.
SQTrack::get_pos_st3
virtual TVector3 get_pos_st3() const =0
Return the track position at Station 3.
SQTrack::get_pos_st1
virtual TVector3 get_pos_st1() const =0
Return the track position at Station 1.
SQTrack::get_mom_st3
virtual TLorentzVector get_mom_st3() const =0
Return the track momentum at Station 3.
SRecDimuon::vtx
TVector3 vtx
3-vector vertex position
Definition: SRecEvent.h:379
SRecDimuon::p_neg
TLorentzVector p_neg
Definition: SRecEvent.h:366
SRecDimuon::p_pos
TLorentzVector p_pos
4-momentum of the muon tracks after vertex fit
Definition: SRecEvent.h:365
SRecDimuon::mass
Double_t mass
Kinematic variables.
Definition: SRecEvent.h:390
SRecDimuon::vtx_neg
TVector3 vtx_neg
Definition: SRecEvent.h:381
SRecDimuon::vtx_pos
TVector3 vtx_pos
Definition: SRecEvent.h:380
SRecEvent::getDimuon
SRecDimuon & getDimuon(Int_t i)
Definition: SRecEvent.h:463
SRecEvent::getNTracks
Int_t getNTracks()
Get tracks.
Definition: SRecEvent.h:455
SRecEvent::getTrack
SRecTrack & getTrack(Int_t i)
Definition: SRecEvent.h:456
SRecEvent::getNDimuons
Int_t getNDimuons()
Get dimuons.
Definition: SRecEvent.h:462
SRecTrack::getTargetPos
TVector3 getTargetPos()
Definition: SRecEvent.h:188
SRecTrack::getVertexPos
TVector3 getVertexPos()
Definition: SRecEvent.h:196
SRecTrack::getMomentumVecSt1
TVector3 getMomentumVecSt1() const
Definition: SRecEvent.h:129
SRecTrack::getVertexMom
TVector3 getVertexMom()
Definition: SRecEvent.h:197
SRecTrack::getTargetMom
TVector3 getTargetMom()
Definition: SRecEvent.h:193