ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/CmsHi/JetAnalysis/src/RecHitTreeProducer.cc

Comparing UserCode/CmsHi/JetAnalysis/src/RecHitTreeProducer.cc (file contents):
Revision 1.4 by yilmaz, Fri Oct 22 14:06:15 2010 UTC vs.
Revision 1.16 by yjlee, Tue Oct 25 13:41:08 2011 UTC

#	Line 12 | Line 12
12		*/
13		//
14		// Original Author:  Yetkin Yilmaz
15	+	// Modified: Frank Ma
16		//         Created:  Tue Sep  7 11:38:19 EDT 2010
17		// $Id$
18		//
19		//
20
21	<
21	>	#define versionTag "v1"
22		// system include files
23		#include <memory>
24		#include <vector>
#	Line 57 | Line 58
58		#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
59		#include "DataFormats/HcalDetId/interface/HcalDetId.h"
60		#include "DataFormats/EcalDetId/interface/EcalDetIdCollections.h"
61	+	#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
62		#include "DataFormats/DetId/interface/DetId.h"
63
64	+	#include "DataFormats/VertexReco/interface/Vertex.h"
65	+	#include "DataFormats/VertexReco/interface/VertexFwd.h"
66	+
67
68		#include "TNtuple.h"
69
#	Line 67 | Line 72 | using namespace std;
72		#define MAXHITS 1000000
73
74		struct MyRecHit{
75	<
75	>	int depth[MAXHITS];
76		int n;
77
78		float e[MAXHITS];
#	Line 75 | Line 80 | struct MyRecHit{
80		float eta[MAXHITS];
81		float phi[MAXHITS];
82		bool isjet[MAXHITS];
83	+	float etvtx[MAXHITS];
84	+	float etavtx[MAXHITS];
85	+	float emEtVtx[MAXHITS];
86	+	float hadEtVtx[MAXHITS];
87	+
88	+	float jtpt;
89	+	float jteta;
90	+	float jtphi;
91
92		};
93
94
95	+	struct MyBkg{
96	+	int n;
97	+	float rho[50];
98	+	float sigma[50];
99	+	};
100	+
101
102		//
103		// class declaration
#	Line 91 | Line 110 | class RecHitTreeProducer : public edm::E
110		double       getEt(const DetId &id, double energy);
111		double       getEta(const DetId &id);
112		double       getPhi(const DetId &id);
113	+	reco::Vertex::Point getVtx(const edm::Event& ev);
114	+
115
116
117		private:
#	Line 110 | Line 131 | class RecHitTreeProducer : public edm::E
131		edm::Handle<HFRecHitCollection> hfHits;
132		edm::Handle<HBHERecHitCollection> hbheHits;
133
134	+	edm::Handle<reco::BasicClusterCollection> bClusters;
135	+	edm::Handle<CaloTowerCollection> towers;
136	+	edm::Handle<reco::VertexCollection> vtxs;
137	+
138		typedef vector<EcalRecHit>::const_iterator EcalIterator;
139
140		edm::Handle<reco::CaloJetCollection> jets;
141	+
142	+	edm::Handle<std::vector<double> > rhos;
143	+	edm::Handle<std::vector<double> > sigmas;
144
145		MyRecHit hbheRecHit;
146		MyRecHit hfRecHit;
147		MyRecHit ebRecHit;
148		MyRecHit eeRecHit;
149	+	MyRecHit myBC;
150	+	MyRecHit myTowers;
151	+	MyBkg bkg;
152
153	+	TNtuple* nt;
154		TTree* hbheTree;
155		TTree* hfTree;
156		TTree* ebTree;
157		TTree* eeTree;
158	<	double cone;
158	>	TTree* bcTree;
159	>	TTree* towerTree;
160	>	TTree* bkgTree;
161
162	+	double cone;
163	+	double hfTowerThreshold_;
164	+	double hfLongThreshold_;
165	+	double hfShortThreshold_;
166	+	double hbheThreshold_;
167	+	double ebThreshold_;
168	+	double eeThreshold_;
169	+
170	+	double hbhePtMin_;
171	+	double hfPtMin_;
172	+	double ebPtMin_;
173	+	double eePtMin_;
174	+	double towerPtMin_;
175	+
176		edm::Service<TFileService> fs;
177		const CentralityBins * cbins_;
178		const CaloGeometry *geo;
#	Line 133 | Line 181 | class RecHitTreeProducer : public edm::E
181		edm::InputTag HcalRecHitHBHESrc_;
182		edm::InputTag EBSrc_;
183		edm::InputTag EESrc_;
184	+	edm::InputTag BCSrc_;
185	+	edm::InputTag TowerSrc_;
186	+	edm::InputTag VtxSrc_;
187	+
188		edm::InputTag JetSrc_;
137	–	bool excludeJets_;
189
190	+	edm::InputTag FastJetTag_;
191	+
192	+	bool useJets_;
193	+	bool doBasicClusters_;
194	+	bool doTowers_;
195	+	bool doEcal_;
196	+	bool doHcal_;
197	+	bool hasVtx_;
198	+
199	+	bool doFastJet_;
200	+
201	+	bool doEbyEonly_;
202		};
203
204		//
#	Line 155 | Line 218 | RecHitTreeProducer::RecHitTreeProducer(c
218		cone(0.5)
219		{
220		//now do what ever initialization is needed
158	–
221		EBSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("EBRecHitSrc",edm::InputTag("ecalRecHit","EcalRecHitsEB"));
222		EESrc_ = iConfig.getUntrackedParameter<edm::InputTag>("EERecHitSrc",edm::InputTag("ecalRecHit","EcalRecHitsEE"));
223		HcalRecHitHFSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHFRecHitSrc",edm::InputTag("hfreco"));
224		HcalRecHitHBHESrc_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHBHERecHitSrc",edm::InputTag("hbhereco"));
225	<	JetSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("JetSrc",edm::InputTag("iterativeCone5CaloJets"));
226	<	excludeJets_ = iConfig.getUntrackedParameter<bool>("excludeJets",false);
225	>	BCSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("BasicClusterSrc1",edm::InputTag("ecalRecHit","EcalRecHitsEB","RECO"));
226	>	TowerSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("towersSrc",edm::InputTag("towerMaker"));
227	>	VtxSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("vtxSrc",edm::InputTag("hiSelectedVertex"));
228	>	JetSrc_ = iConfig.getUntrackedParameter<edm::InputTag>("JetSrc",edm::InputTag("iterativeConePu5CaloJets"));
229	>	useJets_ = iConfig.getUntrackedParameter<bool>("useJets",true);
230	>	doBasicClusters_ = iConfig.getUntrackedParameter<bool>("doBasicClusters",false);
231	>	doTowers_ = iConfig.getUntrackedParameter<bool>("doTowers",true);
232	>	doEcal_ = iConfig.getUntrackedParameter<bool>("doEcal",true);
233	>	doHcal_ = iConfig.getUntrackedParameter<bool>("doHcal",true);
234	>	hasVtx_ = iConfig.getUntrackedParameter<bool>("hasVtx",false);
235	>	doFastJet_ = iConfig.getUntrackedParameter<bool>("doFastJet",true);
236	>	FastJetTag_ = iConfig.getUntrackedParameter<edm::InputTag>("FastJetTag",edm::InputTag("kt4CaloJets"));
237	>	doEbyEonly_ = iConfig.getUntrackedParameter<bool>("doEbyEonly",false);
238	>	hfTowerThreshold_ = iConfig.getUntrackedParameter<double>("HFtowerMin",3.);
239	>	hfLongThreshold_ = iConfig.getUntrackedParameter<double>("HFlongMin",0.5);
240	>	hfShortThreshold_ = iConfig.getUntrackedParameter<double>("HFshortMin",0.85);
241	>	hbhePtMin_ = iConfig.getUntrackedParameter<double>("HBHETreePtMin",0);
242	>	hfPtMin_ = iConfig.getUntrackedParameter<double>("HFTreePtMin",0);
243	>	ebPtMin_ = iConfig.getUntrackedParameter<double>("EBTreePtMin",0);
244	>	eePtMin_ = iConfig.getUntrackedParameter<double>("EETreePtMin",0.);
245	>	towerPtMin_ = iConfig.getUntrackedParameter<double>("TowerTreePtMin",0.);
246		}
247
248
#	Line 187 | Line 268 | RecHitTreeProducer::analyze(const edm::E
268		hbheRecHit.n = 0;
269		ebRecHit.n = 0;
270		eeRecHit.n = 0;
271	+	myBC.n = 0;
272	+	myTowers.n = 0;
273	+	bkg.n = 0;
274	+
275	+	// get vertex
276	+	reco::Vertex::Point vtx;
277	+	if (hasVtx_) vtx = getVtx(ev);
278
279	+	if(doEcal_){
280		ev.getByLabel(EBSrc_,ebHits);
281		ev.getByLabel(EESrc_,eeHits);
282	<
282	>	}
283	>	if(doHcal_){
284		ev.getByLabel(HcalRecHitHFSrc_,hfHits);
285		ev.getByLabel(HcalRecHitHBHESrc_,hbheHits);
286	+	}
287	+	if(useJets_) {
288	+	ev.getByLabel(JetSrc_,jets);
289	+	}
290
291	+	if(doBasicClusters_){
292	+	ev.getByLabel(BCSrc_,bClusters);
293	+	}
294
295	<	if(!excludeJets_) {
296	<	ev.getByLabel(JetSrc_,jets);
295	>	if(doTowers_){
296	>	ev.getByLabel(TowerSrc_,towers);
297	>	}
298	>
299	>	if(doFastJet_){
300	>	ev.getByLabel(edm::InputTag(FastJetTag_.label(),"rhos",FastJetTag_.process()),rhos);
301	>	ev.getByLabel(edm::InputTag(FastJetTag_.label(),"sigmas",FastJetTag_.process()),sigmas);
302	>	bkg.n = rhos->size();
303	>	for(unsigned int i = 0; i < rhos->size(); ++i){
304	>	bkg.rho[i] = (*rhos)[i];
305	>	bkg.sigma[i] = (*sigmas)[i];
306	>	}
307		}
308
309	<	if(0 && !cbins_) cbins_ = getCentralityBinsFromDB(iSetup);
309	>	if(0 && !cbins_) cbins_ = getCentralityBinsFromDB(iSetup);
310
311		if(!geo){
312		edm::ESHandle<CaloGeometry> pGeo;
313		iSetup.get<CaloGeometryRecord>().get(pGeo);
314		geo = pGeo.product();
315		}
316	+
317	+	int nHFlongPlus = 0;
318	+	int nHFshortPlus = 0;
319	+	int nHFtowerPlus = 0;
320	+	int nHFlongMinus = 0;
321	+	int nHFshortMinus = 0;
322	+	int nHFtowerMinus = 0;
323	+
324	+
325
326	+	if(doHcal_){
327		for(unsigned int i = 0; i < hfHits->size(); ++i){
328		const HFRecHit & hit= (*hfHits)[i];
329		hfRecHit.e[hfRecHit.n] = hit.energy();
#	Line 214 | Line 331 | RecHitTreeProducer::analyze(const edm::E
331		hfRecHit.eta[hfRecHit.n] = getEta(hit.id());
332		hfRecHit.phi[hfRecHit.n] = getPhi(hit.id());
333		hfRecHit.isjet[hfRecHit.n] = false;
334	<	if(!excludeJets_){
334	>	hfRecHit.depth[hfRecHit.n] = hit.id().depth();
335	>
336	>	if(hit.id().ieta() > 0){
337	>	if(hit.energy() > hfShortThreshold_ && hit.id().depth() != 1) nHFshortPlus++;
338	>	if(hit.energy() > hfLongThreshold_ && hit.id().depth() == 1) nHFlongPlus++;
339	>	}else{
340	>	if(hit.energy() > hfShortThreshold_ && hit.id().depth() != 1) nHFshortMinus++;
341	>	if(hit.energy() > hfLongThreshold_ && hit.id().depth() == 1) nHFlongMinus++;
342	>	}
343	>
344	>	if(useJets_){
345		for(unsigned int j = 0 ; j < jets->size(); ++j){
346		const reco::Jet& jet = (*jets)[j];
347		double dr = reco::deltaR(hfRecHit.eta[hfRecHit.n],hfRecHit.phi[hfRecHit.n],jet.eta(),jet.phi());
348		if(dr < cone){ hfRecHit.isjet[hfRecHit.n] = true; }
349		}
350		}
351	<	hfRecHit.n++;
351	>	if (hfRecHit.et[hfRecHit.n]>=hfPtMin_) hfRecHit.n++;
352		}
353	<
353	>	if(!doEbyEonly_){
354		for(unsigned int i = 0; i < hbheHits->size(); ++i){
355		const HBHERecHit & hit= (*hbheHits)[i];
356	+	if (getEt(hit.id(),hit.energy())<hbhePtMin_) continue;
357		hbheRecHit.e[hbheRecHit.n] = hit.energy();
358		hbheRecHit.et[hbheRecHit.n] = getEt(hit.id(),hit.energy());
359		hbheRecHit.eta[hbheRecHit.n] = getEta(hit.id());
360		hbheRecHit.phi[hbheRecHit.n] = getPhi(hit.id());
361		hbheRecHit.isjet[hbheRecHit.n] = false;
362	<	if(!excludeJets_){
362	>	if(useJets_){
363		for(unsigned int j = 0 ; j < jets->size(); ++j){
364		const reco::Jet& jet = (*jets)[j];
365		double dr = reco::deltaR(hbheRecHit.eta[hbheRecHit.n],hbheRecHit.phi[hbheRecHit.n],jet.eta(),jet.phi());
#	Line 240 | Line 368 | RecHitTreeProducer::analyze(const edm::E
368		}
369		hbheRecHit.n++;
370		}
371	<
371	>	}
372	>	}
373	>	if(doEcal_ && !doEbyEonly_){
374		for(unsigned int i = 0; i < ebHits->size(); ++i){
375		const EcalRecHit & hit= (*ebHits)[i];
376	+	if (getEt(hit.id(),hit.energy())<ebPtMin_) continue;
377		ebRecHit.e[ebRecHit.n] = hit.energy();
378		ebRecHit.et[ebRecHit.n] = getEt(hit.id(),hit.energy());
379		ebRecHit.eta[ebRecHit.n] = getEta(hit.id());
380		ebRecHit.phi[ebRecHit.n] = getPhi(hit.id());
381		ebRecHit.isjet[ebRecHit.n] = false;
382	<	if(!excludeJets_){
382	>	if(useJets_){
383		for(unsigned int j = 0 ; j < jets->size(); ++j){
384		const reco::Jet& jet = (*jets)[j];
385		double dr = reco::deltaR(ebRecHit.eta[ebRecHit.n],ebRecHit.phi[ebRecHit.n],jet.eta(),jet.phi());
#	Line 260 | Line 391 | RecHitTreeProducer::analyze(const edm::E
391
392		for(unsigned int i = 0; i < eeHits->size(); ++i){
393		const EcalRecHit & hit= (*eeHits)[i];
394	+	if (getEt(hit.id(),hit.energy())<eePtMin_) continue;
395		eeRecHit.e[eeRecHit.n] = hit.energy();
396		eeRecHit.et[eeRecHit.n] = getEt(hit.id(),hit.energy());
397		eeRecHit.eta[eeRecHit.n] = getEta(hit.id());
398		eeRecHit.phi[eeRecHit.n] = getPhi(hit.id());
399		eeRecHit.isjet[eeRecHit.n] = false;
400	<	if(!excludeJets_){
400	>	if(useJets_){
401		for(unsigned int j = 0 ; j < jets->size(); ++j){
402		const reco::Jet& jet = (*jets)[j];
403		double dr = reco::deltaR(eeRecHit.eta[eeRecHit.n],eeRecHit.phi[eeRecHit.n],jet.eta(),jet.phi());
#	Line 274 | Line 406 | RecHitTreeProducer::analyze(const edm::E
406		}
407		eeRecHit.n++;
408		}
409	<
410	<	eeTree->Fill();
411	<	ebTree->Fill();
412	<
413	<	hbheTree->Fill();
414	<	hfTree->Fill();
409	>	}
410	>
411	>	if(doTowers_){
412	>
413	>	for(unsigned int i = 0; i < towers->size(); ++i){
414	>	const CaloTower & hit= (*towers)[i];
415	>	if (getEt(hit.id(),hit.energy())<towerPtMin_) continue;
416	>	myTowers.e[myTowers.n] = hit.energy();
417	>	myTowers.et[myTowers.n] = getEt(hit.id(),hit.energy());
418	>	myTowers.eta[myTowers.n] = getEta(hit.id());
419	>	myTowers.phi[myTowers.n] = getPhi(hit.id());
420	>	myTowers.isjet[myTowers.n] = false;
421	>	if (hasVtx_) {
422	>	myTowers.etvtx[myTowers.n] = hit.p4(vtx).Et();
423	>	myTowers.etavtx[myTowers.n] = hit.p4(vtx).Eta();
424	>	myTowers.emEtVtx[myTowers.n] = hit.emEt(vtx);
425	>	myTowers.hadEtVtx[myTowers.n] = hit.hadEt(vtx);
426	>	}
427	>
428	>	if(hit.ieta() > 29 && hit.energy() > hfTowerThreshold_) nHFtowerPlus++;
429	>	if(hit.ieta() < -29 && hit.energy() > hfTowerThreshold_) nHFtowerMinus++;
430	>
431	>	if(useJets_){
432	>	for(unsigned int j = 0 ; j < jets->size(); ++j){
433	>	const reco::Jet& jet = (*jets)[j];
434	>	double dr = reco::deltaR(myTowers.eta[myTowers.n],myTowers.phi[myTowers.n],jet.eta(),jet.phi());
435	>	if(dr < cone){ myTowers.isjet[myTowers.n] = true; }
436	>	}
437	>	}
438	>	myTowers.n++;
439	>	}
440	>
441	>	}
442	>
443	>	if(doBasicClusters_ && !doEbyEonly_){
444	>	for(unsigned int j = 0 ; j < jets->size(); ++j){
445	>	const reco::Jet& jet = (*jets)[j];
446	>	myBC.n = 0;
447	>	myBC.jtpt = jet.pt();
448	>	myBC.jteta = jet.eta();
449	>	myBC.jtphi = jet.phi();
450	>
451	>	for(unsigned int i = 0; i < bClusters->size(); ++i){
452	>	const reco::BasicCluster & bc= (*bClusters)[i];
453	>	double dr = reco::deltaR(bc.eta(),bc.phi(),jet.eta(),jet.phi());
454	>	if(dr < cone){
455	>	myBC.e[myBC.n] = bc.energy();
456	>	myBC.et[myBC.n] = bc.energy()*sin(bc.position().theta());
457	>	myBC.eta[myBC.n] = bc.eta();
458	>	myBC.phi[myBC.n] = bc.phi();
459	>	myBC.n++;
460	>	}
461	>	}
462	>	bcTree->Fill();
463	>	}
464	>	}
465	>
466	>	if(!doEbyEonly_){
467	>	towerTree->Fill();
468	>
469	>	eeTree->Fill();
470	>	ebTree->Fill();
471	>
472	>	hbheTree->Fill();
473	>	hfTree->Fill();
474	>
475	>	if (doFastJet_) {
476	>	bkgTree->Fill();
477	>	}
478	>	}
479	>
480	>	nt->Fill(nHFtowerPlus,nHFtowerMinus,nHFlongPlus,nHFlongMinus,nHFshortPlus,nHFshortMinus);
481
482		}
483
#	Line 289 | Line 487 | void
487		RecHitTreeProducer::beginJob()
488		{
489
490	<	hbheTree = fs->make<TTree>("hbhe","");
490	>	hbheTree = fs->make<TTree>("hbhe",versionTag);
491		hbheTree->Branch("n",&hbheRecHit.n,"n/I");
492		hbheTree->Branch("e",hbheRecHit.e,"e[n]/F");
493		hbheTree->Branch("et",hbheRecHit.et,"et[n]/F");
494		hbheTree->Branch("eta",hbheRecHit.eta,"eta[n]/F");
495		hbheTree->Branch("phi",hbheRecHit.phi,"phi[n]/F");
496	<	hbheTree->Branch("isjet",hbheRecHit.isjet,"isjet[n]/I");
496	>	hbheTree->Branch("isjet",hbheRecHit.isjet,"isjet[n]/O");
497
498	<	hfTree = fs->make<TTree>("hf","");
498	>	hfTree = fs->make<TTree>("hf",versionTag);
499		hfTree->Branch("n",&hfRecHit.n,"n/I");
500		hfTree->Branch("e",hfRecHit.e,"e[n]/F");
501		hfTree->Branch("et",hfRecHit.et,"et[n]/F");
502		hfTree->Branch("eta",hfRecHit.eta,"eta[n]/F");
503		hfTree->Branch("phi",hfRecHit.phi,"phi[n]/F");
504	<	hfTree->Branch("isjet",hfRecHit.isjet,"isjet[n]/I");
504	>	hfTree->Branch("depth",hfRecHit.depth,"depth[n]/I");
505	>	hfTree->Branch("isjet",hfRecHit.isjet,"isjet[n]/O");
506
507	<	eeTree = fs->make<TTree>("ee","");
507	>	eeTree = fs->make<TTree>("ee",versionTag);
508		eeTree->Branch("n",&eeRecHit.n,"n/I");
509		eeTree->Branch("e",eeRecHit.e,"e[n]/F");
510		eeTree->Branch("et",eeRecHit.et,"et[n]/F");
511		eeTree->Branch("eta",eeRecHit.eta,"eta[n]/F");
512		eeTree->Branch("phi",eeRecHit.phi,"phi[n]/F");
513	<	eeTree->Branch("isjet",eeRecHit.isjet,"isjet[n]/I");
513	>	eeTree->Branch("isjet",eeRecHit.isjet,"isjet[n]/O");
514
515	<	ebTree = fs->make<TTree>("eb","");
515	>	ebTree = fs->make<TTree>("eb",versionTag);
516		ebTree->Branch("n",&ebRecHit.n,"n/I");
517		ebTree->Branch("e",ebRecHit.e,"e[n]/F");
518		ebTree->Branch("et",ebRecHit.et,"et[n]/F");
#	Line 321 | Line 520 | RecHitTreeProducer::beginJob()
520		ebTree->Branch("phi",ebRecHit.phi,"phi[n]/F");
521		ebTree->Branch("isjet",ebRecHit.isjet,"isjet[n]/O");
522
523	+	towerTree = fs->make<TTree>("tower",versionTag);
524	+	towerTree->Branch("n",&myTowers.n,"n/I");
525	+	towerTree->Branch("e",myTowers.e,"e[n]/F");
526	+	towerTree->Branch("et",myTowers.et,"et[n]/F");
527	+	towerTree->Branch("eta",myTowers.eta,"eta[n]/F");
528	+	towerTree->Branch("phi",myTowers.phi,"phi[n]/F");
529	+	towerTree->Branch("isjet",myTowers.isjet,"isjet[n]/O");
530	+	if (hasVtx_) {
531	+	towerTree->Branch("etvtx",myTowers.etvtx,"etvtx[n]/F");
532	+	towerTree->Branch("etavtx",myTowers.etavtx,"etavtx[n]/F");
533	+	towerTree->Branch("emEtVtx",myTowers.emEtVtx,"emEtVtx[n]/F");
534	+	towerTree->Branch("hadEtVtx",myTowers.hadEtVtx,"hadEtVtx[n]/F");
535	+	}
536	+
537	+
538	+	if(doBasicClusters_){
539	+	bcTree = fs->make<TTree>("bc",versionTag);
540	+	bcTree->Branch("n",&myBC.n,"n/I");
541	+	bcTree->Branch("e",myBC.e,"e[n]/F");
542	+	bcTree->Branch("et",myBC.et,"et[n]/F");
543	+	bcTree->Branch("eta",myBC.eta,"eta[n]/F");
544	+	bcTree->Branch("phi",myBC.phi,"phi[n]/F");
545	+	bcTree->Branch("jtpt",&myBC.jtpt,"jtpt/F");
546	+	bcTree->Branch("jteta",&myBC.jteta,"jteta/F");
547	+	bcTree->Branch("jtphi",&myBC.jtphi,"jtphi/F");
548	+	//     bcTree->Branch("isjet",bcRecHit.isjet,"isjet[n]/O");
549	+	}
550	+
551	+	if(doFastJet_){
552	+	bkgTree = fs->make<TTree>("bkg",versionTag);
553	+	bkgTree->Branch("n",&bkg.n,"n/I");
554	+	bkgTree->Branch("rho",bkg.rho,"rho[n]/F");
555	+	bkgTree->Branch("sigma",bkg.sigma,"sigma[n]/F");
556	+	}
557	+
558	+	nt = fs->make<TNtuple>("ntEvent","","nHFplus:nHFminus:nHFlongPlus:nHFlongMinus:nHFshortPlus:nHFshortMinus");
559	+
560		}
561
562		// ------------ method called once each job just after ending the event loop  ------------
#	Line 346 | Line 582 | double RecHitTreeProducer::getPhi(const
582		return et;
583		}
584
585	<
585	>	reco::Vertex::Point RecHitTreeProducer::getVtx(const edm::Event& ev)
586	>	{
587	>	ev.getByLabel(VtxSrc_,vtxs);
588	>	int greatestvtx = 0;
589	>	int nVertex = vtxs->size();
590	>
591	>	for (unsigned int i = 0 ; i< vtxs->size(); ++i){
592	>	unsigned int daughter = (*vtxs)[i].tracksSize();
593	>	if( daughter > (*vtxs)[greatestvtx].tracksSize()) greatestvtx = i;
594	>	//cout <<"Vertex: "<< (*vtxs)[i].position().z()<<" "<<daughter<<endl;
595	>	}
596	>
597	>	if(nVertex<=0){
598	>	return reco::Vertex::Point(-999,-999,-999);
599	>	}
600	>	return (*vtxs)[greatestvtx].position();
601	>	}
602
603		//define this as a plug-in
604		DEFINE_FWK_MODULE(RecHitTreeProducer);

Diff Legend

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