Definition at line 634 of file myFastSimVal.cc.
References funct::abs(), edm::SortedCollection< T, SORT >::begin(), CaloJetAlgorithm1, CaloJetAlgorithm2, CaloJetAlgorithm3, CaloJetAlgorithm4, isolatedTracks_cfi::caloJets, eleHcalExtractorBlocks_cff::caloTowers, funct::cos(), gather_cfg::cout, reco::Candidate::daughter(), SiPixelRawToDigiRegional_cfi::deltaPhi, PbPb_ZMuSkimMuonDPG_cff::deltaR, dEta1, dEta12, dEta2, dEta3, dEta4, dEtaPar1, dEtaPar2, dEtaPar3, dijetMass1, dijetMass101, dijetMass102, dijetMass103, dijetMass12, dijetMass13, dijetMass131, dijetMass132, dijetMass133, dijetMass2, dijetMass22, dijetMass23, dijetMass3, dijetMass32, dijetMass33, dijetMass4, dijetMass401, dijetMass402, dijetMass403, dijetMass42, dijetMass43, dijetMass_2000_101, dijetMass_2000_131, dijetMass_2000_401, dijetMass_5000_101, dijetMass_5000_131, dijetMass_5000_401, dijetMass_700_101, dijetMass_700_131, dijetMass_700_401, dijetMassCor1, dijetMassCor101, dijetMassCor131, dijetMassCor401, dijetMassCor_2000_1, dijetMassCor_2000_101, dijetMassCor_2000_131, dijetMassCor_2000_401, dijetMassCor_5000_1, dijetMassCor_5000_101, dijetMassCor_5000_131, dijetMassCor_5000_401, dijetMassCor_700_1, dijetMassCor_700_101, dijetMassCor_700_131, dijetMassCor_700_401, dijetMassP1, dijetMassP101, dijetMassP131, dijetMassP2, dijetMassP3, dijetMassP401, dijetMassP_2000_101, dijetMassP_2000_131, dijetMassP_5000_101, dijetMassP_5000_131, dijetMassP_5000_401, dijetMassP_700_101, dijetMassP_700_131, dPhi1, dPhi12, dPhi2, dPhi3, dPhi4, dPhiPar1, dPhiPar2, dPhiPar3, dPt1, dPt100Frac1, dPt100Frac2, dPt100Frac3, dPt100Frac4, dPt12, dPt2, dPt20Frac1, dPt20Frac2, dPt20Frac3, dPt20Frac4, dPt3, dPt4, dPt40Frac1, dPt40Frac2, dPt40Frac3, dPt40Frac4, dPt80Frac1, dPt80Frac2, dPt80Frac3, dPt80Frac4, dPtFrac1, dPtFrac2, dPtFrac3, dPtFrac4, dPtPar1, dPtPar2, dPtPar3, dR1, dR12, dR2, dR3, dR4, dRPar1, dRPar2, dRPar3, dRParton, dRPartonMin, emEneLeadJet1, emEneLeadJet12, emEneLeadJet13, emEneLeadJetEta1_1, emEneLeadJetEta1_2, emEneLeadJetEta2_1, emEneLeadJetEta2_2, emEneLeadJetEta3_1, emEneLeadJetEta3_2, edm::SortedCollection< T, SORT >::end(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, relval_steps::gen(), GenJetAlgorithm1, GenJetAlgorithm2, GenJetAlgorithm3, GenJetAlgorithm4, ttbarCategorization_cff::genJets, genParticles2HepMC_cfi::genParticles, edm::Event::getByLabel(), h_etaCal1, h_etaCal2, h_etaCal3, h_etaCal4, h_etaCalL1, h_etaCalL2, h_etaCalL3, h_etaCalL4, h_etaGen1, h_etaGen2, h_etaGen3, h_etaGenL1, h_etaGenL2, h_etaGenL3, h_jet1Pt1, h_jet1Pt2, h_jet1Pt3, h_jet1Pt4, h_jet2Pt1, h_jet2Pt2, h_jet2Pt3, h_jet2Pt4, h_jet3Pt1, h_jet3Pt2, h_jet3Pt3, h_jet3Pt4, h_jet4Pt1, h_jet4Pt2, h_jet4Pt3, h_jet4Pt4, h_jet5Pt1, h_jet5Pt2, h_jet5Pt3, h_jet5Pt4, h_jet6Pt1, h_jet6Pt2, h_jet6Pt3, h_jet6Pt4, h_jet7Pt1, h_jet7Pt2, h_jet7Pt3, h_jet7Pt4, h_jetEt1, h_jetEt2, h_jetEt3, h_lowPtCal11, h_lowPtCal12, h_lowPtCal13, h_lowPtCal1c11, h_lowPtCal1c12, h_lowPtCal1c13, h_lowPtCal21, h_lowPtCal22, h_lowPtCal23, h_lowPtCal2c11, h_lowPtCal2c12, h_lowPtCal2c13, h_lowPtCal31, h_lowPtCal32, h_lowPtCal33, h_lowPtCal3c11, h_lowPtCal3c12, h_lowPtCal3c13, h_lowPtCal41, h_lowPtCal42, h_lowPtCal43, h_lowPtCal4c11, h_lowPtCal4c12, h_lowPtCal4c13, h_missEt1, h_missEt1s, h_missEt2, h_missEt2s, h_missEt3, h_missEt3s, h_nCalJets1, h_nCalJets2, h_nCalJets3, h_nCalJets4, h_nGenJets1, h_nGenJets2, h_nGenJets3, h_phiCal1, h_phiCal2, h_phiCal3, h_phiCal4, h_phiCalL1, h_phiCalL2, h_phiCalL3, h_phiCalL4, h_phiGen1, h_phiGen2, h_phiGen3, h_phiGenL1, h_phiGenL2, h_phiGenL3, h_ptCal1, h_ptCal12, h_ptCal13, h_ptCal2, h_ptCal22, h_ptCal23, h_ptCal3, h_ptCal32, h_ptCal33, h_ptCal4, h_ptCal42, h_ptCal43, h_ptCalL1, h_ptCalL12, h_ptCalL13, h_ptCalL2, h_ptCalL22, h_ptCalL23, h_ptCalL3, h_ptCalL32, h_ptCalL33, h_ptCalL4, h_ptCalL42, h_ptCalL43, h_ptGen1, h_ptGen12, h_ptGen13, h_ptGen2, h_ptGen22, h_ptGen23, h_ptGen3, h_ptGen32, h_ptGen33, h_ptGenL1, h_ptGenL12, h_ptGenL13, h_ptGenL2, h_ptGenL22, h_ptGenL23, h_ptGenL3, h_ptGenL32, h_ptGenL33, h_totMissEt1, h_totMissEt2, h_totMissEt3, hadEneLeadJet1, hadEneLeadJet12, hadEneLeadJet13, hadEneLeadJetEta1_1, hadEneLeadJetEta1_2, hadEneLeadJetEta2_1, hadEneLeadJetEta2_2, hadEneLeadJetEta3_1, hadEneLeadJetEta3_2, hadFracEta11, hadFracEta21, hadFracEta31, hadFracLeadJet1, hf_nJet1, hf_nJet11, hf_nJet1s, hf_nJet2, hf_nJet21, hf_nJet2s, hf_nJet3, hf_nJet31, hf_nJet3s, hf_nJet4, hf_nJet41, hf_nJet4s, hf_PtResponse1, hf_PtResponse2, hf_PtResponse3, hf_PtResponse4, hf_sumTowerAllEx, hf_sumTowerAllEy, hf_TowerDelR1, hf_TowerDelR12, hf_TowerJetEt1, mps_fire::i, globals_cff::id1, globals_cff::id2, cuy::ii, dqmiolumiharvest::j, reco::btau::jetPhi, reco::btau::jetPt, singleTopDQM_cfi::jets, findQualityFiles::jj, EgHLTOffHistBins_cfi::mass, reco::Candidate::mass(), matchedAllPt11, matchedAllPt12, matchedAllPt13, matchedAllPt21, matchedAllPt22, matchedAllPt23, matchedAllPt31, matchedAllPt32, matchedAllPt33, matchedAllPt41, matchedAllPt42, matchedAllPt43, matchedPt11, matchedPt12, matchedPt13, matchedPt21, matchedPt22, matchedPt23, matchedPt31, matchedPt32, matchedPt33, matchedPt41, matchedPt42, matchedPt43, MAXJETS, MET1, MET12, MET13, ALCARECOTkAlWMuNu_cff::minJetPt, reco::Candidate::mother(), jets_cff::nConstituents, nEvent, njet, nPart(), nTowers1, nTowers2, nTowers3, nTowers4, nTowersLeadJet1, nTowersLeadJetPt1, nTowersLeadJetPt2, nTowersLeadJetPt3, nTowersLeadJetPt4, nTowersSecondJet1, reco::Candidate::numberOfDaughters(), AlCaHLTBitMon_ParallelJobs::p, reco::Candidate::p4(), ParMatch1, ParMatch2, ParMatch3, EgammaValidation_cff::pdgId, reco::Candidate::pdgId(), phi, DiDispStaMuonMonitor_cfi::pt, Scenarios_cff::scale, funct::sin(), edm::OwnVector< T, P >::size(), mathSSE::sqrt(), SumEt1, SumEt12, SumEt13, OrderedSet::t, tbarMass, tbarMassGen, tMass, tMassGen, topMass1, topMass2, topMass3, topMassParton, totEneLeadJetEta1_1, totEneLeadJetEta1_2, totEneLeadJetEta2_1, totEneLeadJetEta2_2, totEneLeadJetEta3_1, totEneLeadJetEta3_2, hgcalTowerProducer_cfi::tower, TowerEtLeadJet1, TowerEtLeadJet12, TowerEtLeadJet13, bphysicsOniaDQM_cfi::vertex, reco::Candidate::vx(), reco::Candidate::vy(), reco::Candidate::vz(), ZpMass, ZpMass_2000_10, ZpMass_2000_13, ZpMass_2000_40, ZpMass_5000_10, ZpMass_5000_13, ZpMass_5000_40, ZpMass_700_10, ZpMass_700_13, ZpMass_700_40, ZpMassGen, ZpMassGen10, ZpMassGen13, ZpMassGen40, ZpMassGen_2000_10, ZpMassGen_2000_13, ZpMassGen_2000_40, ZpMassGen_5000_10, ZpMassGen_5000_13, ZpMassGen_5000_40, ZpMassGen_700_10, ZpMassGen_700_13, ZpMassGen_700_40, ZpMassMatched1, ZpMassMatched2, ZpMassMatched3, ZpMassRes101, ZpMassRes131, ZpMassRes401, ZpMassRes402, ZpMassRes403, ZpMassResL101, ZpMassResL102, ZpMassResL103, ZpMassResL131, ZpMassResL132, ZpMassResL133, ZpMassResL401, ZpMassResL402, ZpMassResL403, ZpMassResPRL101, ZpMassResPRL102, ZpMassResPRL103, ZpMassResPRL131, ZpMassResPRL132, ZpMassResPRL133, ZpMassResPRL401, ZpMassResPRL402, ZpMassResPRL403, ZpMassResRL101, ZpMassResRL102, ZpMassResRL103, ZpMassResRL131, ZpMassResRL132, ZpMassResRL133, ZpMassResRL401, ZpMassResRL402, ZpMassResRL403, ZpMassResRLoP101, ZpMassResRLoP102, ZpMassResRLoP103, ZpMassResRLoP131, ZpMassResRLoP132, ZpMassResRLoP133, ZpMassResRLoP401, ZpMassResRLoP402, and ZpMassResRLoP403.
635 int EtaOk10, EtaOk13, EtaOk40;
637 double ZpM, ZpMG, ZpMM;
638 double LeadMass1, LeadMass2, LeadMass3, LeadMass4;
639 double LeadMassP1, LeadMassP2, LeadMassP3;
642 float minJetPt10 = 10.;
643 int jetInd, allJetInd;
646 double matchedDelR = 0.3;
664 for (
int istep = 0; istep < 100; ++istep) {
666 float ptStep = (istep * (5000. / 100.));
667 for (CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++cal) {
668 if (cal->pt() > ptStep)
676 for (
int istep = 0; istep < 100; ++istep) {
678 float ptStep = (istep * (200. / 100.));
679 for (CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++cal) {
680 if (cal->pt() > ptStep)
688 for (
int istep = 0; istep < 100; ++istep) {
690 float ptStep = (istep * (3000. / 100.));
691 for (CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++cal) {
692 if (cal->pt() > ptStep)
715 for (CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++cal) {
718 double corPt = scale * cal->pt();
721 if (corPt > highestPt) {
723 p4cortmp[1] = p4cortmp[0];
725 p4cortmp[0] = scale * cal->p4();
726 }
else if (corPt > nextPt) {
728 p4cortmp[1] = scale * cal->p4();
741 if (allJetInd == 1) {
743 p4tmp[0] = cal->p4();
744 if (fabs(cal->eta()) < 1.0)
746 if (fabs(cal->eta()) < 1.3)
748 if (fabs(cal->eta()) < 4.0)
751 if (allJetInd == 2) {
753 p4tmp[1] = cal->p4();
754 if (fabs(cal->eta()) < 1.0)
756 if (fabs(cal->eta()) < 1.3)
758 if (fabs(cal->eta()) < 4.0)
761 if ((allJetInd == 1) || (allJetInd == 2)) {
781 if (fabs(cal->eta()) < 1.3) {
783 if (cal->pt() > 10.) {
787 if ((fabs(cal->eta()) > 1.3) && (fabs(cal->eta()) < 3.)) {
789 if (cal->pt() > 10.) {
793 if (fabs(cal->eta()) > 3.0) {
795 if (cal->pt() > 10.) {
815 LeadMass1 = (p4tmp[0] + p4tmp[1]).
mass();
844 LeadMass1 = (p4cortmp[0] + p4cortmp[1]).
mass();
897 for (
int istep = 0; istep < 100; ++istep) {
899 float ptStep = (istep * (5000. / 100.));
901 for (CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++cal)
902 if (cal->pt() > ptStep)
908 for (
int istep = 0; istep < 100; ++istep) {
910 float ptStep = (istep * (200. / 100.));
912 for (CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++cal)
913 if (cal->pt() > ptStep)
919 for (
int istep = 0; istep < 100; ++istep) {
921 float ptStep = (istep * (3000. / 100.));
923 for (CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++cal)
924 if (cal->pt() > ptStep)
933 for (CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++cal) {
935 if (allJetInd == 1) {
937 p4tmp[0] = cal->p4();
939 if (allJetInd == 2) {
941 p4tmp[1] = cal->p4();
943 if ((allJetInd == 1) || (allJetInd == 2)) {
962 if (fabs(cal->eta()) < 1.3) {
964 if (cal->pt() > 10.) {
968 if ((fabs(cal->eta()) > 1.3) && (fabs(cal->eta()) < 3.)) {
970 if (cal->pt() > 10.) {
974 if (fabs(cal->eta()) > 3.0) {
976 if (cal->pt() > 10.) {
994 LeadMass2 = (p4tmp[0] + p4tmp[1]).
mass();
1035 for (
int istep = 0; istep < 100; ++istep) {
1037 float ptStep = (istep * (5000. / 100.));
1039 for (CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++cal)
1040 if (cal->pt() > ptStep)
1046 for (
int istep = 0; istep < 100; ++istep) {
1048 float ptStep = (istep * (200. / 100.));
1050 for (CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++cal)
1051 if (cal->pt() > ptStep)
1057 for (
int istep = 0; istep < 100; ++istep) {
1059 float ptStep = (istep * (3000. / 100.));
1061 for (CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++cal)
1062 if (cal->pt() > ptStep)
1068 for (CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++cal) {
1070 if (allJetInd == 1) {
1072 p4tmp[0] = cal->p4();
1074 if (allJetInd == 2) {
1076 p4tmp[1] = cal->p4();
1078 if ((allJetInd == 1) || (allJetInd == 2)) {
1097 if (fabs(cal->eta()) < 1.3) {
1099 if (cal->pt() > 10.) {
1103 if ((fabs(cal->eta()) > 1.3) && (fabs(cal->eta()) < 3.)) {
1105 if (cal->pt() > 10.) {
1109 if (fabs(cal->eta()) > 3.0) {
1111 if (cal->pt() > 10.) {
1130 LeadMass3 = (p4tmp[0] + p4tmp[1]).
mass();
1148 for (
int istep = 0; istep < 100; ++istep) {
1150 float ptStep = (istep * (5000. / 100.));
1152 for (CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++cal)
1153 if (cal->pt() > ptStep)
1159 for (
int istep = 0; istep < 100; ++istep) {
1161 float ptStep = (istep * (200. / 100.));
1163 for (CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++cal)
1164 if (cal->pt() > ptStep)
1170 for (
int istep = 0; istep < 100; ++istep) {
1172 float ptStep = (istep * (3000. / 100.));
1174 for (CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++cal)
1175 if (cal->pt() > ptStep)
1181 for (CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++cal) {
1183 if (allJetInd == 1) {
1185 p4tmp[0] = cal->p4();
1187 if (allJetInd == 2) {
1189 p4tmp[1] = cal->p4();
1191 if ((allJetInd == 1) || (allJetInd == 2)) {
1210 if (fabs(cal->eta()) < 1.3) {
1212 if (cal->pt() > 10.) {
1216 if ((fabs(cal->eta()) > 1.3) && (fabs(cal->eta()) < 3.)) {
1218 if (cal->pt() > 10.) {
1222 if (fabs(cal->eta()) > 3.0) {
1224 if (cal->pt() > 10.) {
1243 LeadMass4 = (p4tmp[0] + p4tmp[1]).
mass();
1259 for (GenJetCollection::const_iterator
gen = genJets1->begin();
gen != genJets1->end(); ++
gen) {
1261 if (allJetInd == 1) {
1262 p4tmp[0] =
gen->p4();
1264 if (allJetInd == 2) {
1265 p4tmp[1] =
gen->p4();
1268 if ((allJetInd == 1) || (allJetInd == 2)) {
1289 LeadMassP1 = (p4tmp[0] + p4tmp[1]).
mass();
1320 for (GenJetCollection::const_iterator
gen = genJets2->begin();
gen != genJets2->end(); ++
gen) {
1322 if (allJetInd == 1) {
1323 p4tmp[0] =
gen->p4();
1325 if (allJetInd == 2) {
1326 p4tmp[1] =
gen->p4();
1328 if ((allJetInd == 1) || (allJetInd == 2)) {
1349 LeadMassP2 = (p4tmp[0] + p4tmp[1]).
mass();
1362 for (GenJetCollection::const_iterator
gen = genJets3->begin();
gen != genJets3->end(); ++
gen) {
1364 if (allJetInd == 1) {
1365 p4tmp[0] =
gen->p4();
1367 if (allJetInd == 2) {
1368 p4tmp[1] =
gen->p4();
1370 if ((allJetInd == 1) || (allJetInd == 2)) {
1391 LeadMassP3 = (p4tmp[0] + p4tmp[1]).
mass();
1417 for (
int i = 0;
i < maxJets; ++
i)
1422 for (GenJetCollection::const_iterator
gen = genJets->begin();
gen != genJets->end() && jetInd < maxJets; ++
gen) {
1423 p4gen[jetInd] =
gen->p4();
1424 dRmin[jetInd] = 1000.0;
1429 for (CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++cal) {
1430 double delR =
deltaR(cal->eta(), cal->phi(),
gen->eta(),
gen->phi());
1432 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
1433 dRmin[jetInd] = delR;
1434 p4cal[jetInd] = cal->p4();
1441 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1444 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1447 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1451 if (usedInd != -1) {
1454 if (p4cal[jetInd].
pt() > minJetPt10)
1457 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1460 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1463 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1478 double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1479 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1480 if ((p4gen[jetInd].
pt() > 20) && (p4gen[jetInd].
pt() < 40)) {
1483 if ((p4gen[jetInd].
pt() > 40) && (p4gen[jetInd].
pt() < 60)) {
1486 if ((p4gen[jetInd].
pt() > 80) && (p4gen[jetInd].
pt() < 100)) {
1489 if ((p4gen[jetInd].
pt() > 100) && (p4gen[jetInd].
pt() < 120)) {
1493 if ((p4gen[jetInd].
pt() > minJetPt10) && (fabs(p4gen[jetInd].
eta()) < 1.3)) {
1494 dR1.Fill(dRmin[jetInd]);
1495 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4gen[jetInd].
phi());
1497 double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1502 if (((dpt / p4gen[jetInd].
pt()) < -0.5) && (fabs(dpt) > 90.)) {
1503 cout <<
" deltaR min = " << dRmin[jetInd] <<
" Ind = " << jetInd <<
" / " << usedInd <<
" / " << used[nj]
1504 <<
" Del pt / frac = " << dpt <<
" / " << dpt / p4gen[jetInd].pt()
1505 <<
" cal/gen pt = " << p4cal[jetInd].pt() <<
" / " << p4gen[jetInd].pt()
1506 <<
" cal/gen eta = " << p4cal[jetInd].eta() <<
" / " << p4gen[jetInd].eta()
1507 <<
" cal/gen phi = " << p4cal[jetInd].phi() <<
" / " << p4gen[jetInd].phi() << endl;
1531 for (
int i = 0;
i < maxJets; ++
i)
1533 for (GenJetCollection::const_iterator
gen = genJets->begin();
gen != genJets->end() && jetInd < maxJets; ++
gen) {
1534 p4gen[jetInd] =
gen->p4();
1535 dRmin[jetInd] = 1000.0;
1540 for (CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++cal) {
1541 double delR =
deltaR(cal->eta(), cal->phi(),
gen->eta(),
gen->phi());
1543 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
1544 dRmin[jetInd] = delR;
1545 p4cal[jetInd] = cal->p4();
1551 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1554 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1557 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1561 if (usedInd != -1) {
1564 if (p4cal[jetInd].
pt() > minJetPt10)
1567 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1570 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1573 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1588 double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1589 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1590 if ((p4gen[jetInd].
pt() > 20) && (p4gen[jetInd].
pt() < 40)) {
1593 if ((p4gen[jetInd].
pt() > 40) && (p4gen[jetInd].
pt() < 60)) {
1596 if ((p4gen[jetInd].
pt() > 80) && (p4gen[jetInd].
pt() < 100)) {
1599 if ((p4gen[jetInd].
pt() > 100) && (p4gen[jetInd].
pt() < 120)) {
1603 if ((p4gen[jetInd].
pt() > minJetPt10) && (fabs(p4gen[jetInd].
eta()) < 1.3)) {
1604 dR2.Fill(dRmin[jetInd]);
1605 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4gen[jetInd].
phi());
1607 double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1633 for (
int i = 0;
i < maxJets; ++
i)
1635 for (GenJetCollection::const_iterator
gen = genJets->begin();
gen != genJets->end() && jetInd < maxJets; ++
gen) {
1636 p4gen[jetInd] =
gen->p4();
1637 dRmin[jetInd] = 1000.0;
1642 for (CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++cal) {
1643 double delR =
deltaR(cal->eta(), cal->phi(),
gen->eta(),
gen->phi());
1645 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
1646 dRmin[jetInd] = delR;
1647 p4cal[jetInd] = cal->p4();
1653 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1656 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1659 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1663 if (usedInd != -1) {
1666 if (p4cal[jetInd].
pt() > minJetPt10)
1669 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1672 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1675 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1690 double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1691 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1692 if ((p4gen[jetInd].
pt() > 20) && (p4gen[jetInd].
pt() < 40)) {
1695 if ((p4gen[jetInd].
pt() > 40) && (p4gen[jetInd].
pt() < 60)) {
1698 if ((p4gen[jetInd].
pt() > 80) && (p4gen[jetInd].
pt() < 100)) {
1701 if ((p4gen[jetInd].
pt() > 100) && (p4gen[jetInd].
pt() < 120)) {
1706 if ((p4gen[jetInd].
pt() > minJetPt10) && (fabs(p4gen[jetInd].
eta()) < 1.3)) {
1707 dR3.Fill(dRmin[jetInd]);
1708 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4gen[jetInd].
phi());
1710 double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1731 for (
int i = 0;
i < maxJets; ++
i)
1733 for (GenJetCollection::const_iterator
gen = genJets->begin();
gen != genJets->end() && jetInd < maxJets; ++
gen) {
1734 p4gen[jetInd] =
gen->p4();
1735 dRmin[jetInd] = 1000.0;
1740 for (CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++cal) {
1741 double delR =
deltaR(cal->eta(), cal->phi(),
gen->eta(),
gen->phi());
1743 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
1744 dRmin[jetInd] = delR;
1745 p4cal[jetInd] = cal->p4();
1751 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1754 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1757 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1761 if (usedInd != -1) {
1764 if (p4cal[jetInd].
pt() > minJetPt10)
1767 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1770 if ((fabs(p4gen[jetInd].
eta()) > 1.3) && (fabs(p4gen[jetInd].
eta()) < 3.)) {
1773 if (fabs(p4gen[jetInd].
eta()) > 3.) {
1789 double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1790 if (fabs(p4gen[jetInd].
eta()) < 1.3) {
1791 if ((p4gen[jetInd].
pt() > 20) && (p4gen[jetInd].
pt() < 40)) {
1794 if ((p4gen[jetInd].
pt() > 40) && (p4gen[jetInd].
pt() < 60)) {
1797 if ((p4gen[jetInd].
pt() > 80) && (p4gen[jetInd].
pt() < 100)) {
1800 if ((p4gen[jetInd].
pt() > 100) && (p4gen[jetInd].
pt() < 120)) {
1805 if ((p4gen[jetInd].
pt() > minJetPt10) && (fabs(p4gen[jetInd].
eta()) < 1.3)) {
1806 dR4.Fill(dRmin[jetInd]);
1807 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4gen[jetInd].
phi());
1809 double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1832 for (
int i = 0;
i < maxJets; ++
i)
1834 for (CaloJetCollection::const_iterator cal1 = calo1Jets->begin(); cal1 != calo1Jets->end() && jetInd < maxJets;
1836 p4gen[jetInd] = cal1->p4();
1837 dRmin[jetInd] = 1000.0;
1840 for (CaloJetCollection::const_iterator cal2 = calo2Jets->begin(); cal2 != calo2Jets->end(); ++cal2) {
1841 double delR =
deltaR(cal1->eta(), cal1->phi(), cal2->eta(), cal2->phi());
1842 if ((delR < dRmin[jetInd]) && (used[nj] == 0)) {
1843 dRmin[jetInd] = delR;
1844 p4cal[jetInd] = cal2->p4();
1851 if (p4gen[jetInd].
pt() > minJetPt) {
1852 dR12.Fill(dRmin[jetInd]);
1853 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4gen[jetInd].
phi());
1855 double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1857 double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1868 evt.
getByLabel(
"genParticleCandidates", genParticles);
1874 for (
size_t i = 0;
i < genParticles->
size();
i++) {
1915 p4Zp[
nPart] = momentum = d1->
p4();
1927 ZpM = (p4Zp[0] + p4Zp[1]).
mass();
1953 for (
int i = 0;
i < maxJets; ++
i)
1955 for (
int i = 0;
i < 2; ++
i) {
1956 dRmin[jetInd] = 1000.0;
1959 for (CaloJetCollection::const_iterator cal1 = calo1Jets->begin(); cal1 != calo1Jets->end() && jetInd < maxJets;
1961 double delR =
deltaR(cal1->eta(), cal1->phi(), p4Zp[
i].eta(), p4Zp[
i].phi());
1964 if ((delR < dRmin[jetInd]) && (used[nj] == 0)) {
1965 dRmin[jetInd] = delR;
1966 p4cal[jetInd] = cal1->p4();
1990 if (usedInd != -1) {
1996 ZpMM = (p4cal[0] + p4cal[1]).
mass();
1999 if ((ZpMG != 0) && (EtaOk40 == 2)) {
2019 if ((ZpMG != 0) && (EtaOk10 == 2)) {
2039 if ((ZpMG != 0) && (EtaOk13 == 2)) {
2064 for (
int i = 0;
i < maxJets; ++
i)
2066 for (
int i = 0;
i < 2; ++
i) {
2067 dRmin[jetInd] = 1000.0;
2070 for (CaloJetCollection::const_iterator cal2 = calo2Jets->begin(); cal2 != calo2Jets->end() && jetInd < maxJets;
2072 double delR =
deltaR(cal2->eta(), cal2->phi(), p4Zp[
i].eta(), p4Zp[
i].phi());
2074 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
2075 dRmin[jetInd] = delR;
2076 p4cal[jetInd] = cal2->p4();
2095 if (usedInd != -1) {
2101 ZpMM = (p4cal[0] + p4cal[1]).
mass();
2110 for (
int i = 0;
i < maxJets; ++
i)
2112 for (
int i = 0;
i < 2; ++
i) {
2113 dRmin[jetInd] = 1000.0;
2116 for (CaloJetCollection::const_iterator cal3 = calo3Jets->begin(); cal3 != calo3Jets->end() && jetInd < maxJets;
2118 double delR =
deltaR(cal3->eta(), cal3->phi(), p4Zp[
i].eta(), p4Zp[
i].phi());
2120 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
2121 dRmin[jetInd] = delR;
2122 p4cal[jetInd] = cal3->p4();
2141 if (usedInd != -1) {
2147 ZpMM = (p4cal[0] + p4cal[1]).
mass();
2152 cout <<
"Z' (3): nPart = " << nPart << endl;
2164 for (
size_t i = 0;
i < genParticles->
size();
i++) {
2172 cout <<
"TOP: id = " <<
id <<
" mass = " << p.
mass() << endl;
2187 p4par[nJet++] = momentum = d1->
p4();
2190 <<
" mother 1 = " << (d1->
mother())->
pdgId() <<
" Momentum " << momentum <<
" GeV/c" << endl;
2193 p4part[0] = momentum = d1->
p4();
2194 cout <<
">>> part0 = " << p4part[0] << endl;
2197 p4part[3] = momentum = d1->
p4();
2198 cout <<
">>> part3 = " << p4part[3] << endl;
2212 p4par[nJet++] = momentum = d2->
p4();
2215 p4part[
ii] = momentum = d2->
p4();
2216 cout <<
">>> part" << ii <<
" = " << p4part[
ii] << endl;
2220 p4part[
jj] = momentum = d2->
p4();
2221 cout <<
">>> part" << jj <<
" = " << p4part[
jj] << endl;
2226 <<
" mother 2 = " << (d2->
mother())->
pdgId() <<
" Momentum " << momentum <<
" GeV/c" << endl;
2239 double tmass = (p4part[0] + p4part[1] + p4part[2]).
mass();
2240 double tbarmass = (p4part[3] + p4part[4] + p4part[5]).
mass();
2245 cout <<
">>> T Mass = " << tmass <<
" / " << tbarmass << endl;
2247 double mindR = 1000.;
2248 for (
size_t i = 0;
i < 6; ++
i) {
2249 for (
size_t j = 0;
j < 6; ++
j) {
2264 for (
int i = 0;
i < maxJets; ++
i)
2266 for (
int i = 0;
i < 6; ++
i) {
2267 dRmin[jetInd] = 1000.0;
2270 for (CaloJetCollection::const_iterator cal1 = calo1Jets->begin(); cal1 != calo1Jets->end() && jetInd < maxJets;
2272 double delR =
deltaR(cal1->eta(), cal1->phi(), p4par[
i].eta(), p4par[
i].phi());
2274 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
2275 dRmin[jetInd] = delR;
2276 p4cal[jetInd] = cal1->p4();
2296 if (usedInd != -1) {
2299 dRPar1.Fill(dRmin[jetInd]);
2300 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4par[genInd].
phi());
2302 double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2304 double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2311 topMass1.Fill((p4cal[0] + p4cal[1] + p4cal[2]).
mass());
2312 topMass1.Fill((p4cal[3] + p4cal[4] + p4cal[5]).
mass());
2325 for (
int i = 0;
i < maxJets; ++
i)
2327 for (
int i = 0;
i < 6; ++
i) {
2328 dRmin[jetInd] = 1000.0;
2331 for (CaloJetCollection::const_iterator cal2 = calo2Jets->begin(); cal2 != calo2Jets->end() && jetInd < maxJets;
2333 double delR =
deltaR(cal2->eta(), cal2->phi(), p4par[
i].eta(), p4par[
i].phi());
2335 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
2336 dRmin[jetInd] = delR;
2337 p4cal[jetInd] = cal2->p4();
2353 if (usedInd != -1) {
2356 dRPar2.Fill(dRmin[jetInd]);
2357 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4par[genInd].
phi());
2359 double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2361 double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2369 topMass2.Fill((p4cal[0] + p4cal[1] + p4cal[2]).
mass());
2370 topMass2.Fill((p4cal[3] + p4cal[4] + p4cal[5]).
mass());
2383 for (
int i = 0;
i < maxJets; ++
i)
2385 for (
int i = 0;
i < 6; ++
i) {
2386 dRmin[jetInd] = 1000.0;
2389 for (CaloJetCollection::const_iterator cal3 = calo3Jets->begin(); cal3 != calo3Jets->end() && jetInd < maxJets;
2391 double delR =
deltaR(cal3->eta(), cal3->phi(), p4par[
i].eta(), p4par[
i].phi());
2393 if ((delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0)) {
2394 dRmin[jetInd] = delR;
2395 p4cal[jetInd] = cal3->p4();
2410 if (usedInd != -1) {
2413 dRPar3.Fill(dRmin[jetInd]);
2414 double dphi =
deltaPhi(p4cal[jetInd].
phi(), p4par[genInd].
phi());
2416 double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2418 double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2426 topMass3.Fill((p4cal[0] + p4cal[1] + p4cal[2]).
mass());
2427 topMass3.Fill((p4cal[3] + p4cal[4] + p4cal[5]).
mass());
2437 int nTow1, nTow2, nTow3, nTow4;
2448 for (CaloJetCollection::const_iterator ijet = jets->begin(); ijet != jets->end(); ijet++) {
2451 float hadEne = ijet->hadEnergyInHB() + ijet->hadEnergyInHO() + ijet->hadEnergyInHE() + ijet->hadEnergyInHF();
2452 float emEne = ijet->emEnergyInEB() + ijet->emEnergyInEE() + ijet->emEnergyInHF();
2453 float had = ijet->energyFractionHadronic();
2455 float j_et = ijet->et();
2457 if (fabs(ijet->eta()) < 1.3) {
2466 if (ijet->pt() > minJetPt10)
2469 if ((fabs(ijet->eta()) > 1.3) && (fabs(ijet->eta()) < 3.)) {
2478 if (ijet->pt() > minJetPt10)
2481 if (fabs(ijet->eta()) > 3.) {
2490 if (ijet->pt() > minJetPt10)
2504 const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2508 nTow1 = nTow2 = nTow3 = nTow4 = 0;
2510 float et = jetCaloRefs[
i]->et();
2530 if ((jjet == 1) && (fabs(ijet->eta()) < 1.3)) {
2534 float t_et = jetCaloRefs[
i]->et();
2535 double delR =
deltaR(ijet->eta(), ijet->phi(), jetCaloRefs[
i]->eta(), jetCaloRefs[
i]->phi());
2544 if ((jjet == 2) && (fabs(ijet->eta()) < 1.3)) {
2555 double SumEtNotJets(0);
2556 double SumEtJets(0);
2557 double SumEtTowers(0);
2562 double sumTowerAllPx(0);
2563 double sumTowerAllPy(0);
2565 double sumTowerAllEx(0);
2566 double sumTowerAllEy(0);
2568 std::vector<CaloTowerPtr> UsedTowerList;
2569 std::vector<CaloTower> TowerUsedInJets;
2570 std::vector<CaloTower> TowerNotUsedInJets;
2578 nTow1 = nTow2 = nTow3 = nTow4 = 0;
2580 double sum_et = 0.0;
2581 double sum_ex = 0.0;
2582 double sum_ey = 0.0;
2587 Double_t et =
tower->et();
2600 double phix =
tower->phi();
2605 double et =
tower->et();
2609 sum_ex += et *
cos(phix);
2610 sum_ey += et *
sin(phix);
2614 SumEtTowers +=
tower->et();
2616 sumTowerAllEx += et *
cos(phi);
2617 sumTowerAllEy += et *
sin(phi);
2625 MET1.Fill(
sqrt(sum_ex * sum_ex + sum_ey * sum_ey));
2626 MET12.Fill(
sqrt(sum_ex * sum_ex + sum_ey * sum_ey));
2627 MET13.Fill(
sqrt(sum_ex * sum_ex + sum_ey * sum_ey));
2648 UsedTowerList.clear();
2649 TowerUsedInJets.clear();
2650 TowerNotUsedInJets.clear();
2654 for (CaloJetCollection::const_iterator ijet = jets->begin(); ijet != jets->end(); ijet++) {
2655 Double_t
jetPt = ijet->pt();
2656 Double_t
jetPhi = ijet->phi();
2660 Double_t jetPx = jetPt *
cos(jetPhi);
2661 Double_t jetPy = jetPt *
sin(jetPhi);
2666 const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2667 int nConstituents = jetCaloRefs.size();
2669 UsedTowerList.push_back(jetCaloRefs[
i]);
2676 int NTowersUsed = UsedTowerList.size();
2681 Double_t et =
tower->et();
2684 Double_t phi =
tower->phi();
2685 SumEtTowers +=
tower->et();
2687 sumTowerAllPx += et *
cos(phi);
2688 sumTowerAllPy += et *
sin(phi);
2692 for (
int i = 0;
i < NTowersUsed;
i++) {
2693 if (
tower->id() == UsedTowerList[
i]->id()) {
2700 TowerUsedInJets.push_back(t);
2702 TowerNotUsedInJets.push_back(t);
2707 int nUsed = TowerUsedInJets.size();
2708 int nNotUsed = TowerNotUsedInJets.size();
2713 for (
int i = 0;
i < nUsed;
i++) {
2714 SumEtJets += TowerUsedInJets[
i].et();
2718 for (
int i = 0;
i < nNotUsed;
i++) {
2719 if (TowerNotUsedInJets[
i].
et() > 0.5)
2720 SumEtNotJets += TowerNotUsedInJets[
i].et();
2729 UsedTowerList.clear();
2730 TowerUsedInJets.clear();
2731 TowerNotUsedInJets.clear();
2735 for (CaloJetCollection::const_iterator ijet = jets->begin(); ijet != jets->end(); ijet++) {
2736 Double_t jetPt = ijet->pt();
2737 Double_t jetPhi = ijet->phi();
2741 Double_t jetPx = jetPt *
cos(jetPhi);
2742 Double_t jetPy = jetPt *
sin(jetPhi);
2747 const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2748 int nConstituents = jetCaloRefs.size();
2750 UsedTowerList.push_back(jetCaloRefs[
i]);
2760 NTowersUsed = UsedTowerList.size();
2765 Double_t et =
tower->et();
2768 Double_t phi =
tower->phi();
2770 SumEtTowers +=
tower->et();
2772 sumTowerAllPx += et *
cos(phi);
2773 sumTowerAllPy += et *
sin(phi);
2777 for (
int i = 0;
i < NTowersUsed;
i++) {
2778 if (
tower->id() == UsedTowerList[
i]->id()) {
2785 TowerUsedInJets.push_back(t);
2787 TowerNotUsedInJets.push_back(t);
2792 nUsed = TowerUsedInJets.size();
2793 nNotUsed = TowerNotUsedInJets.size();
2798 for (
int i = 0;
i < nUsed;
i++) {
2799 SumEtJets += TowerUsedInJets[
i].et();
2803 for (
int i = 0;
i < nNotUsed;
i++) {
2804 if (TowerNotUsedInJets[
i].
et() > 0.5)
2805 SumEtNotJets += TowerNotUsedInJets[
i].et();
2814 UsedTowerList.clear();
2815 TowerUsedInJets.clear();
2816 TowerNotUsedInJets.clear();
2820 for (CaloJetCollection::const_iterator ijet = jets->begin(); ijet != jets->end(); ijet++) {
2821 Double_t jetPt = ijet->pt();
2822 Double_t jetPhi = ijet->phi();
2826 Double_t jetPx = jetPt *
cos(jetPhi);
2827 Double_t jetPy = jetPt *
sin(jetPhi);
2832 const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2833 int nConstituents = jetCaloRefs.size();
2835 UsedTowerList.push_back(jetCaloRefs[
i]);
2845 NTowersUsed = UsedTowerList.size();
2850 Double_t et =
tower->et();
2861 for (
int i = 0;
i < NTowersUsed;
i++) {
2862 if (
tower->id() == UsedTowerList[
i]->id()) {
2869 TowerUsedInJets.push_back(t);
2871 TowerNotUsedInJets.push_back(t);
2876 nUsed = TowerUsedInJets.size();
2877 nNotUsed = TowerNotUsedInJets.size();
2882 for (
int i = 0;
i < nUsed;
i++) {
2883 SumEtJets += TowerUsedInJets[
i].et();
2887 for (
int i = 0;
i < nNotUsed;
i++) {
2888 if (TowerNotUsedInJets[
i].
et() > 0.5)
2889 SumEtNotJets += TowerNotUsedInJets[
i].et();
std::string GenJetAlgorithm3
TH1F dijetMassCor_700_101
virtual double vx() const =0
x coordinate of vertex position
std::string CaloJetAlgorithm3
std::string GenJetAlgorithm1
std::string CaloJetAlgorithm4
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
Sin< T >::type sin(const T &t)
std::vector< CaloTower >::const_iterator const_iterator
TH1F dijetMassCor_5000_101
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
TH1F dijetMassCor_2000_131
virtual double vy() const =0
y coordinate of vertex position
TH1F dijetMassCor_5000_131
TH1F dijetMassCor_2000_101
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
TH1F dijetMassCor_2000_401
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
virtual int pdgId() const =0
PDG identifier.
Cos< T >::type cos(const T &t)
Abs< T >::type abs(const T &t)
TH1F dijetMassCor_700_131
def gen(fragment, howMuch)
Production test section ####.
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
const_iterator end() const
TH1F dijetMassCor_5000_401
std::string GenJetAlgorithm4
TString nPart(Int_t part, TString string, TString delimit=";", Bool_t removerest=true)
TH1F dijetMassCor_700_401
XYZVectorD XYZVector
spatial vector with cartesian internal representation
virtual double mass() const =0
mass
std::string CaloJetAlgorithm1
virtual double vz() const =0
z coordinate of vertex position
std::string GenJetAlgorithm2
std::string CaloJetAlgorithm2
virtual size_type numberOfDaughters() const =0
number of daughters
const_iterator begin() const
Public Types inherited from edm::EDAnalyzer
1.8.11
