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

Comparing UserCode/CmsHi/JetAnalysis/src/RecHitComparison.cc (file contents):
Revision 1.3 by yilmaz, Wed Oct 20 14:11:46 2010 UTC vs.
Revision 1.5 by yilmaz, Fri Oct 22 13:34:09 2010 UTC

# Line 106 | Line 106 | class RecHitComparison : public edm::EDA
106    edm::InputTag EESrc1_;
107    edm::InputTag EESrc2_;
108     edm::InputTag signalTag_;
109 +   edm::InputTag centTag_;
110  
111     TNtuple* ntEB;
112    TNtuple* ntEE;
# Line 114 | Line 115 | class RecHitComparison : public edm::EDA
115     TNtuple* ntjet;
116  
117     double cone;
118 +   bool jetsOnly_;
119  
120     edm::Service<TFileService> fs;
121     const CentralityBins * cbins_;
# Line 137 | Line 139 | RecHitComparison::RecHitComparison(const
139     cone(0.5)
140   {
141     //now do what ever initialization is needed
142 +   centTag_ =  iConfig.getUntrackedParameter<edm::InputTag>("centrality",edm::InputTag("hiCentrality","","RECO"));
143 +
144 +   jetsOnly_ = iConfig.getUntrackedParameter<bool>("jetsOnly",false);
145     signalTag_ = iConfig.getUntrackedParameter<edm::InputTag>("signalJets",edm::InputTag("iterativeCone5CaloJets","","SIGNAL"));
146  
147    HcalRecHitHFSrc1_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHFRecHitSrc1",edm::InputTag("hfreco"));
# Line 176 | Line 181 | RecHitComparison::analyze(const edm::Eve
181        geo = pGeo.product();
182     }
183  
179 –   edm::InputTag centTag("hiCentrality","","RECO");
180 –
184     using namespace edm;
185  
186 <   ev.getByLabel(centTag,cent);
186 >   ev.getByLabel(centTag_,cent);
187     ev.getByLabel(EBSrc1_,ebHits1);
188     ev.getByLabel(EBSrc1_,ebHits2);
189     ev.getByLabel(signalTag_,signalJets);
# Line 287 | Line 290 | RecHitComparison::analyze(const edm::Eve
290        double eta2 = pos2.eta();
291        double phi2 = pos2.eta();
292        double et2 = e2*sin(pos2.theta());
293 <      if(isjet) ntEB->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
293 >      if(!jetsOnly_ ||  isjet) ntEB->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
294     }
295  
296     for(unsigned int i = 0; i < eeHits1->size(); ++i){
# Line 324 | Line 327 | RecHitComparison::analyze(const edm::Eve
327       double eta2 = pos2.eta();
328       double phi2 = pos2.eta();
329       double et2 = e2*sin(pos2.theta());
330 <     if(isjet) ntEE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
330 >     if(!jetsOnly_ || isjet) ntEE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
331     }
332    
333     for(unsigned int i = 0; i < hbheHits1->size(); ++i){
# Line 361 | Line 364 | RecHitComparison::analyze(const edm::Eve
364       double eta2 = pos2.eta();
365       double phi2 = pos2.eta();
366       double et2 = e2*sin(pos2.theta());
367 <     if(isjet) ntHBHE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
367 >     if(!jetsOnly_ || isjet) ntHBHE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
368     }
369  
370     for(unsigned int i = 0; i < hfHits1->size(); ++i){
# Line 398 | Line 401 | RecHitComparison::analyze(const edm::Eve
401       double eta2 = pos2.eta();
402       double phi2 = pos2.eta();
403       double et2 = e2*sin(pos2.theta());
404 <     if(isjet) ntHF->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
404 >     if(!jetsOnly_ || isjet) ntHF->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet);
405     }
406  
407     for(unsigned int j1 = 0 ; j1 < signalJets->size(); ++j1){

Diff Legend

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