Changeset 38 in svn for trunk/Delphes.cpp

Timestamp:

Nov 18, 2008, 8:54:30 AM (16 years ago)

Author:

severine ovyn

Message:

remove MY bug jets

File:

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

trunk/Delphes.cpp

@@ -164 +164 @@
   
   TLorentzVector genMomentum(0,0,0,0);
+  TLorentzVector genMomentumCalo(0,0,0,0);
   LorentzVector jetMomentum;
   vector<TLorentzVector> TrackCentral;  
@@ -339 +340 @@
               PTmis = PTmis + genMomentum;//ptmis
               if(pid !=pMU) {
-                
-                towers.push_back(PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())));
+                PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
+                towers.push_back(CaloTower);
                 // create a fastjet::PseudoJet with these components and put it onto
                 // back of the input_particles vector
                 input_particles.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())); 
+                genMomentumCalo.SetPxPyPzE(CaloTower.fourVector.px,CaloTower.fourVector.py,CaloTower.fourVector.pz,CaloTower.fourVector.E);
                 elementCalo = (TRootCalo*) branchCalo->NewEntry();
-                elementCalo->Set(genMomentum);
+                elementCalo->Set(genMomentumCalo);
               }
             }
@@ -492 +494 @@
         }
       for (unsigned int i = 0; i < sorted_jets.size(); i++) {
-        elementJet = (TRootJet*) branchJet->NewEntry();
         TLorentzVector JET;
         JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
@@ -510 +511 @@
         if(JET.Pt() > DET->JET_pt)
           {
+            elementJet = (TRootJet*) branchJet->NewEntry();
             elementJet->Set(JET);
             // b-jets