125 addLine(TString::Format(
" %d items in selection", (
int) sted.size()));
127 addLine(
"--------------------------------------------------+--------------");
128 addLine(
" px py pz pT | Collection");
129 addLine(
"--------------------------------------------------+--------------");
136 double sum_len_xy = 0;
140 for (std::set<FWModelId>::const_iterator
i = sted.begin();
i != sted.end(); ++
i, ++
n)
144 TClass *model_class =
const_cast<TClass*
>(
i->item()->modelType());
145 void *model_data =
const_cast<void*
> (
i->item()->modelData(
i->index()));
151 (model_class->DynamicCast(rc_class, model_data));
156 v.SetXYZ(rc->
px(), rc->
py(), rc->
pz());
161 (model_class->DynamicCast(rtb_class, model_data));
166 v.SetXYZ(rtb->
px(), rtb->
py(), rtb->
pz());
173 sum_len += TMath::Sqrt(v.mag2());
174 sum_len_xy += TMath::Sqrt(v.perp2());
176 line = TString::Format(
" %+10.3f %+10.3f %+10.3f %10.3f", v.x(), v.y(), v.z(), TMath::Sqrt(v.perp2()));
181 line = TString::Format(
" -------- not a Candidate or TrackBase --------");
183 line += TString::Format(
" | %s[%d]",
i->item()->name().c_str(),
i->index());
187 if (n == 0) first =
v;
else if (n == 1) second =
v;
190 addLine(
"--------------------------------------------------+--------------");
191 addLine(TString::Format(
" %+10.3f %+10.3f %+10.3f %10.3f | Sum", sum.x(), sum.y(), sum.z(), TMath::Sqrt(sum.perp2())));
193 addLine(TString::Format(
"m = %10.3f", TMath::Sqrt(
TMath::Max(0.0, sum_len * sum_len - sum.mag2()))));
194 addLine(TString::Format(
"mT = %10.3f", TMath::Sqrt(
TMath::Max(0.0, sum_len_xy * sum_len_xy - sum.perp2()))));
195 addLine(TString::Format(
"HT = %10.3f", sum_len_xy));
198 addLine(TString::Format(
"deltaPhi = %+6.4f",
deltaPhi(first.Phi(), second.Phi())));
199 addLine(TString::Format(
"deltaEta = %+6.4f", first.Eta()- second.Eta()));
200 addLine(TString::Format(
"deltaR = % 6.4f",
deltaR(first.Eta(), first.Phi(), second.Eta(), second.Phi())));
const std::set< FWModelId > & selected() const
virtual double pz() const =0
z coordinate of momentum vector
double px() const
x coordinate of momentum vector
U second(std::pair< T, U > const &p)
virtual double py() const =0
y coordinate of momentum vector
void addLine(const TString &line)
double pz() const
z coordinate of momentum vector
double deltaR(double eta1, double eta2, double phi1, double phi2)
XYZVectorD XYZVector
spatial vector with cartesian internal representation
virtual double px() const =0
x coordinate of momentum vector
FWSelectionManager * m_selectionMgr
double py() const
y coordinate of momentum vector