|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementBrem.h"
#include "DataFormats/ParticleFlowReco/interface/PFRecTrack.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/ParticleFlowReco/interface/PFTrajectoryPoint.h"
#include <iomanip>
using namespace reco;
using namespace std;
PFBlockElementBrem::PFBlockElementBrem(const GsfPFRecTrackRef& gsfref,
const double DeltaP,
const double SigmaDeltaP,
const unsigned int indTrajPoint)
: PFBlockElement(BREM),
GsftrackRefPF_(gsfref),
GsftrackRef_(gsfref->gsfTrackRef()),
deltaP_(DeltaP),
sigmadeltaP_(SigmaDeltaP),
indPoint_(indTrajPoint) {
const reco::PFTrajectoryPoint& atECAL =
((*GsftrackRefPF()).PFRecBrem()[(indPoint_ - 2)]).extrapolatedPoint(reco::PFTrajectoryPoint::ECALEntrance);
if (atECAL.isValid())
positionAtECALEntrance_.SetCoordinates(atECAL.position().x(), atECAL.position().y(), atECAL.position().z());
}
void PFBlockElementBrem::Dump(ostream& out, const char* tab) const {
if (!out)
return;
if (!GsftrackRefPF_.isNull()) {
double charge = 0.;
double dp = deltaP_;
double sigmadp = sigmadeltaP_;
int indextrj = (indPoint_ - 2);
out << setprecision(0);
out << tab << setw(7) << "charge=" << setw(3) << charge;
out << setprecision(3);
out << setiosflags(ios::right);
out << setiosflags(ios::fixed);
out << ", DeltaP= " << dp;
out << ", SigmaDeltaP= " << sigmadp;
out << ", Traj Point= " << indextrj;
out << resetiosflags(ios::right | ios::fixed);
}
}
|