Changeset 319 in svn for trunk/src

Timestamp:

Mar 10, 2009, 11:59:08 PM (16 years ago)

Author:

Xavier Rouby

Message:

from TRootGenParticle to GenParticle in B-tag

File:

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

trunk/src/SmearUtil.cc

@@ -1034 +1034 @@
   for(unsigned int i=0; i < tracks.size(); i++) { 
       if(tracks[i].PT < pt_track ) continue;
+      //float dr = DeltaR(phi,eta,tracks[i].PhiOuter,tracks[i].EtaOuter);
       float dr = DeltaR(phi,eta,tracks[i].Phi,tracks[i].Eta);
       if (dr > TAU_track_scone) continue; 
@@ -1041 +1042 @@
   return numbtrack;
 }
-/*
-unsigned int RESOLution::NumTracks(unsigned int& NumTracksR, float& charge, const vector<TRootTracks> &tracks, const float pt_track, const float eta, const float phi, const float R) {
-  unsigned int numbtrack=0; // number of track in the tau-jet cone, which is smaller than R;
-  charge=0;      NumTracksR=0; // total number of tracks
-  if (R<TAU_track_scone) cout << "Warning: jet radius smaller than tau-jet radius in RESOLution::NumTracks\n";
-  for(unsigned int i=0; i < tracks.size(); i++) {
-      if(tracks[i].PT < pt_track ) continue;
-      float dr = DeltaR(phi,eta,tracks[i].Phi,tracks[i].Eta);
-      if (dr > R) continue;
-
-      NumTracksR++;
-      if(dr < TAU_track_scone) {// previously R==TAU_track_scone, for tau-jets only
-        numbtrack++;
-        charge += tracks[i].Charge; // total charge in the cone for Tau-jet
-      }
-
-  }
-  return numbtrack;
-}
-*/
-
 
 //*** Returns the PID of the particle with the highest energy, in a cone with a radius CONERADIUS and an axis (eta,phi)  *********
 //used by Btaggedjet
 ///// Attention : bug removed => CONERADIUS/2 -> CONERADIUS !!
-int RESOLution::Bjets(const TSimpleArray<TRootGenParticle> &subarray, const float eta, const float phi) {
+int RESOLution::Bjets(const TSimpleArray<GenParticle> &subarray, const float& eta, const float& phi) {
   float emax=0;
   int Ppid=0;
@@ -1083 +1063 @@
 
 //******************** Simulates the b-tagging efficiency for real bjet, or the misendentification for other jets****************
-bool RESOLution::Btaggedjet(const TLorentzVector &JET, const TSimpleArray<TRootGenParticle> &subarray)  {
+bool RESOLution::Btaggedjet(const TLorentzVector &JET, const TSimpleArray<GenParticle> &subarray) {
   if(      rand()%100 < (BTAG_b+1)    && Bjets(subarray,JET.Eta(),JET.Phi())==pB ) return true; // b-tag of b-jets is 40%
   else if( rand()%100 < (BTAG_mistag_c+1) && Bjets(subarray,JET.Eta(),JET.Phi())==pC ) return true; // b-tag of c-jets is 10%
@@ -1092 +1072 @@
 //***********************Isolation criteria***********************
 //****************************************************************
-//bool RESOLution::Isolation(const float phi, const float eta,const vector<TLorentzVector> &tracks, const float pt_second_track)
-bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float pt_second_track, const float isolCone)
+bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float& pt_second_track, const float& isolCone)
 {
    bool isolated = false;
@@ -1121 +1100 @@
                 iEta = (eta>0) ? TOWER_eta_edges[i-1] : -TOWER_eta_edges[i];
                 index = i-1;
-                //cout << setw(15) << left << eta << "\t" << iEta << endl;
                 break;
         }