ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc
(Generate patch)

Comparing UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc (file contents):
Revision 1.62 by degrutto, Fri Mar 16 09:12:46 2012 UTC vs.
Revision 1.65 by arizzi, Fri Mar 30 12:26:25 2012 UTC

# Line 268 | Line 268 | HbbAnalyzerNew::produce(edm::Event& iEve
268          }
269          auxInfo->mcW.push_back(wtemp);
270        }
271 +
272 +      if(abs(id)==15) {
273 +        VHbbEventAuxInfo::ParticleMCInfo tautemp;
274 +        tautemp.status=st;
275 +        tautemp.charge=p.charge();
276 +        if(p.mother(0)!=0) tautemp.momid=p.mother(0)->pdgId();
277 +        if(p.mother(0)!=0 && p.mother(0)->mother(0)!=0) tautemp.gmomid=p.mother(0)->mother(0)->pdgId();
278 +        tautemp.p4=GENPTOLOR(p);
279 +
280 +        int ndau = p.numberOfDaughters();
281 +        for(int j = 0; j < ndau; ++ j) {
282 +          const Candidate * Taudau = p.daughter( j );
283 +          tautemp.dauid.push_back(Taudau->pdgId());
284 +          tautemp.dauFourMomentum.push_back(GENPTOLORP(Taudau));
285 +        }
286 +        auxInfo->mcTau.push_back(tautemp);
287 +      }
288          
289        if(abs(id)==23){
290  
# Line 1493 | Line 1510 | BTagSFContainer btagSFs;
1510      Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
1511      double acop = deltaphi.value();
1512      tf.acop=acop;
1513 <    if (tau->isTauIDAvailable("againstElectronLoose")) tf.idagainstElectronLoose=tau->tauID("againstElectronLoose");
1514 <    if (tau->isTauIDAvailable("againstElectronMedium")) tf.idagainstElectronMedium=tau->tauID("againstElectronMedium");
1515 <    if (tau->isTauIDAvailable("againstElectronTight")) tf.idagainstElectronTight=tau->tauID("againstElectronTight");
1516 <    if (tau->isTauIDAvailable("againstMuonLoose")) tf.idagainstMuonLoose=tau->tauID("againstMuonLoose");
1517 <    if (tau->isTauIDAvailable("againstMuonTight")) tf.idagainstMuonTight=tau->tauID("againstMuonTight");
1518 <    if (tau->isTauIDAvailable("byLooseIsolation")) tf.idbyLooseIsolation=tau->tauID("byLooseIsolation");
1519 <    if (tau->isTauIDAvailable("byMediumIsolation")) tf.idbyMediumIsolation=tau->tauID("byMediumIsolation");
1520 <    if (tau->isTauIDAvailable("byTightIsolation")) tf.idbyTightIsolation=tau->tauID("byTightIsolation");
1521 <    if (tau->isTauIDAvailable("byVLooseIsolation")) tf.idbyVLooseIsolation=tau->tauID("byVLooseIsolation");
1522 <    if (tau->isTauIDAvailable("decayModeFinding")) tf.iddecayModeFinding=tau->tauID("decayModeFinding");
1523 <    if (tau->isTauIDAvailable("byIsolation")) tf.idbyIso=tau->tauID("byIsolation");
1524 <    if (tau->isTauIDAvailable("trackIsolation")) tf.idbyTrackIso=tau->tauID("trackIsolation");
1525 <    if (tau->isTauIDAvailable("byTaNCfrOnePercent")) tf.idbyTaNCfrOnePercent=tau->tauID("byTaNCfrOnePercent");
1526 <    if (tau->isTauIDAvailable("byTaNCfrHalfPercent")) tf.idbyTaNCfrHalfPercent=tau->tauID("byTaNCfrHalfPercent");
1527 <    if (tau->isTauIDAvailable("byTaNCfrQuarterPercent")) tf.idbyTaNCfrQuarterPercent=tau->tauID("byTaNCfrQuarterPercent");
1528 <    if (tau->isTauIDAvailable("byTaNCfrTenthPercent")) tf.idbyTaNCfrTenthPercent=tau->tauID("byTaNCfrTenthPercent");
1529 <    if (tau->isTauIDAvailable("byTaNC")) tf.idbyTaNC=tau->tauID("byTaNC");
1513 >    if (tau->isTauIDAvailable("againstElectronLoose")) tf.againstElectronLoose=tau->tauID("againstElectronLoose");
1514 >    if (tau->isTauIDAvailable("againstElectronMedium")) tf.againstElectronMedium=tau->tauID("againstElectronMedium");
1515 >    if (tau->isTauIDAvailable("againstElectronTight")) tf.againstElectronTight=tau->tauID("againstElectronTight");
1516 >    if (tau->isTauIDAvailable("againstMuonLoose")) tf.againstMuonLoose=tau->tauID("againstMuonLoose");
1517 >    if (tau->isTauIDAvailable("againstMuonTight")) tf.againstMuonTight=tau->tauID("againstMuonTight");
1518 >    if (tau->isTauIDAvailable("byLooseIsolation")) tf.byLooseIsolation=tau->tauID("byLooseIsolation");
1519 >    if (tau->isTauIDAvailable("byMediumIsolation")) tf.byMediumIsolation=tau->tauID("byMediumIsolation");
1520 >    if (tau->isTauIDAvailable("byTightIsolation")) tf.byTightIsolation=tau->tauID("byTightIsolation");
1521 >    if (tau->isTauIDAvailable("byVLooseIsolation")) tf.byVLooseIsolation=tau->tauID("byVLooseIsolation");
1522 >    if (tau->isTauIDAvailable("decayModeFinding")) tf.decayModeFinding=tau->tauID("decayModeFinding");
1523 >    if (tau->isTauIDAvailable("byIsolation")) tf.byIsolation=tau->tauID("byIsolation");
1524 >    if (tau->isTauIDAvailable("trackIsolation")) tf.trackIsolation=tau->tauID("trackIsolation");
1525 >    if (tau->isTauIDAvailable("byTaNCfrOnePercent")) tf.byTaNCfrOnePercent=tau->tauID("byTaNCfrOnePercent");
1526 >    if (tau->isTauIDAvailable("byTaNCfrHalfPercent")) tf.byTaNCfrHalfPercent=tau->tauID("byTaNCfrHalfPercent");
1527 >    if (tau->isTauIDAvailable("byTaNCfrQuarterPercent")) tf.byTaNCfrQuarterPercent=tau->tauID("byTaNCfrQuarterPercent");
1528 >    if (tau->isTauIDAvailable("byTaNCfrTenthPercent")) tf.byTaNCfrTenthPercent=tau->tauID("byTaNCfrTenthPercent");
1529 >    if (tau->isTauIDAvailable("byTaNC")) tf.byTaNC=tau->tauID("byTaNC");
1530 >    if (tau->isTauIDAvailable("byLooseCombinedIsolationDeltaBetaCorr")) tf.byLooseCombinedIsolationDeltaBetaCorr=tau->tauID("byLooseCombinedIsolationDeltaBetaCorr");
1531 >    if (tau->isTauIDAvailable("againstElectronMVA")) tf.againstElectronMVA=tau->tauID("againstElectronMVA");
1532      if (tau->isPFTau()) {
1533        tf.isolationPFChargedHadrCandsPtSum = tau->isolationPFChargedHadrCandsPtSum();
1534        tf.isolationPFGammaCandsEtSum = tau->isolationPFGammaCandsEtSum();
1535 <      if (tau->leadPFChargedHadrCand().isAvailable()) tf.leadPFChargedHadrCandPt = tau->leadPFChargedHadrCand()->pt();
1535 >      if (tau->leadPFChargedHadrCand().isAvailable()) tf.leadPFChargedHadrCandPt = tau->leadPFChargedHadrCand()->pt();
1536        tf.NsignalPFChargedHadrCands = tau->signalPFChargedHadrCands().size();
1537        tf.NsignalPFGammaCands = tau->signalPFGammaCands().size();
1538      }
1539      hbbInfo->tauInfo.push_back(tf);
1540      if (verbose_) {
1541 <      std::cout << "SCZ DEBUG: againstElectronLoose is " << tf.idagainstElectronLoose << std::endl;
1542 <      std::cout << "SCZ DEBUG: againstElectronMedium is " << tf.idagainstElectronMedium << std::endl;
1543 <      std::cout << "SCZ DEBUG: againstElectronTight is " << tf.idagainstElectronTight << std::endl;
1544 <      std::cout << "SCZ DEBUG: againstMuonLoose is " << tf.idagainstMuonLoose << std::endl;
1545 <      std::cout << "SCZ DEBUG: againstMuonTight is " << tf.idagainstMuonTight << std::endl;
1546 <      std::cout << "SCZ DEBUG: byLooseIsolation is " << tf.idbyLooseIsolation << std::endl;
1547 <      std::cout << "SCZ DEBUG: byMediumIsolation is " << tf.idbyMediumIsolation << std::endl;
1548 <      std::cout << "SCZ DEBUG: byTightIsolation is " << tf.idbyTightIsolation << std::endl;
1549 <      std::cout << "SCZ DEBUG: byVLooseIsolation is " << tf.idbyVLooseIsolation << std::endl;
1550 <      std::cout << "SCZ DEBUG: decayModeFinding is " << tf.iddecayModeFinding << std::endl;
1551 <      std::cout << "SCZ DEBUG: byIsolation is " << tf.idbyIso<< std::endl;
1552 <      std::cout << "SCZ DEBUG: trackIsolation is " << tf.idbyTrackIso << std::endl;
1553 <      std::cout << "SCZ DEBUG: byTaNCfrOnePercent is " << tf.idbyTaNCfrOnePercent << std::endl;
1554 <      std::cout << "SCZ DEBUG: byTaNCfrHalfPercent is " << tf.idbyTaNCfrHalfPercent << std::endl;
1555 <      std::cout << "SCZ DEBUG: byTaNCfrQuarterPercent is " << tf.idbyTaNCfrQuarterPercent << std::endl;
1556 <      std::cout << "SCZ DEBUG: byTaNCfrTenthPercent is " << tf.idbyTaNCfrTenthPercent << std::endl;
1557 <      std::cout << "SCZ DEBUG: byTaNC is " << tf.idbyTaNC << std::endl;
1541 >      std::cout << "SCZ DEBUG: againstElectronLoose is " << tf.againstElectronLoose << std::endl;
1542 >      std::cout << "SCZ DEBUG: againstElectronMedium is " << tf.againstElectronMedium << std::endl;
1543 >      std::cout << "SCZ DEBUG: againstElectronTight is " << tf.againstElectronTight << std::endl;
1544 >      std::cout << "SCZ DEBUG: againstMuonLoose is " << tf.againstMuonLoose << std::endl;
1545 >      std::cout << "SCZ DEBUG: againstMuonTight is " << tf.againstMuonTight << std::endl;
1546 >      std::cout << "SCZ DEBUG: byLooseIsolation is " << tf.byLooseIsolation << std::endl;
1547 >      std::cout << "SCZ DEBUG: byMediumIsolation is " << tf.byMediumIsolation << std::endl;
1548 >      std::cout << "SCZ DEBUG: byTightIsolation is " << tf.byTightIsolation << std::endl;
1549 >      std::cout << "SCZ DEBUG: byVLooseIsolation is " << tf.byVLooseIsolation << std::endl;
1550 >      std::cout << "SCZ DEBUG: decayModeFinding is " << tf.decayModeFinding << std::endl;
1551 >      std::cout << "SCZ DEBUG: byIsolation is " << tf.byIsolation<< std::endl;
1552 >      std::cout << "SCZ DEBUG: trackIsolation is " << tf.trackIsolation << std::endl;
1553 >      std::cout << "SCZ DEBUG: byTaNCfrOnePercent is " << tf.byTaNCfrOnePercent << std::endl;
1554 >      std::cout << "SCZ DEBUG: byTaNCfrHalfPercent is " << tf.byTaNCfrHalfPercent << std::endl;
1555 >      std::cout << "SCZ DEBUG: byTaNCfrQuarterPercent is " << tf.byTaNCfrQuarterPercent << std::endl;
1556 >      std::cout << "SCZ DEBUG: byTaNCfrTenthPercent is " << tf.byTaNCfrTenthPercent << std::endl;
1557 >      std::cout << "SCZ DEBUG: byTaNC is " << tf.byTaNC << std::endl;
1558        std::cout << "SCZ DEBUG: isolationPFChargedHadrCandsPtSum is " << tf.isolationPFChargedHadrCandsPtSum << std::endl;
1559        std::cout << "SCZ DEBUG: isolationPFGammaCandsEtSum is " << tf.isolationPFGammaCandsEtSum << std::endl;
1560        std::cout << "SCZ DEBUG: isolationPFGammaCandsEtSum is " << tf.leadPFChargedHadrCandPt << std::endl;
1561        std::cout << "SCZ DEBUG: NsignalPFChargedHadrCands is " << tf.NsignalPFChargedHadrCands << std::endl;
1562        std::cout << "SCZ DEBUG: NsignalPFGammaCands is " << tf.NsignalPFGammaCands << std::endl;
1563 +      std::cout << "SCZ DEBUG: byLooseCombinedIsolationDeltaBetaCorr is " << tf.byLooseCombinedIsolationDeltaBetaCorr << std::endl;
1564 +      std::cout << "SCZ DEBUG: againstElectronMVA is " << tf.againstElectronMVA << std::endl;
1565      }
1566    }
1567  
# Line 1784 | Line 1805 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1805      sj.jpb=jet_iter->bDiscriminator("jetBProbabilityBJetTags");
1806      sj.ssvhe=jet_iter->bDiscriminator("simpleSecondaryVertexHighEffBJetTags");
1807      sj.csv=jet_iter->bDiscriminator("combinedSecondaryVertexBJetTags");
1808 <   sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1808 >    sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1809 >    sj.csvivf=jet_iter->bDiscriminator("combinedInclusiveSecondaryVertexBJetTags");
1810 >    sj.cmva=jet_iter->bDiscriminator("combinedMVABJetTags");
1811      sj.charge=jet_iter->jetCharge();
1812      sj.ntracks=jet_iter->associatedTracks().size();
1813      sj.p4=GENPTOLORP(jet_iter);

Diff Legend

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