ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/Morgan/src/ZeeVertexAnalyzer.cc

Comparing UserCode/Morgan/src/ZeeVertexAnalyzer.cc (file contents):
Revision 1.1 by lethuill, Sat Oct 10 20:35:19 2009 UTC vs.
Revision 1.3 by lethuill, Fri Oct 23 14:25:45 2009 UTC

#	Line 19 | Line 19 | ZeeVertexAnalyzer::~ZeeVertexAnalyzer()
19		}
20
21
22	<	bool ZeeVertexAnalyzer::getVertices(const edm::Event& iEvent, const edm::EventSetup& iSetup, TClonesArray* rootVertices)
22	>	bool ZeeVertexAnalyzer::getVertices(const edm::Event& iEvent, const edm::EventSetup& iSetup, TClonesArray* rootVertices, TRootBardak* rootBardak)
23		{
24		using namespace edm;
25		int iRootVertex = 0;
#	Line 74 | Line 74 | bool ZeeVertexAnalyzer::getVertices(cons
74		if(verbosity_>4) std::cout << "     Electron 2 - (Et,eta,phi)=(" << electrons->at(iElectron2).pt() << "," << electrons->at(iElectron2).eta() << "," << electrons->at(iElectron2).phi() << ")" << std::endl;
75		if(verbosity_>4) std::cout << "     Mee=" << p4Z.mass() << std::endl;
76
77	<	if( p4Z.mass()<60.0 || p4Z.mass()>120.0) return false;
77	>	//if( p4Z.mass()<60.0 || p4Z.mass()>120.0) return false;
78
79		// Find CTF track best matching the GSF track assoicated to the electron
80		// We will remove these tracks in the track collection for the primary vertex reconstruction
81		reco::TrackRef tk1 = (electrons->at(iElectron1)).closestCtfTrackRef();
82		reco::TrackRef tk2 = (electrons->at(iElectron2)).closestCtfTrackRef();
83	–	if(verbosity_>4) std::cout << "     CTF Track 1 - (Et,eta,phi)=(" << tk1->pt() << "," << tk1->eta() << "," << tk1->phi() << ")  -  Fraction of common hits with GSF =" <<  (electrons->at(iElectron1)).shFracInnerHits() << std::endl;
84	–	if(verbosity_>4) std::cout << "     CTF Track 2 - (Et,eta,phi)=(" << tk2->pt() << "," << tk2->eta() << "," << tk2->phi() << ")  -  Fraction of common hits with GSF =" <<  (electrons->at(iElectron2)).shFracInnerHits()  << std::endl;
83
84
85		///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
#	Line 115 | Line 113 | bool ZeeVertexAnalyzer::getVertices(cons
113		}
114
115		if(verbosity_>4) std::cout  << std::endl << "   Private Primary Vertex reconstruction with all tracks - Number of vertices = " << zeeVertices.size() << std::endl;
116	<
116	>
117		for (unsigned int j=0; j<zeeVertices.size(); j++)
118		{
119		reco::Vertex* vertex = & (zeeVertices.at(j));
#	Line 169 | Line 167 | bool ZeeVertexAnalyzer::getVertices(cons
167		//    Private Primary Vertex reconstruction with all CTF tracks except the 2 two tracks matching selected electrons
168		///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
169
170	<	// Remove the two CTF tracks matching the two electron GSF tracks from the track collection
170	>	// Remove the two CTF tracks matching the two electron GSF tracks from the transientTracks collection
171	>	// Keep the two CTF electron tracks in a separate vector (eeTransientTracks)
172		unsigned int itk1 = 999999;
173		unsigned int itk2 = 999999;
174	+
175		for (unsigned int j=0; j<transientTracks.size(); j++)
176		{
177		const reco::Track* tk = & (transientTracks.at(j).track());
#	Line 192 | Line 192 | bool ZeeVertexAnalyzer::getVertices(cons
192		std::cout << "   ***** NO CTF TRACKS MATCHING THE ELECTRON *****" << std::endl;
193		return false;
194		}
195	+
196	+	if(verbosity_>4)
197	+	{
198	+	std::cout << "     CTF Track 1 - (Et,eta,phi)=(" << transientTracks.at(itk1).track().pt() << "," << transientTracks.at(itk1).track().eta() << "," << transientTracks.at(itk1).track().phi() << "  -  Fraction of common hits with GSF =" <<  (electrons->at(iElectron1)).shFracInnerHits();
199	+	if (!transientTracks.at(itk1).stateAtBeamLine().isValid()) std::cout << " IPSig=-1";
200	+	else std::cout << " IPSig=" << transientTracks.at(itk1).stateAtBeamLine().transverseImpactParameter().significance();
201	+	std::cout << " normChi2=" << transientTracks.at(itk1).normalizedChi2() << "  nSiHits=" <<  transientTracks.at(itk1).hitPattern().numberOfValidHits() << "  nPxlHits=" <<  transientTracks.at(itk1).hitPattern().numberOfValidPixelHits() << std::endl;
202	+
203	+	std::cout << "     CTF Track 2 - (Et,eta,phi)=(" << transientTracks.at(itk2).track().pt() << "," << transientTracks.at(itk2).track().eta() << "," << transientTracks.at(itk2).track().phi() << "  -  Fraction of common hits with GSF =" <<  (electrons->at(iElectron2)).shFracInnerHits();
204	+	if (!transientTracks.at(itk2).stateAtBeamLine().isValid()) std::cout << " IPSig=-1";
205	+	else std::cout << " IPSig=" << transientTracks.at(itk2).stateAtBeamLine().transverseImpactParameter().significance();
206	+	std::cout << " normChi2=" << transientTracks.at(itk2).normalizedChi2() << "  nSiHits=" <<  transientTracks.at(itk2).hitPattern().numberOfValidHits() << "  nPxlHits=" << transientTracks.at(itk2).hitPattern().numberOfValidPixelHits()<< std::endl;
207	+	}
208	+
209
210	+	rootBardak->setEle1(
211	+	TRootTrack(
212	+	transientTracks.at(itk1).track().px()
213	+	,transientTracks.at(itk1).track().py()
214	+	,transientTracks.at(itk1).track().pz()
215	+	,transientTracks.at(itk1).track().p()
216	+	,transientTracks.at(itk1).track().vx()
217	+	,transientTracks.at(itk1).track().vy()
218	+	,transientTracks.at(itk1).track().vz()
219	+	,0
220	+	,transientTracks.at(itk1).track().charge()
221	+	,transientTracks.at(itk1).hitPattern().pixelLayersWithMeasurement()
222	+	,transientTracks.at(itk1).hitPattern().stripLayersWithMeasurement()
223	+	,transientTracks.at(itk1).track().chi2()
224	+	,transientTracks.at(itk1).track().d0()
225	+	,transientTracks.at(itk1).track().d0Error()
226	+	,transientTracks.at(itk1).track().dz()
227	+	,transientTracks.at(itk1).track().dzError()
228	+	)
229	+	,transientTracks.at(itk1).hitPattern().numberOfValidHits()
230	+	,transientTracks.at(itk1).hitPattern().numberOfValidPixelHits()
231	+	,transientTracks.at(itk1).stateAtBeamLine().transverseImpactParameter().significance()
232	+	,transientTracks.at(itk1).normalizedChi2()
233	+	);
234	+
235	+	rootBardak->setEle2(
236	+	TRootTrack(
237	+	transientTracks.at(itk2).track().px()
238	+	,transientTracks.at(itk2).track().py()
239	+	,transientTracks.at(itk2).track().pz()
240	+	,transientTracks.at(itk2).track().p()
241	+	,transientTracks.at(itk2).track().vx()
242	+	,transientTracks.at(itk2).track().vy()
243	+	,transientTracks.at(itk2).track().vz()
244	+	,0
245	+	,transientTracks.at(itk2).track().charge()
246	+	,transientTracks.at(itk2).hitPattern().pixelLayersWithMeasurement()
247	+	,transientTracks.at(itk2).hitPattern().stripLayersWithMeasurement()
248	+	,transientTracks.at(itk2).track().chi2()
249	+	,transientTracks.at(itk2).track().d0()
250	+	,transientTracks.at(itk2).track().d0Error()
251	+	,transientTracks.at(itk2).track().dz()
252	+	,transientTracks.at(itk2).track().dzError()
253	+	)
254	+	,transientTracks.at(itk2).hitPattern().numberOfValidHits()
255	+	,transientTracks.at(itk2).hitPattern().numberOfValidPixelHits()
256	+	,transientTracks.at(itk2).stateAtBeamLine().transverseImpactParameter().significance()
257	+	,transientTracks.at(itk2).normalizedChi2()
258	+	);
259	+
260	+	vector<reco::TransientTrack> eeTransientTracks;
261	+	eeTransientTracks.push_back(transientTracks.at(itk1));
262	+	eeTransientTracks.push_back(transientTracks.at(itk2));
263		if (itk1<itk2) itk2--;
264		transientTracks.erase (transientTracks.begin()+itk1);
265		transientTracks.erase (transientTracks.begin()+itk2);
#	Line 225 | Line 292 | bool ZeeVertexAnalyzer::getVertices(cons
292		// Put your vertex selection here....
293		if (! vertex->isValid() ) continue;
294		if ( vertex->isFake() ) continue;
295	<
295	>
296		Int_t ntracks = 0;
297		Float_t higherPt = 0.;
298		Float_t scalarSumPt = 0.;
299		Float_t vectorSumPt = 0.;
300		math::XYZVector vectorSum(0.,0.,0.);
301	<
301	>
302		for( std::vector< reco::TrackBaseRef >::const_iterator it = vertex->tracks_begin(); it != vertex->tracks_end(); it++)
303		{
304		scalarSumPt += (**it).pt();
#	Line 240 | Line 307 | bool ZeeVertexAnalyzer::getVertices(cons
307		ntracks++;
308		}
309		vectorSumPt = sqrt(vectorSum.Perp2());
310	<
310	>
311		// No refitted tracks embeded in reco::Vertex....
312		//cout  << "vertex->refittedTracks().size()=" << vertex->refittedTracks().size() << endl;
313	<
313	>
314		TRootVertex localVertex(
315		vertex->x()
316		,vertex->y()
#	Line 252 | Line 319 | bool ZeeVertexAnalyzer::getVertices(cons
319		,vertex->yError()
320		,vertex->zError()
321		);
322	<
322	>
323		localVertex.setAlgoName("NoEE");
324		localVertex.setChi2( vertex->chi2() );
325		localVertex.setNdof( vertex->ndof() );
#	Line 260 | Line 327 | bool ZeeVertexAnalyzer::getVertices(cons
327		localVertex.setHigherTrackPt( higherPt );
328		localVertex.setScalarSumPt( scalarSumPt );
329		localVertex.setVectorSumPt( vectorSumPt );
330	+
331	+	new( (*rootVertices)[iRootVertex] ) TRootVertex(localVertex);
332	+	if(verbosity_>2) cout << "   ["<< setw(3) << iRootVertex << "] " << localVertex << endl;
333	+	iRootVertex++;
334	+	}
335
336	+
337	+	///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
338	+	//    Private Primary Vertex reconstruction with only the two CTF tracks matched to the selected electrons
339	+	///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
340	+
341	+	vector<TransientVertex> transientOnlyEEVertices = theVertexProducer.vertices(eeTransientTracks, vertexBeamSpot);
342	+	reco::VertexCollection onlyEEVertices;
343	+
344	+	for (vector<TransientVertex>::const_iterator vertex = transientOnlyEEVertices.begin(); vertex != transientOnlyEEVertices.end(); vertex++)
345	+	{
346	+	reco::Vertex v = *vertex;
347	+	onlyEEVertices.push_back(v);
348	+	}
349	+
350	+	if(verbosity_>4) std::cout  << std::endl << "   Private Primary Vertex reconstruction using only the two electron CTF tracks - Number of vertices = " << onlyEEVertices.size() << std::endl;
351	+
352	+	for (unsigned int j=0; j<onlyEEVertices.size(); j++)
353	+	{
354	+	reco::Vertex* vertex = & (onlyEEVertices.at(j));
355	+
356	+	// Put your vertex selection here....
357	+	if (! vertex->isValid() ) continue;
358	+	if ( vertex->isFake() ) continue;
359	+
360	+	Int_t ntracks = 0;
361	+	Float_t higherPt = 0.;
362	+	Float_t scalarSumPt = 0.;
363	+	Float_t vectorSumPt = 0.;
364	+	math::XYZVector vectorSum(0.,0.,0.);
365	+
366	+	for( std::vector< reco::TrackBaseRef >::const_iterator it = vertex->tracks_begin(); it != vertex->tracks_end(); it++)
367	+	{
368	+	scalarSumPt += (**it).pt();
369	+	vectorSum += (**it).momentum();
370	+	if( (**it).pt()>higherPt ) higherPt=(**it).pt();
371	+	ntracks++;
372	+	}
373	+	vectorSumPt = sqrt(vectorSum.Perp2());
374	+
375	+	// No refitted tracks embeded in reco::Vertex....
376	+	//cout  << "vertex->refittedTracks().size()=" << vertex->refittedTracks().size() << endl;
377	+
378	+	TRootVertex localVertex(
379	+	vertex->x()
380	+	,vertex->y()
381	+	,vertex->z()
382	+	,vertex->xError()
383	+	,vertex->yError()
384	+	,vertex->zError()
385	+	);
386	+
387	+	localVertex.setAlgoName("OnlyEE");
388	+	localVertex.setChi2( vertex->chi2() );
389	+	localVertex.setNdof( vertex->ndof() );
390	+	localVertex.setNtracks( ntracks );
391	+	localVertex.setHigherTrackPt( higherPt );
392	+	localVertex.setScalarSumPt( scalarSumPt );
393	+	localVertex.setVectorSumPt( vectorSumPt );
394	+
395		new( (*rootVertices)[iRootVertex] ) TRootVertex(localVertex);
396		if(verbosity_>2) cout << "   ["<< setw(3) << iRootVertex << "] " << localVertex << endl;
397		iRootVertex++;
398		}
399
400	<	return true;
400	>	return true;
401		}
271	–

Diff Legend

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