Changeset 40 in svn for trunk

Timestamp:

Nov 18, 2008, 10:31:32 AM (16 years ago)

Author:

severine ovyn

Message:

ETmis using CaloTowers

File:

• trunk/Delphes.cpp (modified) (3 diffs)

trunk/Delphes.cpp

@@ -338 +338 @@
             // for calorimetric towers and mission PT
             if(genMomentum.E()!=0) {
-              PTmis = PTmis + genMomentum;//ptmis
               if(pid !=pMU) {
                 PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
@@ -345 +344 @@
                 // back of the input_particles vector
                 input_particles.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())); 
+                //genMomentumCalo.SetPtEtaPhiE(CaloTower.Et(),CaloTower.iEta(),CaloTower.iPhi(),CaloTower.fourVector.E);
                 genMomentumCalo.SetPxPyPzE(CaloTower.fourVector.px,CaloTower.fourVector.py,CaloTower.fourVector.pz,CaloTower.fourVector.E);
                 elementCalo = (TRootCalo*) branchCalo->NewEntry();
@@ -472 +472 @@
       
       // computes the Missing Transverse Momentum
+      TLorentzVector Att(0.,0.,0.,0.);
+      for(unsigned int i=0; i < towers.size(); i++)
+        {
+          Att.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E);
+          PTmis = PTmis + Att;
+        }
       elementEtmis = (TRootETmis*) branchETmis->NewEntry();
-      elementEtmis->ET = (-PTmis).Pt();
+      elementEtmis->ET = (PTmis).Pt();
       elementEtmis->Phi = (-PTmis).Phi();
       elementEtmis->Px = (-PTmis).Px();