Changeset 40 in svn
- Timestamp:
- Nov 18, 2008, 10:31:32 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Delphes.cpp
r38 r40 338 338 // for calorimetric towers and mission PT 339 339 if(genMomentum.E()!=0) { 340 PTmis = PTmis + genMomentum;//ptmis341 340 if(pid !=pMU) { 342 341 PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())); … … 345 344 // back of the input_particles vector 346 345 input_particles.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())); 346 //genMomentumCalo.SetPtEtaPhiE(CaloTower.Et(),CaloTower.iEta(),CaloTower.iPhi(),CaloTower.fourVector.E); 347 347 genMomentumCalo.SetPxPyPzE(CaloTower.fourVector.px,CaloTower.fourVector.py,CaloTower.fourVector.pz,CaloTower.fourVector.E); 348 348 elementCalo = (TRootCalo*) branchCalo->NewEntry(); … … 472 472 473 473 // computes the Missing Transverse Momentum 474 TLorentzVector Att(0.,0.,0.,0.); 475 for(unsigned int i=0; i < towers.size(); i++) 476 { 477 Att.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E); 478 PTmis = PTmis + Att; 479 } 474 480 elementEtmis = (TRootETmis*) branchETmis->NewEntry(); 475 elementEtmis->ET = ( -PTmis).Pt();481 elementEtmis->ET = (PTmis).Pt(); 476 482 elementEtmis->Phi = (-PTmis).Phi(); 477 483 elementEtmis->Px = (-PTmis).Px();
Note:
See TracChangeset
for help on using the changeset viewer.