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

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs.
Revision 1.38 by fabstoec, Tue Apr 24 11:45:53 2012 UTC

+#include "MitPhysics/Utils/interface/ElectronTools.h"
+#include "MitAna/DataTree/interface/StableData.h"
+#include <TFile.h>
-+
+#include "MitPhysics/ElectronLikelihood/interface/ElectronLikelihood.h"
-+
+#include "MitPhysics/ElectronLikelihood/interface/LikelihoodSwitches.h"
-+
+#include "MitPhysics/ElectronLikelihood/interface/LikelihoodMeasurements.h"
+ClassImp(mithep::ElectronTools)
+  };
+  Double_t VBTFWorkingPoint80[6][8] = {
-<
+    {0.04,  0.04,  0.04,  0.04,  0.025,  0.025,  0.025,  0.025}, //hovere
->
+    {0.04,  0.04,  0.04,  0.04,  0.10,   0.10,   0.10,   0.10 }, //hovere
+    {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
+    {0.06,  0.06,  0.06,  0.06,  0.03,   0.03,   0.03,   0.03 }, //deltaphiin
+    {0.004, 0.004, 0.004, 0.004, 0.007,  0.007,  0.007,  0.007}, //deltaetain
+  if(eSeedOverPin<fCuts[4][cat+4*eb])
+    return kFALSE;
-+
+  Bool_t isoCut = kTRUE;
-+
+  if(idType == kVBTFWorkingPointFakeableId){
-+
+    double isoEcal = ele->EcalRecHitIsoDr03();
-+
+    if(ele->IsEB()) isoEcal = isoEcal - 1.0;
-+
+    isoCut = (ele->TrackIsolationDr03() < ele->Pt()*0.2) &&
-+
+             (isoEcal                   < ele->Pt()*0.2) &&
-+
+             (ele->HcalTowerSumEtDr03() < ele->Pt()*0.2);
-+
+  }
-+
+  if(isoCut == kFALSE) return kFALSE;
-+
+  // Cuts only for pt<20 region and kVBTFWorkingPointLowPtId
+  if(ele->Pt() < 20 && idType == kVBTFWorkingPointLowPtId) {
+    Bool_t isGoodLowPtEl = fBrem > 0.15 ||
+  Bool_t d0cut = kFALSE;
+  Double_t d0_real = 1e30;
-+
-+
+  if( nVertex >= (int) vertices->GetEntries() )
-+
+    nVertex = vertices->GetEntries() - 1;
-+
+  if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
+  else            {
+    Double_t distVtx = 999.0;
+  Bool_t dzcut = kFALSE;
+  Double_t distVtx = 999.0;
-+
-+
+  if( nVertex >= (int) vertices->GetEntries() )
-+
+    nVertex = vertices->GetEntries()-1;
-+
+  if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
+  else {
+    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
-<
->
+  if(distVtx < fDZCut) dzcut = kTRUE;
+  return dzcut;
+  return accept;
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+Double_t ElectronTools::Likelihood(ElectronLikelihood *LH, const Electron *ele)
-+
+{
-+
+  if (!LH) {
-+
+    std::cout << "Error: Likelihood not properly initialized\n";
-+
+    return -9999;
-+
+  }
-+
-+
+  LikelihoodMeasurements measurements;
-+
+  measurements.pt = ele->Pt();
-+
+  if (ele->IsEB() && ele->AbsEta()<1.0) measurements.subdet = 0;
-+
+  else if (ele->IsEB())                 measurements.subdet = 1;
-+
+  else                                  measurements.subdet = 2;
-+
+  measurements.deltaPhi = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
-+
+  measurements.deltaEta = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
-+
+  measurements.eSeedClusterOverPout = ele->ESeedClusterOverPout();
-+
+  measurements.eSuperClusterOverP = ele->ESuperClusterOverP();
-+
+  measurements.hadronicOverEm = ele->HadronicOverEm();
-+
+  measurements.sigmaIEtaIEta = ele->CoviEtaiEta();
-+
+  measurements.sigmaIPhiIPhi = TMath::Sqrt(ele->SCluster()->Seed()->CoviPhiiPhi());
-+
+  measurements.fBrem = ele->FBrem();
-+
+  measurements.nBremClusters = ele->NumberOfClusters() - 1;
-+
+  //measurements.OneOverEMinusOneOverP = (1.0 / ele->SCluster()->Energy()) - (1.0 / ele->BestTrk()->P());
-+
+  measurements.OneOverEMinusOneOverP = (1.0 / ele->ESuperClusterOverP() / ele->BestTrk()->P()) - (1.0 / ele->BestTrk()->P());
-+
+  double likelihood = LH->result(measurements);
-+
-+
+  double newLik = 0.0;
-+
+  if     (likelihood<=0) newLik = -20.0;
-+
+  else if(likelihood>=1) newLik =  20.0;
-+
+  else                   newLik = log(likelihood/(1.0-likelihood));
-+
-+
+  Bool_t isDebug = kFALSE;
-+
+  if(isDebug == kTRUE){
-+
+    printf("LIKELIHOOD: %f %d %f %f %f %f %f %f %f %f %d %f %f %f - %f %f\n",measurements.pt,measurements.subdet,
-+
+    measurements.deltaPhi          ,measurements.deltaEta      ,measurements.eSeedClusterOverPout,
-+
+    measurements.eSuperClusterOverP,measurements.hadronicOverEm,measurements.sigmaIEtaIEta,
-+
+    measurements.sigmaIPhiIPhi     ,measurements.fBrem         ,measurements.nBremClusters,
-+
+    measurements.OneOverEMinusOneOverP,ele->SCluster()->Energy(),ele->BestTrk()->P(),
-+
+    likelihood,newLik);
-+
+  }
-+
-+
+  return newLik;
-+
-+
+}
-+
-+
+Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
-+
-+
+  Double_t EffectiveArea = 0;
-+
-+
+  if (fabs(SCEta) < 1.0) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.045;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.034;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.079;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00016;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
-+
+  } else if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.025;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.052;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.050;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.073;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00022;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
-+
+  } else if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.170;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.060;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.187;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00030;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;
-+
+  } else if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.25 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.022;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.623;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.055;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.659;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00054;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
-+
+  } else if (fabs(SCEta) >= 2.25 && fabs(SCEta) < 2.5 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.018;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 1.198;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.073;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 1.258;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00082;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
-+
+  }
-+
-+
+  //NoCorrections
-+
+  if (EffectiveAreaTarget == kEleEANoCorr) {
-+
+    return 0.0;
-+
+  }
-+
+  //2011 Data Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEAData2011) {
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.041;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.039;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.048;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.054;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
-+
+     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.003;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.019;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.025;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.084;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.107;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.228;
-+
+    }
-+
+  }
-+
-+
+  //Summer11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEASummer11MC) {
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.023;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.046;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.046;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.058;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.073;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.083;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.053;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.048;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.062;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.085;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.118;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.004;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.040;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.086;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.105;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.169;
-+
+    }
-+
+  }
-+
-+
+  //Fall11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEAFall11MC) {
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.012;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.068;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.067;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.080;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.040;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.083;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.123;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.095;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.133;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.007;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.110;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.220;
-+
+    }
-+
+  }
-+
-+
-+
-+
+  return EffectiveArea;
-+
+}
-+
-+
-+
+Bool_t ElectronTools::PassHggLeptonTagID(const Electron* ele) {
-+
-+
+  float dist = ( ele->ConvPartnerDist()      == -9999.? 9999:TMath::Abs(ele->ConvPartnerDist()));
-+
+  float dcot = ( ele->ConvPartnerDCotTheta() == -9999.? 9999:TMath::Abs(ele->ConvPartnerDCotTheta()));
-+
-+
+  if (dist < 0.02) return false;
-+
+  if (dcot < 0.02) return false;
-+
-+
+  int numInnerHits = ele->Trk()->NExpectedHitsInner();
-+
+  if( numInnerHits > 0 ) return false;
-+
-+
+  float coviEtaiEta = ele->CoviEtaiEta();
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
-+
+  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.031 ) return false;
-+
-+
+  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
-+
+  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
-+
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.039 || deltaEtaIn > 0.005 ) ) return false;
-+
+  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.028 || deltaEtaIn > 0.007 ) ) return false;
-+
-+
+  return true;
-+
+}

Diff Legend

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

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents): Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs. Revision 1.38 by fabstoec, Tue Apr 24 11:45:53 2012 UTC

Diff Legend

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs.
Revision 1.38 by fabstoec, Tue Apr 24 11:45:53 2012 UTC