CommonTools/AnalysisJet.cc at master · pandolf/CommonTools · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
#include "AnalysisJet.h"

#include <map>
#include <iostream>

bool AnalysisJet::btag_loose() const {

  return trackCountingHighEffBJetTag>1.7;

}


bool AnalysisJet::btag_medium() const {

  return trackCountingHighEffBJetTag>3.3;

}
Bool_t AnalysisJet::passedJetID( const std::string& strength ) const {

 Bool_t returnBool = false;

 std::map<string,long> numbers;
   	numbers["minimal"]=1;

switch ( numbers[strength] ) {
	case 1: //minimal
	 	returnBool = ( this->nConstituents() > 1) &&
			( (ePhotons + eHFEM)/Energy() < 0.99 ) &&
			( (eNeutralHadrons)/Energy() < 0.99 );

		if( fabs(Eta())<2.4 ) {returnBool = returnBool && ( eChargedHadrons > 0. );}
		break;
	default:
	std::cout << "Jet ID '" << strength << "' not implemented yet. Returning FALSE." << std::endl;
	}

 return returnBool;

}