| 1 |
grimes |
1.1 |
#ifndef SimDataFormats_TrackingDataPrint_h
|
| 2 |
|
|
#define SimDataFormats_TrackingDataPrint_h
|
| 3 |
|
|
/** Concrete TrackingParticle.
|
| 4 |
|
|
* All track parameters are passed in the constructor and stored internally.
|
| 5 |
|
|
*/
|
| 6 |
|
|
|
| 7 |
|
|
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
|
| 8 |
|
|
#include "SimDataFormats/TrackingAnalysis/interface/TrackingVertex.h"
|
| 9 |
|
|
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticleFwd.h"
|
| 10 |
|
|
#include <iostream>
|
| 11 |
|
|
|
| 12 |
|
|
std::ostream& operator<< (std::ostream& s, const TrackingParticle & tp)
|
| 13 |
|
|
{
|
| 14 |
|
|
|
| 15 |
|
|
// Compare momenta from sources
|
| 16 |
|
|
s << "T.P. Track Momentum, q , ID, & Event # "
|
| 17 |
|
|
<< tp.p4() << " " << tp.charge() << " "
|
| 18 |
|
|
<< tp.pdgId() << " "
|
| 19 |
|
|
<< tp.eventId().bunchCrossing() << "." << tp.eventId().event()
|
| 20 |
|
|
<< std::endl;
|
| 21 |
|
|
s << " Hits for this track: " << tp.trackPSimHit().size()
|
| 22 |
|
|
<< std::endl;
|
| 23 |
|
|
|
| 24 |
|
|
for (TrackingParticle::genp_iterator hepT = tp.genParticle_begin();
|
| 25 |
|
|
hepT != tp.genParticle_end(); ++hepT)
|
| 26 |
|
|
{
|
| 27 |
|
|
s << " HepMC Track Momentum " << (*hepT)->momentum().mag() << std::endl;
|
| 28 |
|
|
}
|
| 29 |
|
|
for (TrackingParticle::g4t_iterator g4T = tp.g4Track_begin();
|
| 30 |
|
|
g4T != tp.g4Track_end(); ++g4T)
|
| 31 |
|
|
{
|
| 32 |
|
|
s << " Geant Track Momentum " << g4T->momentum() << std::endl;
|
| 33 |
|
|
s << " Geant Track ID & type " << g4T->trackId() << " "
|
| 34 |
|
|
<< g4T->type() << std::endl;
|
| 35 |
|
|
if (g4T->type() != tp.pdgId())
|
| 36 |
|
|
{
|
| 37 |
|
|
s << " Mismatch b/t TrackingParticle and Geant types"
|
| 38 |
|
|
<< std::endl;
|
| 39 |
|
|
}
|
| 40 |
|
|
}
|
| 41 |
|
|
return s;
|
| 42 |
|
|
}
|
| 43 |
|
|
|
| 44 |
|
|
#endif // SimDataFormats_TrackingDataPrint_H
|
