ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/Morgan/src/JetAnalyzer.cc
(Generate patch)

Comparing UserCode/Morgan/src/JetAnalyzer.cc (file contents):
Revision 1.6 by lethuill, Thu Jan 8 15:43:30 2009 UTC vs.
Revision 1.9 by lethuill, Fri Apr 10 08:27:40 2009 UTC

# Line 156 | Line 156 | void JetAnalyzer::Process(const edm::Eve
156                          localJet.setBtag_trackCountingHighEff(patJet->bDiscriminator("trackCountingHighEffBJetTags"));
157                          localJet.setBtag_trackCountingHighPur(patJet->bDiscriminator("trackCountingHighPurBJetTags"));
158                          localJet.setBtag_jetProbability(patJet->bDiscriminator("jetProbabilityBJetTags"));
159 <                        localJet.setBCorrection(patJet->correctionFactor(pat::Jet::bCorrection));
160 <                        localJet.setCCorrection(patJet->correctionFactor(pat::Jet::cCorrection));
161 <                        localJet.setUDSCorrection(patJet->correctionFactor(pat::Jet::udsCorrection));
162 <                        localJet.setGCorrection(patJet->correctionFactor(pat::Jet::gCorrection));
159 >
160 >                        // see DataFormats/PatCandidates/interface/JetCorrFactors.h
161 >                        localJet.setBCorrection(patJet->jetCorrFactors().correction(pat::JetCorrFactors::L7b,patJet->jetCorrStep() ));
162 >                        localJet.setCCorrection(patJet->jetCorrFactors().correction(pat::JetCorrFactors::L7c,patJet->jetCorrStep() ));
163 >                        localJet.setUDSCorrection(patJet->jetCorrFactors().correction(pat::JetCorrFactors::L7uds,patJet->jetCorrStep() ));
164 >                        localJet.setGCorrection(patJet->jetCorrFactors().correction(pat::JetCorrFactors::L7g,patJet->jetCorrStep() ));
165  
166                          // Use  associated tracks to calculate charged broadness of the jet
167                          // FIXME - Check generalTracks collection is present
# Line 241 | Line 243 | void JetAnalyzer::Process(const edm::Eve
243                                          localJet.setGenParticleIndex(-1);
244                                  }
245  
244 –                                // Matched generated jet
245 –                                // FIXME - Keep genJet index or UID instead
246 –                                // genJet stocked in MCParticles branch
247 –                                /*
248 –                                if ( patJet->genJet() != NULL )
249 –                                {
250 –                                        localJet.setMomentumMCJet(TLorentzVector(patJet->genJet()->px(),patJet->genJet()->py(),patJet->genJet()->pz(),patJet->genJet()->energy()));
251 –                                        localJet.setVertexMCJet(TVector3(patJet->genJet()->vx(),patJet->genJet()->vy(),patJet->genJet()->vz() ));
252 –                                        localJet.setPdgIdMCJet(patJet->genJet()->pdgId());
253 –                                }
254 –                                */
255 –
246                                  // check if jet comes from a top
247                                  bool IsTopJet =  false;
248                                  if(patJet->genParton())
# Line 272 | Line 262 | void JetAnalyzer::Process(const edm::Eve
262                                          }
263                                  }
264                                  localJet.setIsTopJet(IsTopJet);
265 +
266 +                                // Matched generated jet
267 +                                if ( patJet->genJet() != NULL )
268 +                                {
269 +                                        localJet.setGenJetEnergy( patJet->genJet()->energy() );
270 +                                }
271 +
272                          }
273                  }
274  

Diff Legend

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