ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc

Comparing UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc (file contents):
Revision 1.37 by degrutto, Thu Sep 15 14:00:01 2011 UTC vs.
Revision 1.44 by tboccali, Mon Sep 19 14:08:08 2011 UTC

#	Line 420 | Line 420 | BTagSFContainer btagSFs;
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
#	Line 435 | Line 436 | BTagSFContainer btagSFs;
436		//genJet
437		const reco::GenJet *gJ = jet_iter->genJet();
438		//physical parton for mother info ONLY
439	<	if( (jet_iter->genParton())
440	<	and (jet_iter->genParton()->mother()) )
441	<	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
439	>	if( (jet_iter->genParton()) ){
440	>	sj.bestMCid = jet_iter->genParton()->pdgId();
441	>	if( (jet_iter->genParton()->mother()) )
442	>	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
443	>	}
444		TLorentzVector gJp4;
445		if(gJ){
446		gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
447	<	sj. bestMCp4mom = gJp4;
447	>	sj.bestMCp4 = gJp4;
448		if(verbose_){
449		std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
450		std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
#	Line 460 | Line 463 | BTagSFContainer btagSFs;
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);
#	Line 475 | Line 479 | BTagSFContainer btagSFs;
479		//genJet
480		const reco::GenJet *gJ = jet_iter->genJet();
481		//physical parton for mother info ONLY
482	<	if( (jet_iter->genParton())
483	<	and (jet_iter->genParton()->mother()) )
484	<	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
482	>	if( (jet_iter->genParton()) ){
483	>	sj.bestMCid = jet_iter->genParton()->pdgId();
484	>	if( (jet_iter->genParton()->mother()) )
485	>	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
486	>	}
487		TLorentzVector gJp4;
488		if(gJ){
489		gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
490	<	sj. bestMCp4mom = gJp4;
490	>	sj.bestMCp4 = gJp4;
491		if(verbose_){
492		std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
493		std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
#	Line 500 | Line 506 | BTagSFContainer btagSFs;
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);
#	Line 515 | Line 522 | BTagSFContainer btagSFs;
522		//genJet
523		const reco::GenJet *gJ = jet_iter->genJet();
524		//physical parton for mother info ONLY
525	<	if( (jet_iter->genParton())
526	<	and (jet_iter->genParton()->mother()) )
527	<	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
525	>	if( (jet_iter->genParton()) ){
526	>	sj.bestMCid = jet_iter->genParton()->pdgId();
527	>	if( (jet_iter->genParton()->mother()) )
528	>	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
529	>	}
530		TLorentzVector gJp4;
531		if(gJ){
532		gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
533	<	sj. bestMCp4mom = gJp4;
533	>	sj.bestMCp4 = gJp4;
534		if(verbose_){
535		std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
536		std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
#	Line 540 | Line 549 | BTagSFContainer btagSFs;
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);
#	Line 595 | Line 605 | BTagSFContainer btagSFs;
605		//genJet
606		const reco::GenJet *gJ = jet_iter->genJet();
607		//physical parton for mother info ONLY
608	<	if( (jet_iter->genParton())
609	<	and (jet_iter->genParton()->mother()) )
610	<	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
608	>	if( (jet_iter->genParton()) ){
609	>	sj.bestMCid = jet_iter->genParton()->pdgId();
610	>	if( (jet_iter->genParton()->mother()) )
611	>	sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
612	>	}
613		TLorentzVector gJp4;
614		if(gJ){
615		gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
616	<	sj. bestMCp4mom = gJp4;
616	>	sj.bestMCp4 = gJp4;
617		if(verbose_){
618		std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
619		std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
#	Line 681 | Line 693 | BTagSFContainer btagSFs;
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);
#	Line 732 | Line 745 | BTagSFContainer btagSFs;
745		if (verbose_)     std::cout <<" METTC "<<     hbbInfo->tcmet.metSig <<" " <<     hbbInfo->tcmet.sumEt<<std::endl;
746		}
747
748	+	edm::Handle<edm::View<reco::MET> > pfMETNoPUHandle;
749	+	iEvent.getByLabel("pfMETNoPU",pfMETNoPUHandle);
750	+	edm::View<reco::MET> metspfMETNoPU = *pfMETNoPUHandle;
751	+	if(metspfMETNoPU.size()){
752	+	hbbInfo->metNoPU.sumEt=(metspfMETNoPU[0]).sumEt();
753	+	hbbInfo->metNoPU.metSig=(metspfMETNoPU[0]).significance();
754	+	hbbInfo->metNoPU.eLong=(metspfMETNoPU[0]).e_longitudinal();
755	+	hbbInfo->metNoPU.p4=GENPTOLOR((metspfMETNoPU[0]));
756	+	if (verbose_)     std::cout <<" pfMETNoPU "<<     hbbInfo->metNoPU.metSig <<" " <<     hbbInfo->metNoPU.sumEt<<std::endl;
757	+	}
758
759		edm::Handle<edm::View<reco::MET> > mHTHandle;
760		iEvent.getByLabel("patMETsHT",mHTHandle);
#	Line 744 | Line 767 | BTagSFContainer btagSFs;
767		if (verbose_)     std::cout <<" METHT "<<     hbbInfo->mht.metSig <<" " <<     hbbInfo->mht.sumEt<<std::endl;
768		}
769
770	<	edm::Handle<edm::View<pat::MET> > metHandle;
770	>	edm::Handle<edm::View<reco::MET> > metHandle;
771		iEvent.getByLabel(metLabel_,metHandle);
772	<	edm::View<pat::MET> mets = *metHandle;
772	>	edm::View<reco::MET> mets = *metHandle;
773
774		if(mets.size()){
775		hbbInfo->calomet.sumEt=(mets[0]).sumEt();
#	Line 813 | Line 836 | BTagSFContainer btagSFs;
836		const reco::HitPattern& q = gTrack->hitPattern();
837		mf.globChi2=gTrack.get()->normalizedChi2();
838		mf.globNHits=q.numberOfValidMuonHits();
839	<	mf. validMuStations = q. muonStationsWithValidHits();
839	>	mf.validMuStations = q.muonStationsWithValidHits();
840		}else{
841		mf.globChi2=-99;
842		mf.globNHits=-99;
#	Line 876 | Line 899 | BTagSFContainer btagSFs;
899		ef.pfPhoIso=elec->photonIso();
900		ef.pfNeuIso=elec->neutralHadronIso();
901
902	<	Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
902	>	Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
903		ef.acop = deltaphi.value();
904		//
905		// fill eleids
#	Line 935 | Line 958 | BTagSFContainer btagSFs;
958		tf.tIso=tau->trackIso();
959		tf.eIso=tau->ecalIso();
960		tf.hIso=tau->hcalIso();
961	<	Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
961	>	Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
962		double acop = deltaphi.value();
963		tf.acop=acop;
964		tf.idbyIso=tau->tauID("byIsolation");
#	Line 1185 | Line 1208 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1208		sj.jpb=jet_iter->bDiscriminator("jetBProbabilityBJetTags");
1209		sj.ssvhe=jet_iter->bDiscriminator("simpleSecondaryVertexHighEffBJetTags");
1210		sj.csv=jet_iter->bDiscriminator("combinedSecondaryVertexBJetTags");
1211	<	sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1211	>	sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1212		sj.charge=jet_iter->jetCharge();
1213		sj.ntracks=jet_iter->associatedTracks().size();
1214		sj.p4=GENPTOLORP(jet_iter);
1215	+	//    std::cout << " ECCO "<<sj.csv<< " "<< sj.p4.Pt()<<std::endl;
1216		sj.chargedTracksFourMomentum=(getChargedTracksMomentum(&*(jet_iter)));
1217		sj.SF_CSVL=1;
1218		sj.SF_CSVM=1;
#	Line 1204 | Line 1228 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1228		sj.neutralHadronEFraction = jet_iter-> neutralHadronEnergyFraction ();
1229		sj.chargedEmEFraction = jet_iter-> chargedEmEnergyFraction ();
1230		sj.neutralEmEFraction = jet_iter-> neutralEmEnergyFraction ();
1231	<	sj. nConstituents = jet_iter->getPFConstituents().size();
1231	>	sj.nConstituents = jet_iter->getPFConstituents().size();
1232
1233		}
1234		//

Diff Legend

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