1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 //_________________________________________________________________________
19 // A Particle modified by PHOS response and produced by AliPHOSvFast
20 // To become a general class of AliRoot ?
22 // This is also a base class for AliPHOSRecParticle produced by AliPHOSPIDv1
23 // The rec.particle type is to be defined by AliPHOSvFast or AliPHOSPIDv1
25 //*-- Author: Yves Schutz (SUBATECH)
27 // --- ROOT system ---
29 // --- Standard library ---
31 // --- AliRoot header files ---
33 #include "AliPHOSFastRecParticle.h"
35 #include "TPaveText.h"
37 ClassImp(AliPHOSFastRecParticle) ;
39 //____________________________________________________________________________
40 AliPHOSFastRecParticle::AliPHOSFastRecParticle() : TParticle()
46 //____________________________________________________________________________
47 AliPHOSFastRecParticle::AliPHOSFastRecParticle(const AliPHOSFastRecParticle & rp)
52 fPdgCode = rp.fPdgCode;
53 fStatusCode = rp.fStatusCode;
54 fMother[0] = rp.fMother[0];
55 fMother[1] = rp.fMother[1];
56 fDaughter[0] = rp.fDaughter[0];
57 fDaughter[1] = rp.fDaughter[1];
59 fCalcMass = rp.fCalcMass;
68 fPolarTheta = rp.fPolarTheta;
69 fPolarPhi = rp.fPolarPhi;
70 fParticlePDG = rp.fParticlePDG;
73 //____________________________________________________________________________
74 AliPHOSFastRecParticle::AliPHOSFastRecParticle(const TParticle & pp)
76 // ctor from a TParticle (crummy?!)
78 TParticle & pnoconst = (TParticle &)(pp) ;
79 AliPHOSFastRecParticle & p = (AliPHOSFastRecParticle &)(pnoconst) ;
81 fPdgCode = p.fPdgCode;
82 fStatusCode = p.fStatusCode;
83 fMother[0] = p.fMother[0];
84 fMother[1] = p.fMother[1];
85 fDaughter[0] = p.fDaughter[0];
86 fDaughter[1] = p.fDaughter[1];
88 fCalcMass = p.fCalcMass;
97 fPolarTheta = p.fPolarTheta;
98 fPolarPhi = p.fPolarPhi;
99 fParticlePDG = p.fParticlePDG;
103 //____________________________________________________________________________
104 Int_t AliPHOSFastRecParticle::DistancetoPrimitive(Int_t px, Int_t py)
106 // Compute distance from point px,py to a AliPHOSFastRecParticle considered as a Tmarker
107 // Compute the closest distance of approach from point px,py to this marker.
108 // The distance is computed in pixels units.
110 Double_t kRADDEG = 180. / TMath::Pi() ;
111 Coord_t x = Phi() * kRADDEG ;
112 Coord_t y = Theta() * kRADDEG ;
113 const Int_t kMaxDiff = 10;
114 Int_t pxm = gPad->XtoAbsPixel(x);
115 Int_t pym = gPad->YtoAbsPixel(y);
116 Int_t dist = (px-pxm)*(px-pxm) + (py-pym)*(py-pym);
118 if (dist > kMaxDiff) return 9999;
122 //___________________________________________________________________________
123 void AliPHOSFastRecParticle::Draw(Option_t *option)
125 // Draw this AliPHOSFastRecParticle with its current attributes
130 //______________________________________________________________________________
131 void AliPHOSFastRecParticle::ExecuteEvent(Int_t event, Int_t px, Int_t py)
133 // Execute action corresponding to one event
134 // This member function is called when a AliPHOSFastRecParticle is clicked with the locator
136 if (!gPad->IsEditable())
139 static TPaveText * clustertext = 0 ;
144 Double_t kRADDEG = 180. / TMath::Pi() ;
145 Coord_t x = Phi() * kRADDEG ;
146 Coord_t y = Theta() * kRADDEG ;
147 clustertext = new TPaveText(x-1, y+1, x+5, y+3, "") ;
150 sprintf( line1, "PID: %s ", (const char*)Name() ) ;
151 sprintf( line2, "ENERGY: %f ", Energy() ) ;
152 clustertext ->AddText(line1) ;
153 clustertext ->AddText(line2) ;
154 clustertext ->Draw("");
169 //____________________________________________________________________________
170 Bool_t AliPHOSFastRecParticle::IsPhoton(const char* purity) const
172 // Rec.Particle is a photon if it has a photon-like shape, fast and neutral
173 // photon-like shape is defined with a purity "low", "medium" or "high"
175 Bool_t photonLike = kFALSE;
176 if (strcmp(purity,"low") ==0) photonLike = TestPIDBit(6);
177 else if (strcmp(purity,"medium")==0) photonLike = TestPIDBit(7);
178 else if (strcmp(purity,"high") ==0) photonLike = TestPIDBit(8);
179 else Error("IsPhoton","Wrong purity type: \'%s\'",purity);
180 if (photonLike && // photon by PCA
181 (TestPIDBit(5)||TestPIDBit(4)||TestPIDBit(3))&& // fast by TOF
182 (TestPIDBit(2)||TestPIDBit(1)||TestPIDBit(0))) // neutral by CPV
188 //____________________________________________________________________________
189 Bool_t AliPHOSFastRecParticle::IsPi0(const char* purity) const
191 // Rec.Particle is a pi0 if it has a pi0-like shape, fast and neutral
192 // pi0-like shape is defined with a purity "low", "medium" or "high"
194 Bool_t pi0Like = kFALSE;
195 if (strcmp(purity,"low") ==0) pi0Like = TestPIDBit(9);
196 else if (strcmp(purity,"medium")==0) pi0Like = TestPIDBit(10);
197 else if (strcmp(purity,"high") ==0) pi0Like = TestPIDBit(11);
198 else Error("IsPi0","Wrong purity type: %s",purity);
199 if (pi0Like && // pi0 by PCA
200 (TestPIDBit(5)||TestPIDBit(4)||TestPIDBit(3))&& // fast by TOF
201 (TestPIDBit(2)||TestPIDBit(1)||TestPIDBit(0))) // neutral by CPV
207 //____________________________________________________________________________
208 Bool_t AliPHOSFastRecParticle::IsElectron(const char* purity) const
210 // Rec.Particle is an electron if it has a photon-like shape, fast and charged
211 // photon-like shape is defined with a purity "low", "medium" or "high"
213 Bool_t photonLike = kFALSE;
214 if (strcmp(purity,"low") ==0) photonLike = TestPIDBit(6);
215 else if (strcmp(purity,"medium")==0) photonLike = TestPIDBit(7);
216 else if (strcmp(purity,"high") ==0) photonLike = TestPIDBit(8);
217 else Error("IsElectron","Wrong purity type: %s",purity);
218 if (photonLike && // photon by PCA
219 (TestPIDBit(5)||TestPIDBit(4)||TestPIDBit(3))&& // fast by TOF
220 !(TestPIDBit(2)||TestPIDBit(1)||TestPIDBit(0))) // charged by CPV
226 //____________________________________________________________________________
227 Bool_t AliPHOSFastRecParticle::IsHadron() const
229 // Rec.Particle is an hadron if it does not look like
230 // a low-purity photon nor low-purity pi0
232 if ( !TestPIDBit(6) && !TestPIDBit(9) ) // not photon nor pi0
238 //____________________________________________________________________________
239 Bool_t AliPHOSFastRecParticle::IsChargedHadron() const
241 // Rec.Particle is a charged hadron if it does not look like
242 // a low-purity photon nor low-purity pi0 and is low-purity charged
244 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
245 !TestPIDBit(2)) // charged by CPV
251 //____________________________________________________________________________
252 Bool_t AliPHOSFastRecParticle::IsNeutralHadron() const
254 // Rec.Particle is a neutral hadron if it does not look like
255 // a low-purity photon nor low-purity pi0 and is high-purity neutral
257 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
258 TestPIDBit(2)) // neutral by CPV
264 //____________________________________________________________________________
265 Bool_t AliPHOSFastRecParticle::IsFastChargedHadron() const
267 // Rec.Particle is a fast charged hadron if it does not look like
268 // a low-purity photon nor low-purity pi0, is low-purity charged
269 // and is high-purity fast
271 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
272 !TestPIDBit(2) && // charged by CPV
273 TestPIDBit(5)) // fast by TOF
279 //____________________________________________________________________________
280 Bool_t AliPHOSFastRecParticle::IsSlowChargedHadron() const
282 // Rec.Particle is a slow neutral hadron if it does not look like
283 // a low-purity photon nor low-purity pi0, is high-purity neutral
284 // and is not high-purity fast
286 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
287 !TestPIDBit(2) && // charged by CPV
288 !TestPIDBit(5)) // slow by TOF
295 //____________________________________________________________________________
296 Bool_t AliPHOSFastRecParticle::IsFastNeutralHadron() const
298 // Rec.Particle is a fast neutral hadron if it does not look like
299 // a low-purity photon nor low-purity pi0, is high-purity neutral
300 // and is high-purity fast
302 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
303 TestPIDBit(2) && // neutral by CPV
304 TestPIDBit(5)) // fast by TOF
310 //____________________________________________________________________________
311 Bool_t AliPHOSFastRecParticle::IsSlowNeutralHadron() const
313 // Rec.Particle is a slow neutral hadron if it does not look like
314 // a low-purity photon nor low-purity pi0, is high-purity neutral
315 // and is not high-purity fast
317 if ( !TestPIDBit(6) && !TestPIDBit(9) && // not photon nor pi0
318 TestPIDBit(2) && // neutral by CPV
319 !TestPIDBit(5)) // slow by TOF
325 //____________________________________________________________________________
326 TString AliPHOSFastRecParticle::Name() const
328 // Returns the name of the particle type (only valid if PIDv1 is employed)
332 name = "Undefined particle" ;
335 name = "Photon low purity ";
336 else if (IsPhoton("medium"))
337 name = "Photon medium purity";
338 else if (IsPhoton("high"))
339 name = "Photon high purity ";
342 name += "Pi0 low purity ";
343 else if (IsPi0("medium"))
344 name += "Pi0 medium purity ";
345 else if (IsPi0("high"))
346 name += "Pi0 high purity ";
348 if (IsElectron("low"))
349 name += "Electron low purity ";
350 else if (IsElectron("medium"))
351 name += "Electron medium purity ";
352 else if (IsElectron("high"))
353 name += "Electron high purity ";
357 if (IsChargedHadron()) {
358 name.Prepend("charged ");
359 if (IsFastChargedHadron())
360 name.Prepend("fast ");
361 else if (IsSlowChargedHadron())
362 name.Prepend("slow ");
364 else if (IsNeutralHadron()) {
365 name.Prepend("neutral ");
366 if (IsFastNeutralHadron())
367 name.Prepend("fast ");
368 else if (IsSlowNeutralHadron())
369 name.Prepend("slow ");
377 //______________________________________________________________________________
378 void AliPHOSFastRecParticle::SetType(Int_t type) {
379 // sets the particle type
380 // bit-mask of the particle type means the following:
381 // bits 0,1,2 - neutral particle with low, medium and high purity
382 // bits 3.4,5 - fast particle with low, medium and high purity
383 // bits 6.7,8 - photon shower with low, medium and high purity
384 // bits 9,10,11 - hard-pi0 shower with low, medium and high purity
388 if((type == 127) || (fType == 511) || (fType == 255) ||(fType == 383)||(fType == 447)){
393 if ((fType == 63)|| ((fType < 8)&&(fType > 0)) ){
397 if ( ((fType == 504) || (fType == 505) ||(fType == 248)||(fType == 249)||(fType == 120)||(fType == 121)) ){
401 if ((fType == 448) || (fType == 449) ||(fType == 192)||(fType == 193)||(fType == 64)||(fType == 64)){
405 if((fType == 56)||(fType == 57)){
416 //______________________________________________________________________________
417 void AliPHOSFastRecParticle::Paint(Option_t *)
419 // Paint this ALiRecParticle in theta,phi coordinate as a TMarker with its current attributes
421 Double_t kRADDEG = 180. / TMath::Pi() ;
422 Coord_t x = Phi() * kRADDEG ;
423 Coord_t y = Theta() * kRADDEG ;
424 Color_t markercolor = 1 ;
425 Size_t markersize = 1. ;
426 Style_t markerstyle = 5 ;
428 if (!gPad->IsBatch()) {
429 gVirtualX->SetMarkerColor(markercolor) ;
430 gVirtualX->SetMarkerSize (markersize) ;
431 gVirtualX->SetMarkerStyle(markerstyle) ;
433 gPad->SetAttMarkerPS(markercolor,markerstyle,markersize) ;
434 gPad->PaintPolyMarker(1,&x,&y,"") ;
437 //____________________________________________________________________________
438 void AliPHOSFastRecParticle::Print(Option_t * opt)const
440 // Print the type, energy and momentum of the reconstructed particle
443 message = "\n PID bits are %d%d%d %d%d%d %d%d%d %d%d%d" ;
444 message += ", type is \"%s\"\n" ;
445 message += " (E,Px,Py,Pz) = (% .3e, % .3e, % .3e, % .3e) GeV\n" ;
446 Info("Print", message.Data(),
447 TestPIDBit(0),TestPIDBit(1),
448 TestPIDBit(2),TestPIDBit(3),
449 TestPIDBit(4),TestPIDBit(5),
450 TestPIDBit(6),TestPIDBit(7),
451 TestPIDBit(8),TestPIDBit(9),
452 TestPIDBit(10),TestPIDBit(11),