| 420 |
|
// if(jet_iter->pt()>50) |
| 421 |
|
// njetscounter++; |
| 422 |
|
VHbbEvent::SimpleJet sj; |
| 423 |
+ |
// std::cout <<" sj1"<<std::endl; |
| 424 |
|
fillSimpleJet(sj,jet_iter); |
| 425 |
|
// if(!runOnMC_) |
| 426 |
|
|
| 463 |
|
// if(jet_iter->pt()>50) |
| 464 |
|
// njetscounter++; |
| 465 |
|
VHbbEvent::SimpleJet sj; |
| 466 |
+ |
// std::cout <<" sj3"<<std::endl; |
| 467 |
|
fillSimpleJet(sj,jet_iter); |
| 468 |
|
// if(!runOnMC_) |
| 469 |
|
setJecUnc(sj,jecUnc); |
| 506 |
|
// if(jet_iter->pt()>50) |
| 507 |
|
// njetscounter++; |
| 508 |
|
VHbbEvent::SimpleJet sj; |
| 509 |
+ |
// std::cout <<" sj4"<<std::endl; |
| 510 |
|
fillSimpleJet(sj,jet_iter); |
| 511 |
|
// if(!runOnMC_) |
| 512 |
|
setJecUnc(sj,jecUnc); |
| 549 |
|
for(edm::View<pat::Jet>::const_iterator jet_iter = simplejets2.begin(); jet_iter!=simplejets2.end(); ++jet_iter){ |
| 550 |
|
|
| 551 |
|
VHbbEvent::SimpleJet sj; |
| 552 |
+ |
// std::cout <<" sj2"<<std::endl; |
| 553 |
|
fillSimpleJet(sj,jet_iter); |
| 554 |
|
// if(!runOnMC_) |
| 555 |
|
setJecUnc(sj,jecUnc); |
| 693 |
|
<< "," << subjet_iter->bDiscriminator("combinedSecondaryVertexBJetTags") << "\n";} |
| 694 |
|
|
| 695 |
|
VHbbEvent::SimpleJet sj; |
| 696 |
+ |
// std::cout <<" sub jet "<<std::endl; |
| 697 |
|
fillSimpleJet(sj,subjet_iter); |
| 698 |
|
// if(!runOnMC_) |
| 699 |
|
setJecUnc(sj,jecUnc); |
| 814 |
|
double acop = deltaphi.value(); |
| 815 |
|
mf.acop=acop; |
| 816 |
|
|
| 817 |
+ |
mf.nMatches = mu->numberOfMatches(); |
| 818 |
+ |
|
| 819 |
|
mf.ipDb=mu->dB(); |
| 820 |
|
mf.ipErrDb=mu->edB(); |
| 821 |
|
mf.cat=0; |
| 832 |
|
TrackRef iTrack1 = mu->innerTrack(); |
| 833 |
|
const reco::HitPattern& p1 = iTrack1->hitPattern(); |
| 834 |
|
mf.nPixelHits=p1.pixelLayersWithMeasurement(); |
| 835 |
< |
} |
| 835 |
> |
|
| 836 |
> |
mf.nValidTracker = p1.numberOfValidTrackerHits(); |
| 837 |
> |
mf.nValidPixel = p1.numberOfValidPixelHits(); |
| 838 |
> |
|
| 839 |
> |
|
| 840 |
> |
|
| 841 |
> |
} |
| 842 |
|
if(mu->isGlobalMuon()){ |
| 843 |
|
TrackRef gTrack = mu->globalTrack(); |
| 844 |
|
const reco::HitPattern& q = gTrack->hitPattern(); |
| 910 |
|
Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px())); |
| 911 |
|
ef.acop = deltaphi.value(); |
| 912 |
|
// |
| 913 |
+ |
ef.sihih = elec->sigmaIetaIeta(); |
| 914 |
+ |
ef.Dphi = elec->deltaPhiSuperClusterTrackAtVtx(); |
| 915 |
+ |
ef.Deta = elec->deltaEtaSuperClusterTrackAtVtx(); |
| 916 |
+ |
ef.HoE = elec->hadronicOverEm(); |
| 917 |
+ |
ef.convDist = elec->convDist(); |
| 918 |
+ |
ef.convDcot = elec->convDcot(); |
| 919 |
+ |
if(elec->gsfTrack().isNonnull()) ef.innerHits = elec->gsfTrack()->trackerExpectedHitsInner().numberOfHits(); |
| 920 |
+ |
ef.isEB = elec->isEB(); |
| 921 |
+ |
ef.isEE = elec->isEE(); |
| 922 |
+ |
// |
| 923 |
|
// fill eleids |
| 924 |
|
// |
| 925 |
|
/* ef.id95 = elec->electronID("simpleEleId95cIso"); |
| 1206 |
|
// test |
| 1207 |
|
// |
| 1208 |
|
|
| 1209 |
< |
return; |
| 1209 |
> |
// return; |
| 1210 |
|
double eta = sj.p4.Eta(); |
| 1211 |
|
double pt = sj.p4.Pt(); |
| 1212 |
|
|
| 1226 |
|
sj.jpb=jet_iter->bDiscriminator("jetBProbabilityBJetTags"); |
| 1227 |
|
sj.ssvhe=jet_iter->bDiscriminator("simpleSecondaryVertexHighEffBJetTags"); |
| 1228 |
|
sj.csv=jet_iter->bDiscriminator("combinedSecondaryVertexBJetTags"); |
| 1229 |
< |
sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags"); |
| 1229 |
> |
sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags"); |
| 1230 |
|
sj.charge=jet_iter->jetCharge(); |
| 1231 |
|
sj.ntracks=jet_iter->associatedTracks().size(); |
| 1232 |
|
sj.p4=GENPTOLORP(jet_iter); |
| 1233 |
+ |
// std::cout << " ECCO "<<sj.csv<< " "<< sj.p4.Pt()<<std::endl; |
| 1234 |
|
sj.chargedTracksFourMomentum=(getChargedTracksMomentum(&*(jet_iter))); |
| 1235 |
|
sj.SF_CSVL=1; |
| 1236 |
|
sj.SF_CSVM=1; |
