Changeset 551 in svn for trunk/src/SmearUtil.cc

Timestamp:

Feb 23, 2010, 2:20:23 AM (15 years ago)

Author:

Xavier Rouby

Message:

EtRatio for muon and electron; apparently bug free

File:

• trunk/src/SmearUtil.cc (modified) (5 diffs)

trunk/src/SmearUtil.cc

@@ -1445 +1445 @@
   // iTower is the index of the tower, in [0, n_tower]. iTower points only towers in positive range
 
-  // etrat, which is a percentage between 00 and 99. It is the ratio of the transverse energy 
-  // in a N×N grid surrounding the muon to the pT of the muon. For well-isolated muons, both ptiso and etrat will be small.
-  // available parameters: ISOL_Calo_ET ,  ISOL_Calo_Grid
-
    if(ISOL_Calo_ET>1E10) return UNDEFINED; // avoid doing anything unreasonable...
    float et_sum=0;
@@ -1457 +1453 @@
                       // !! TOWER_eta_edges is only with eta>0
 
-cout << "iEta= " << iEta << "\tiPhi= " << iPhi << endl;
-   //cout << "la particule est dans la " << index << "eme cellule, entre " << TOWER_eta_edges[index] << " et " << TOWER_eta_edges[index+1] << endl;
-   //float centre = (iEta>0) ? ((TOWER_eta_edges[index]+TOWER_eta_edges[index+1])/2.0) : -((TOWER_eta_edges[index]+TOWER_eta_edges[index+1])/2.0);
-   //cout << "le centre de cette tour est " << centre << endl; 
-
+   //cout << "iEta= " << iEta << "\tiPhi= " << iPhi << endl;
    if(index != iUNDEFINED) {
       int range = (int) (N-1)/2;           // the (int) conversion is needed, as N is "unsigned int"
@@ -1467 +1459 @@
       // goal : to identify the center of each cell in the NxN grid.
       // the eta/phi coordinates of each center will be in (eta_ith_tower,phi_ith_tower)
-
          unsigned int i_eta = (index+i>=0) ? index + i : (int) -1*(index + i +1); // i_eta is the index in [0; index_eta_max]
-         //cout << i_eta << "\t";
+         if(i_eta>=TOWER_number) continue; // avoid going too far: TOWER_dphi[TOWER_number], TOWER_eta_edges[TOWER_number+1]
+         //cout << "i_eta" << i_eta << "\t TOWER_number = " << TOWER_number << endl;
+
          float dphi = TOWER_dphi[i_eta]*pi/180.; // in rad // size in phi of the cells for this eta
          float eta_ith_tower = (TOWER_eta_edges[i_eta]+TOWER_eta_edges[i_eta+1])/2.0;
          if(iEta<0) eta_ith_tower *= -1;     // needed if the central tower is in negative region
-         if(index+i<0) eta_ith_tower *= -1;  // needed if the central part is crossed by the grid
+         if(index+i<0) eta_ith_tower *= -1;  // needed if the "eta=0"-axis is crossed by the grid
          //cout << "pour eta_ith_tower=" << eta_ith_tower << ", on va dphi = " << dphi << endl;
          // !!! at this point, eta_ith_tower is the value in eta of the center of the current calo cell
@@ -1484 +1477 @@
             
             D_CaloTower calMuon(towers.getElement(eta_ith_tower,phi_ith_tower));
-            cout << "eta_ith_tower= "   << eta_ith_tower     << "\tphi_ith_tower= " << phi_ith_tower <<"\t"
-                 << "calMuon.getEta= " << calMuon.getEta() << "\tcalMuon.getPhi()= " << calMuon.getPhi() <<"\t";
-
-            if(calMuon.getEta() != UNDEFINED )
-             if(calMuon.getET() > ISOL_Calo_ET) {
-              et_sum += calMuon.getET();
-              //cout << "eta/phi/Et = " << eta_ith_tower << "\t" << phi_ith_tower << "\t" << calMuon.getET() << " GeV" << endl;
-              cout << calMuon.getET() << " GeV";
+            //cout << "eta_ith_tower= "   << eta_ith_tower     << "\tphi_ith_tower= " << phi_ith_tower <<"\t"
+            //     << "calMuon.getEta= " << calMuon.getEta() << "\tcalMuon.getPhi()= " << calMuon.getPhi() <<"\t";
+
+            if(calMuon.getEta() != UNDEFINED && calMuon.getET() > ISOL_Calo_ET) {
+               et_sum += calMuon.getET();
+               //cout << calMuon.getET() << " GeV";
             }
-            //else cout << "eta/phi = " << eta_ith_tower << "\t" << phi_ith_tower << "\tnot active\n";
-            else cout << phi_ith_tower << "\tnot active";
-            cout << endl;
-         }
-      }
+            //else cout << " not active";
+         } // j-loop (phi indices)
+      } // i-loop (eta indices)
    }  // undefined index
    else {
@@ -1503 +1492 @@
      if (CEN_max_mu < CEN_max_calo_fwd) 
         cout << "**    ERROR in RESOLution::CaloIsolation: 'muon'-tower not found!  **" << endl; 
-   }
-        // should never happen ! this would be a bug
-  //cout << "etrat = " << et_sum << "\t Pt=" << part.Pt() << endl;
-  
-  // should return a number between 0 and 99 (due to LHCO definitions)
-  // which is Pt(muon) / sum(ET)
-  float etrat = 0.;
-/*  if(et_sum==0) etrat = 99.;
-  //else if(et_sum>0) etrat = 100*part.Pt()/et_sum;
-  if(etrat<0) cout << "Error: negative etrat in CaloIsolation (" << etrat <<")\n";
-  //else if(etrat>99) cout << "Error: etrat should be in [0;99] in CaloIsolation (" << etrat <<")\n";
-  if(etrat>99) etrat = 99;
-*/
-  etrat = et_sum/part.Pt();
-  return etrat; 
+   }      // should never happen ! this would be a bug
+  
+  return (part.Pt()==0)? UNDEFINED: et_sum/part.Pt();
 }