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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.7 by loizides, Thu Dec 4 13:53:34 2008 UTC vs.
Revision 1.29 by ceballos, Mon May 28 06:43:33 2012 UTC

+#include "MitPhysics/Utils/interface/MuonTools.h"
+#include <TFile.h>
-+
+ClassImp(mithep::MuonTools)
-+
+using namespace mithep;
+//--------------------------------------------------------------------------------------------------
+Double_t MuonTools::GetCaloCompatability(const Muon *iMuon,
+                                         Bool_t iEMSpecial, Bool_t iCorrectedHCAL) const
-<
+  // todo
->
+  // Get calo compatibility value for given muon based on calorimeter templates.
->
+  // If iEMSpecial is true, then a use different arrangement of ECAL for compatibility.
+  Double_t lEta = iMuon->Eta();
+  Double_t aEta = TMath::Abs(lEta);
+      return iMuon->PromptTight(Muon::kAny);
+      break;
+    case kTMOneStationLoose:
-<
+      return iMuon->TMOneStation(99999,999999);
->
+      return iMuon->TMOneStation(999999,999999);
+      break;
+    case kTMOneStationTight:
+      return iMuon->TMOneStation();
+      break;
-<
+  Double_t lVal = 1.2*GetSegmentCompatability(iMuon);
-<
+  if (lVal/1.2 == 0.5)
->
+  Double_t lVal = GetSegmentCompatability(iMuon);
->
+  if (lVal == 0.5) // exclude this border case
+    return kFALSE;
-+
+  lVal *= 1.2;
+  lVal += 0.8*GetCaloCompatability(iMuon,kTRUE,kTRUE);
+  if (lVal > tm2dcut)
+    return kTRUE;
-+
+  return kFALSE;
+//--------------------------------------------------------------------------------------------------
+Double_t MuonTools::GetSegmentCompatability(const mithep::Muon *iMuon) const
-<
+  // todo
->
+  // Get segment compatability for given muon based on likelihood of well defined
->
+  // track through chambers.
+  Int_t lNStationsCrossed = 0;
+  Int_t lNStationsSegment = 0;
+  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)
+  ret->SetDirectory(0);
+  return ret;
-+
+//--------------------------------------------------------------------------------------------------
-+
+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;
-+
+  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;
-+
-+
+  return d0cut;
-+
+}
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+Bool_t MuonTools::PassD0Cut(const Muon *mu, const BeamSpotCol *beamspots, Double_t fD0Cut)
-+
+{
-+
+  Bool_t d0cut = kFALSE;
-+
+  const Track *mt = mu->BestTrk();
-+
+  if (!mt) return kFALSE;
-+
-+
+  // d0 cut
-+
+  Double_t d0_real = 99999;
-+
+  for(UInt_t i0 = 0; i0 < beamspots->GetEntries(); i0++) {
-+
+    Double_t pD0 = mt->D0Corrected(*beamspots->At(i0));
-+
+    if(TMath::Abs(pD0) < TMath::Abs(d0_real)) d0_real = TMath::Abs(pD0);
-+
+  }
-+
+  if(d0_real < fD0Cut) d0cut = kTRUE;
-+
-+
+  return d0cut;
-+
+}
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+Bool_t MuonTools::PassDZCut(const Muon *mu, const VertexCol *vertices, Double_t fDZCut, Int_t nVertex)
-+
+{
-+
+  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;
-+
+  }
-+
-+
+  //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 (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;
-+
+    }
-+
+  }
-+
-+
+  //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;
-+
-+
+}
-+

Diff Legend

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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents): Revision 1.7 by loizides, Thu Dec 4 13:53:34 2008 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.7 by loizides, Thu Dec 4 13:53:34 2008 UTC vs.
Revision 1.29 by ceballos, Mon May 28 06:43:33 2012 UTC