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

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.37 by ceballos, Fri Apr 20 20:25:58 2012 UTC vs.
Revision 1.51 by pharris, Wed Jan 9 18:22:07 2013 UTC

+  if(eSeedOverPin<fCuts[4][cat+4*eb])
+    return kFALSE;
-+
+  // Apply detector isolation at high pt only
+  Bool_t isoCut = kTRUE;
+  if(idType == kVBTFWorkingPointFakeableId){
+    double isoEcal = ele->EcalRecHitIsoDr03();
+//--------------------------------------------------------------------------------------------------
-+
+Bool_t ElectronTools::PassConversionFilterPFAOD(const Electron *ele,
-+
+                                           const DecayParticleCol *conversions,
-+
+                                           const BaseVertex *vtx,
-+
+                                           UInt_t nWrongHitsMax,
-+
+                                           Double_t probMin,
-+
+                                           Double_t lxyMin,
-+
+                                           Bool_t matchCkf,
-+
+                                           Bool_t requireArbitratedMerged,
-+
+                                           Double_t trkptMin)
-+
+{
-+
-+
+  Bool_t isGoodConversion              = false;
-+
+  for(unsigned int ifc=0; ifc<conversions->GetEntries(); ++ifc){
-+
+    if(!(conversions->At(ifc)->Prob() > probMin) && (!requireArbitratedMerged || conversions->At(ifc)->Quality().Quality(ConversionQuality::arbitratedMerged)) && (conversions->At(ifc)->LxyCorrected((BaseVertex*)vtx) > lxyMin)) continue;
-+
+    bool conversionMatchFound = false;
-+
+    for(unsigned int d=0; d<conversions->At(ifc)->NDaughters(); ++d){
-+
+      const ChargedParticle *pParticle = 0;
-+
+      pParticle = dynamic_cast<const ChargedParticle*>(conversions->At(ifc)->Daughter(d));
-+
+      if(pParticle == 0) continue;
-+
+      const Track* trk = 0;
-+
+      trk = pParticle->Trk();
-+
+      if(trk == 0) continue;
-+
+      if( ele->GsfTrk() == trk || (matchCkf && ele->TrackerTrk()==trk) ){ conversionMatchFound = true; break; }
-+
+      const StableData* sd = dynamic_cast<const StableData*> (conversions->At(ifc)->DaughterDat(d));
-+
+      isGoodConversion = true;
-+
+      if( sd->NWrongHits() > nWrongHitsMax ){ isGoodConversion = false; }
-+
+    }
-+
+    if(isGoodConversion) break;
-+
+  }
-+
+  return !isGoodConversion;
-+
+}
-+
-+
+//--------------------------------------------------------------------------------------------------
+Bool_t ElectronTools::PassD0Cut(const Electron *ele, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
+  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;
+  Double_t EffectiveArea = 0;
-+
+  if(type == ElectronTools::kEleNeutralIso04){
-+
+    if     (fabs(SCEta) < 1.0                           ) EffectiveArea = 0.19;
-+
+    else if(fabs(SCEta) >= 1.0   && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
-+
+    else if(fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.000 ) EffectiveArea = 0.12;
-+
+    else if(fabs(SCEta) >= 2.000 && fabs(SCEta) < 2.200 ) EffectiveArea = 0.21;
-+
+    else if(fabs(SCEta) >= 2.200 && fabs(SCEta) < 2.300 ) EffectiveArea = 0.27;
-+
+    else if(fabs(SCEta) >= 2.300 && fabs(SCEta) < 2.400 ) EffectiveArea = 0.44;
-+
+    else                                                  EffectiveArea = 0.52;
-+
+    return EffectiveArea;
-+
+  }
-+
+  if (fabs(SCEta) < 1.0) {
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
+  if (EffectiveAreaTarget == kEleEANoCorr) {
+    return 0.0;
-+
+  else if (EffectiveAreaTarget == kEleEAData2012) {
-+
+    if (type == kEleGammaIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.122;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.147;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.106;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.138;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.221;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.211;
-+
+    }
-+
+    if (type == kEleGammaIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.176;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.206;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.094;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.172;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.244;
-+
+      if (fabs(SCEta) >= 2.4 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.333;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.348;
-+
+    }
-+
+    if (type == kEleNeutralHadronIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.019;
-+
+    }
-+
+    if (type == kEleNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.022;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.052;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.063;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.028;
-+
+    }
-+
+    if (type == kEleGammaAndNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.12;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.21;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.27;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.44;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.52;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.013;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.102;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.034;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.042;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.039;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.053;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.059;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.069;
-+
+    }
-+
+    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.045;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
-+
+     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.056;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.065;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.074;
-+
+    }
-+
+    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.006;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.008;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.010;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.007;
-+
+      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.009;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.028;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.063;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.143;
-+
+    }
-+
+  }
-+
+  //2011 Data Effective Areas
+  else if (EffectiveAreaTarget == kEleEAData2011) {
-+
+    if (type == kEleGammaAndNeutralHadronIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.100;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.085;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.110;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.130;
-+
+    }
-+
+    if (type == kEleGammaAndNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.180;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.200;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.150;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.190;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.210;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.220;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.290;
-+
+    }
+    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;
-+
+  return EffectiveArea;
-+
+}
-<
+  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 > 1 ) 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;
->
+}
->
->
+Bool_t ElectronTools::PassHggLeptonTagID2012(const Electron* ele) {
->
->
+  if (TMath::Abs(1./ele->E()-1./ele->Pt())>0.05) return false;
->
->
+  int numInnerHits = ele->Trk()->NExpectedHitsInner();
->
+  if( numInnerHits > 1 ) 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.03 ) return false; // h
->
->
+  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
->
+  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
->
->
+  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.15 || deltaEtaIn > 0.007 || ele->HadronicOverEm()>0.12) ) return false;   // h
->
+  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.10 || deltaEtaIn > 0.009 || ele->HadronicOverEm()>0.10 ) ) return false;   // h
->
->
+  return true;
-+
+std::pair<Double_t,Double_t> ElectronTools::ComputeEPCombination( const Electron * ele, const float regression_energy,
-+
+                                                  const float regression_energy_error) {
-+
-+
+  enum Classification { GSF_ELECTRON_UNKNOWN=-1,
-+
+                        GSF_ELECTRON_GOLDEN=0,
-+
+                        GSF_ELECTRON_BIGBREM=1,
-+
+                        GSF_ELECTRON_BADTRACK=2,
-+
+                        GSF_ELECTRON_SHOWERING=3,
-+
+                        GSF_ELECTRON_GAP=4 } ;
-+
+  float newEnergyError_ = regression_energy_error;
-+
+  float scEnergy = regression_energy;
-+
-+
+  int elClass = ele->Classification();
-+
-+
+  float trackMomentum  = ele->PIn() ;
-+
+  float errorTrackMomentum_ = 999. ;
-+
-+
+  // the electron's track momentum error was not available in bambu versions less than 029
-+
+  if(ele->TrackMomentumError() > 0)
-+
+    errorTrackMomentum_ = ele->TrackMomentumError();
-+
+  else if ( ele->GsfTrk()->PtErr() > 0)
-+
+    errorTrackMomentum_ = ele->GsfTrk()->PtErr()*cosh(ele->GsfTrk()->Eta());
-+
+  else
-+
+    assert(0);
-+
-+
+  float finalMomentum = ele->E(); // initial
-+
+  float finalMomentumError = 999.;
-+
-+
+  // first check for large errors
-+
-+
+  if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy <= 0.5) {
-+
+    finalMomentum = regression_energy;    finalMomentumError = regression_energy_error;
-+
+  }
-+
+  else if (errorTrackMomentum_/trackMomentum <= 0.5 && regression_energy_error/regression_energy > 0.5){
-+
+    finalMomentum = trackMomentum;  finalMomentumError = errorTrackMomentum_;
-+
+  }
-+
+  else if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy > 0.5){
-+
+    if (errorTrackMomentum_/trackMomentum < regression_energy_error/regression_energy) {
-+
+      finalMomentum = trackMomentum; finalMomentumError = errorTrackMomentum_;
-+
+    }
-+
+    else{
-+
+      finalMomentum = regression_energy; finalMomentumError = regression_energy_error;
-+
+    }
-+
+  }
-+
+  // then apply the combination algorithm
-+
+  else {
-+
+     // calculate E/p and corresponding error
-+
+    float eOverP = regression_energy / trackMomentum;
-+
+    float errorEOverP = sqrt(
-+
+                             (regression_energy_error/trackMomentum)*(regression_energy_error/trackMomentum) +
-+
+                             (regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum)*
-+
+                             (regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum));
-+
-+
-+
+    bool eleIsNotInCombination = false ;
-+
+    if ( (eOverP  > 1 + 2.5*errorEOverP) || (eOverP  < 1 - 2.5*errorEOverP) || (eOverP < 0.8) || (eOverP > 1.3) )
-+
+      { eleIsNotInCombination = true ; }
-+
+     if (eleIsNotInCombination)
-+
+       {
-+
+         if (eOverP > 1)
-+
+           { finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
-+
+         else
-+
+           {
-+
+             if (elClass == GSF_ELECTRON_GOLDEN)
-+
+               { finalMomentum = regression_energy; finalMomentumError = regression_energy_error; }
-+
+             if (elClass == GSF_ELECTRON_BIGBREM)
-+
+               {
-+
+                 if (regression_energy<36)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
-+
+               }
-+
+             if (elClass == GSF_ELECTRON_BADTRACK)
-+
+               { finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
-+
+             if (elClass == GSF_ELECTRON_SHOWERING)
-+
+               {
-+
+                 if (regression_energy<30)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy; finalMomentumError = regression_energy_error;}
-+
+               }
-+
+             if (elClass == GSF_ELECTRON_GAP)
-+
+               {
-+
+                 if (regression_energy<60)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
-+
+               }
-+
+           }
-+
+       }
-+
+     else
-+
+       {
-+
+         // combination
-+
+         finalMomentum = (regression_energy/regression_energy_error/regression_energy_error + trackMomentum/errorTrackMomentum_/errorTrackMomentum_) /
-+
+           (1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
-+
+         float finalMomentumVariance = 1 / (1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
-+
+         finalMomentumError = sqrt(finalMomentumVariance);
-+
+       }
-+
+  }
-+
-+
+  std::pair<Double_t,Double_t> value;
-+
+  value.first = finalMomentum;
-+
+  value.second = finalMomentumError;
-+
+  return value;
-+
-+
+}

Diff Legend

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

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents): Revision 1.37 by ceballos, Fri Apr 20 20:25:58 2012 UTC vs. Revision 1.51 by pharris, Wed Jan 9 18:22:07 2013 UTC

Diff Legend

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.37 by ceballos, Fri Apr 20 20:25:58 2012 UTC vs.
Revision 1.51 by pharris, Wed Jan 9 18:22:07 2013 UTC