ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/Morgan/interface/TRootPhoton.h

Comparing UserCode/Morgan/interface/TRootPhoton.h (file contents):
Revision 1.7 by lethuill, Mon Dec 15 19:08:52 2008 UTC vs.
Revision 1.10 by lethuill, Mon Feb 9 16:11:43 2009 UTC

#	Line 13 | Line 13 | using namespace std;
13
14		class TRootPhoton : public TRootParticle
15		{
16	<
16	>
17		public:
18	<
18	>
19		TRootPhoton() :
20	<	TRootParticle()
21	<	,caloPosition_()
22	<	,hoe_(0.)
23	<	,hasPixelSeed_(0.)
24	<	,eMax_(0.)
25	<	,e3x3_(0.)
26	<	,e5x5_(0.)
27	<	,scRawEnergy_(0.)
28	<	,pi0nn_(-1.)
29	<	,isoEcalIsland_(-1.)
30	<	,isoEcalDoubleCone_(-1.)
31	<	,isoHcal_(-1.)
32	<	,isoTracks_(-1.)
33	<	,isoNTracks_(-1)
34	<	,scIndex_()
35	<	,scRef_()
36	<	,convNTracks_(0)
37	<	,convEoverP_(-1.)
38	<	,convMass_(-1.)
39	<	,convCotanTheta_(-999.)
40	<	,convLikely_(-1.)
41	<	,convVertex_()
42	<	,convEcalImpactPosition1_()
43	<	,convEcalImpactPosition2_()
44	<	,convIndexTrack1_(-1)
45	<	,convIndexTrack2_(-1)
46	<	,convTrack1_()
47	<	,convTrack2_()
48	<	,photonMCIndex_(-1)
49	<	,momentumMCPhoton_()
50	<	,vertexMCPhoton_()
51	<	,pdgIdMCPhoton_(0)
52	<	{;}
53	<
20	>	TRootParticle()
21	>	,caloPosition_()
22	>	,hoe_(0.)
23	>	,hasPixelSeed_(0.)
24	>	,eMax_(0.)
25	>	,e3x3_(0.)
26	>	,e5x5_(0.)
27	>	,scRawEnergy_(0.)
28	>	,pi0nn_(-1.)
29	>	,isoEcalIsland_(-1.)
30	>	,isoEcalDoubleCone_(-1.)
31	>	,isoHcal_(-1.)
32	>	,isoTracks_(-1.)
33	>	,isoNTracks_(-1)
34	>	,scIndex_()
35	>	,scRef_()
36	>	,convNTracks_(0)
37	>	,convEoverP_(-1.)
38	>	,convMass_(-1.)
39	>	,convCotanTheta_(-999.)
40	>	,convLikely_(-1.)
41	>	,convVertex_()
42	>	,convEcalImpactPosition1_()
43	>	,convEcalImpactPosition2_()
44	>	,convIndexTrack1_(-1)
45	>	,convIndexTrack2_(-1)
46	>	,convTrack1_()
47	>	,convTrack2_()
48	>	{;}
49	>
50		TRootPhoton(const TRootPhoton& photon) :
51	<	TRootParticle(photon)
52	<	,caloPosition_(photon.caloPosition_)
53	<	,hoe_(photon.hoe_)
54	<	,hasPixelSeed_(photon.hasPixelSeed_)
55	<	,eMax_(photon.eMax_)
56	<	,e3x3_(photon.e3x3_)
57	<	,e5x5_(photon.e5x5_)
58	<	,scRawEnergy_(photon.scRawEnergy_)
59	<	,pi0nn_(photon.pi0nn_)
60	<	,isoEcalIsland_(photon.isoEcalIsland_)
61	<	,isoEcalDoubleCone_(photon.isoEcalDoubleCone_)
62	<	,isoHcal_(photon.isoHcal_)
63	<	,isoTracks_(photon.isoTracks_)
64	<	,isoNTracks_(photon.isoNTracks_)
65	<	,scIndex_(photon.scIndex_)
66	<	,scRef_(photon.scRef_)
67	<	,convNTracks_(photon.convNTracks_)
68	<	,convEoverP_(photon.convEoverP_)
69	<	,convMass_(photon.convMass_)
70	<	,convCotanTheta_(photon.convCotanTheta_)
71	<	,convLikely_(photon.convLikely_)
72	<	,convVertex_(photon.convVertex_)
73	<	,convEcalImpactPosition1_(photon.convEcalImpactPosition1_)
74	<	,convEcalImpactPosition2_(photon.convEcalImpactPosition2_)
75	<	,convIndexTrack1_(photon.convIndexTrack1_)
76	<	,convIndexTrack2_(photon.convIndexTrack2_)
77	<	,convTrack1_(photon.convTrack1_)
78	<	,convTrack2_(photon.convTrack2_)
79	<	,photonMCIndex_(photon.photonMCIndex_)
80	<	,momentumMCPhoton_(photon.momentumMCPhoton_)
85	<	,vertexMCPhoton_(photon.vertexMCPhoton_)
86	<	,pdgIdMCPhoton_(photon.pdgIdMCPhoton_)
87	<	{;}
88	<
51	>	TRootParticle(photon)
52	>	,caloPosition_(photon.caloPosition_)
53	>	,hoe_(photon.hoe_)
54	>	,hasPixelSeed_(photon.hasPixelSeed_)
55	>	,eMax_(photon.eMax_)
56	>	,e3x3_(photon.e3x3_)
57	>	,e5x5_(photon.e5x5_)
58	>	,scRawEnergy_(photon.scRawEnergy_)
59	>	,pi0nn_(photon.pi0nn_)
60	>	,isoEcalIsland_(photon.isoEcalIsland_)
61	>	,isoEcalDoubleCone_(photon.isoEcalDoubleCone_)
62	>	,isoHcal_(photon.isoHcal_)
63	>	,isoTracks_(photon.isoTracks_)
64	>	,isoNTracks_(photon.isoNTracks_)
65	>	,scIndex_(photon.scIndex_)
66	>	,scRef_(photon.scRef_)
67	>	,convNTracks_(photon.convNTracks_)
68	>	,convEoverP_(photon.convEoverP_)
69	>	,convMass_(photon.convMass_)
70	>	,convCotanTheta_(photon.convCotanTheta_)
71	>	,convLikely_(photon.convLikely_)
72	>	,convVertex_(photon.convVertex_)
73	>	,convEcalImpactPosition1_(photon.convEcalImpactPosition1_)
74	>	,convEcalImpactPosition2_(photon.convEcalImpactPosition2_)
75	>	,convIndexTrack1_(photon.convIndexTrack1_)
76	>	,convIndexTrack2_(photon.convIndexTrack2_)
77	>	,convTrack1_(photon.convTrack1_)
78	>	,convTrack2_(photon.convTrack2_)
79	>	{;}
80	>
81		TRootPhoton(Double_t px, Double_t py, Double_t pz, Double_t e) :
82	<	TRootParticle(px,py,pz,e)
83	<	,caloPosition_()
84	<	,hoe_(0.)
85	<	,hasPixelSeed_(0.)
86	<	,eMax_(0.)
87	<	,e3x3_(0.)
88	<	,e5x5_(0.)
89	<	,scRawEnergy_(0.)
90	<	,pi0nn_(-1.)
91	<	,isoEcalIsland_(-1.)
92	<	,isoEcalDoubleCone_(-1.)
93	<	,isoHcal_(-1.)
94	<	,isoTracks_(-1.)
95	<	,isoNTracks_(-1)
96	<	,scIndex_()
97	<	,scRef_()
98	<	,convNTracks_(0)
99	<	,convEoverP_(-1.)
100	<	,convMass_(-1.)
101	<	,convCotanTheta_(-999.)
102	<	,convLikely_(-1.)
103	<	,convVertex_()
104	<	,convEcalImpactPosition1_()
105	<	,convEcalImpactPosition2_()
106	<	,convIndexTrack1_(-1)
107	<	,convIndexTrack2_(-1)
108	<	,convTrack1_()
109	<	,convTrack2_()
110	<	,photonMCIndex_(-1)
119	<	,momentumMCPhoton_()
120	<	,vertexMCPhoton_()
121	<	,pdgIdMCPhoton_(0)
122	<	{;}
82	>	TRootParticle(px,py,pz,e)
83	>	,caloPosition_()
84	>	,hoe_(0.)
85	>	,hasPixelSeed_(0.)
86	>	,eMax_(0.)
87	>	,e3x3_(0.)
88	>	,e5x5_(0.)
89	>	,scRawEnergy_(0.)
90	>	,pi0nn_(-1.)
91	>	,isoEcalIsland_(-1.)
92	>	,isoEcalDoubleCone_(-1.)
93	>	,isoHcal_(-1.)
94	>	,isoTracks_(-1.)
95	>	,isoNTracks_(-1)
96	>	,scIndex_()
97	>	,scRef_()
98	>	,convNTracks_(0)
99	>	,convEoverP_(-1.)
100	>	,convMass_(-1.)
101	>	,convCotanTheta_(-999.)
102	>	,convLikely_(-1.)
103	>	,convVertex_()
104	>	,convEcalImpactPosition1_()
105	>	,convEcalImpactPosition2_()
106	>	,convIndexTrack1_(-1)
107	>	,convIndexTrack2_(-1)
108	>	,convTrack1_()
109	>	,convTrack2_()
110	>	{;}
111
112		TRootPhoton(Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vtx_x, Double_t vtx_y, Double_t vtx_z) :
113	<	TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z)
114	<	,caloPosition_()
115	<	,hoe_(0.)
116	<	,hasPixelSeed_(0.)
117	<	,eMax_(0.)
118	<	,e3x3_(0.)
119	<	,e5x5_(0.)
120	<	,scRawEnergy_(0.)
121	<	,pi0nn_(-1.)
122	<	,isoEcalIsland_(-1.)
123	<	,isoEcalDoubleCone_(-1.)
124	<	,isoHcal_(-1.)
125	<	,isoTracks_(-1.)
126	<	,isoNTracks_(-1)
127	<	,scIndex_()
128	<	,scRef_()
129	<	,convNTracks_(0)
130	<	,convEoverP_(-1.)
131	<	,convMass_(-1.)
132	<	,convCotanTheta_(-999.)
133	<	,convLikely_(-1.)
134	<	,convVertex_()
135	<	,convEcalImpactPosition1_()
136	<	,convEcalImpactPosition2_()
137	<	,convIndexTrack1_(-1)
138	<	,convIndexTrack2_(-1)
139	<	,photonMCIndex_(-1)
140	<	,momentumMCPhoton_()
153	<	,vertexMCPhoton_()
154	<	,pdgIdMCPhoton_(0)
155	<	{;}
156	<
113	>	TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z)
114	>	,caloPosition_()
115	>	,hoe_(0.)
116	>	,hasPixelSeed_(0.)
117	>	,eMax_(0.)
118	>	,e3x3_(0.)
119	>	,e5x5_(0.)
120	>	,scRawEnergy_(0.)
121	>	,pi0nn_(-1.)
122	>	,isoEcalIsland_(-1.)
123	>	,isoEcalDoubleCone_(-1.)
124	>	,isoHcal_(-1.)
125	>	,isoTracks_(-1.)
126	>	,isoNTracks_(-1)
127	>	,scIndex_()
128	>	,scRef_()
129	>	,convNTracks_(0)
130	>	,convEoverP_(-1.)
131	>	,convMass_(-1.)
132	>	,convCotanTheta_(-999.)
133	>	,convLikely_(-1.)
134	>	,convVertex_()
135	>	,convEcalImpactPosition1_()
136	>	,convEcalImpactPosition2_()
137	>	,convIndexTrack1_(-1)
138	>	,convIndexTrack2_(-1)
139	>	{;}
140	>
141		TRootPhoton(Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vtx_x, Double_t vtx_y, Double_t vtx_z, Int_t type, Float_t charge) :
142	<	TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z,type,charge)
143	<	,caloPosition_()
144	<	,hoe_(0.)
145	<	,hasPixelSeed_(0.)
146	<	,eMax_(0.)
147	<	,e3x3_(0.)
148	<	,e5x5_(0.)
149	<	,scRawEnergy_(0.)
150	<	,pi0nn_(-1.)
151	<	,isoEcalIsland_(-1.)
152	<	,isoEcalDoubleCone_(-1.)
153	<	,isoHcal_(-1.)
154	<	,isoTracks_(-1.)
155	<	,isoNTracks_(-1)
156	<	,scIndex_()
157	<	,scRef_()
158	<	,convNTracks_(0)
159	<	,convEoverP_(-1.)
160	<	,convMass_(-1.)
161	<	,convCotanTheta_(-999.)
162	<	,convLikely_(-1.)
163	<	,convVertex_()
164	<	,convEcalImpactPosition1_()
165	<	,convEcalImpactPosition2_()
166	<	,convIndexTrack1_(-1)
167	<	,convIndexTrack2_(-1)
168	<	,convTrack1_()
169	<	,convTrack2_()
170	<	,photonMCIndex_(-1)
171	<	,momentumMCPhoton_()
188	<	,vertexMCPhoton_()
189	<	,pdgIdMCPhoton_(0)
190	<	{;}
191	<
142	>	TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z,type,charge)
143	>	,caloPosition_()
144	>	,hoe_(0.)
145	>	,hasPixelSeed_(0.)
146	>	,eMax_(0.)
147	>	,e3x3_(0.)
148	>	,e5x5_(0.)
149	>	,scRawEnergy_(0.)
150	>	,pi0nn_(-1.)
151	>	,isoEcalIsland_(-1.)
152	>	,isoEcalDoubleCone_(-1.)
153	>	,isoHcal_(-1.)
154	>	,isoTracks_(-1.)
155	>	,isoNTracks_(-1)
156	>	,scIndex_()
157	>	,scRef_()
158	>	,convNTracks_(0)
159	>	,convEoverP_(-1.)
160	>	,convMass_(-1.)
161	>	,convCotanTheta_(-999.)
162	>	,convLikely_(-1.)
163	>	,convVertex_()
164	>	,convEcalImpactPosition1_()
165	>	,convEcalImpactPosition2_()
166	>	,convIndexTrack1_(-1)
167	>	,convIndexTrack2_(-1)
168	>	,convTrack1_()
169	>	,convTrack2_()
170	>	{;}
171	>
172		TRootPhoton(const TLorentzVector &momentum) :
173	<	TRootParticle(momentum)
174	<	,caloPosition_()
175	<	,hoe_(0.)
176	<	,hasPixelSeed_(0.)
177	<	,eMax_(0.)
178	<	,e3x3_(0.)
179	<	,e5x5_(0.)
180	<	,scRawEnergy_(0.)
181	<	,pi0nn_(-1.)
182	<	,isoEcalIsland_(-1.)
183	<	,isoEcalDoubleCone_(-1.)
184	<	,isoHcal_(-1.)
185	<	,isoTracks_(-1.)
186	<	,isoNTracks_(-1)
187	<	,scIndex_()
188	<	,scRef_()
189	<	,convNTracks_(0)
190	<	,convEoverP_(-1.)
191	<	,convMass_(-1.)
192	<	,convCotanTheta_(-999.)
193	<	,convLikely_(-1.)
194	<	,convVertex_()
195	<	,convEcalImpactPosition1_()
196	<	,convEcalImpactPosition2_()
197	<	,convIndexTrack1_(-1)
198	<	,convIndexTrack2_(-1)
199	<	,convTrack1_()
200	<	,convTrack2_()
201	<	,photonMCIndex_(-1)
202	<	,momentumMCPhoton_()
223	<	,vertexMCPhoton_()
224	<	,pdgIdMCPhoton_(0)
225	<	{;}
226	<
173	>	TRootParticle(momentum)
174	>	,caloPosition_()
175	>	,hoe_(0.)
176	>	,hasPixelSeed_(0.)
177	>	,eMax_(0.)
178	>	,e3x3_(0.)
179	>	,e5x5_(0.)
180	>	,scRawEnergy_(0.)
181	>	,pi0nn_(-1.)
182	>	,isoEcalIsland_(-1.)
183	>	,isoEcalDoubleCone_(-1.)
184	>	,isoHcal_(-1.)
185	>	,isoTracks_(-1.)
186	>	,isoNTracks_(-1)
187	>	,scIndex_()
188	>	,scRef_()
189	>	,convNTracks_(0)
190	>	,convEoverP_(-1.)
191	>	,convMass_(-1.)
192	>	,convCotanTheta_(-999.)
193	>	,convLikely_(-1.)
194	>	,convVertex_()
195	>	,convEcalImpactPosition1_()
196	>	,convEcalImpactPosition2_()
197	>	,convIndexTrack1_(-1)
198	>	,convIndexTrack2_(-1)
199	>	,convTrack1_()
200	>	,convTrack2_()
201	>	{;}
202	>
203		TRootPhoton(const TLorentzVector &momentum, const TVector3 &vertex, Int_t type, Float_t charge) :
204	<	TRootParticle(momentum, vertex, type, charge)
205	<	,caloPosition_()
206	<	,hoe_(0.)
207	<	,hasPixelSeed_(0.)
208	<	,eMax_(0.)
209	<	,e3x3_(0.)
210	<	,e5x5_(0.)
211	<	,scRawEnergy_(0.)
212	<	,pi0nn_(-1.)
213	<	,isoEcalIsland_(-1.)
214	<	,isoEcalDoubleCone_(-1.)
215	<	,isoHcal_(-1.)
216	<	,isoTracks_(-1.)
217	<	,isoNTracks_(-1)
218	<	,scIndex_()
219	<	,scRef_()
220	<	,convNTracks_(0)
221	<	,convEoverP_(-1.)
222	<	,convMass_(-1.)
223	<	,convCotanTheta_(-999.)
224	<	,convLikely_(-1.)
225	<	,convVertex_()
226	<	,convEcalImpactPosition1_()
227	<	,convEcalImpactPosition2_()
228	<	,convIndexTrack1_(-1)
229	<	,convIndexTrack2_(-1)
230	<	,convTrack1_()
231	<	,convTrack2_()
232	<	,photonMCIndex_(-1)
233	<	,momentumMCPhoton_()
258	<	,vertexMCPhoton_()
259	<	,pdgIdMCPhoton_(0)
260	<	{;}
261	<
204	>	TRootParticle(momentum, vertex, type, charge)
205	>	,caloPosition_()
206	>	,hoe_(0.)
207	>	,hasPixelSeed_(0.)
208	>	,eMax_(0.)
209	>	,e3x3_(0.)
210	>	,e5x5_(0.)
211	>	,scRawEnergy_(0.)
212	>	,pi0nn_(-1.)
213	>	,isoEcalIsland_(-1.)
214	>	,isoEcalDoubleCone_(-1.)
215	>	,isoHcal_(-1.)
216	>	,isoTracks_(-1.)
217	>	,isoNTracks_(-1)
218	>	,scIndex_()
219	>	,scRef_()
220	>	,convNTracks_(0)
221	>	,convEoverP_(-1.)
222	>	,convMass_(-1.)
223	>	,convCotanTheta_(-999.)
224	>	,convLikely_(-1.)
225	>	,convVertex_()
226	>	,convEcalImpactPosition1_()
227	>	,convEcalImpactPosition2_()
228	>	,convIndexTrack1_(-1)
229	>	,convIndexTrack2_(-1)
230	>	,convTrack1_()
231	>	,convTrack2_()
232	>	{;}
233	>
234		~TRootPhoton() {;}
235
236
237	<	TVector3 caloPosition() const  { return caloPosition_;}
237	>	TVector3 caloPosition() const { return caloPosition_; }
238		Float_t hoe() const { return hoe_; }
239	<	Bool_t hasPixelSeed() const  { return hasPixelSeed_; }
240	<	Float_t eMax() const   { return eMax_; }
241	<	Float_t e3x3() const   { return e3x3_; }
242	<	Float_t e5x5() const   { return e5x5_; }
243	<	Float_t scRawEnergy() const   { return scRawEnergy_; }
244	<	Float_t pi0nn() const  { return pi0nn_; }
245	<	Float_t isoEcalIsland() const  { return isoEcalIsland_; }
246	<	Float_t isoEcalDoubleCone() const  { return isoEcalDoubleCone_; }
247	<	Float_t isoHcal() const  { return isoHcal_; }
248	<	Float_t isoTracks() const  { return isoTracks_; }
249	<	Int_t isoNTracks() const  { return isoNTracks_; }
250	<
239	>	Bool_t hasPixelSeed() const { return hasPixelSeed_; }
240	>	Float_t eMax() const  { return eMax_; }
241	>	Float_t e3x3() const  { return e3x3_; }
242	>	Float_t e5x5() const  { return e5x5_; }
243	>	Float_t scRawEnergy() const { return scRawEnergy_; }
244	>	Float_t pi0nn() const { return pi0nn_; }
245	>	Float_t isoEcalIsland() const { return isoEcalIsland_; }
246	>	Float_t isoEcalDoubleCone() const { return isoEcalDoubleCone_; }
247	>	Float_t isoHcal() const { return isoHcal_; }
248	>	Float_t isoTracks() const { return isoTracks_; }
249	>	Int_t isoNTracks() const { return isoNTracks_; }
250	>
251		// Index in TCloneArray of the SC seed of type "type"
252		Int_t scIndex(Int_t type)
253		{
254		map<Int_t,Int_t>::iterator it=scIndex_.find(type);
255		return ( it ==scIndex_.end() ? -1 : (*it).second );
256		}
257	<	map<Int_t,Int_t> scIndexMap() const { return scIndex_;}
257	>	map<Int_t,Int_t> scIndexMap() const { return scIndex_; }
258
259		// return first SC seed
260		TObject* superCluster()
#	Line 297 | Line 269 | public:
269		map<Int_t,TRef>::iterator it=scRef_.find(type);
270		return ( it ==scRef_.end() ? 0 : (*it).second.GetObject() );
271		}
272	<	map<Int_t,TRef> scRefMap() const { return scRef_;}
272	>	map<Int_t,TRef> scRefMap() const { return scRef_; }
273
274		// Converted photon infos
275	<	Int_t convNTracks() const  { return convNTracks_;}
276	<	Float_t convEoverP() const  { return convEoverP_;}
277	<	Float_t convMass() const  { return convMass_;}
278	<	Float_t convCotanTheta() const  { return convCotanTheta_;}
279	<	Float_t convLikely() const  { return convLikely_;}
280	<	TVector3 convVertex() const  { return convVertex_;}
281	<	TVector3 convEcalImpactPosition1() const  { return convEcalImpactPosition1_;}
282	<	TVector3 convEcalImpactPosition2() const  { return convEcalImpactPosition2_;}
283	<	TObject* convTrack1() const  { return convTrack1_.GetObject();}
284	<	TObject* convTrack2() const  { return convTrack2_.GetObject();}
275	>	Int_t convNTracks() const { return convNTracks_; }
276	>	Float_t convEoverP() const { return convEoverP_; }
277	>	Float_t convMass() const { return convMass_; }
278	>	Float_t convCotanTheta() const { return convCotanTheta_; }
279	>	Float_t convLikely() const { return convLikely_; }
280	>	TVector3 convVertex() const { return convVertex_; }
281	>	TVector3 convEcalImpactPosition1() const { return convEcalImpactPosition1_; }
282	>	TVector3 convEcalImpactPosition2() const { return convEcalImpactPosition2_; }
283	>	TObject* convTrack1() const { return convTrack1_.GetObject(); }
284	>	TObject* convTrack2() const { return convTrack2_.GetObject(); }
285		// FIXME - remove convIndexTrack1,2() (kept for backward compatibility)
286	<	Int_t convIndexTrack1() const  { return convIndexTrack1_;}
287	<	Int_t convIndexTrack2() const  { return convIndexTrack2_;}
288	<
286	>	Int_t convIndexTrack1() const { return convIndexTrack1_; }
287	>	Int_t convIndexTrack2() const { return convIndexTrack2_; }
288	>
289		// Photon ID
290	<	Bool_t isLooseEM() const{return isLooseEM_;}
291	<	Bool_t isLoosePhoton() const{return isLoosePhoton_;}
292	<	Bool_t isTightPhoton() const{return isTightPhoton_;}
293	<	Bool_t isEBPho() const{return isEBPho_;}        // if photon is in ECAL barrel
294	<	Bool_t isEEPho() const{return isEEPho_;}        // if photon is in ECAL endcap
295	<	Bool_t isEBGap() const{return isEBGap_;}        // if photon is in EB, and inside the boundaries in super crystals/modules
296	<	Bool_t isEEGap() const{return isEEGap_;}        // if photon is in EE, and inside the boundaries in supercrystal/D
297	<	Bool_t isEBEEGap() const{return isEBEEGap_;}    // if photon is in boundary between EB and EE
298	<	Bool_t isAlsoElectron() const{return isAlsoElectron_;}  // if this is also a GsfElectron
290	>	Bool_t isLooseEM() const { return isLooseEM_; }
291	>	Bool_t isLoosePhoton() const { return isLoosePhoton_; }
292	>	Bool_t isTightPhoton() const { return isTightPhoton_; }
293	>	Bool_t isEBPho() const { return isEBPho_; }     // if photon is in ECAL barrel
294	>	Bool_t isEEPho() const { return isEEPho_; }     // if photon is in ECAL endcap
295	>	Bool_t isEBGap() const { return isEBGap_; }     // if photon is in EB, and inside the boundaries in super crystals/modules
296	>	Bool_t isEEGap() const { return isEEGap_; }     // if photon is in EE, and inside the boundaries in supercrystal/D
297	>	Bool_t isEBEEGap() const { return isEBEEGap_; } // if photon is in boundary between EB and EE
298	>	Bool_t isAlsoElectron() const { return isAlsoElectron_; }       // if this is also a GsfElectron
299
300		// ratio of Emax/E(3x3)
301		Float_t r19() const { return ( e3x3_ == 0. ? -1. : eMax_/e3x3_ ); }
#	Line 332 | Line 304 | public:
304		// ratio of E(3x3)/ESC
305		Float_t r9() const { return ( scRawEnergy_ == 0. ? -1. : e3x3_/scRawEnergy_ ); }
306
307	<	Int_t photonMCIndex() const { return photonMCIndex_; }
308	<	TLorentzVector momentumMCPhoton() const  { return momentumMCPhoton_; }
309	<	TVector3 vertexMCPhoton() const { return vertexMCPhoton_; }
338	<	Int_t pdgIdMCPhoton() const { return pdgIdMCPhoton_; }
307	>	//TObject* genPhoton() const { return genPhoton_.GetObject() ;}
308	>
309	>	virtual TString typeName() const { return "TRootPhoton"; }
310
311
312		void setCaloPosition(TVector3 caloPosition) { caloPosition_ = caloPosition; }
#	Line 355 | Line 326 | public:
326		void setIsoNTracks(Int_t isoNTracks) { isoNTracks_ = isoNTracks; }
327		void setSCIndex(Int_t type, Int_t index) { scIndex_[type]=index; }
328		void setSCRef(Int_t type, TObject* superCluster) { scRef_[type]=superCluster; }
329	<
329	>
330		void setIsolation(Float_t isoEcalIsland, Float_t isoEcalDoubleCone, Float_t isoHcal, Float_t isoTracks, Int_t isoNTracks)
331		{
332		isoEcalIsland_ = isoEcalIsland;
#	Line 365 | Line 336 | public:
336		isoNTracks_ = isoNTracks;
337		}
338
339	<
339	>
340		// Attach new vertex position to the photon - Et, Eta, Phi corrected accordingly
341		void setVertex(TVector3 vertex) {
342		TVector3 direction = caloPosition() - vertex;
#	Line 374 | Line 345 | public:
345		this->SetXYZT(momentum.x(), momentum.y(), momentum.z(), energy );
346		vertex_ = vertex;
347		}
348	<
348	>
349
350		void setConvNTracks(Int_t convNTracks) { convNTracks_=convNTracks; }
351		void setConvEoverP(Float_t convEoverP) { convEoverP_=convEoverP; }
#	Line 415 | Line 386 | public:
386		isEBEEGap_ = isEBEEGap;
387		isAlsoElectron_ = isAlsoElectron;
388		}
418	–
419	–	void setphotonMCIndex(Int_t photonMCIndex) { photonMCIndex_ = photonMCIndex; }
420	–	void setMomentumMCPhoton(Float_t px, Float_t py, Float_t pz , Float_t energy) { momentumMCPhoton_ =  TLorentzVector(px,py,pz,energy); }
421	–	void setVertexMCPhoton(TVector3 vertexMCPhoton) { vertexMCPhoton_ = vertexMCPhoton; }
422	–	void setPdgIdMCPhoton(Int_t pdgIdMCPhoton) { pdgIdMCPhoton_ = pdgIdMCPhoton; }
389
390
391		// TODO - Add photonID bits and conversion infos in output
392		friend std::ostream& operator<< (std::ostream& stream, const TRootPhoton& photon) {
393	<	stream << "TRootPhoton  (E,Et,eta,phi)=("<< photon.Mag() <<","<< photon.Et() <<","<< photon.Eta() <<","<< photon.Phi() << ")"
393	>	stream << "TRootPhoton  (E,Et,eta,phi)=("<< photon.Energy() <<","<< photon.Et() <<","<< photon.Eta() <<","<< photon.Phi() << ")"
394		<< " hasPxlSeed=" << photon. hasPixelSeed() <<"  E5x5=" << photon.e5x5() <<"  R19=" << photon.r19()
395		<<"  R9=" << photon.r9() <<"  hoe=" << photon.hoe()  <<"  Pi0=" << photon.pi0nn()
396		<< " vertex=("<< photon.vx() <<","<< photon.vy() <<","<< photon.vz() << ")"
#	Line 432 | Line 398 | public:
398		<< " iso(ecal1,ecal2,hcal,track,ntracks)=("<< photon.isoEcalIsland() << "," << photon.isoEcalDoubleCone() << "," << photon.isoHcal() << "," << photon.isoTracks() << "," << photon.isoNTracks() << ")";
399		return stream;
400		};
401	<
401	>
402		void Print()
403		{
404	<	cout << "TRootPhoton  (E,Et,eta,phi)=("<< this->Mag() <<","<< this->Et() <<","<< this->Eta() <<","<< this->Phi() << ")"
404	>	cout << "TRootPhoton  (E,Et,eta,phi)=("<< this->Energy() <<","<< this->Et() <<","<< this->Eta() <<","<< this->Phi() << ")"
405		<< " hasPxlSeed=" << this-> hasPixelSeed() <<"  E5x5=" << this->e5x5() <<"  R19=" << this->r19()
406		<<"  R9=" << this->r9() <<"  hoe=" << this->hoe()  <<"  Pi0=" << this->pi0nn()
407		<< " vertex=("<< this->vx() <<","<< this->vy() <<","<< this->vz() << ")"
408		<< " caloPosition=(" << this->caloPosition().X() << "," << this->caloPosition().Y() << "," << this->caloPosition().Z() << ")"
409		<< " iso(ecal1,ecal2,hcal,track,ntracks)=("<< this->isoEcalIsland() << "," << this->isoEcalDoubleCone() << "," << this->isoHcal() << "," << this->isoTracks() << "," << this->isoNTracks() << ")";
410		};
411	<
411	>
412		private:
413
414		TVector3 caloPosition_;         // SuperCluster centroid position
#	Line 460 | Line 426 | private:
426		Int_t isoNTracks_;                      // Tracks multiplicity in a cone around the photon
427
428		map<Int_t,Int_t> scIndex_;      // indexes of the seed SuperClusters (one per SC tytpe)
429	<	map<Int_t,TRef> scRef_;     // references to the seed SuperClusters (one per SC tytpe)
429	>	map<Int_t,TRef> scRef_;       // references to the seed SuperClusters (one per SC tytpe)
430
431		// Conversion Infos
432		Int_t convNTracks_;                     // Number of reconstructed tracks associated to the converted photon candidate
#	Line 476 | Line 442 | private:
442		Int_t convIndexTrack2_;         // Index of second conversion track (in tracks TCloneArray)
443		TRef convTrack1_;
444		TRef convTrack2_;
445	<
445	>
446		// Photon ID:
447		//              LooseEM:  ECalRecHit isolation < 20 GeV   AND   HCalRecHit isolation < 10 GeV
448		//              LoosePhoton:  LooseEM   AND   Hollow Cone track isolation < 30 GeV
#	Line 490 | Line 456 | private:
456		Bool_t isLooseEM_;
457		Bool_t isLoosePhoton_;
458		Bool_t isTightPhoton_;
459	<
459	>
460		Bool_t isEBPho_;
461		Bool_t isEEPho_;
462		Bool_t isEBGap_;
#	Line 498 | Line 464 | private:
464		Bool_t isEBEEGap_;
465		Bool_t isAlsoElectron_;
466
467	<	Int_t photonMCIndex_;
502	<	TLorentzVector momentumMCPhoton_;
503	<	TVector3 vertexMCPhoton_;
504	<	Int_t pdgIdMCPhoton_;
505	<
506	<	ClassDef (TRootPhoton,1);
467	>	ClassDef (TRootPhoton,2);
468		};
469
470		#endif

Diff Legend

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