9 #include <daq/dqm/modules/TrackAna.h>
11 #include <framework/dataobjects/EventMetaData.h>
12 #include <framework/datastore/StoreObjPtr.h>
13 #include <framework/datastore/StoreArray.h>
14 #include <mdst/dataobjects/Track.h>
15 #include <mdst/dataobjects/TrackFitResult.h>
32 setDescription(
"The simplest physics analysis");
33 setPropertyFlags(c_ParallelProcessingCertified);
36 B2INFO(
"TrackAna: Constructor done.");
40 TrackAnaModule::~TrackAnaModule()
44 void TrackAnaModule::initialize()
49 void TrackAnaModule::defineHisto()
51 h_multi =
new TH1F(
"Multi",
"Particle Multiplicity", 50, 0.0, 50.0);
52 h_p[0] =
new TH1F(
"Px",
"Particle Momentum X", 100, -5.0, 5.0);
53 h_p[1] =
new TH1F(
"Py",
"Particle Momentum Y", 100, -5.0, 5.0);
54 h_p[2] =
new TH1F(
"Pz",
"Particle Momentum Z", 100, -5.0, 5.0);
55 h_p[3] =
new TH1F(
"E",
"Particle Energy", 100, 0.0, 10.0);
58 void TrackAnaModule::beginRun()
60 B2INFO(
"TrackAna: started to measure elapsed time.");
64 void TrackAnaModule::event()
68 int expno = evt->getExperiment();
69 int runno = evt->getRun();
70 int subrno = evt->getSubrun();
71 int evtno = evt->getEvent();
77 h_multi->Fill((
float)ntrk);
79 for (
int i = 0; i < ntrk; i++) {
80 Track* trk = trklist[i];
82 TLorentzVector p4 = fit->get4Momentum();
83 for (
int j = 0; j < 4; j++) {
89 void TrackAnaModule::endRun()
94 void TrackAnaModule::terminate()
96 B2INFO(
"TrackAna: terminate called");