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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.6 by loizides, Fri Nov 28 09:13:35 2008 UTC vs.
Revision 1.25 by sixie, Sun Apr 15 12:07:00 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);
+  if (aEta < 1.1) {
+    if(iCorrectedHCAL)
-<
+      lHad *= TMath::Sin(2*TMath::ATan(TMath::Exp(lEta))); //todo ask!
->
+      lHad *= TMath::Sin(2*TMath::ATan(TMath::Exp(lEta)));
+    lTPionHad  = fpion_had_etaB;
+    lTMuonHad  = fmuon_had_etaB;
+      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;
+    } else
+      lStCrossed[i0]  = 0;
-<
+    if(iMuon->GetDX(i0) < 999999.) { //Use iMuon->GetSegmentX--> CHECK
->
+    if(iMuon->GetDX(i0) < 999999.) {
+      lNStationsSegment++;
+      lStSegmentmatch[i0] = 1;
+    } else
+  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;
-+
-+
+  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;
-+
-+
+  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->BestTrk()->NHits() <= 10) 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;
-+
+  }
-+
-+
-+
+  //NoCorrections
-+
+  if (EffectiveAreaTarget == kMuEANoCorr) {
-+
+    return 0.0;
-+
+  }
-+
+  //2011 Data Effective Areas
-+
+  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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.005;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.011;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.021;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.011;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.032;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.030;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.052;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.054;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.007;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.014;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.008;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.024;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.022;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.038;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.077;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.114;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.195;
-+
+    }
-+
+  }
-+
-+
+  //Summer11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kMuEASummer11MC) {
-+
+    //Fall11 MC Effective Areas
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.006;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.017;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.015;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.007;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.036;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.030;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.062;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.058;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.093;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.074;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.003;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.011;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.013;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.012;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.024;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.017;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.044;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.082;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.118;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.168;
-+
+    }
-+
+  }
-+
-+
+  //Fall11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kMuEAFall11MC) {
-+
+    //Fall11 MC Effective Areas
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.011;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.003;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.024;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.034;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.055;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.068;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.083;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.115;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.128;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.014;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.011;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.017;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.014;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.026;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.042;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.086;
-+
+    }
-+
+    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 && fabs(Eta) < 2.4 ) EffectiveArea = 0.135;
-+
+      if (fabs(Eta) >= 2.4) EffectiveArea = 0.200;
-+
+    }
-+
+  }
-+
-+
+  return EffectiveArea;
-+
+}
-+

Diff Legend

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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents): Revision 1.6 by loizides, Fri Nov 28 09:13:35 2008 UTC vs. Revision 1.25 by sixie, Sun Apr 15 12:07:00 2012 UTC

Diff Legend

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.6 by loizides, Fri Nov 28 09:13:35 2008 UTC vs.
Revision 1.25 by sixie, Sun Apr 15 12:07:00 2012 UTC