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

Comparing UserCode/MitPhysics/Validation/src/JetValidationMod.cc (file contents):
Revision 1.4 by loizides, Mon Jun 15 15:00:23 2009 UTC vs.
Revision 1.5 by loizides, Fri Jul 10 14:04:44 2009 UTC

# Line 119 | Line 119 | void JetValidationMod::Process()
119        
120        fIC5GenJetRecoJetDeltaR->Fill(dR);
121        fIC5GenJetRecoJetDeltaEta->Fill(fabs(GoodIC5GenJets->At(i)->Eta() - iter->second->Eta()));
122 <      fIC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodIC5GenJets->At(i)->Phi(),iter->second->Phi()));
122 >      fIC5GenJetRecoJetDeltaPhi->Fill(MathUtils::DeltaPhi(GoodIC5GenJets->At(i)->Phi(),
123 >                                                          iter->second->Phi()));
124        
125        if (fabs(iter->second->Eta()) <= 1.4)
126 <        fIC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());
126 >        fIC5JetResponseVsGenJetPtInBarrel->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
127 >                                                GoodIC5GenJets->At(i)->Pt());
128        else if (fabs(iter->second->Eta()) <= 3.0)
129 <        fIC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());
129 >        fIC5JetResponseVsGenJetPtInEndcap->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
130 >                                                GoodIC5GenJets->At(i)->Pt());
131        else if (fabs(iter->second->Eta()) <= 5.0)
132 <        fIC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),GoodIC5GenJets->At(i)->Pt());                
132 >        fIC5JetResponseVsGenJetPtForward->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
133 >                                               GoodIC5GenJets->At(i)->Pt());                
134        
135        //avoid the low pt region
136        if (iter->second->Pt() > 30)
137 <        fIC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),iter->second->Eta());
137 >        fIC5JetResponseVsCaloJetEta->Fill(iter->second->Et()/GoodIC5GenJets->At(i)->Et(),
138 >                                          iter->second->Eta());
139        
140        if (fabs(iter->second->Eta()) <= 3.0) {
141          fIC5CentralGenJetRecoJetDeltaR->Fill(dR);
# Line 148 | Line 153 | void JetValidationMod::Process()
153        if (GoodIC5GenJets->At(i)->Pt() > 60.0 && GoodIC5GenJets->At(i)->Pt() < 80.0)
154          fIC5NMatchedCaloJetsVsGenJetEta_Pt60To80->Fill(GoodIC5GenJets->At(i)->Eta());
155      
156 <      fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() * iter->second->L2RelativeCorrectionScale() * iter->second->L3AbsoluteCorrectionScale() / GoodIC5GenJets->At(i)->Pt(), GoodIC5GenJets->At(i)->Pt());
156 >      fIC5CorrPtCaloJetsOverGenJetsPtVsGenJetPt->Fill(iter->second->Pt() *
157 >                                                      iter->second->L2RelativeCorrectionScale() *
158 >                                                      iter->second->L3AbsoluteCorrectionScale() /
159 >                                                      GoodIC5GenJets->At(i)->Pt(),
160 >                                                      GoodIC5GenJets->At(i)->Pt());
161      }            
162    }
163  
# Line 157 | Line 166 | void JetValidationMod::Process()
166        IC5JetToIC5GenJetMap.find(GoodIC5Jets->At(i));      
167      //unmatched ones
168      if (iter == IC5JetToIC5GenJetMap.end()) {          
169 <      fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(GoodIC5Jets->At(i)->Pt() * GoodIC5Jets->At(i)->L2RelativeCorrectionScale() * GoodIC5Jets->At(i)->L3AbsoluteCorrectionScale());
169 >      fIC5NUnmatchedCaloJetsVsCorrectedCaloJetPt->Fill(
170 >        GoodIC5Jets->At(i)->Pt() * GoodIC5Jets->At(i)->L2RelativeCorrectionScale() *
171 >        GoodIC5Jets->At(i)->L3AbsoluteCorrectionScale());
172        fIC5NUnmatchedCalojetsVsCorrectedCaloJetEta->Fill(GoodIC5Jets->At(i)->Eta());
173      }
174    }

Diff Legend

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