Skip to content

Commit

Permalink
---
Browse files Browse the repository at this point in the history 
yaml
---
r: 62327
b: refs/heads/l1tmuon-upgrade-dev
c: 651d0a5
h: refs/heads/l1tmuon-upgrade-dev
i:
  62325: 55e38a4
  62323: 55cf25e
  62319: 0558075
  • Loading branch information
Sebastian Mirko Naumann-Emme committed Mar 19, 2009
1 parent bd62bc4 commit 773b77f
Show file tree
Hide file tree
Showing 2 changed files with 47 additions and 47 deletions.
2 changes: 1 addition & 1 deletion [refs]
View file
Original file line number Diff line number Diff line change
@@ -1,2 +1,2 @@
---
refs/heads/l1tmuon-upgrade-dev: ba65feb8faa5fb3650799c8a63d5569fa0a3c4c0
refs/heads/l1tmuon-upgrade-dev: 651d0a560a96e09d20e6b769db9046955d164a9d
92 changes: 46 additions & 46 deletions trunk/TopQuarkAnalysis/TopKinFitter/src/TtSemiLepKinFitter.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -205,50 +205,50 @@ int TtSemiLepKinFitter::fit(const std::vector<pat::Jet>& jets, const pat::Lepton
res::HelperJet jetRes;
switch(jetParam_){
case kEMom :
m1b(0,0) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(1,1) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(2,2) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(3,3) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m2b(0,0) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(1,1) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(2,2) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(3,3) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m3b(0,0) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(1,1) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(2,2) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(3,3) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m4b(0,0) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(1,1) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(2,2) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(3,3) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m1b(0,0) = pow(jetRes.a (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(1,1) = pow(jetRes.b (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(2,2) = pow(jetRes.c (q1pt, q1eta, res::HelperJet::kUds), 2);
m1b(3,3) = pow(jetRes.d (q1pt, q1eta, res::HelperJet::kUds), 2);
m2b(0,0) = pow(jetRes.a (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(1,1) = pow(jetRes.b (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(2,2) = pow(jetRes.c (q2pt, q2eta, res::HelperJet::kUds), 2);
m2b(3,3) = pow(jetRes.d (q2pt, q2eta, res::HelperJet::kUds), 2);
m3b(0,0) = pow(jetRes.a (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(1,1) = pow(jetRes.b (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(2,2) = pow(jetRes.c (b1pt, b1eta, res::HelperJet::kB ), 2);
m3b(3,3) = pow(jetRes.d (b1pt, b1eta, res::HelperJet::kB ), 2);
m4b(0,0) = pow(jetRes.a (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(1,1) = pow(jetRes.b (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(2,2) = pow(jetRes.c (b2pt, b2eta, res::HelperJet::kB ), 2);
m4b(3,3) = pow(jetRes.d (b2pt, b2eta, res::HelperJet::kB ), 2);
break;
case kEtEtaPhi :
m1 (0,0) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (0,0) = pow(jetRes.et (q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (1,1) = pow(jetRes.eta(q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (2,2) = pow(jetRes.phi(q1pt, q1eta, res::HelperJet::kUds), 2);
m2 (0,0) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (0,0) = pow(jetRes.et (q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (1,1) = pow(jetRes.eta(q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (2,2) = pow(jetRes.phi(q2pt, q2eta, res::HelperJet::kUds), 2);
m3 (0,0) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (0,0) = pow(jetRes.et (b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (1,1) = pow(jetRes.eta(b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (2,2) = pow(jetRes.phi(b1pt, b1eta, res::HelperJet::kB ), 2);
m4 (0,0) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (0,0) = pow(jetRes.et (b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (1,1) = pow(jetRes.eta(b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (2,2) = pow(jetRes.phi(b2pt, b2eta, res::HelperJet::kB ), 2);
break;
case kEtThetaPhi :
m1 (0,0) = pow(jetRes.pt (q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (1,1) = pow(jetRes.eta(q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (2,2) = pow(jetRes.phi(q1pt, q1eta, res::HelperJet::kUds), 2);
m2 (0,0) = pow(jetRes.pt (q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (1,1) = pow(jetRes.eta(q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (2,2) = pow(jetRes.phi(q2pt, q2eta, res::HelperJet::kUds), 2);
m3 (0,0) = pow(jetRes.pt (b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (1,1) = pow(jetRes.eta(b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (2,2) = pow(jetRes.phi(b1pt, b1eta, res::HelperJet::kB ), 2);
m4 (0,0) = pow(jetRes.pt (b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (1,1) = pow(jetRes.eta(b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (2,2) = pow(jetRes.phi(b2pt, b2eta, res::HelperJet::kB ), 2);
m1 (0,0) = pow(jetRes.et (q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (1,1) = pow(jetRes.theta(q1pt, q1eta, res::HelperJet::kUds), 2);
m1 (2,2) = pow(jetRes.phi (q1pt, q1eta, res::HelperJet::kUds), 2);
m2 (0,0) = pow(jetRes.et (q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (1,1) = pow(jetRes.theta(q2pt, q2eta, res::HelperJet::kUds), 2);
m2 (2,2) = pow(jetRes.phi (q2pt, q2eta, res::HelperJet::kUds), 2);
m3 (0,0) = pow(jetRes.et (b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (1,1) = pow(jetRes.theta(b1pt, b1eta, res::HelperJet::kB ), 2);
m3 (2,2) = pow(jetRes.phi (b1pt, b1eta, res::HelperJet::kB ), 2);
m4 (0,0) = pow(jetRes.et (b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (1,1) = pow(jetRes.theta(b2pt, b2eta, res::HelperJet::kB ), 2);
m4 (2,2) = pow(jetRes.phi (b2pt, b2eta, res::HelperJet::kB ), 2);
break;
}
}
Expand All @@ -263,43 +263,43 @@ int TtSemiLepKinFitter::fit(const std::vector<pat::Jet>& jets, const pat::Lepton
res::HelperElectron elecRes;
switch(lepParam_){
case kEMom :
m5(0,0) = pow(elecRes.pt (pt, eta), 2);
m5(1,1) = pow(elecRes.pt (pt, eta), 2);
m5(2,2) = pow(elecRes.pt (pt, eta), 2);
m5(0,0) = pow(elecRes.a (pt, eta), 2);
m5(1,1) = pow(elecRes.b (pt, eta), 2);
m5(2,2) = pow(elecRes.c (pt, eta), 2);
break;
case kEtEtaPhi :
m5(0,0) = pow(elecRes.pt (pt, eta), 2);
m5(0,0) = pow(elecRes.et (pt, eta), 2);
m5(1,1) = pow(elecRes.eta(pt, eta), 2);
m5(2,2) = pow(elecRes.phi(pt, eta), 2);
break;
case kEtThetaPhi :
m5(0,0) = pow(elecRes.pt (pt, eta), 2);
m5(1,1) = pow(elecRes.eta(pt, eta), 2);
m5(0,0) = pow(elecRes.et (pt, eta), 2);
m5(1,1) = pow(elecRes.theta(pt, eta), 2);
m5(2,2) = pow(elecRes.phi(pt, eta), 2);
break;
}
}
// add neutrino resolutions
{
//FIXME this dirty hack needs a clean solution soon!
double met = neutrino.pt();
double pt = neutrino.pt();

res::HelperMET metRes;
switch(metParam_){
case kEMom :
m6(0,0) = pow(metRes.met(met), 2);
m6(1,1) = pow( 9999., 2);
m6(2,2) = pow(metRes.phi(met), 2);
m6(0,0) = pow(metRes.a(pt), 2);
m6(1,1) = pow(metRes.b(pt), 2);
m6(2,2) = pow(metRes.c(pt), 2);
break;
case kEtEtaPhi :
m6(0,0) = pow(metRes.met(met), 2);
m6(0,0) = pow(metRes.et(pt), 2);
m6(1,1) = pow( 9999., 2);
m6(2,2) = pow(metRes.phi(met), 2);
m6(2,2) = pow(metRes.phi(pt), 2);
break;
case kEtThetaPhi :
m6(0,0) = pow(metRes.met(met), 2);
m6(0,0) = pow(metRes.et(pt), 2);
m6(1,1) = pow( 9999., 2);
m6(2,2) = pow(metRes.phi(met), 2);
m6(2,2) = pow(metRes.phi(pt), 2);
break;
}
}
Expand Down

0 comments on commit 773b77f

Please sign in to comment.