[ViewVC] Diff of: cvsroot/UserCode/MitHzz4l/LeptonSelection/src/IsolationSelection.cc

Comparing UserCode/MitHzz4l/LeptonSelection/src/IsolationSelection.cc (file contents):
Revision 1.14 by khahn, Thu May 10 11:14:33 2012 UTC vs.
Revision 1.16 by khahn, Thu May 10 22:53:21 2012 UTC

+mithep::ElectronIDMVA * eleIsoMVA;
+mithep::ElectronTools   eleT;
-+
+// global hack to sync
-+
+double gChargedIso;
-+
+double gGammaIso;
-+
+double gNeutralIso;
-+
-+
+//--------------------------------------------------------------------------------------------------
+Float_t computePFMuonIso(const mithep::Muon *muon,
+                         const mithep::Vertex & vtx,
+          IsLeptonFootprint = kTRUE;
+      } // loop over electrons
-+
+      // KH, comment to sync
-+
+      /*
-+
+      //
-+
+      // Check for muons
-+
+      //
-+
+      for (Int_t q=0; q < muonsToVeto.size(); ++q) {
-+
+        const mithep::Muon *tmpmu = muonsToVeto[q];
-+
+        // 4l muon
-+
+        if( pf->HasTrackerTrk() ) {
-+
+          if( pf->TrackerTrk() == tmpmu->TrackerTrk() )
-+
+            IsLeptonFootprint = kTRUE;
-+
+        }
-+
+        // PF charged
-+
+        if (pf->Charge() != 0 && mithep::MathUtils::DeltaR(tmpmu->Phi(),tmpmu->Eta(), pf->Phi(), pf->Eta()) < 0.01)
-+
+          IsLeptonFootprint = kTRUE;
-+
+      } // loop over muons
-+
+      */
-+
-+
+    if (IsLeptonFootprint)
-+
+      continue;
-+
-+
+    //
-+
+    // Charged Iso
-+
+    //
-+
+    if (pf->Charge() != 0 ) {
-+
-+
+      //if( dr < 0.01 ) continue; // only for muon iso mva?
-+
+      if (abs(pf->PFType()) == PFCandidate::eElectron || abs(pf->PFType()) == PFCandidate::eMuon) continue;
-+
-+
+      if( pf->HasTrackerTrk() ) {
-+
+        if (abs(pf->TrackerTrk()->DzCorrected(vtx)) > 0.2) continue;
-+
+        if( ctrl.debug ) cout << "charged:: " << pf->PFType() << " " << pf->Pt() << " "
-+
+                              << abs(pf->TrackerTrk()->DzCorrected(vtx)) << " "
-+
+                              << dr << endl;
-+
+      }
-+
+      if( pf->HasGsfTrk() ) {
-+
+        if (abs(pf->GsfTrk()->DzCorrected(vtx)) > 0.2) continue;
-+
+        if( ctrl.debug ) cout << "charged:: " << pf->PFType() << " " << pf->Pt() << " "
-+
+                              << abs(pf->GsfTrk()->DzCorrected(vtx)) << " "
-+
+                              << dr << endl;
-+
+      }
-+
-+
-+
+      fChargedIso += pf->Pt();
-+
+    }
-+
-+
+    //
-+
+    // Gamma Iso
-+
+    //
-+
+    else if (abs(pf->PFType()) == PFCandidate::eGamma) {
-+
+      // KH, add to sync
-+
+      if( pf->Pt() > 0.5 )
-+
+      fGammaIso += pf->Pt();
-+
+    }
-+
-+
+    //
-+
+    // Other Neutrals
-+
+    //
-+
+    else {
-+
+      // KH, add to sync
-+
+      if( pf->Pt() > 0.5 )
-+
+        fNeutralHadronIso += pf->Pt();
-+
+    }
-+
-+
+    }
-+
-+
+  }
-+
-+
+  double rho = 0;
-+
+//   if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
-+
+//     rho = fPUEnergyDensity->At(0)->Rho();
-+
+  if (!(isnan(fPUEnergyDensity->At(0)->RhoLowEta()) || isinf(fPUEnergyDensity->At(0)->RhoLowEta())))
-+
+    rho = fPUEnergyDensity->At(0)->RhoLowEta();
-+
-+
+  // WARNING!!!!
-+
+  // hardcode for sync ...
-+
+  EffectiveAreaVersion = muT.kMuEAData2011;
-+
+  // WARNING!!!!
-+
-+
-+
+  double pfIso = fChargedIso + max(0.0,(fGammaIso + fNeutralHadronIso
-+
+                                        -rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
-+
+                                                                   mu->Eta(),EffectiveAreaVersion)));
-+
-+
+  gChargedIso = fChargedIso;
-+
+  gGammaIso = fGammaIso;
-+
+  gNeutralIso = fNeutralHadronIso;
-+
+  return pfIso;
-+
+}
-+
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+// hacked version
-+
+double  muonPFIso04(ControlFlags &ctrl,
-+
+                    const mithep::Muon * mu,
-+
+                    const mithep::Vertex & vtx,
-+
+                    const mithep::Array<mithep::PFCandidate> * fPFCandidates,
-+
+                    float rho,
-+
+                    mithep::MuonTools::EMuonEffectiveAreaTarget EffectiveAreaVersion,
-+
+                    vector<const mithep::Muon*> muonsToVeto,
-+
+                    vector<const mithep::Electron*> electronsToVeto)
-+
+//--------------------------------------------------------------------------------------------------
-+
+{
-+
-+
+  extern double gChargedIso;
-+
+  extern double  gGammaIso;
-+
+  extern double  gNeutralIso;
-+
-+
+  if( ctrl.debug ) {
-+
+    cout << "muonIsoMVASelection :: muons to veto " << endl;
-+
+    for( int i=0; i<muonsToVeto.size(); i++ ) {
-+
+      const mithep::Muon * vmu = muonsToVeto[i];
-+
+      cout << "\tpt: " << vmu->Pt()
-+
+           << "\teta: " << vmu->Eta()
-+
+           << "\tphi: " << vmu->Phi()
-+
+           << endl;
-+
+    }
-+
+    cout << "muonIsoMVASelection :: electrson to veto " << endl;
-+
+    for( int i=0; i<electronsToVeto.size(); i++ ) {
-+
+      const mithep::Electron * vel = electronsToVeto[i];
-+
+      cout << "\tpt: " << vel->Pt()
-+
+           << "\teta: " << vel->Eta()
-+
+           << "\tphi: " << vel->Phi()
-+
+           << endl;
-+
+    }
-+
+  }
-+
-+
+  //
-+
+  // final iso
-+
+  //
-+
+  Double_t fChargedIso  = 0.0;
-+
+  Double_t fGammaIso  = 0.0;
-+
+  Double_t fNeutralHadronIso  = 0.0;
-+
-+
+  //
-+
+  //Loop over PF Candidates
-+
+  //
-+
+  for(int k=0; k<fPFCandidates->GetEntries(); ++k) {
-+
+    const mithep::PFCandidate *pf = (mithep::PFCandidate*)((*fPFCandidates)[k]);
-+
-+
+    Double_t deta = (mu->Eta() - pf->Eta());
-+
+    Double_t dphi = mithep::MathUtils::DeltaPhi(Double_t(mu->Phi()),Double_t(pf->Phi()));
-+
+    Double_t dr = mithep::MathUtils::DeltaR(mu->Phi(),mu->Eta(), pf->Phi(), pf->Eta());
-+
+    if (dr > 0.4) continue;
-+
-+
+    if (pf->HasTrackerTrk() && (pf->TrackerTrk() == mu->TrackerTrk()) ) continue;
-+
-+
+    //
-+
+    // Lepton Footprint Removal
-+
+    //
-+
+    Bool_t IsLeptonFootprint = kFALSE;
-+
+    if (dr < 1.0) {
-+
-+
+      //
-+
+      // Check for electrons
-+
+      //
-+
+      for (Int_t q=0; q < electronsToVeto.size(); ++q) {
-+
+        const mithep::Electron *tmpele = electronsToVeto[q];
-+
+        // 4l electron
-+
+        if( pf->HasTrackerTrk() ) {
-+
+          if( pf->TrackerTrk() == tmpele->TrackerTrk() )
-+
+            IsLeptonFootprint = kTRUE;
-+
+        }
-+
+        if( pf->HasGsfTrk() ) {
-+
+          if( pf->GsfTrk() == tmpele->GsfTrk() )
-+
+            IsLeptonFootprint = kTRUE;
-+
+        }
-+
+        // PF charged
-+
+        if (pf->Charge() != 0 && fabs(tmpele->SCluster()->Eta()) > 1.479
-+
+            && mithep::MathUtils::DeltaR(tmpele->Phi(),tmpele->Eta(), pf->Phi(), pf->Eta()) < 0.015)
-+
+          IsLeptonFootprint = kTRUE;
-+
+        // PF gamma
-+
+        if (abs(pf->PFType()) == PFCandidate::eGamma && fabs(tmpele->SCluster()->Eta()) > 1.479
-+
+            && mithep::MathUtils::DeltaR(tmpele->Phi(),tmpele->Eta(), pf->Phi(), pf->Eta()) < 0.08)
-+
+          IsLeptonFootprint = kTRUE;
-+
+      } // loop over electrons
-+
+      /* KH - comment for sync
+      //
+      // Check for muons
-<
+  double rho = 0;
-<
+  if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
-<
+    rho = fPUEnergyDensity->At(0)->Rho();
->
+//   double rho = 0;
->
+//   if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
->
+//     rho = fPUEnergyDensity->At(0)->Rho();
+  // WARNING!!!!
+  // hardcode for sync ...
+  double pfIso = fChargedIso + max(0.0,(fGammaIso + fNeutralHadronIso
+                                        -rho*muT.MuonEffectiveArea(muT.kMuGammaAndNeutralHadronIso04,
+                                                                   mu->Eta(),EffectiveAreaVersion)));
-+
+  gChargedIso = fChargedIso;
-+
+  gGammaIso   = fGammaIso;
-+
+  gNeutralIso = fNeutralHadronIso;
+  return pfIso;
-+
+//--------------------------------------------------------------------------------------------------
+SelectionStatus muonReferenceIsoSelection(ControlFlags &ctrl,
+                                          const mithep::Muon * mu,
+  SelectionStatus status;
-<
->
+  double pfIso = muonPFIso04( ctrl, mu, vtx, fPFCandidates, fPUEnergyDensity,
+                              EffectiveAreaVersion, muonsToVeto ,electronsToVeto );
-+
+  cout << "--------------> setting muon isoPF04 to" << pfIso << endl;
-+
+  status.isoPF04 = pfIso;
-+
+  status.chisoPF04 = gChargedIso;
-+
+  status.gaisoPF04 = gGammaIso;
-+
+  status.neisoPF04 = gNeutralIso;
-+
-+
+  bool pass = false;
-+
+  if( (pfIso/mu->Pt()) < MUON_REFERENCE_PFISO_CUT ) pass = true;
-+
-+
+  if( pass ) {
-+
+    status.orStatus(SelectionStatus::LOOSEISO);
-+
+    status.orStatus(SelectionStatus::TIGHTISO);
-+
+  }
-+
+  if(ctrl.debug) cout << "returning status : " << hex << status.getStatus() << dec << endl;
-+
+  return status;
-+
-+
+}
-+
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+// hacked version
-+
+SelectionStatus muonReferenceIsoSelection(ControlFlags &ctrl,
-+
+                                          const mithep::Muon * mu,
-+
+                                          const mithep::Vertex & vtx,
-+
+                                          const mithep::Array<mithep::PFCandidate> * fPFCandidates,
-+
+                                          float rho,
-+
+                                          mithep::MuonTools::EMuonEffectiveAreaTarget EffectiveAreaVersion,
-+
+                                          vector<const mithep::Muon*> muonsToVeto,
-+
+                                          vector<const mithep::Electron*> electronsToVeto)
-+
+//--------------------------------------------------------------------------------------------------
-+
+{
-+
-+
+  SelectionStatus status;
-+
-+
+  double pfIso = muonPFIso04( ctrl, mu, vtx, fPFCandidates, rho,
-+
+                              EffectiveAreaVersion, muonsToVeto ,electronsToVeto );
+  bool pass = false;
+  if( (pfIso/mu->Pt()) < MUON_REFERENCE_PFISO_CUT ) pass = true;
+      } // loop over electrons
-+
-+
+      //
-+
+      // Check for muons
-+
+      //
-+
+      for (Int_t q=0; q < muonsToVeto.size(); ++q) {
-+
+        const mithep::Muon *tmpmu = muonsToVeto[q];
-+
+        // 4l muon
-+
+        if( pf->HasTrackerTrk() ) {
-+
+          if (pf->TrackerTrk() == tmpmu->TrackerTrk() ){
-+
+            if( ctrl.debug) cout << "\tmatches 4L mu ..." << endl;
-+
+            IsLeptonFootprint = kTRUE;
-+
+          }
-+
+        }
-+
+        // PF charged
-+
+        if (pf->Charge() != 0 && mithep::MathUtils::DeltaR(tmpmu->Phi(),tmpmu->Eta(), pf->Phi(), pf->Eta()) < 0.01) {
-+
+          if( ctrl.debug) cout << "\tcharged trk, dR matches 4L mu ..." << endl;
-+
+          IsLeptonFootprint = kTRUE;
-+
+        }
-+
+      } // loop over muons
-+
-+
-+
+    if (IsLeptonFootprint)
-+
+      continue;
-+
-+
+    //
-+
+    // Charged Iso
-+
+    //
-+
+    if (pf->Charge() != 0 ) {
-+
-+
+      if( pf->HasTrackerTrk() )
-+
+        if (abs(pf->TrackerTrk()->DzCorrected(vtx)) > 0.2) continue;
-+
+      if( pf->HasGsfTrk() )
-+
+        if (abs(pf->GsfTrk()->DzCorrected(vtx)) > 0.2) continue;
-+
-+
+      // Veto any PFmuon, or PFEle
-+
+      if (abs(pf->PFType()) == PFCandidate::eElectron || abs(pf->PFType()) == PFCandidate::eMuon) continue;
-+
-+
+      // Footprint Veto
-+
+      if (fabs(ele->SCluster()->Eta()) > 1.479 && dr < 0.015) continue;
-+
-+
+      if( ctrl.debug) cout << "charged:: pt: " << pf->Pt()
-+
+                           << "\ttype: " << pf->PFType()
-+
+                           << "\ttrk: " << pf->TrackerTrk() << endl;
-+
-+
+      fChargedIso += pf->Pt();
-+
+    }
-+
-+
+    //
-+
+    // Gamma Iso
-+
+    //
-+
+    else if (abs(pf->PFType()) == PFCandidate::eGamma) {
-+
-+
+      if (fabs(ele->SCluster()->Eta()) > 1.479) {
-+
+        if (mithep::MathUtils::DeltaR(ele->Phi(),ele->Eta(), pf->Phi(), pf->Eta()) < 0.08) continue;
-+
+      }
-+
+      if( ctrl.debug) cout << "gamma:: " << pf->Pt() << " "
-+
+                           << dr << endl;
-+
+      fGammaIso += pf->Pt();
-+
+    }
-+
-+
+    //
-+
+    // Neutral Iso
-+
+    //
-+
+    else {
-+
+      if( ctrl.debug) cout << "neutral:: " << pf->Pt() << " "
-+
+                           << dr << endl;
-+
+        fNeutralHadronIso += pf->Pt();
-+
+    }
-+
-+
+    }
-+
-+
+  }
-+
-+
+  double rho = 0;
-+
+//   if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
-+
+//     rho = fPUEnergyDensity->At(0)->Rho();
-+
+  if (!(isnan(fPUEnergyDensity->At(0)->RhoLowEta()) || isinf(fPUEnergyDensity->At(0)->RhoLowEta())))
-+
+    rho = fPUEnergyDensity->At(0)->RhoLowEta();
-+
-+
+  // WARNING!!!!
-+
+  // hardcode for sync ...
-+
+  EffectiveAreaVersion = eleT.kEleEAData2011;
-+
+  // WARNING!!!!
-+
-+
-+
+  double pfIso = fChargedIso +
-+
+    max(0.0,fGammaIso -rho*eleT.ElectronEffectiveArea(eleT.kEleGammaIso04,
-+
+                                                ele->Eta(),EffectiveAreaVersion)) +
-+
+    max(0.0,fNeutralHadronIso -rho*eleT.ElectronEffectiveArea(eleT.kEleNeutralHadronIso04,
-+
+                                                        ele->Eta(),EffectiveAreaVersion)) ;
-+
-+
+  gChargedIso = fChargedIso;
-+
+  gGammaIso = fGammaIso;
-+
+  gNeutralIso = fNeutralHadronIso;
-+
-+
+  return pfIso;
-+
+}
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+// hacked version
-+
+float electronPFIso04(ControlFlags &ctrl,
-+
+                      const mithep::Electron * ele,
-+
+                      const mithep::Vertex & vtx,
-+
+                      const mithep::Array<mithep::PFCandidate> * fPFCandidates,
-+
+                      float rho,
-+
+                      mithep::ElectronTools::EElectronEffectiveAreaTarget EffectiveAreaVersion,
-+
+                      vector<const mithep::Muon*> muonsToVeto,
-+
+                      vector<const mithep::Electron*> electronsToVeto)
-+
+//--------------------------------------------------------------------------------------------------
-+
+{
-+
-+
+  if( ctrl.debug ) {
-+
+    cout << "electronIsoMVASelection :: muons to veto " << endl;
-+
+    for( int i=0; i<muonsToVeto.size(); i++ ) {
-+
+      const mithep::Muon * vmu = muonsToVeto[i];
-+
+      cout << "\tpt: " << vmu->Pt()
-+
+           << "\teta: " << vmu->Eta()
-+
+           << "\tphi: " << vmu->Phi()
-+
+           << endl;
-+
+    }
-+
+    cout << "electronIsoMVASelection :: electrson to veto " << endl;
-+
+    for( int i=0; i<electronsToVeto.size(); i++ ) {
-+
+      const mithep::Electron * vel = electronsToVeto[i];
-+
+      cout << "\tpt: " << vel->Pt()
-+
+           << "\teta: " << vel->Eta()
-+
+           << "\tphi: " << vel->Phi()
-+
+           << "\ttrk: " << vel->TrackerTrk()
-+
+           << endl;
-+
+    }
-+
+  }
-+
-+
-+
+  //
-+
+  // final iso
-+
+  //
-+
+  Double_t fChargedIso = 0.0;
-+
+  Double_t fGammaIso = 0.0;
-+
+  Double_t fNeutralHadronIso = 0.0;
-+
-+
-+
+  //
-+
+  //Loop over PF Candidates
-+
+  //
-+
+  for(int k=0; k<fPFCandidates->GetEntries(); ++k) {
-+
+    const mithep::PFCandidate *pf = (mithep::PFCandidate*)((*fPFCandidates)[k]);
-+
+    Double_t deta = (ele->Eta() - pf->Eta());
-+
+    Double_t dphi = mithep::MathUtils::DeltaPhi(Double_t(ele->Phi()),Double_t(pf->Phi()));
-+
+    Double_t dr = mithep::MathUtils::DeltaR(ele->Phi(),ele->Eta(), pf->Phi(), pf->Eta());
-+
+    if (dr >= 0.4) continue;
-+
+    if(ctrl.debug) {
-+
+      cout << "pf :: type: " << pf->PFType() << "\tpt: " << pf->Pt();
-+
+      if( pf->HasTrackerTrk() ) cout << "\tdZ: " << pf->TrackerTrk()->DzCorrected(vtx);
-+
+      cout << endl;
-+
+    }
-+
-+
-+
+    if ( (pf->HasTrackerTrk() && (pf->TrackerTrk() == ele->TrackerTrk())) ||
-+
+         (pf->HasGsfTrk() && (pf->GsfTrk() == ele->GsfTrk()))) continue;
-+
-+
-+
+    //
-+
+    // Lepton Footprint Removal
-+
+    //
-+
+    Bool_t IsLeptonFootprint = kFALSE;
-+
+    if (dr < 1.0) {
-+
-+
+      //
-+
+      // Check for electrons
-+
+      //
-+
+      for (Int_t q=0; q < electronsToVeto.size(); ++q) {
-+
+        const mithep::Electron *tmpele = electronsToVeto[q];
-+
+        // 4l electron
-+
+        if( pf->HasTrackerTrk()  ) {
-+
+          if( pf->TrackerTrk() == tmpele->TrackerTrk() ) {
-+
+            if( ctrl.debug) cout << "\tcharged tktrk, matches 4L ele ..." << endl;
-+
+            IsLeptonFootprint = kTRUE;
-+
+          }
-+
+        }
-+
+        if( pf->HasGsfTrk()  ) {
-+
+          if( pf->GsfTrk() == tmpele->GsfTrk() ) {
-+
+            if( ctrl.debug) cout << "\tcharged gsftrk, matches 4L ele ..." << endl;
-+
+            IsLeptonFootprint = kTRUE;
-+
+          }
-+
+        }
-+
+        // PF charged
-+
+        if (pf->Charge() != 0 && fabs(tmpele->SCluster()->Eta()) > 1.479
-+
+            && mithep::MathUtils::DeltaR(tmpele->Phi(),tmpele->Eta(), pf->Phi(), pf->Eta()) < 0.015) {
-+
+          if( ctrl.debug) cout << "\tcharged trk, dR matches 4L ele ..." << endl;
-+
+          IsLeptonFootprint = kTRUE;
-+
+        }
-+
+        // PF gamma
-+
+        if (abs(pf->PFType()) == PFCandidate::eGamma && fabs(tmpele->SCluster()->Eta()) > 1.479
-+
+            && mithep::MathUtils::DeltaR(tmpele->Phi(),tmpele->Eta(), pf->Phi(), pf->Eta()) < 0.08) {
-+
+          if( ctrl.debug) cout << "\tPF gamma, matches 4L ele ..." << endl;
-+
+          IsLeptonFootprint = kTRUE;
-+
+        }
-+
+      } // loop over electrons
-+
+      /* KH - comment for sync
+      //
+      // Check for muons
-<
+  double rho = 0;
-<
+  if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
-<
+    rho = fPUEnergyDensity->At(0)->Rho();
->
+//   double rho = 0;
->
+//   if (!(isnan(fPUEnergyDensity->At(0)->Rho()) || isinf(fPUEnergyDensity->At(0)->Rho())))
->
+//     rho = fPUEnergyDensity->At(0)->Rho();
+  // WARNING!!!!
+  // hardcode for sync ...
+  return pfIso;
-+
+//--------------------------------------------------------------------------------------------------
+SelectionStatus electronReferenceIsoSelection(ControlFlags &ctrl,
+                                              const mithep::Electron * ele,
+  double pfIso = electronPFIso04( ctrl, ele, vtx, fPFCandidates, fPUEnergyDensity,
+                                  EffectiveAreaVersion, muonsToVeto ,electronsToVeto );
-+
+  cout << "--------------> setting electron isoPF04 to " << pfIso << endl;
-+
+  status.isoPF04 = pfIso;
-+
+  status.chisoPF04 = gChargedIso;
-+
+  status.gaisoPF04 = gGammaIso;
-+
+  status.neisoPF04 = gNeutralIso;
-+
-+
+  bool pass = false;
-+
+  if( (pfIso/ele->Pt()) < ELECTRON_REFERENCE_PFISO_CUT ) pass = true;
-+
-+
+  if( pass ) {
-+
+    status.orStatus(SelectionStatus::LOOSEISO);
-+
+    status.orStatus(SelectionStatus::TIGHTISO);
-+
+  }
-+
+  if(ctrl.debug) cout << "returning status : " << hex << status.getStatus() << dec << endl;
-+
+  return status;
-+
-+
+}
-+
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+// hacked version
-+
+SelectionStatus electronReferenceIsoSelection(ControlFlags &ctrl,
-+
+                                              const mithep::Electron * ele,
-+
+                                              const mithep::Vertex & vtx,
-+
+                                              const mithep::Array<mithep::PFCandidate> * fPFCandidates,
-+
+                                              float rho,
-+
+                                              mithep::ElectronTools::EElectronEffectiveAreaTarget EffectiveAreaVersion,
-+
+                                              vector<const mithep::Muon*> muonsToVeto,
-+
+                                              vector<const mithep::Electron*> electronsToVeto)
-+
+//--------------------------------------------------------------------------------------------------
-+
+{
-+
-+
+  SelectionStatus status;
-+
-+
+  double pfIso = electronPFIso04( ctrl, ele, vtx, fPFCandidates, rho,
-+
+                                  EffectiveAreaVersion, muonsToVeto ,electronsToVeto );
+  bool pass = false;
+  if( (pfIso/ele->Pt()) < ELECTRON_REFERENCE_PFISO_CUT ) pass = true;

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/MitHzz4l/LeptonSelection/src/IsolationSelection.cc (file contents): Revision 1.14 by khahn, Thu May 10 11:14:33 2012 UTC vs. Revision 1.16 by khahn, Thu May 10 22:53:21 2012 UTC

Diff Legend

Comparing UserCode/MitHzz4l/LeptonSelection/src/IsolationSelection.cc (file contents):
Revision 1.14 by khahn, Thu May 10 11:14:33 2012 UTC vs.
Revision 1.16 by khahn, Thu May 10 22:53:21 2012 UTC