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.74 by arizzi, Thu May 17 10:58:34 2012 UTC vs.
Revision 1.75 by arizzi, Thu May 24 15:55:16 2012 UTC

#	Line 221 | Line 221 | HbbAnalyzerNew::produce(edm::Event& iEve
221		std::vector< PileupSummaryInfo> pu = (*puHandle);
222		for (std::vector<PileupSummaryInfo>::const_iterator it= pu.begin(); it!=pu.end(); ++it){
223		int bx = (*it).getBunchCrossing();
224	+	if(bx == 0) { auxInfo->puInfo.truePU = (*it).getTrueNumInteractions();}
225		unsigned int num = (*it).getPU_NumInteractions();
226		//      std::cout <<" PU PUSHING "<<bx<<" " <<num<<std::endl;
227		auxInfo->puInfo.pus[bx]  =num;
#	Line 1361 | Line 1362 | BTagSFContainer btagSFs;
1362		}
1363
1364		edm::Handle<edm::View<pat::MET> > metPFHandle;
1365	<	iEvent.getByLabel("patMETs",metPFHandle);
1365	>	iEvent.getByLabel("patMETsPFlow",metPFHandle);
1366		edm::View<pat::MET> metsPF = *metPFHandle;
1367
1368		if(metsPF.size()){
#	Line 1503 | Line 1504 | BTagSFContainer btagSFs;
1504
1505		/* Check if there are photons sharing the super cluster*/
1506		for(size_t k=0;k<photonsForIso.size();k++) {
1506	–	// std::cout << "PHO: eta " <<   photonsForIso[k].eta()  << " phi " <<  photonsForIso[k].phi() <<  " pt " <<  photonsForIso[k].pt() << " " <<  ef.pfPhoIso<< " refs " <<
1507	–	//photonsForIso[k].superClusterRef().key() << " vs " << elec->superCluster().key() <<     std::endl;
1508	–	//std::cout << "dR , asym " <<  deltaR(elec->eta(),elec->phi(),photonsForIso[k].eta(),photonsForIso[k].phi())  << " "  << abs(photonsForIso[k].pt()-elec->pt())/(photonsForIso[k].pt()+elec->pt()) << std::endl;
1507		if(deltaR(elec->eta(),elec->phi(),photonsForIso[k].eta(),photonsForIso[k].phi()) < 0.05 && abs(photonsForIso[k].pt()-elec->pt())/(photonsForIso[k].pt()+elec->pt()) < 0.05 ) {
1508		std::cout << "Double counting of supercluster!" << std::endl;
1509		ef.pfPhoIsoDoubleCounted+=photonsForIso[k].pt();
#	Line 1531 | Line 1529 | BTagSFContainer btagSFs;
1529		ef.HoE = elec->hadronicOverEm();
1530		ef.convDist = elec->convDist();
1531		ef.convDcot = elec->convDcot();
1532	<	if(elec->gsfTrack().isNonnull()) ef.innerHits = elec->gsfTrack()->trackerExpectedHitsInner().numberOfHits();
1532	>	if(elec->gsfTrack().isNonnull())
1533	>	{
1534	>	ef.innerHits = elec->gsfTrack()->trackerExpectedHitsInner().numberOfHits();
1535	>	}
1536		ef.isEB = elec->isEB();
1537		ef.isEE = elec->isEE();
1538		/* 2012 ELEID*/

Diff Legend

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