*/
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
-#include <TRandom.h>
+#include <TPDGCode.h>
#include <TSystem.h>
#include <TVirtualMC.h>
#include <TGeant3.h>
#include "PMD/AliPMDv1.h"
#include "T0/AliT0v1.h"
#include "EMCAL/AliEMCALv1.h"
-#include "CRT/AliCRTv0.h"
+#include "ACORDE/AliACORDEv0.h"
#include "VZERO/AliVZEROv2.h"
#endif
// This part for configuration
// EG_t eg = test50;
// EG_t eg = kParam_fmd;
- // EG_t eg = kParam_2000; // kPythia;
- EG_t eg = kFMDFlat;
+ EG_t eg = kParam_2000; // kPythia;
+ // EG_t eg = kFMDFlat;
Geo_t geo = kNoHoles;
Rad_t rad = kGluonRadiation;
Mag_t mag = k5kG;
//
// Libraries needed by GEANT 3.21
//
+ gSystem->Load("EGPythia6.so");
gSystem->Load("libgeant321");
//
gMC->SetProcess("MUNU",1);
gMC->SetProcess("CKOV",1);
gMC->SetProcess("HADR",1);
- gMC->SetProcess("LOSS",2);
+ gMC->SetProcess("LOSS",2); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
gMC->SetProcess("MULS",1);
gMC->SetProcess("RAYL",1);
// Used detectors
//
Bool_t useABSO = kFALSE;
- Bool_t useCRT = kFALSE;
+ Bool_t useACORDE = kFALSE;
Bool_t useDIPO = kFALSE;
Bool_t useFMD = kTRUE;
Bool_t useFRAME = kFALSE;
AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCAL_55_25");
}
- if (useCRT) {
- // =================== CRT parameters ============================
- AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+ if (useACORDE) {
+ // =================== ACORDE parameters ============================
+ AliACORDE *ACORDE = new AliACORDEv0("ACORDE", "normal ACORDE");
}
if (useVZERO) {
{
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;
}