ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/HbbAnalysis/src/JetFlavour.cc

Comparing UserCode/HbbAnalysis/src/JetFlavour.cc (file contents):
Revision 1.2 by amagnan, Tue Jun 2 13:24:26 2009 UTC vs.
Revision 1.5 by amagnan, Thu Sep 17 19:39:52 2009 UTC

#	Line 158 | Line 158 | namespace HbbAnalysis {
158		unsigned int *eid  = new unsigned int[nPartons];
159		unsigned int *muid = new unsigned int[nPartons];
160
161	–	assert (nPartons <= 2);
162	–
161		for (unsigned int ip(0); ip<nPartons; ip++){//loop on partons
162		nDaugh[ip] = 0;
163		flav[ip] = 1;
#	Line 266 | Line 264 | namespace HbbAnalysis {
264		}//for b and c jets
265
266
267	<	if (parents[ip].size() > 250) std::cout << "WARNING ! Initial default size of parents not sufficient = " << parents[ip].size() << std::endl;
268	<	if (daughters[ip].size() > 250) std::cout << "WARNING ! Initial default size of daughters not sufficient = " << daughters[ip].size() << std::endl;
267	>	if (parents[ip].size() > 500) std::cout << "WARNING ! Initial default size of parents not sufficient = " << parents[ip].size() << std::endl;
268	>	if (daughters[ip].size() > 500) std::cout << "WARNING ! Initial default size of daughters not sufficient = " << daughters[ip].size() << std::endl;
269
270
271		}//loop on partons
#	Line 319 | Line 317 | namespace HbbAnalysis {
317		}//fillPartons
318
319
320	<	std::pair<int,unsigned int> JetFlavour::partonMatchingGenJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, double aDR){//parton match with gen jet
320	>	const std::pair<int,unsigned int> JetFlavour::partonMatchingGenJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, const double aDR) const{//parton match with gen jet
321
322		int aIndex = -1;
323		double deltaRMin = 1000;
#	Line 339 | Line 337 | namespace HbbAnalysis {
337
338		}//parton match with gen jet
339
340	<	unsigned int JetFlavour::leptonMatchingGenJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, double aDR){//lepton match with gen jet
340	>	const unsigned int JetFlavour::leptonMatchingGenJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, const double aDR) const{//lepton match with gen jet
341
342		int aIndex = -1;
343		double deltaRMin = 1000;
#	Line 375 | Line 373 | namespace HbbAnalysis {
373
374		}//lepton match with gen jet
375
376	<	unsigned int JetFlavour::partonMatchingRecoJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, double aDR){//parton match with reco jet
376	>	const unsigned int JetFlavour::partonMatchingRecoJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, const double aDR) const{//parton match with reco jet
377
378		int aIndex = -1;
379		double deltaRMin = 1000;
#	Line 393 | Line 391 | namespace HbbAnalysis {
391
392		}//parton match with reco jet
393
394	<	unsigned int JetFlavour::leptonMatchingRecoJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, double aDR){//lepton match with reco jet
394	>	const unsigned int JetFlavour::leptonMatchingRecoJet(const pat::Jet & aJet, const edm::Handle<reco::GenParticleCollection> & aGenParticles, const double aDR) const{//lepton match with reco jet
395
396		int aIndex = -1;
397		double deltaRMin = 1000;
#	Line 457 | Line 455 | namespace HbbAnalysis {
455
456		}
457
458	<	void JetFlavour::printSummary() {
458	>	void JetFlavour::printSummary() const{
459
460		std::cout <<
461		"=====================================================================" << std::endl <<

Diff Legend

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