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.6 by yilmaz, Sun Oct 24 14:39:32 2010 UTC vs.
Revision 1.7 by yilmaz, Mon Nov 1 21:48:31 2010 UTC

# Line 55 | Line 55
55   #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
56   #include "DataFormats/HcalDetId/interface/HcalDetId.h"
57  
58 + #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
59 +
60   #include "TNtuple.h"
61  
62   using namespace std;
# Line 88 | Line 90 | class RecHitComparison : public edm::EDA
90    edm::Handle<HBHERecHitCollection> hbheHits1;
91    edm::Handle<HBHERecHitCollection> hbheHits2;
92  
93 +   edm::Handle<reco::BasicClusterCollection> bClusters1;
94 +   edm::Handle<reco::BasicClusterCollection> bClusters2;
95 +
96 +
97     typedef vector<EcalRecHit>::const_iterator EcalIterator;
98    typedef vector<HFRecHit>::const_iterator HFIterator;
99    typedef vector<HBHERecHit>::const_iterator HBHEIterator;
# Line 102 | Line 108 | class RecHitComparison : public edm::EDA
108    edm::InputTag EBSrc2_;
109    edm::InputTag EESrc1_;
110    edm::InputTag EESrc2_;
111 +
112 +   edm::InputTag BCSrc1_;
113 +   edm::InputTag BCSrc2_;
114 +
115     edm::InputTag signalTag_;
116  
117 +   TNtuple* ntBC;
118     TNtuple* ntEB;
119    TNtuple* ntEE;
120    TNtuple* ntHBHE;
# Line 112 | Line 123 | class RecHitComparison : public edm::EDA
123  
124     double cone;
125     bool jetsOnly_;
126 +   bool doBasicClusters_;
127  
128     edm::Service<TFileService> fs;
129     CentralityProvider* centrality_;
# Line 136 | Line 148 | RecHitComparison::RecHitComparison(const
148   {
149     //now do what ever initialization is needed
150     jetsOnly_ = iConfig.getUntrackedParameter<bool>("jetsOnly",false);
151 +   doBasicClusters_ = iConfig.getUntrackedParameter<bool>("doBasicClusters",false);
152     signalTag_ = iConfig.getUntrackedParameter<edm::InputTag>("signalJets",edm::InputTag("iterativeCone5CaloJets","","SIGNAL"));
153  
154    HcalRecHitHFSrc1_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHFRecHitSrc1",edm::InputTag("hfreco"));
# Line 146 | Line 159 | RecHitComparison::RecHitComparison(const
159    EBSrc2_ = iConfig.getUntrackedParameter<edm::InputTag>("EBRecHitSrc2",edm::InputTag("ecalRecHit","EcalRecHitsEB","SIGNALONLY"));
160    EESrc1_ = iConfig.getUntrackedParameter<edm::InputTag>("EERecHitSrc1",edm::InputTag("ecalRecHit","EcalRecHitsEE","RECO"));
161    EESrc2_ = iConfig.getUntrackedParameter<edm::InputTag>("EERecHitSrc2",edm::InputTag("ecalRecHit","EcalRecHitsEE","SIGNALONLY"));
162 <
162 >  BCSrc1_ = iConfig.getUntrackedParameter<edm::InputTag>("BasicClusterSrc1",edm::InputTag("ecalRecHit","EcalRecHitsEB","RECO"));
163 >  BCSrc2_ = iConfig.getUntrackedParameter<edm::InputTag>("BasicClusterSrc2",edm::InputTag("ecalRecHit","EcalRecHitsEB","SIGNALONLY"));
164   }
165  
166  
# Line 187 | Line 201 | RecHitComparison::analyze(const edm::Eve
201     ev.getByLabel(EESrc1_,eeHits1);
202     ev.getByLabel(EESrc2_,eeHits2);
203  
204 +   if(doBasicClusters_){
205 +      ev.getByLabel(BCSrc1_,bClusters1);
206 +      ev.getByLabel(BCSrc2_,bClusters2);
207 +   }
208 +
209     centrality_->newEvent(ev,iSetup);
210     double hf = centrality_->centralityValue(ev);
211     int bin = centrality_->getBin(ev);
# Line 412 | Line 431 | RecHitComparison::beginJob()
431     ntEE = fs->make<TNtuple>("ntEE","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet");
432     ntHBHE = fs->make<TNtuple>("ntHBHE","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet");
433     ntHF = fs->make<TNtuple>("ntHF","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet");
434 +
435 +   ntBC = fs->make<TNtuple>("ntBC","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet");
436 +
437     ntjet = fs->make<TNtuple>("ntjet","","bin:pt:eta:ethit:f05:f1:f15:f2:f25:f3:em:emf");
438    
439   }

Diff Legend

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