Changes in modules/MomentumSmearing.cc [84a097e:341014c] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/MomentumSmearing.cc
r84a097e r341014c 95 95 { 96 96 Candidate *candidate, *mother; 97 Double_t pt, eta, phi, e, m,res;97 Double_t pt, eta, phi, e, res; 98 98 99 99 fItInputArray->Reset(); … … 106 106 pt = candidateMomentum.Pt(); 107 107 e = candidateMomentum.E(); 108 m = candidateMomentum.M();108 res = fFormula->Eval(pt, eta, phi, e); 109 109 110 res = fFormula->Eval(pt, eta, phi, e); 110 // apply smearing formula 111 //pt = gRandom->Gaus(pt, fFormula->Eval(pt, eta, phi, e) * pt); 112 111 113 res = (res > 1.0) ? 1.0 : res; 112 114 113 115 pt = LogNormal(pt, res * pt); 116 117 //if(pt <= 0.0) continue; 114 118 115 119 mother = candidate; … … 117 121 eta = candidateMomentum.Eta(); 118 122 phi = candidateMomentum.Phi(); 119 candidate->Momentum.SetPtEtaPhiM(pt, eta, phi, m); 120 candidate->ErrorPT = res*pt; 121 candidate->PT = pt; 122 123 candidate->Momentum.SetPtEtaPhiE(pt, eta, phi, pt * TMath::CosH(eta)); 124 //candidate->TrackResolution = fFormula->Eval(pt, eta, phi, e); 123 125 candidate->TrackResolution = res; 124 126 candidate->AddCandidate(mother);
Note:
See TracChangeset
for help on using the changeset viewer.