*/
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
-#include <TRandom.h>
+#include <TPDGCode.h>
#include <TSystem.h>
#include <TVirtualMC.h>
#include <TGeant3.h>
{
comment = comment.Append(" Flat in FMD1 range");
AliGenBox* gener = new AliGenBox(2000);
- gener->SetPart(211);
+ gener->SetPart(kPiPlus);
gener->SetMomentumRange(3,4);
gener->SetPhiRange(0, 360);
gener->SetThetaRange(0.77, 3.08);
{
comment = comment.Append(" Flat in FMD2 range");
AliGenBox* gener = new AliGenBox(2000);
- gener->SetPart(211);
+ gener->SetPart(kPiPlus);
gener->SetMomentumRange(3,4);
gener->SetPhiRange(0, 360);
gener->SetThetaRange(2.95, 20.42);
{
comment = comment.Append(" Flat in FMD3 range");
AliGenBox* gener = new AliGenBox(2000);
- gener->SetPart(211);
+ gener->SetPart(kPiPlus);
gener->SetMomentumRange(3,4);
gener->SetPhiRange(0, 360);
gener->SetThetaRange(155.97, 176.73);
gener->SetPhiRange(0, 360);
AliGenBox* gener3 = new AliGenBox(2000);
gener3->SetThetaRange(155.97, 176.73);
- gener3->SetPart(211);
+ gener3->SetPart(kPiPlus);
gener->AddGenerator(gener3, "FMD3", .33);
AliGenBox* gener2 = new AliGenBox(2000);
gener2->SetThetaRange(2.95, 20.42);
- gener2->SetPart(211);
+ gener2->SetPart(kPiPlus);
gener->AddGenerator(gener2, "FMD2", .33);
AliGenBox* gener1 = new AliGenBox(2000);
gener1->SetThetaRange(0.77, 3.08);
- gener1->SetPart(211);
+ gener1->SetPart(kPiPlus);
gener->AddGenerator(gener1, "FMD1", .34);
gGener = gener;
}
//
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
gener->SetOrigin(0.,0.,-800.);
//vertex position
gener->SetSigma(0.1,0.1,0.0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(211); //GEANT particle type
+ gener->SetPart(kPiPlus);
gener->Init();
// Activate this line if you want the vertex smearing to happen
// track by track
//
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
gener->SetOrigin(0.,0.,-50.);
//vertex position
gener->SetSigma(0.1,0.1,0.0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(211); //GEANT particle type
+ gener->SetPart(kPiPlus);
gener->Init();
// Activate this line if you want the vertex smearing to happen
// track by track
//vertex position
gener->SetSigma(1.,1.,0.0); //Sigma in (X,Y,Z) (cm) on IP position
- //gener->SetPart(211); //GEANT particle type
+ //gener->SetPart(kPiPlus);
//cout<<"kPiPlus="<<kPiPlus<<endl;
- gener->SetPart(kPiPlus); //GEANT particle type
+ gener->SetPart(kPiPlus);
gener->Init();
// Activate this line if you want the vertex smearing to happen
// track by track
+#include <TPDGCode.h>
+
void Config(){
// Set Random Number seed
// gRandom->SetSeed(12345);
gener->SetOrigin(0.,0.,-50.);
//vertex position
gener->SetSigma(0.1,0.1,0.0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(211); //GEANT particle type
+ gener->SetPart(kPiPlus);
gener->Init();
// Activate this line if you want the vertex smearing to happen
// track by track
+#include <TPDGCode.h>
+
void Config(){
// 7-DEC-2000 09:00
// Switch on Transition Radiation simulation. 6/12/00 18:00
gener->SetPhi(180.);
gener->SetTheta(95.);
gener->SetOrigin(0,0,0); //vertex position
- gener->SetPart(13); //GEANT particle type
+ gener->SetPart(kMuonMinus);
gener->Init();
//
// Activate this line if you want the vertex smearing to happen
gener->SetMomentumRange(20.,20.1);
gener->SetPhiRange(0., 360.);
gener->SetThetaRange(171.000,178.001);
- gener->SetPart(13); // Muons
+ gener->SetPart(kMuonMinus); // Muons
gener->SetOrigin(0.,0., 0.); //vertex position
gener->SetSigma(0.0, 0.0, 0.0); //Sigma in (X,Y,Z) (cm) on IP position
}
gener->SetPhiRange(0.);
gener->SetThetaRange(0.);
gener->SetOrigin(30,30,-1200);//vertex position
- gener->SetPart(13); //GEANT particle type 13 is muons
+ gener->SetPart(kMuonMinus); //GEANT particle type 13 is muons
}
if (!strcmp(option,"scan")) {
AliGenScan *gener = new AliGenScan(-1);
+#include <TPDGCode.h>
+
static Int_t eventsPerRun = 100;
enum PprGeo_t
{
// gener->AddGenerator(hijingparam,"HIJING PARAM",1);
// AliGenBox *genbox = new AliGenBox(nParticles);
-// genbox->SetPart(22);
+// genbox->SetPart(kGamma);
// genbox->SetPtRange(0.3, 10.00);
// gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
// gener->Init();
gener->SetOrigin(0,0,0); //vertex position
gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(22);
+ gener->SetPart(kGamma);
gener->Init();
//
#include <TSystem.h>
#include <TVirtualMC.h>
#include <TGeant3TGeo.h>
+#include <TPDGCode.h>
#include "STEER/AliRunLoader.h"
#include "STEER/AliRun.h"
#include "STEER/AliConfig.h"
gener->SetSigma(0, 0, 0); //Sigma in (X,Y,Z) (cm) on IP position
AliGenBox *genbox = new AliGenBox(nParticles);
- genbox->SetPart(11);
+ genbox->SetPart(kElectron);
genbox->SetPtRange(0.3, 10.00);
gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
gener->Init();
+#include <TPDGCode.h>
+
static Int_t eventsPerRun = 100;
enum PprGeo_t
{
void Config()
{
+
cout << "==> Config.C..." << endl;
// Set Random Number seed
gener->SetMomentumRange(0.5, 5.);
gener->SetPhiRange(260., 280.);
gener->SetThetaRange(82.,98.);
- gener->SetPart(22);
+ gener->SetPart(kGamma);
gener->SetOrigin(0, 0, 0); //vertex position
gener->SetSigma(0, 0, 0); //Sigma in (X,Y,Z) (cm) on IP position
+#include <TPGCode.h>
+
static Int_t eventsPerRun = 100;
enum PprGeo_t
{
void Config()
{
+
// 7-DEC-2000 09:00
// Switch on Transition Radiation simulation. 6/12/00 18:00
// iZDC=1 7/12/00 09:00
// gener->AddGenerator(hijingparam,"HIJING PARAM",1);
// AliGenBox *genbox = new AliGenBox(nParticles);
-// genbox->SetPart(22);
+// genbox->SetPart(kGamma);
// genbox->SetPtRange(0.3, 10.00);
// gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
// gener->Init();
gener->SetOrigin(0,0,0); //vertex position
gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(22);
+ gener->SetPart(kGamma);
gener->Init();
//
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
gener->SetPhi(180.);
gener->SetTheta(5.);
gener->SetOrigin(0,0,0); //vertex position
- gener->SetPart(13); //GEANT particle type
+ gener->SetPart(kNeutron);
gGener = gener;
}
break;
gener->SetOrigin(0,0,0);
//vertex position
gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(5); //GEANT particle type
+ gener->SetPart(kMuonPlus);
gGener = gener;
}
break;
gener->SetThetaRange(9,9);
//vertex position
gener->SetSigma(6,6,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(5);
+ gener->SetPart(kMuonPlus);
gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
gGener = gener;
}
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TDatime.h>
#include <TSystem.h>
// electron generator
AliGenBox *eminus= new AliGenBox(Npart);
- eminus->SetPart(11);
+ eminus->SetPart(kElectron);
eminus->SetMomentumRange(momen,momen);
eminus->SetPhiRange(0,360);
eminus->SetThetaRange(thmin,thmax);
// positron generator
AliGenBox *eplus= new AliGenBox(Npart);
- eplus->SetPart(-11);
+ eplus->SetPart(kPositron);
eplus->SetMomentumRange(momen,momen);
eplus->SetPhiRange(0,360);
eplus->SetThetaRange(thmin,thmax);
// mu- generator
AliGenBox *muminus= new AliGenBox(Npart);
- muminus->SetPart(13);
+ muminus->SetPart(kMuonMinus);
muminus->SetMomentumRange(momen,momen);
muminus->SetPhiRange(0,360);
muminus->SetThetaRange(thmin,thmax);
// mu+ generator
AliGenBox *muplus= new AliGenBox(Npart);
- muplus->SetPart(-13);
+ muplus->SetPart(kMuonPlus);
muplus->SetMomentumRange(momen,momen);
muplus->SetPhiRange(0,360);
muplus->SetThetaRange(thmin,thmax);
// pi- generator
AliGenBox *piminus= new AliGenBox(Npart);
- piminus->SetPart(-211);
+ piminus->SetPart(kPiMinus);
piminus->SetMomentumRange(momen,momen);
piminus->SetPhiRange(0,360);
piminus->SetThetaRange(thmin,thmax);
// pi+ generator
AliGenBox *piplus= new AliGenBox(Npart);
- piplus->SetPart(211);
+ piplus->SetPart(kPiPlus);
piplus->SetMomentumRange(momen,momen);
piplus->SetPhiRange(0,360);
piplus->SetThetaRange(thmin,thmax);
// K- generator
AliGenBox *kminus= new AliGenBox(Npart);
- kminus->SetPart(-321);
+ kminus->SetPart(kKMinus);
kminus->SetMomentumRange(momen,momen);
kminus->SetPhiRange(0,360);
kminus->SetThetaRange(thmin,thmax);
// K+ generator
AliGenBox *kplus= new AliGenBox(Npart);
- kplus->SetPart(321);
+ kplus->SetPart(kKPlus);
kplus->SetMomentumRange(momen,momen);
kplus->SetPhiRange(0,360);
kplus->SetThetaRange(thmin,thmax);
// p generator
AliGenBox *proton= new AliGenBox(Npart);
- proton->SetPart(2212);
+ proton->SetPart(kProton);
proton->SetMomentumRange(momen,momen);
proton->SetPhiRange(0,360);
proton->SetThetaRange(thmin,thmax);
// anti-p generator
AliGenBox *aproton= new AliGenBox(Npart);
- aproton->SetPart(-2212);
+ aproton->SetPart(kProtonBar);
aproton->SetMomentumRange(momen,momen);
aproton->SetPhiRange(0,360);
aproton->SetThetaRange(thmin,thmax);
+#include <TPDGCode.h>
+
void Config()
{
-
Int_t iField = 0;
// libraries required by geant321
AliGenBox *genEl = new AliGenBox(100);
genEl->SetOrigin(0,0,0); // Vertex position
genEl->SetSigma(0,0,0); // Sigma in (X,Y,Z) (cm) on IP position
- genEl->SetPart(11); // Only electrons
+ genEl->SetPart(kElectron); // Only electrons
AliGenBox *genPi = new AliGenBox(100);
genPi->SetOrigin(0,0,0); // Vertex position
genPi->SetSigma(0,0,0); // Sigma in (X,Y,Z) (cm) on IP position
- genPi->SetPart(-211); // Only pions
+ genPi->SetPart(kPiMinus); // Only pions
gener->AddGenerator(genEl,"Electrons",1);
gener->AddGenerator(genPi,"Pions" ,1);
genbox->SetMomentumRange(10,11.);
genbox->SetPhiRange(270.5,270.7);
genbox->SetThetaRange(90.5,90.7);
- genbox->SetPart(22);
+ genbox->SetPart(kGamma);
gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
gener->Init(); // Initialization of the coctail generator
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
gener->AddGenerator(hijingparam,"HIJING PARAM",1);
// AliGenBox *genbox = new AliGenBox(nParticles);
-// genbox->SetPart(22);
+// genbox->SetPart(kGamma);
// genbox->SetPtRange(0.3, 10.00);
// gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
gener->Init();
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
// g1->SetMomentumRange(0, 3);
// g1->SetOrigin(0,0,0); //vertex position
// g1->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
- // g1->SetPart(310); //K0short(310), Lambda(3122)
+ // g1->SetPart(kK0Short); //K0short(310), Lambda(3122)
/***************** M E V S I M ********************************/
gener->SetPhi(180.);
gener->SetTheta(5.);
gener->SetOrigin(0,0,0); //vertex position
- gener->SetPart(13); //GEANT particle type
+ gener->SetPart(kNeutron);
gGener = gener;
}
break;
gener->SetOrigin(0,0,0);
//vertex position
gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(5); //GEANT particle type
+ gener->SetPart(kMuonPlus);
gGener = gener;
}
break;
gener->SetThetaRange(9,9);
//vertex position
gener->SetSigma(6,6,0); //Sigma in (X,Y,Z) (cm) on IP position
- gener->SetPart(5);
+ gener->SetPart(kMuonPlus);
gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
gGener = gener;
}
//gener->SetPhiRange(0,0);
//gener->SetThetaRange(5., 5.);
//gener->SetOrigin(0,0,0); //vertex position
-//gener->SetPart(14) //GEANT particle type
+//gener->SetPart(kProton);
//*********************************************
// Example for Moving Particle Gun *
gener->SetOrigin(0,0,0);
//vertex position
gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
-gener->SetPart(14) //GEANT particle type
+gener->SetPart(kProton);
*/
//**************************************
// Example for HIJING Parameterisation *
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
// Particle guns for the barrel part (taken from RichConfig)
AliGenFixed *pG1=new AliGenFixed(1);
- pG1->SetPart(2212);
+ pG1->SetPart(kProton);
pG1->SetMomentum(2.5);
pG1->SetTheta(109.5-3);
pG1->SetPhi(10);
gener->AddGenerator(pG1,"g1",1);
AliGenFixed *pG2=new AliGenFixed(1);
- pG2->SetPart(211);
+ pG2->SetPart(kPiPlus);
pG2->SetMomentum(1.0);
pG2->SetTheta( 90.0-3);
pG2->SetPhi(10);
gener->AddGenerator(pG2,"g2",1);
AliGenFixed *pG3=new AliGenFixed(1);
- pG3->SetPart(-211);
+ pG3->SetPart(kPiMinus);
pG3->SetMomentum(1.5);
pG3->SetTheta(109.5-3);
pG3->SetPhi(30);
gener->AddGenerator(pG3,"g3",1);
AliGenFixed *pG4=new AliGenFixed(1);
- pG4->SetPart(321);
+ pG4->SetPart(kKPlus);
pG4->SetMomentum(0.7);
pG4->SetTheta( 90.0-3);
pG4->SetPhi(30);
gener->AddGenerator(pG4,"g4",1);
AliGenFixed *pG5=new AliGenFixed(1);
- pG5->SetPart(-321);
+ pG5->SetPart(kKMinus);
pG5->SetMomentum(1.0);
pG5->SetTheta( 70.0-3);
pG5->SetPhi(30);
gener->AddGenerator(pG5,"g5",1);
AliGenFixed *pG6=new AliGenFixed(1);
- pG6->SetPart(-2212);
+ pG6->SetPart(kProtonBar);
pG6->SetMomentum(2.5);
pG6->SetTheta( 90.0-3);
pG6->SetPhi(50);
gener->AddGenerator(pG6,"g6",1);
AliGenFixed *pG7=new AliGenFixed(1);
- pG7->SetPart(-211);
+ pG7->SetPart(kPiMinus);
pG7->SetMomentum(0.7);
pG7->SetTheta( 70.0-3);
pG7->SetPhi(50);
// Electrons for TRD
AliGenFixed *pG8=new AliGenFixed(1);
- pG8->SetPart(11);
+ pG8->SetPart(kElectron);
pG8->SetMomentum(1.2);
pG8->SetTheta( 95.0);
pG8->SetPhi(190);
gener->AddGenerator(pG8,"g8",1);
AliGenFixed *pG9=new AliGenFixed(1);
- pG9->SetPart(-11);
+ pG9->SetPart(kPositron);
pG9->SetMomentum(1.2);
pG9->SetTheta( 85.0);
pG9->SetPhi(190);
gphos->SetMomentumRange(10,11.);
gphos->SetPhiRange(270.5,270.7);
gphos->SetThetaRange(90.5,90.7);
- gphos->SetPart(22);
+ gphos->SetPart(kGamma);
gener->AddGenerator(gphos,"GENBOX GAMMA for PHOS",1);
// EMCAL
gemcal->SetMomentumRange(10,11.);
gemcal->SetPhiRange(90.5,199.5);
gemcal->SetThetaRange(90.5,90.7);
- gemcal->SetPart(22);
+ gemcal->SetPart(kGamma);
gener->AddGenerator(gemcal,"GENBOX GAMMA for EMCAL",1);
// MUON
gmuon1->SetMomentumRange(20.,20.1);
gmuon1->SetPhiRange(0., 360.);
gmuon1->SetThetaRange(171.000,178.001);
- gmuon1->SetPart(13); // Muons
+ gmuon1->SetPart(kMuonMinus); // Muons
gener->AddGenerator(gmuon1,"GENBOX MUON1",1);
AliGenBox * gmuon2 = new AliGenBox(1);
gmuon2->SetMomentumRange(20.,20.1);
gmuon2->SetPhiRange(0., 360.);
gmuon2->SetThetaRange(171.000,178.001);
- gmuon2->SetPart(-13); // Muons
+ gmuon2->SetPart(kMuonPlus); // Muons
gener->AddGenerator(gmuon2,"GENBOX MUON1",1);
//TOF
AliGenFixed *gtof=new AliGenFixed(1);
- gtof->SetPart(2212);
+ gtof->SetPart(kProton);
gtof->SetMomentum(2.5);
gtof->SetTheta(95);
pG1->SetPhi(340);