- // Set External decayer
- AliDecayer* decayer = new AliDecayerPythia();
- decayer->SetForceDecay(kAll);
- decayer->Init();
- gMC->SetExternalDecayer(decayer);
-
- // Physics process control
- gMC ->SetProcess("DCAY",1);
- gMC ->SetProcess("PAIR",1);
- gMC ->SetProcess("COMP",1);
- gMC ->SetProcess("PHOT",1);
- gMC ->SetProcess("PFIS",0);
- gMC ->SetProcess("DRAY",0);
- gMC ->SetProcess("ANNI",1);
- gMC ->SetProcess("BREM",1);
- gMC ->SetProcess("MUNU",1);
- //xx gMC ->SetProcess("CKOV",1);
- gMC ->SetProcess("HADR",1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
- gMC ->SetProcess("LOSS",2);
- gMC ->SetProcess("MULS",1);
- //xx gMC ->SetProcess("RAYL",1);
-
- // Energy cuts
- // (in development)
- Float_t cut = 1.e-3; // 1MeV cut by default
- Float_t tofmax = 1.e10;
-
- gMC ->SetCut("CUTGAM",cut);
- gMC ->SetCut("CUTELE",cut);
- gMC ->SetCut("CUTNEU",cut);
- gMC ->SetCut("CUTHAD",cut);
- gMC ->SetCut("CUTMUO",cut);
- gMC ->SetCut("BCUTE",cut);
- gMC ->SetCut("BCUTM",cut);
- gMC ->SetCut("DCUTE",cut);
- gMC ->SetCut("DCUTM",cut);
- gMC ->SetCut("PPCUTM",cut);
- gMC ->SetCut("TOFMAX",tofmax);
-
- // =============================
- // Event generator
- // =============================
-
- // --- Specify event type to be tracked through the ALICE setup
- // --- All positions are in cm, angles in degrees, and P and E in GeV
-
- Int_t nParticles = 100;
- if (gSystem->Getenv("CONFIG_NPARTICLES"))
- nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
-
- AliGenCocktail *gener = new AliGenCocktail();
- gener->SetPhiRange(0, 360);
- // Set pseudorapidity range from -8 to 8.
- Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
- Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
- gener->SetThetaRange(thmin,thmax);
- gener->SetOrigin(0, 0, 0); //vertex position
- gener->SetSigma(0, 0, 0); //Sigma in (X,Y,Z) (cm) on IP position
-
- AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(nParticles);
- hijingparam->SetMomentumRange(0.2, 999);
- gener->AddGenerator(hijingparam,"HIJING PARAM",1);
-
- // AliGenBox *genbox = new AliGenBox(nParticles);
- // genbox->SetPart(22);
- // genbox->SetPtRange(0.3, 10.00);
- // gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
- gener->Init();