Fork me on GitHub

Changeset 15 in svn for trunk/Delphes.cpp


Ignore:
Timestamp:
Nov 6, 2008, 4:44:48 PM (16 years ago)
Author:
severine ovyn
Message:

remove PT sorted bug

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Delphes.cpp

    r11 r15  
    216216      input_particles.clear();
    217217      TSimpleArray<TRootGenParticle> NFCentralQ;
    218 
     218      inclusive_jets.clear();
     219      sorted_jets.clear();
    219220
    220221      // Loop over all particles in event
     
    471472         
    472473          // sort jets into increasing pt
    473           sorted_by_pt(inclusive_jets);
     474          sorted_jets = sorted_by_pt(inclusive_jets);
    474475        }
    475      
    476       for (unsigned int i = 0; i < inclusive_jets.size(); i++) {
     476      cout<<"**********"<<endl;
     477      for (unsigned int i = 0; i < sorted_jets.size(); i++) {
    477478        elementJet = (TRootJet*) branchJet->NewEntry();
    478479        TLorentzVector JET;
    479         JET.SetPxPyPzE(inclusive_jets[i].px(),inclusive_jets[i].py(),inclusive_jets[i].pz(),inclusive_jets[i].E());
     480        JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
     481        cout<<"JET.Pt() "<<JET.Pt()<<endl;
    480482        elementJet->Set(JET);
    481483        // b-jets
Note: See TracChangeset for help on using the changeset viewer.