cout << " primary # " << index << " = " << listofprimaries[index] << endl ;
switch(recparticle->GetType())
{
- case kGAMMA:
+ case AliPHOSFastRecParticle::kGAMMA:
fhPhotonEnergy->Fill(recparticle->Energy() ) ;
//fhPhotonPositionX->Fill(recpart. ) ;
//fhPhotonPositionY->Fill(recpart. ) ;
cout << "PHOTON" << endl;
break;
- case kELECTRON:
+ case AliPHOSFastRecParticle::kELECTRON:
fhElectronEnergy->Fill(recparticle->Energy() ) ;
//fhElectronPositionX->Fill(recpart. ) ;
//fhElectronPositionY->Fill(recpart. ) ;
cout << "ELECTRON" << endl;
break;
- case kNEUTRALHA:
+ case AliPHOSFastRecParticle::kNEUTRALHA:
fhNeutralHadronEnergy->Fill(recparticle->Energy() ) ;
//fhNeutralHadronPositionX->Fill(recpart. ) ;
//fhNeutralHadronPositionY->Fill(recpart. ) ;
cout << "NEUTRAl HADRON" << endl;
break ;
- case kNEUTRALEM:
+ case AliPHOSFastRecParticle::kNEUTRALEM:
fhNeutralEMEnergy->Fill(recparticle->Energy() ) ;
//fhNeutralEMPositionX->Fill(recpart. ) ;
//fhNeutralEMPositionY->Fill(recpart. ) ;
//cout << "NEUTRAL EM" << endl;
break ;
- case kCHARGEDHA:
+ case AliPHOSFastRecParticle::kCHARGEDHA:
fhChargedHadronEnergy->Fill(recparticle->Energy() ) ;
//fhChargedHadronPositionX->Fill(recpart. ) ;
//fhChargedHadronPositionY->Fill(recpart. ) ;
cout << "CHARGED HADRON" << endl;
break ;
- case kGAMMAHA:
+ case AliPHOSFastRecParticle::kGAMMAHA:
fhPhotonHadronEnergy->Fill(recparticle->Energy() ) ;
//fhPhotonHadronPositionX->Fill(recpart. ) ;
//fhPhotonHadronPositionY->Fill(recpart. ) ;
sprintf(histoname,"Event %d: Incident particles in module %d", evt, module) ;
Double_t tm, tM, pm, pM ; // min and Max theta and phi covered by module
- fGeom->EmcModuleCoverage(module, tm, tM, pm, pM, kDegre) ;
+ fGeom->EmcModuleCoverage(module, tm, tM, pm, pM, AliPHOSGeometry::kDegre) ;
Double_t theta, phi ;
- fGeom->EmcXtalCoverage(theta, phi, kDegre) ;
+ fGeom->EmcXtalCoverage(theta, phi, AliPHOSGeometry::kDegre) ;
Int_t tdim = (Int_t)( (tM - tm) / theta ) ;
Int_t pdim = (Int_t)( (pM - pm) / phi ) ;
Text_t histoname[80] ;
sprintf(histoname,"Event %d: Reconstructed particles in module %d", fEvt, module) ;
Double_t tm, tM, pm, pM ; // min and Max theta and phi covered by module
- fGeom->EmcModuleCoverage(module, tm, tM, pm, pM, kDegre) ;
+ fGeom->EmcModuleCoverage(module, tm, tM, pm, pM, AliPHOSGeometry::kDegre) ;
Double_t theta, phi ;
- fGeom->EmcXtalCoverage(theta, phi, kDegre) ;
+ fGeom->EmcXtalCoverage(theta, phi, AliPHOSGeometry::kDegre) ;
Int_t tdim = (Int_t)( (tM - tm) / theta ) ;
Int_t pdim = (Int_t)( (pM - pm) / phi ) ;
tm -= theta ;
#include "TPad.h"
#include "TPaveText.h"
-ClassImp(AliPHOSFastRecParticle)
+ClassImp(AliPHOSFastRecParticle) ;
//____________________________________________________________________________
AliPHOSFastRecParticle::AliPHOSFastRecParticle(const AliPHOSFastRecParticle & rp)
// --- AliRoot header files ---
-const static Int_t kUNDEFINED =-1 ;
-const static Int_t kNEUTRALEM = 0 ;
-const static Int_t kNEUTRALHA = 1 ;
-const static Int_t kGAMMA = 2 ;
-const static Int_t kGAMMAHA = 3 ;
-const static Int_t kABSURDEM = 4 ;
-const static Int_t kABSURDHA = 5 ;
-const static Int_t kELECTRON = 6 ;
-const static Int_t kCHARGEDHA = 7 ;
-
-
class AliPHOSFastRecParticle : public TParticle {
public:
fIndexInList = val ;
}
+ enum EParticleType { kUNDEFINED=-1, kNEUTRALEM, kNEUTRALHA, kGAMMA , kGAMMAHA ,
+ kABSURDEM, kABSURDHA , kELECTRON, kCHARGEDHA } ;
+
protected:
Int_t fIndexInList ; // the index of this RecParticle in the list stored in TreeR (to be set by analysis)
#include "AliPHOSPpsdRecPoint.h"
#include "AliConst.h"
-ClassImp(AliPHOSGeometry)
+ClassImp(AliPHOSGeometry) ;
- AliPHOSGeometry * AliPHOSGeometry::fgGeom = 0 ;
+AliPHOSGeometry * AliPHOSGeometry::fgGeom = 0 ;
+
+const TString AliPHOSGeometry::kDegre = "deg" ;
+const TString AliPHOSGeometry::kRadian = "rad" ;
//____________________________________________________________________________
AliPHOSGeometry::~AliPHOSGeometry(void)
{
// Initializes the PHOS parameters
- kDegre = "deg" ;
- kRadian = "rad" ;
-
fRotMatrixArray = new TObjArray(fNModules) ;
cout << "PHOS geometry setup: parameters for option " << fName << " " << fTitle << endl ;
virtual void GetGlobal(const AliRecPoint* RecPoint, TVector3 & gpos, TMatrix & gmat) ;
virtual void GetGlobal(const AliRecPoint* RecPoint, TVector3 & gpos) ;
- static TString kDegre ;
- static TString kRadian ;
+ const static TString kDegre ;
+ const static TString kRadian ;
protected:
#include "AliPHOSReconstructioner.h"
#include "AliPHOSClusterizer.h"
+#include "AliPHOSFastRecParticle.h"
ClassImp(AliPHOSReconstructioner)
primaries = (rp->GetPHOSTrackSegment())->GetPrimariesEmc(nprimaries);
switch(rp->GetType())
{
- case kNEUTRALEM:
+ case AliPHOSFastRecParticle::kNEUTRALEM:
strcpy( particle, "NEUTRAL_EM");
break;
- case kNEUTRALHA:
+ case AliPHOSFastRecParticle::kNEUTRALHA:
strcpy(particle, "NEUTRAL_HA");
break;
- case kGAMMA:
+ case AliPHOSFastRecParticle::kGAMMA:
strcpy(particle, "GAMMA");
break ;
- case kGAMMAHA:
+ case AliPHOSFastRecParticle::kGAMMAHA:
strcpy(particle, "GAMMA_H");
break ;
- case kABSURDEM:
+ case AliPHOSFastRecParticle::kABSURDEM:
strcpy(particle, "ABSURD_EM") ;
break ;
- case kABSURDHA:
+ case AliPHOSFastRecParticle::kABSURDHA:
strcpy(particle, "ABSURD_HA") ;
break ;
- case kELECTRON:
+ case AliPHOSFastRecParticle::kELECTRON:
strcpy(particle, "ELECTRON") ;
break ;
- case kCHARGEDHA:
+ case AliPHOSFastRecParticle::kCHARGEDHA:
strcpy(particle, "CHARGED_HA") ;
break ;
}
{
// Generate a particle type using the performance of the EMC+PPSD setup
- Int_t rv = kUNDEFINED ;
+ Int_t rv = AliPHOSFastRecParticle::kUNDEFINED ;
Int_t charge = (Int_t)rp.GetPDG()->Charge() ;
Int_t test ;
Float_t ran ;
case 22: // it's a photon
ran = fRan.Rndm() ;
if ( ran <= 0.5 ) // 50 %
- rv = kGAMMA ;
+ rv = AliPHOSFastRecParticle::kGAMMA ;
else {
ran = fRan.Rndm() ;
if( ran <= 0.9498 )
- rv = kNEUTRALEM ;
+ rv = AliPHOSFastRecParticle::kNEUTRALEM ;
else
- rv = kNEUTRALHA ;
+ rv = AliPHOSFastRecParticle::kNEUTRALHA ;
}
break ;
case 2112: // it's a neutron
ran = fRan.Rndm() ;
if ( ran <= 0.9998 )
- rv = kNEUTRALHA ;
+ rv = AliPHOSFastRecParticle::kNEUTRALHA ;
else
- rv = kNEUTRALEM ;
+ rv = AliPHOSFastRecParticle::kNEUTRALEM ;
break ;
case -2112: // it's a anti-neutron
ran = fRan.Rndm() ;
if ( ran <= 0.9984 )
- rv = kNEUTRALHA ;
+ rv = AliPHOSFastRecParticle::kNEUTRALHA ;
else
- rv = kNEUTRALEM ;
+ rv = AliPHOSFastRecParticle::kNEUTRALEM ;
break ;
case 11: // it's a electron
ran = fRan.Rndm() ;
if ( ran <= 0.9996 )
- rv = kELECTRON ;
+ rv = AliPHOSFastRecParticle::kELECTRON ;
else
- rv = kCHARGEDHA ;
+ rv = AliPHOSFastRecParticle::kCHARGEDHA ;
break;
case -11: // it's a positon
ran = fRan.Rndm() ;
if ( ran <= 0.9996 )
- rv = kELECTRON ;
+ rv = AliPHOSFastRecParticle::kELECTRON ;
else
- rv = kCHARGEDHA ;
+ rv = AliPHOSFastRecParticle::kCHARGEDHA ;
break;
case -1: // it's a charged
ran = fRan.Rndm() ;
if ( ran <= 0.9996 )
- rv = kCHARGEDHA ;
+ rv = AliPHOSFastRecParticle::kCHARGEDHA ;
else
- rv = kGAMMA ;
+ rv = AliPHOSFastRecParticle::kGAMMA ;
break ;
}