CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetIdSelector.cc
Go to the documentation of this file.
1 
33 
37 
50 
51 #include <memory>
52 #include <vector>
53 #include <sstream>
54 
56 // class definition
58 template<typename T>
60 {
61 public:
62  typedef std::vector<T> JetCollection;
63  // construction/destruction
64  explicit JetIdSelector(const edm::ParameterSet& iConfig);
65  virtual ~JetIdSelector();
66 
67  // member functions
68  void produce(edm::Event& iEvent,const edm::EventSetup& iSetup) override;
69  void endJob() override;
70 
71 private:
72  // member data
77 
78 
79  unsigned int nJetsTot_;
80  unsigned int nJetsPassed_;
85 };
86 
87 
88 using namespace std;
89 
90 
92 // construction/destruction
94 //______________________________________________________________________________
95 template<typename T>
97  : src_ (iConfig.getParameter<edm::InputTag> ("src"))
98  , qualityStr (iConfig.getParameter<string> ("idLevel"))
99  , jetIDMap_(iConfig.getUntrackedParameter<edm::InputTag> ("jetIDMap", edm::InputTag("ak5JetID")))
100  , moduleLabel_(iConfig.getParameter<string> ("@module_label"))
101  , nJetsTot_(0)
102  , nJetsPassed_(0)
103 {
104  produces<JetCollection>();
105 
106  use_pfloose = false;
107  use_pfmedium = false;
108  use_pftight = false;
109 
110  if ( qualityStr == "MINIMAL" ) {
111  jetIDFunctor
113  use_pfloose = true;
114  }
115  else if ( qualityStr == "LOOSE_AOD" ) {
116  jetIDFunctor
118  use_pfloose = true;
119  }
120  else if ( qualityStr == "LOOSE" ) {
121  jetIDFunctor
123  use_pfloose = true;
124  }
125  else if ( qualityStr == "MEDIUM" ) {
126  jetIDFunctor
128  // There is no medium quality for CaloJet !!
129  use_pfmedium = true;
130  }
131  else if ( qualityStr == "TIGHT" ) {
132  jetIDFunctor
134  use_pftight = true;
135  }
136  else
137  throw cms::Exception("InvalidInput")
138  << "Expect quality to be one of MINIMAL, LOOSE_AOD, LOOSE, MEDIUM, TIGHT" << std::endl;
139 }
140 
141 
142 //______________________________________________________________________________
143 template<typename T>
145  if(jetIDFunctor) delete jetIDFunctor;
146 }
147 
148 
149 
151 // implementation of member functions
153 
154 //______________________________________________________________________________
155 template<typename T>
157 {
158  auto_ptr<JetCollection> selectedJets(new JetCollection);
159  edm::Handle<reco::JetView> jets; // uncorrected jets!
160  iEvent.getByLabel(src_,jets);
161 
162 
163  // handle to the jet ID variables
165 
166  if(typeid((*jets)[0]) == typeid(reco::CaloJet))
167  iEvent.getByLabel( jetIDMap_, hJetIDMap );
168 
169  unsigned int idx=0;
170  bool passed = false;
171 
172  for ( edm::View<reco::Jet>::const_iterator ibegin = jets->begin(),
173  iend = jets->end(), iJet = ibegin;
174  iJet != iend; ++iJet ) {
175 
176  // initialize the boolean flag to false
177  passed = false;
178 
179  //calculate the Calo jetID
180  const std::type_info & type = typeid((*jets)[idx]);
181  if( type == typeid(reco::CaloJet) ) {
182  const reco::CaloJet calojet = static_cast<const reco::CaloJet &>((*jets)[idx]);
183  edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);
184  reco::JetID const & jetId = (*hJetIDMap)[ jetRef ];
185  passed = (*jetIDFunctor)( calojet, jetId);
186  }
187 
188  //calculate the PF jetID
189  if ( type == typeid(reco::PFJet) ) {
190  const reco::PFJet pfjet = static_cast<const reco::PFJet &>((*jets)[idx]);
191  bool passingLoose=false;
192  bool passingMedium=false;
193  bool passingTight=false;
194  bool ThisIsClean=true;
195  //apply following only if |eta|<2.4: CHF>0, CEMF<0.99, chargedMultiplicity>0
196  if(( pfjet.chargedHadronEnergy()/ pfjet.energy())<= 0.0
197  && fabs(pfjet.eta())<2.4) ThisIsClean=false;
198  if( (pfjet.chargedEmEnergy()/pfjet.energy())>= 0.99
199  && fabs(pfjet.eta())<2.4 ) ThisIsClean=false;
200  if( pfjet.chargedMultiplicity()<=0 && fabs(pfjet.eta())<2.4 )
201  ThisIsClean=false;
202 
203  // always require #Constituents > 1
204  if( pfjet.nConstituents() <=1 ) ThisIsClean=false;
205 
206  if(ThisIsClean &&
207  (pfjet.neutralHadronEnergy()/pfjet.energy())< 0.99
208  && (pfjet.neutralEmEnergy()/pfjet.energy())<0.99)
209  passingLoose=true;
210 
211  if(ThisIsClean &&
212  (pfjet.neutralHadronEnergy()/pfjet.energy())< 0.95
213  && (pfjet.neutralEmEnergy()/pfjet.energy())<0.95)
214  passingMedium=true;
215 
216  if(ThisIsClean &&
217  (pfjet.neutralHadronEnergy()/pfjet.energy())< 0.90
218  && (pfjet.neutralEmEnergy()/pfjet.energy())<0.90)
219  passingTight=true;
220 
221 
222  if ( use_pfloose && passingLoose) passed = true;
223  if ( use_pfmedium && passingMedium) passed = true;
224  if ( use_pftight && passingTight) passed = true;
225  }
226 
227  if ( type == typeid(reco::GenJet) || type == typeid(reco::JPTJet)) {
228  edm::LogWarning( "JetId" )<<
229  "Criteria for jets other than CaloJets and PFJets are not yet implemented";
230  passed = true;
231  } // close GenJet, JPT jet
232 
233 
234  const T& goodJet = static_cast<const T&>((*jets)[idx]);
235  if(passed) selectedJets->push_back( goodJet );
236 
237  idx++;
238  } // close jet iterator
239 
240 
241 
242  nJetsTot_ +=jets->size();
243  nJetsPassed_+=selectedJets->size();
244  iEvent.put(selectedJets);
245 }
246 
247 
248 
249 //______________________________________________________________________________
250 template<typename T>
252 {
253  stringstream ss;
254  ss<<"nJetsTot="<<nJetsTot_<<" nJetsPassed="<<nJetsPassed_
255  <<" fJetsPassed="<<100.*(nJetsPassed_/(double)nJetsTot_)<<"%\n";
256  cout<<"++++++++++++++++++++++++++++++++++++++++++++++++++"
257  <<"\n"<< moduleLabel_ << "(JetIdSelector) SUMMARY:\n"<<ss.str()
258  <<"++++++++++++++++++++++++++++++++++++++++++++++++++"
259  <<endl;
260 }
261 
262 
264 // plugin definition
270 
275 
276 
type
Definition: HCALResponse.h:21
virtual double energy() const GCC11_FINAL
energy
virtual ~JetIdSelector()
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:81
Jets made from CaloTowers.
Definition: CaloJet.h:29
std::string qualityStr
float chargedEmEnergy() const
chargedEmEnergy
Definition: PFJet.h:138
JetIDSelectionFunctor * jetIDFunctor
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
edm::InputTag src_
JetIdSelector< reco::CaloJet > CaloJetIdSelector
void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) override
JetIdSelector(const edm::ParameterSet &iConfig)
Jet ID object.
Definition: JetID.h:16
int chargedMultiplicity() const
chargedMultiplicity
Definition: PFJet.h:151
unsigned int nJetsPassed_
Jets made from PFObjects.
Definition: PFJet.h:21
float neutralEmEnergy() const
neutralEmEnergy
Definition: PFJet.h:146
std::string moduleLabel_
unsigned int nJetsTot_
int iEvent
Definition: GenABIO.cc:230
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:116
JetIdSelector< reco::JPTJet > JPTJetIdSelector
Jets made from CaloJets corrected for ZSP and tracks.
Definition: JPTJet.h:29
vector< PseudoJet > jets
Jets made from MC generator particles.
Definition: GenJet.h:24
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:390
virtual float eta() const GCC11_FINAL
momentum pseudorapidity
Jet selector for pat::Jets and for CaloJets.
JetIdSelector< reco::GenJet > GenJetIdSelector
tuple idx
DEBUGGING if hasattr(process,&quot;trackMonIterativeTracking2012&quot;): print &quot;trackMonIterativeTracking2012 D...
virtual int nConstituents() const
of constituents
Definition: Jet.h:65
float neutralHadronEnergy() const
neutralHadronEnergy
Definition: PFJet.h:98
tuple cout
Definition: gather_cfg.py:121
std::vector< T > JetCollection
JetIdSelector< reco::PFJet > PFJetIdSelector
void endJob() override
moduleLabel_(iConfig.getParameter< string >("@module_label"))
long double T
edm::InputTag jetIDMap_
float chargedHadronEnergy() const
chargedHadronEnergy
Definition: PFJet.h:94