[ViewVC] Diff of: cvsroot/UserCode/MitPhysics/Utils/src/MuonTools.cc

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.14 by ceballos, Fri Jan 21 11:25:29 2011 UTC vs.
Revision 1.29 by ceballos, Mon May 28 06:43:33 2012 UTC

+  Double_t lStBoundary[8];
+  Double_t lWeight  = 0.;
-–
+  Bool_t lAdjust    = kTRUE;
+  for (Int_t i0 = 0; i0 < 8; ++i0) {
+    lStBoundary[i0] = 0.;
+    if(iMuon->GetTrackDist(i0) < 999999. ) {
+    } else
+      lStSegmentmatch[i0] = 0;
-–
+    if(iMuon->GetDY(i0) < 999999.)
-–
+      lAdjust = kFALSE;
+  if (lNStationsCrossed == 0)
+  return ret;
+//--------------------------------------------------------------------------------------------------
-<
+Bool_t MuonTools::PassD0Cut(const Muon *mu, const VertexCol *vertices, Double_t fD0Cut)
->
+Bool_t MuonTools::PassD0Cut(const Muon *mu, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
+  Bool_t d0cut = kFALSE;
+  const Track *mt = mu->BestTrk();
+  if (!mt) return kFALSE;
-+
+  if( nVertex >= (int) vertices->GetEntries() )
-+
+    nVertex = vertices->GetEntries() - 1;
-+
+  Double_t d0_real = 1e30;
-<
+  for(UInt_t i0 = 0; i0 < vertices->GetEntries(); i0++) {
-<
+    if(vertices->At(i0)->NTracks() > 0){
-<
+      Double_t pD0 = mt->D0Corrected(*vertices->At(i0));
-<
+      d0_real = TMath::Abs(pD0);
-<
+      break;
->
+  if(nVertex >= 0) d0_real = TMath::Abs(mt->D0Corrected(*vertices->At(nVertex)));
->
+  else            {
->
+    Double_t distVtx = 999.0;
->
+    Int_t closestVtx = 0;
->
+    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
->
+      double dz = TMath::Abs(mt->DzCorrected(*vertices->At(nv)));
->
+      if(dz < distVtx) {
->
+        distVtx    = dz;
->
+        closestVtx = nv;
->
+      }
-+
+    d0_real = TMath::Abs(mt->D0Corrected(*vertices->At(closestVtx)));
+  if(d0_real < fD0Cut) d0cut = kTRUE;
+//--------------------------------------------------------------------------------------------------
-<
+Bool_t MuonTools::PassSoftMuonCut(const Muon *mu, const VertexCol *vertices)
->
+Bool_t MuonTools::PassDZCut(const Muon *mu, const VertexCol *vertices, Double_t fDZCut, Int_t nVertex)
-<
+  if(mu->Pt() <= 3.0) return kFALSE;
->
+  Bool_t dzcut = kFALSE;
->
+  const Track *mt = mu->BestTrk();
->
+  if (!mt) return kFALSE;
->
->
+  if( nVertex >= (int) vertices->GetEntries() )
->
+    nVertex = vertices->GetEntries() - 1;
->
->
+  Double_t distVtx = 999.0;
->
+  if(nVertex >= 0) distVtx = TMath::Abs(mt->DzCorrected(*vertices->At(nVertex)));
->
+  else {
->
+    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
->
+      double dz = TMath::Abs(mt->DzCorrected(*vertices->At(nv)));
->
+      if(dz < distVtx) {
->
+        distVtx = dz;
->
+      }
->
+    }
->
+  }
->
->
+  if(distVtx < fDZCut) dzcut = kTRUE;
-+
+  return dzcut;
-+
+}
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+Bool_t MuonTools::PassSoftMuonCut(const Muon *mu, const VertexCol *vertices, const Double_t fDZCut,
-+
+                                  const Bool_t applyIso)
-+
+{
-+
+  if(mu->Pt() <= 3.0) return kFALSE;
-+
+  if(!mu->IsTrackerMuon()) return kFALSE;
+  if(!mu->Quality().Quality(MuonQuality::TMLastStationAngTight)) return kFALSE;
-+
-+
+  if(mu->NTrkLayersHit() <= 5) return kFALSE;
-+
-+
+  if(!PassD0Cut(mu, vertices, 0.2, 0)) return kFALSE;
-+
-+
+  if(!PassDZCut(mu, vertices, fDZCut, 0)) return kFALSE;
-+
-+
+  if(applyIso == kTRUE){
-+
+    Double_t totalIso = 1.0 * mu->IsoR03SumPt() +
-+
+.0 * mu->IsoR03EmEt() +
-+
+.0 * mu->IsoR03HadEt();
-+
+    if (totalIso < (mu->Pt()*0.10) && mu->Pt() > 20.0) return kFALSE;
-+
+  }
-+
-+
+  return kTRUE;
-+
+}
-+
-+
+Double_t MuonTools::MuonEffectiveArea(EMuonEffectiveAreaType type, Double_t Eta,
-+
+                                      EMuonEffectiveAreaTarget EffectiveAreaTarget) {
-+
-+
+  Double_t EffectiveArea = 0;
-+
+  if (fabs(Eta) < 1.0) {
-+
+    if (type == kMuChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso03) EffectiveArea = 0.080;
-+
+    if (type == kMuChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso04) EffectiveArea = 0.163;
-+
+    if (type == kMuHadEnergy)    EffectiveArea = 0.000;
-+
+    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuHadS9Energy)  EffectiveArea = 0.016;
-+
+    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuTrkIso03)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso03)      EffectiveArea = 0.080;
-+
+    if (type == kMuHadIso03)     EffectiveArea = 0.025;
-+
+    if (type == kMuTrkIso05)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso05)      EffectiveArea = 0.290;
-+
+    if (type == kMuHadIso05)     EffectiveArea = 0.091;
-+
+  } else if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) {
-+
+    if (type == kMuChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso03) EffectiveArea = 0.083;
-+
+    if (type == kMuChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso04) EffectiveArea = 0.168;
-+
+    if (type == kMuHadEnergy)    EffectiveArea = 0.005;
-+
+    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuHadS9Energy)  EffectiveArea = 0.041;
-+
+    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuTrkIso03)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso03)      EffectiveArea = 0.043;
-+
+    if (type == kMuHadIso03)     EffectiveArea = 0.028;
-+
+    if (type == kMuTrkIso05)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso05)      EffectiveArea = 0.184;
-+
+    if (type == kMuHadIso05)     EffectiveArea = 0.106;
-+
+  } else if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) {
-+
+    if (type == kMuChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso03) EffectiveArea = 0.060;
-+
+    if (type == kMuChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso04) EffectiveArea = 0.131;
-+
+    if (type == kMuHadEnergy)    EffectiveArea = 0.020;
-+
+    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuHadS9Energy)  EffectiveArea = 0.072;
-+
+    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuTrkIso03)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso03)      EffectiveArea = 0.025;
-+
+    if (type == kMuHadIso03)     EffectiveArea = 0.036;
-+
+    if (type == kMuTrkIso05)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso05)      EffectiveArea = 0.124;
-+
+    if (type == kMuHadIso05)     EffectiveArea = 0.140;
-+
+  } else if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.25 ) {
-+
+    if (type == kMuChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso03) EffectiveArea = 0.066;
-+
+    if (type == kMuChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso04) EffectiveArea = 0.149;
-+
+    if (type == kMuHadEnergy)    EffectiveArea = 0.056;
-+
+    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuHadS9Energy)  EffectiveArea = 0.148;
-+
+    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuTrkIso03)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso03)      EffectiveArea = 0.025;
-+
+    if (type == kMuHadIso03)     EffectiveArea = 0.050;
-+
+    if (type == kMuTrkIso05)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso05)      EffectiveArea = 0.120;
-+
+    if (type == kMuHadIso05)     EffectiveArea = 0.186;
-+
+  } else if (fabs(Eta) >= 2.25 && fabs(Eta) < 2.4 ) {
-+
+    if (type == kMuChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso03) EffectiveArea = 0.098;
-+
+    if (type == kMuChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == kMuNeutralIso04) EffectiveArea = 0.200;
-+
+    if (type == kMuHadEnergy)    EffectiveArea = 0.093;
-+
+    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
-+
+    if (type == kMuHadS9Energy)  EffectiveArea = 0.260;
-+
+    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
-+
+    if (type == kMuTrkIso03)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso03)      EffectiveArea = 0.027;
-+
+    if (type == kMuHadIso03)     EffectiveArea = 0.060;
-+
+    if (type == kMuTrkIso05)     EffectiveArea = 0.000;
-+
+    if (type == kMuEMIso05)      EffectiveArea = 0.139;
-+
+    if (type == kMuHadIso05)     EffectiveArea = 0.228;
-+
+  }
-+
-+
+  if (EffectiveAreaTarget == kMuEANoCorr) {
-+
+    return 0.0;
-+
+  }
-<
+  if(mu->BestTrk()->NHits() <= 10) return kFALSE;
->
+  //2012 Data Effective Areas
->
+  else if (EffectiveAreaTarget == kMuEAData2012) {
->
+            if (type == kMuGammaIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.023;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.3 ) EffectiveArea = 0.000;
->
+    }
->
+    if (type == kMuGammaIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.007;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 2.3 ) EffectiveArea = 0.016;
->
+    }
->
+    if (type == kMuGammaIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.027;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.021;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.018;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.025;
->
+    }
->
+    if (type == kMuGammaIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.044;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.031;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.024;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.021;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.040;
->
+    }
->
+    if (type == kMuGammaIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.060;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.045;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.027;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.037;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.043;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.052;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.003;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.003;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.003;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.007;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.015;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.013;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.014;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.024;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.023;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.077;
->
+    }
->
->
->
->
+    if (type == kMuGammaIso04){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.50419;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.30582;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.19765;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.28723;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.52529;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.48818;
->
+    }
->
+    if (type == kMuNeutralHadronIso04){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.16580;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.25904;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.24695;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.22021;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.34045;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.21592;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso04){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.674;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.565;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.442;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.515;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.821;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.660;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso03){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.382;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.317;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.242;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.326;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.462;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.372;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso04Tight){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.340;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.310;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.315;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.415;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.658;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.405;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso03Tight){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.207;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.183;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.177;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.271;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.348;
->
+      if (fabs(Eta) >= 2.3 )                  EffectiveArea = 0.246;
->
+    }
->
+  }
->
->
+  //2011 Data Effective Areas
->
+  else if (EffectiveAreaTarget == kMuEAData2011) {
->
->
+    if (type == kMuGammaIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.3 ) EffectiveArea = 0.005;
->
+    }
->
+    if (type == kMuGammaIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.011;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 2.3 ) EffectiveArea = 0.011;
->
+    }
->
+    if (type == kMuGammaIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.023;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.014;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.017;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.021;
->
+    }
->
+    if (type == kMuGammaIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.036;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.017;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.023;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.028;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.032;
->
+    }
->
+    if (type == kMuGammaIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.051;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.037;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.028;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.033;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.042;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.052;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.001;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.001;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.001;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.007;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.014;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.017;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.017;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.024;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.021;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.024;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.032;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.038;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.038;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.020;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.033;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.045;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.051;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.114;
->
+    }
->
+    /// BEGIN FROM SLIDE 11 OF  https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=188494
->
+    /// NOTE: to be used with the rho from ALL pf candidates within |eta|<2.5
->
+    if (type == kMuGammaIso03){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.049;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.030;
->
+      if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.022;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.034;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.041;
->
+      if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.048;
->
+    }
->
+    if (type == kMuGammaIso04){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.085;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.052;
->
+      if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.038;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.055;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.070;
->
+      if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.081;
->
+    }
->
+    if (type == kMuNeutralHadronIso03){
->
+        if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.027;
->
+        if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.039;
->
+        if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.044;
->
+        if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.047;
->
+        if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.055;
->
+        if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.065;
->
+    }
->
+    if (type == kMuNeutralHadronIso04){
->
+        if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.046;
->
+        if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.067;
->
+        if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.074;
->
+        if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.083;
->
+        if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.095;
->
+        if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.105;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso03){
->
+        if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.076;
->
+        if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.070;
->
+        if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.067;
->
+        if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.082;
->
+        if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.097;
->
+        if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.115;
->
+    }
->
+    if (type == kMuGammaAndNeutralHadronIso04){
->
+        if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.132;
->
+        if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.5 ) EffectiveArea = 0.120;
->
+        if (fabs(Eta) >= 1.5 && fabs(Eta) < 2.0 ) EffectiveArea = 0.114;
->
+        if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.139;
->
+        if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.168;
->
+        if (fabs(Eta) >= 2.3 )                EffectiveArea = 0.189;
->
+    }
->
+    /// END FROM SLIDE 11 OF  https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=188494
-<
+  if(!PassD0Cut(mu, vertices, 0.2)) return kFALSE;
->
+    if (type == kMuGammaIso05){
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.05317;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.03502;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.03689;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.05221;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.06668;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.0744;
->
+    }
->
+    if (type == kMuNeutralIso05) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.06408;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.07557;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.08864;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.11492;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.13784;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.18745;
->
+    }
->
+  }
-<
+  Double_t totalIso = 1.0 * mu->IsoR03SumPt() +
-<
+.0 * mu->IsoR03EmEt() +
-<
+.0 * mu->IsoR03HadEt();
-<
+  if (totalIso < (mu->Pt()*0.10) && mu->Pt() > 20.0) return kFALSE;
->
+  //Summer11 MC Effective Areas
->
+  else if (EffectiveAreaTarget == kMuEASummer11MC) {
->
+    if (type == kMuGammaIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.006;
->
+    }
->
+    if (type == kMuGammaIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.007;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.015;
->
+    }
->
+    if (type == kMuGammaIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.023;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.018;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.024;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.036;
->
+    }
->
+    if (type == kMuGammaIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.038;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.027;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.033;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.041;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.062;
->
+    }
->
+    if (type == kMuGammaIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.055;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.038;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.032;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.052;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.066;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.093;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.000;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.003;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.013;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.020;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.024;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.021;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.025;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.030;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.044;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.030;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.038;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.048;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.118;
->
+    }
->
+  }
->
->
+  //Fall11 MC Effective Areas
->
+  else if (EffectiveAreaTarget == kMuEAFall11MC) {
->
+    if (type == kMuGammaIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.003;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.003;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.011;
->
+    }
->
+    if (type == kMuGammaIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.008;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.006;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.019;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.024;
->
+    }
->
+    if (type == kMuGammaIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.020;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.012;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.022;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.027;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.034;
->
+    }
->
+    if (type == kMuGammaIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.042;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.033;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.022;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.036;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.059;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.068;
->
+    }
->
+    if (type == kMuGammaIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.060;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.043;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.036;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.055;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.092;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.115;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p0To0p1) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.002;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.004;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.010;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.014;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p1To0p2) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 )   EffectiveArea = 0.005;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.007;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 )   EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 )   EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 2.3  )                 EffectiveArea = 0.017;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p2To0p3) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.009;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.015;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.016;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.018;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.022;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.026;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p3To0p4) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.013;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.021;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.032;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.037;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.042;
->
+    }
->
+    if (type == kMuNeutralHadronIsoDR0p4To0p5) {
->
+      if (fabs(Eta) >= 0.0 && fabs(Eta) < 1.0 ) EffectiveArea = 0.017;
->
+      if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) EffectiveArea = 0.026;
->
+      if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) EffectiveArea = 0.035;
->
+      if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.2 ) EffectiveArea = 0.046;
->
+      if (fabs(Eta) >= 2.2 && fabs(Eta) < 2.3 ) EffectiveArea = 0.063;
->
+      if (fabs(Eta) >= 2.3  ) EffectiveArea = 0.135;
->
+    }
->
+  }
->
->
+  return EffectiveArea;
-–
+  return kTRUE;
-+
-+

Diff Legend

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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents): Revision 1.14 by ceballos, Fri Jan 21 11:25:29 2011 UTC vs. Revision 1.29 by ceballos, Mon May 28 06:43:33 2012 UTC

Diff Legend

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.14 by ceballos, Fri Jan 21 11:25:29 2011 UTC vs.
Revision 1.29 by ceballos, Mon May 28 06:43:33 2012 UTC