Changeset 33 in svn for trunk/Delphes.cpp

Timestamp:

Nov 17, 2008, 9:01:32 PM (16 years ago)

Author:

severine ovyn

Message:

add offline Thresholds

File:

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

trunk/Delphes.cpp

@@ -339 +339 @@
               PTmis = PTmis + genMomentum;//ptmis
               if(pid !=pMU) {
+                
                 towers.push_back(PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())));
                 // create a fastjet::PseudoJet with these components and put it onto
@@ -449 +450 @@
       
       for(unsigned int i=0; i < electron.size(); i++) {
-        if(electron[i].E()!=0 && fabs(electron[i].Eta()) < DET->MAX_TRACKER)
+        if(electron[i].E()!=0 && fabs(electron[i].Eta()) < DET->MAX_TRACKER && electron[i].Pt() > DET->ELEC_pt)
           {
             elementElec = (TRootElectron*) branchElectron->NewEntry();
@@ -458 +459 @@
       }
       for(unsigned int i=0; i < muon.size(); i++) {
-        if(muon[i].E()!=0 && fabs(muon[i].Eta()) < DET->MAX_MU) 
+        if(muon[i].E()!=0 && fabs(muon[i].Eta()) < DET->MAX_MU && muon[i].Pt() > DET->MUON_pt) 
           {
             elementMu = (TRootMuon*) branchMuon->NewEntry();
@@ -494 +495 @@
         TLorentzVector JET;
         JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
-        //cout<<"Jet.Pt() "<<JET.Pt()<<endl;
-        elementJet->Set(JET);
-        // b-jets
-        bool btag=false;
-        if((fabs(JET.Eta()) < DET->MAX_TRACKER && DET->Btaggedjet(JET, NFCentralQ)))btag=true;
-        elementJet->Btag = btag;
-        
-        // Tau jet identification : 1! track and electromagnetic collimation
-        if(fabs(JET.Eta()) < (DET->MAX_TRACKER - DET->TAU_CONE_TRACKS)) {
-          double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi());
-          if(
-             ( Energie_tau_central/JET.E() > DET->TAU_EM_COLLIMATION ) &&
-             ( DET->NumTracks(TrackCentral,DET->PT_TRACK_TAU,JET.Eta(),JET.Phi()) == 1 )
-             ) {
-            elementTauJet = (TRootTauJet*) branchTauJet->NewEntry();
-            elementTauJet->Set(JET);
-          } // if tau jet
-        } // if JET.eta < tracker - tau_cone : Tau jet identification
+        // Tau jet identification : 1! track and electromagnetic collimation
+        if(fabs(JET.Eta()) < (DET->MAX_TRACKER - DET->TAU_CONE_TRACKS)) {
+        double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi());
+        if(
+           ( Energie_tau_central/JET.E() > DET->TAU_EM_COLLIMATION ) &&
+           ( DET->NumTracks(TrackCentral,DET->PT_TRACK_TAU,JET.Eta(),JET.Phi()) == 1 ) &&
+           ( JET.Pt() > DET->TAUJET_pt)
+          ) {
+               elementTauJet = (TRootTauJet*) branchTauJet->NewEntry();
+               elementTauJet->Set(JET);
+            } // if tau jet
+           } // if JET.eta < tracker - tau_cone : Tau jet identification
+  
+        if(JET.Pt() > DET->JET_pt)
+          {
+            elementJet->Set(JET);
+            // b-jets
+            bool btag=false;
+            if((fabs(JET.Eta()) < DET->MAX_TRACKER && DET->Btaggedjet(JET, NFCentralQ)))btag=true;
+            elementJet->Btag = btag;
+          }
       } // for itJet : loop on all jets