Changeset 1345 in svn for trunk/modules/FastJetFinder.cc
- Timestamp:
- Dec 21, 2013, 3:00:11 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/modules/FastJetFinder.cc
r1335 r1345 193 193 TLorentzVector momentum; 194 194 Double_t deta, dphi, detaMax, dphiMax; 195 Double_t time, weightTime, avTime; 195 196 Int_t number; 196 197 Double_t rho = 0; … … 259 260 candidate = factory->NewCandidate(); 260 261 262 time=0; 263 weightTime=0; 264 261 265 inputList.clear(); 262 266 inputList = sequence->constituents(*itOutputList); … … 269 273 if(deta > detaMax) detaMax = deta; 270 274 if(dphi > dphiMax) dphiMax = dphi; 271 275 276 time += TMath::Sqrt(constituent->Momentum.E())*(constituent->Position.T()); 277 weightTime += TMath::Sqrt(constituent->Momentum.E()); 278 272 279 candidate->AddCandidate(constituent); 273 280 } 281 282 avTime = time/weightTime; 274 283 275 284 candidate->Momentum = momentum; 285 candidate->Position.SetT(avTime); 276 286 candidate->Area.SetPxPyPzE(area.px(), area.py(), area.pz(), area.E()); 277 287
Note:
See TracChangeset
for help on using the changeset viewer.