This is an old revision of the document!
PUPPI
git cms-addpkg CommonTools/PileupAlgos
git remote add nhan-remote https://github.com/nhanvtran/cmssw.git
git fetch nhan-remote puppi-bugfix-for-miniaod
ELE ID
git cms-merge-topic 9003
HEEP implementation
float et = ele.energy()!=0. ? ele.et()/ele.energy()*ele.caloEnergy() : 0.;
float eta = ele.superCluster()->eta();
double iso;
double isoCut;
double rho = *(rho_.product());
double dxy = ( vertices_->size() ? ele.gsfTrack()->dxy(vertices_->at(0).position()) : ele.gsfTrack()->dxy() );
if (ele.gsfTrack().isNonnull()){
if( et > 35. ) {
if( fabs(eta) < 1.4442 ){
iso = ele.dr03EcalRecHitSumEt() + ele.dr03HcalDepth1TowerSumEt();
isoCut = 2 + 0.03*et + 0.28*rho;
if( ele.ecalDriven() == 1 && dEtaInSeed( ele ) < 0.004 && ele.deltaPhiSuperClusterTrackAtVtx() < 0.06 &&
ele.hadronicOverEm() < (2./ele.superCluster()->energy()+0.05) &&
(ele.full5x5_e2x5Max()/ele.full5x5_e5x5() > 0.94 || ele.full5x5_e1x5()/ele.full5x5_e5x5() > 0.83) &&
ele.dr03TkSumPt() < 5. && ele.gsfTrack()->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS) <= 1 &&
iso < isoCut && fabs(dxy) < 0.02 ) isHEEP = true;
}
if( fabs(eta) > 1.566 && fabs(eta) < 2.5 ){
iso = ele.dr03EcalRecHitSumEt() + ele.dr03HcalDepth1TowerSumEt();
if( et <= 50 )
isoCut = 2.5 + 0.28*rho;
else
isoCut = 2.5+0.03*(et-50.) + 0.28*rho;
if( ele.ecalDriven() == 1 && dEtaInSeed( ele ) < 0.006 && ele.deltaPhiSuperClusterTrackAtVtx() < 0.06 &&
ele.hadronicOverEm() < (12.5/ele.superCluster()->energy()+0.05) && ele.full5x5_sigmaIetaIeta() < 0.03 &&
ele.dr03TkSumPt() < 5. && ele.gsfTrack()->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS) <= 1 &&
iso < isoCut && fabs(dxy) < 0.05 ) isHEEP = true;
}
}
}