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

Comparing UserCode/CmsHi/JetAnalysis/src/RecHitTreeProducer.cc (file contents):
Revision 1.1 by yilmaz, Sat Oct 16 13:46:13 2010 UTC vs.
Revision 1.2 by yilmaz, Mon Oct 18 16:13:37 2010 UTC

# Line 104 | Line 104 | class RecHitTreeProducer : public edm::E
104     edm::Handle<vector<double> > ktRhos;
105     edm::Handle<vector<double> > akRhos;
106  
107 <   edm::Handle<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > jets1;
108 <   edm::Handle<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > jets2;
107 >  edm::Handle<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > ebHits;
108 >  edm::Handle<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > eeHits;
109 >
110 >  edm::Handle<HFRecHitCollection> hfHits;
111 >  edm::Handle<HBHERecHitCollection> hbheHits;
112  
110 –   //   typedef edm::Handle<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > >::const_iterator EcalIterator;
113     typedef vector<EcalRecHit>::const_iterator EcalIterator;
114  
115     edm::Handle<reco::CaloJetCollection> signalJets;
116  
117    MyRecHit hbheRecHit;
118    MyRecHit hfRecHit;
119 +  MyRecHit ebRecHit;
120 +  MyRecHit eeRecHit;
121  
122    TTree* hbheTree;
123    TTree* hfTree;
124 <
125 <   double cone;
124 >  TTree* ebTree;
125 >  TTree* eeTree;
126 >  double cone;
127  
128     edm::Service<TFileService> fs;
129     const CentralityBins * cbins_;
# Line 126 | Line 131 | class RecHitTreeProducer : public edm::E
131  
132    edm::InputTag HcalRecHitHFSrc_;
133    edm::InputTag HcalRecHitHBHESrc_;
134 +  edm::InputTag EBSrc_;
135 +  edm::InputTag EESrc_;
136  
137  
138   };
# Line 147 | Line 154 | RecHitTreeProducer::RecHitTreeProducer(c
154     cone(0.5)
155   {
156     //now do what ever initialization is needed
157 +
158 +  EBSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("EBRecHitSrc",edm::InputTag("ecalRecHit","EcalRecHitsEB"));
159 +  EESrc_ = iConfig.getUntrackedParameter<edm::InputTag>("EERecHitSrc",edm::InputTag("ecalRecHit","EcalRecHitsEE"));
160    HcalRecHitHFSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHFRecHitSrc",edm::InputTag("hfreco"));
161    HcalRecHitHBHESrc_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHBHERecHitSrc",edm::InputTag("hbhereco"));
162   }
# Line 172 | Line 182 | RecHitTreeProducer::analyze(const edm::E
182  
183    hfRecHit.n = 0;
184    hbheRecHit.n = 0;
185 +  ebRecHit.n = 0;
186 +  eeRecHit.n = 0;
187  
188 <  edm::Handle<HFRecHitCollection> hfHits;
189 <  edm::Handle<HBHERecHitCollection> hbheHits;
188 >  ev.getByLabel(EBSrc_,ebHits);
189 >  ev.getByLabel(EESrc_,eeHits);
190  
191    ev.getByLabel(HcalRecHitHFSrc_,hfHits);
192    ev.getByLabel(HcalRecHitHBHESrc_,hbheHits);
193  
182 –
194     if(0 && !cbins_) cbins_ = getCentralityBinsFromDB(iSetup);
195  
196     if(!geo){
# Line 194 | Line 205 | RecHitTreeProducer::analyze(const edm::E
205       hfRecHit.et[hfRecHit.n] = getEt(hit.id(),hit.energy());
206       hfRecHit.eta[hfRecHit.n] = getEta(hit.id());
207       hfRecHit.phi[hfRecHit.n] = getPhi(hit.id());
197 –     //     hfRecHit.ieta[hfRecHit.n] = hit.id().ieta();
198 –     //     hfRecHit.iphi[hfRecHit.n] = hit.id().iphi();
199 –    
208       hfRecHit.n++;
209     }
210  
# Line 206 | Line 214 | RecHitTreeProducer::analyze(const edm::E
214       hbheRecHit.et[hfRecHit.n] = getEt(hit.id(),hit.energy());
215       hbheRecHit.eta[hfRecHit.n] = getEta(hit.id());
216       hbheRecHit.phi[hfRecHit.n] = getPhi(hit.id());
209 –     //     hbheRecHit.ieta[hfRecHit.n] = hit.id().ieta();
210 –     //     hbheRecHit.iphi[hfRecHit.n] = hit.id().iphi();
211 –
217       hbheRecHit.n++;
218     }
219  
220 +   for(unsigned int i = 0; i < ebHits->size(); ++i){
221 +     const EcalRecHit & hit= (*ebHits)[i];
222 +     ebRecHit.e[ebRecHit.n] = hit.energy();
223 +     ebRecHit.et[ebRecHit.n] = getEt(hit.id(),hit.energy());
224 +     ebRecHit.eta[ebRecHit.n] = getEta(hit.id());
225 +     ebRecHit.phi[ebRecHit.n] = getPhi(hit.id());
226 +     ebRecHit.n++;
227 +   }
228 +
229 +   for(unsigned int i = 0; i < eeHits->size(); ++i){
230 +     const EcalRecHit & hit= (*eeHits)[i];
231 +     eeRecHit.e[eeRecHit.n] = hit.energy();
232 +     eeRecHit.et[eeRecHit.n] = getEt(hit.id(),hit.energy());
233 +     eeRecHit.eta[eeRecHit.n] = getEta(hit.id());
234 +     eeRecHit.phi[eeRecHit.n] = getPhi(hit.id());
235 +     eeRecHit.n++;
236 +   }
237 +
238 +   eeTree->Fill();
239 +   ebTree->Fill();
240 +
241     hbheTree->Fill();
242     hfTree->Fill();
243  
# Line 237 | Line 263 | RecHitTreeProducer::beginJob()
263    hfTree->Branch("eta",hfRecHit.eta,"eta[n]/F");
264    hfTree->Branch("phi",hfRecHit.phi,"phi[n]/F");
265  
266 +  eeTree = fs->make<TTree>("ee","");
267 +  eeTree->Branch("n",&eeRecHit.n,"n/I");
268 +  eeTree->Branch("e",eeRecHit.e,"e[n]/F");
269 +  eeTree->Branch("et",eeRecHit.et,"et[n]/F");
270 +  eeTree->Branch("eta",eeRecHit.eta,"eta[n]/F");
271 +  eeTree->Branch("phi",eeRecHit.phi,"phi[n]/F");
272 +
273 +  ebTree = fs->make<TTree>("eb","");
274 +  ebTree->Branch("n",&ebRecHit.n,"n/I");
275 +  ebTree->Branch("e",ebRecHit.e,"e[n]/F");
276 +  ebTree->Branch("et",ebRecHit.et,"et[n]/F");
277 +  ebTree->Branch("eta",ebRecHit.eta,"eta[n]/F");
278 +  ebTree->Branch("phi",ebRecHit.phi,"phi[n]/F");
279 +
280   }
281  
282   // ------------ method called once each job just after ending the event loop  ------------

Diff Legend

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