ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/src/ElectronTools.cc
(Generate patch)

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.46 by anlevin, Wed Jun 6 15:05:40 2012 UTC vs.
Revision 1.49 by mingyang, Mon Oct 8 17:22:59 2012 UTC

# Line 1165 | Line 1165 | Double_t ElectronTools::ElectronEffectiv
1165      return 0.0;
1166    }
1167    else if (EffectiveAreaTarget == kEleEAData2012) {
1168 +    if (type == kEleGammaIso03) {
1169 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.122;
1170 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.147;
1171 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.055;
1172 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.106;
1173 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.138;
1174 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.221;
1175 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.211;
1176 +    }
1177 +    if (type == kEleGammaIso04) {
1178 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.176;
1179 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.206;
1180 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.094;
1181 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.172;
1182 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.244;
1183 +      if (fabs(SCEta) >= 2.4 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.333;
1184 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.348;
1185 +    }
1186 +    if (type == kEleNeutralHadronIso03) {
1187 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1188 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1189 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1190 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1191 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1192 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1193 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.019;
1194 +    }
1195 +    if (type == kEleNeutralHadronIso04) {
1196 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.022;
1197 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.036;
1198 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1199 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
1200 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.052;
1201 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.063;
1202 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.028;
1203 +    }
1204      if (type == kEleGammaAndNeutralHadronIso04) {
1205        if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
1206        if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
# Line 1577 | Line 1613 | Bool_t ElectronTools::PassHggLeptonTagID
1613    
1614    if (dist < 0.02) return false;
1615    if (dcot < 0.02) return false;
1616 <
1616 >  
1617    int numInnerHits = ele->Trk()->NExpectedHitsInner();
1618 <  if( numInnerHits > 0 ) return false;
1618 >  if( numInnerHits > 1 ) return false;
1619  
1620    float coviEtaiEta = ele->CoviEtaiEta();
1621    if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
# Line 1593 | Line 1629 | Bool_t ElectronTools::PassHggLeptonTagID
1629  
1630    return true;
1631   }
1632 +
1633 + Bool_t ElectronTools::PassHggLeptonTagID2012(const Electron* ele) {
1634 +  
1635 +  if (TMath::Abs(1./ele->E()-1./ele->Pt())>0.05) return false;
1636 +
1637 +  int numInnerHits = ele->Trk()->NExpectedHitsInner();
1638 +  if( numInnerHits > 1 ) return false;
1639 +
1640 +  float coviEtaiEta = ele->CoviEtaiEta();
1641 +  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1642 +  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.03 ) return false; // h
1643 +
1644 +  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1645 +  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1646 +
1647 +  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.15 || deltaEtaIn > 0.007 || ele->HadronicOverEm()>0.12) ) return false;   // h
1648 +  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.10 || deltaEtaIn > 0.009 || ele->HadronicOverEm()>0.10 ) ) return false;   // h
1649 +
1650 +  return true;
1651 + }
1652 +

Diff Legend

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