Changeset 547 in svn for trunk/src/JetsUtil.cc

Timestamp:

Feb 22, 2010, 11:35:46 AM (14 years ago)

Author:

severine ovyn

Message:

bug fix(jet Ntracks), 3-prong taus, vertex coordinates for tracks

File:

• trunk/src/JetsUtil.cc (modified) (2 diffs)

trunk/src/JetsUtil.cc

@@ -165 +165 @@
             for (unsigned int j = 0; j < constituents.size(); j++)
                 {
-                  D_CaloTower calConsti(list_of_active_towers.getElement(constituents[j].eta(),constituents[j].phi()));
+                  //D_CaloTower calConsti(list_of_active_towers.getElement(constituents[j].eta(),constituents[j].phi()));
+                  // bug fix! bad association of phi [-pi;pi] <-> [0 ; 2pi]
+                  D_CaloTower calConsti(list_of_active_towers.getElement(constituents[j].eta(),constituents[j].phi_std()));
                   EmVal  += calConsti.getEem();
                   HadVal += calConsti.getEhad();
@@ -224 +226 @@
   for (unsigned int i = 0; i < sorted_jets.size(); i++) {
     JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
-    // Tau jet identification : 1! track and electromagnetic collimation
+    // Tau jet identification : 1! or 3! track and electromagnetic collimation
     if(fabs(JET.Eta()) < (DET->CEN_max_tracker - DET->TAU_track_scone)) {
       double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi());
+      int NumTrackTau = DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi());
       if(
          ( Energie_tau_central/JET.E() > DET->TAU_energy_frac ) &&
-         ( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi()) == 1 ) &&
+         (( NumTrackTau == 1 ) || ( NumTrackTau == 3 )) &&
          ( JET.Pt() > DET->PTCUT_taujet)
          ) {