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.13 by tboccali, Mon Jul 25 08:55:40 2011 UTC vs.
Revision 1.17 by arizzi, Wed Aug 10 13:35:43 2011 UTC

#	Line 294 | Line 294 | HbbAnalyzerNew::produce(edm::Event& iEve
294		Handle<reco::VertexCollection> recVtxs;
295		iEvent.getByLabel("offlinePrimaryVertices", recVtxs);
296
297	+	auxInfo->pvInfo.nVertices = recVtxs->size();
298	+
299		for(size_t i = 0; i < recVtxs->size(); ++ i) {
300		const Vertex &vtx = (*recVtxs)[i];
301		double RecVtxProb=TMath::Prob(vtx.chi2(),vtx.ndof());
#	Line 311 | Line 313 | HbbAnalyzerNew::produce(edm::Event& iEve
313
314
315		edm::Handle<double> rhoHandle;
316	<	iEvent.getByLabel(edm::InputTag("kt6PFJets", "rho"),rhoHandle); // configure srcRho = cms.InputTag('kt6PFJets")
315	<	auxInfo->puInfo.rho = *rhoHandle;
316	>	iEvent.getByLabel(edm::InputTag("kt6PFJets", "rho"),rhoHandle);   auxInfo->puInfo.rho = *rhoHandle;
317
318		//// real start
319
#	Line 738 | Line 739 | HbbAnalyzerNew::produce(edm::Event& iEve
739		mf.tIso=mu->trackIso();
740		mf.eIso=mu->ecalIso();
741		mf.hIso=mu->hcalIso();
742	+	mf.pfChaIso=mu->chargedHadronIso();
743	+	mf.pfPhoIso=mu->photonIso();
744	+	mf.pfNeuIso=mu->neutralHadronIso();
745		Geom::Phi<double> deltaphi(mu->phi()-atan2(mf.p4.Px(), mf.p4.Py()));
746		double acop = deltaphi.value();
747		mf.acop=acop;
#	Line 801 | Line 805 | HbbAnalyzerNew::produce(edm::Event& iEve
805		ef.tIso=elec->trackIso();
806		ef.eIso=elec->ecalIso();
807		ef.hIso=elec->hcalIso();
808	+	ef.pfChaIso=elec->chargedHadronIso();
809	+	ef.pfPhoIso=elec->photonIso();
810	+	ef.pfNeuIso=elec->neutralHadronIso();
811	+
812		Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
813		ef.acop = deltaphi.value();
814		//
815		// fill eleids
816		//
817	<	ef.id95 = elec->electronID("simpleEleId95cIso");
817	>	/*    ef.id95 = elec->electronID("simpleEleId95cIso");
818		ef.id85 = elec->electronID("simpleEleId85cIso");
819		ef.id70 = elec->electronID("simpleEleId70cIso");
820		ef.id95r = elec->electronID("simpleEleId95relIso");
821		ef.id70r = elec->electronID("simpleEleId70relIso");
822		ef.id85r = elec->electronID("simpleEleId85relIso");
823	+	*/
824	+	ef.id95 =elec->electronID("eidVBTFCom95");
825	+	ef.id95r=elec->electronID("eidVBTFRel95");
826	+	ef.id85 =elec->electronID("eidVBTFCom85");
827	+	ef.id85r=elec->electronID("eidVBTFRel85");
828	+	ef.id80 =elec->electronID("eidVBTFCom80");
829	+	ef.id80r=elec->electronID("eidVBTFRel80");
830	+	ef.id70 =elec->electronID("eidVBTFCom70");
831	+	ef.id70r=elec->electronID("eidVBTFRel70");
832
833		if(runOnMC_){
834		const GenParticle* elecMc = elec->genLepton();
#	Line 956 | Line 973 | HbbAnalyzerNew::produce(edm::Event& iEve
973		df.daughter2.hIso = elecDau1.hcalIso();
974
975		// ids
976	<	df.daughter1.id95 = elecDau0.electronID("simpleEleId95cIso");
976	>	/*df.daughter1.id95 = elecDau0.electronID("simpleEleId95cIso");
977		df.daughter1.id85 = elecDau0.electronID  ("simpleEleId85cIso");
978		df.daughter1.id70 = elecDau0.electronID  ("simpleEleId70cIso");
979		df.daughter1.id95r = elecDau0.electronID ("simpleEleId95relIso");
#	Line 970 | Line 987 | HbbAnalyzerNew::produce(edm::Event& iEve
987		df.daughter2.id95r = elecDau1.electronID ("simpleEleId95relIso");
988		df.daughter2.id85r = elecDau1.electronID ("simpleEleId85relIso");
989		df.daughter2.id70r = elecDau1.electronID ("simpleEleId70relIso");
990	<
990	>	*/
991		hbbInfo->diElectronInfo.push_back(df);
992
993		}
#	Line 989 | Line 1006 | HbbAnalyzerNew::produce(edm::Event& iEve
1006
1007
1008		iEvent.put(hbbInfo);
1009	+	iEvent.put(auxInfo);
1010
1011
1012		}

Diff Legend

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