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.37 by ceballos, Fri Apr 20 20:25:58 2012 UTC vs.
Revision 1.46 by anlevin, Wed Jun 6 15:05:40 2012 UTC

# Line 189 | Line 189 | Bool_t ElectronTools::PassCustomID(const
189    if(eSeedOverPin<fCuts[4][cat+4*eb])
190      return kFALSE;
191  
192 +  // Apply detector isolation at high pt only
193    Bool_t isoCut = kTRUE;
194    if(idType == kVBTFWorkingPointFakeableId){
195      double isoEcal = ele->EcalRecHitIsoDr03();
# Line 359 | Line 360 | Bool_t ElectronTools::PassD0Cut(const El
360    Bool_t d0cut = kFALSE;
361  
362    Double_t d0_real = 1e30;
363 +
364 +  if( nVertex >= (int) vertices->GetEntries() )
365 +    nVertex = vertices->GetEntries() - 1;
366 +  
367    if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
368    else            {
369      Double_t distVtx = 999.0;
# Line 398 | Line 403 | Bool_t ElectronTools::PassDZCut(const El
403    Bool_t dzcut = kFALSE;
404  
405    Double_t distVtx = 999.0;
406 +
407 +  if( nVertex >= (int) vertices->GetEntries() )
408 +    nVertex = vertices->GetEntries()-1;
409 +  
410    if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
411    else {
412      for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
# Line 407 | Line 416 | Bool_t ElectronTools::PassDZCut(const El
416        }
417      }
418    }
419 <
419 >  
420    if(distVtx < fDZCut) dzcut = kTRUE;
421    
422    return dzcut;
# Line 1073 | Line 1082 | Double_t ElectronTools::ElectronEffectiv
1082  
1083    Double_t EffectiveArea = 0;
1084    
1085 +  if(type == ElectronTools::kEleNeutralIso04){
1086 +    if     (fabs(SCEta) < 1.0                           ) EffectiveArea = 0.19;
1087 +    else if(fabs(SCEta) >= 1.0   && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1088 +    else if(fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.000 ) EffectiveArea = 0.12;
1089 +    else if(fabs(SCEta) >= 2.000 && fabs(SCEta) < 2.200 ) EffectiveArea = 0.21;
1090 +    else if(fabs(SCEta) >= 2.200 && fabs(SCEta) < 2.300 ) EffectiveArea = 0.27;
1091 +    else if(fabs(SCEta) >= 2.300 && fabs(SCEta) < 2.400 ) EffectiveArea = 0.44;
1092 +    else                                                  EffectiveArea = 0.52;
1093 +    return EffectiveArea;
1094 +  }
1095 +
1096    if (fabs(SCEta) < 1.0) {
1097      if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1098      if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
# Line 1144 | Line 1164 | Double_t ElectronTools::ElectronEffectiv
1164    if (EffectiveAreaTarget == kEleEANoCorr) {
1165      return 0.0;
1166    }
1167 +  else if (EffectiveAreaTarget == kEleEAData2012) {
1168 +    if (type == kEleGammaAndNeutralHadronIso04) {
1169 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
1170 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1171 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.12;
1172 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.21;
1173 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.27;
1174 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.44;
1175 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.52;
1176 +    }                          
1177 +    if (type == kEleGammaIsoDR0p0To0p1) {
1178 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1179 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1180 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
1181 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1182 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1183 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.013;
1184 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.013;
1185 +    }
1186 +    if (type == kEleGammaIsoDR0p1To0p2) {
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.013;
1189 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1190 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
1191 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1192 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1193 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.102;
1194 +    }
1195 +    if (type == kEleGammaIsoDR0p2To0p3) {
1196 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
1197 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1198 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
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.041;
1201 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.034;
1202 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.042;
1203 +    }
1204 +    if (type == kEleGammaIsoDR0p3To0p4) {
1205 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.039;
1206 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1207 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1208 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1209 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.053;
1210 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.059;
1211 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.069;
1212 +    }
1213 +    if (type == kEleGammaIsoDR0p4To0p5) {
1214 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1215 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.045;
1216 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
1217 +     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.043;
1218 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.056;
1219 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.065;
1220 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.074;
1221 +    }
1222 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1223 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1224 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.006;
1225 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
1226 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1227 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.001;
1228 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.003;
1229 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.008;
1230 +    }
1231 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1232 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1233 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
1234 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1235 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.003;
1236 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1237 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.006;
1238 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.010;
1239 +    }
1240 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1241 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.007;
1242 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1243 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1244 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
1245 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.007;
1246 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.018;
1247 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.028;
1248 +    }
1249 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1250 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1251 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1252 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
1253 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
1254 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.018;
1255 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.026;
1256 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.063;
1257 +    }
1258 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1259 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
1260 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.012;
1261 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1262 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
1263 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1264 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.051;
1265 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.143;
1266 +    }
1267 +  }
1268 +
1269    //2011 Data Effective Areas
1270    else if (EffectiveAreaTarget == kEleEAData2011) {
1271 +    if (type == kEleGammaAndNeutralHadronIso03) {
1272 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.100;
1273 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.120;
1274 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.085;
1275 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.110;
1276 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.120;
1277 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.120;
1278 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.130;
1279 +    }
1280 +    if (type == kEleGammaAndNeutralHadronIso04) {
1281 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.180;
1282 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.200;
1283 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.150;
1284 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.190;
1285 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.210;
1286 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.220;
1287 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.290;
1288 +    }
1289      if (type == kEleGammaIsoDR0p0To0p1) {
1290        if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1291        if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
# Line 1426 | Line 1566 | Double_t ElectronTools::ElectronEffectiv
1566      }
1567    }
1568  
1429 –
1430 –
1569    return EffectiveArea;  
1570   }
1571  
1572  
1573 + Bool_t ElectronTools::PassHggLeptonTagID(const Electron* ele) {
1574 +  
1575 +  float dist = ( ele->ConvPartnerDist()      == -9999.? 9999:TMath::Abs(ele->ConvPartnerDist()));
1576 +  float dcot = ( ele->ConvPartnerDCotTheta() == -9999.? 9999:TMath::Abs(ele->ConvPartnerDCotTheta()));  
1577 +  
1578 +  if (dist < 0.02) return false;
1579 +  if (dcot < 0.02) return false;
1580 +
1581 +  int numInnerHits = ele->Trk()->NExpectedHitsInner();
1582 +  if( numInnerHits > 0 ) return false;
1583 +
1584 +  float coviEtaiEta = ele->CoviEtaiEta();
1585 +  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1586 +  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.031 ) return false;
1587 +
1588 +  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1589 +  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1590 +
1591 +  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.039 || deltaEtaIn > 0.005 ) ) return false;
1592 +  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.028 || deltaEtaIn > 0.007 ) ) return false;
1593 +
1594 +  return true;
1595 + }

Diff Legend

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