#include "EVGEN/AliGenSlowNucleons.h"
#include "EVGEN/AliSlowNucleonModelExp.h"
#include "PYTHIA6/AliGenPythia.h"
-#include "STEER/AliMagFMaps.h"
+#include "STEER/AliMagF.h"
#include "STRUCT/AliBODY.h"
#include "STRUCT/AliMAG.h"
#include "STRUCT/AliABSOv0.h"
#include "ITS/AliITSvPPRasymmFMD.h"
#include "TPC/AliTPCv2.h"
#include "TOF/AliTOFv4T0.h"
-#include "RICH/AliRICHv1.h"
+#include "HMPID/AliHMPIDv1.h"
#include "ZDC/AliZDCv2.h"
#include "TRD/AliTRDv1.h"
#include "FMD/AliFMDv1.h"
#include "MUON/AliMUONv1.h"
#include "PHOS/AliPHOSv1.h"
#include "PMD/AliPMDv1.h"
-#include "START/AliSTARTv1.h"
+#include "T0/AliT0v1.h"
#include "EMCAL/AliEMCALv1.h"
-#include "CRT/AliCRTv0.h"
+#include "ACORDE/AliACORDEv1.h"
#include "VZERO/AliVZEROv3.h"
#endif
kGluonRadiation, kNoGluonRadiation
};
-enum PprMag_t
-{
- k2kG, k4kG, k5kG
-};
-
-
// This part for configuration
static PprRun_t srun = test50;
//static PprRun_t srun = kPythia6;
static PprGeo_t sgeo = kNoHoles;
static PprRad_t srad = kGluonRadiation;
-static PprMag_t smag = k5kG;
+//static AliMagF::BMap_t smag = AliMagF::k5kG;
// Comment line
static TString comment;
gener->SetVertexSmear(kPerEvent);
gener->SetTrackingFlag(1);
gener->Init();
-
- if (smag == k2kG) {
+
+ /*
+ if (smag == AliMagF::k2kG) {
comment = comment.Append(" | L3 field 0.2 T");
- } else if (smag == k4kG) {
- comment = comment.Append(" | L3 field 0.4 T");
- } else if (smag == k5kG) {
+ } else if (smag == AliMagF::k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
}
-
+ */
if (srad == kGluonRadiation)
{
if (sgeo == kHoles)
{
- comment = comment.Append(" | Holes for PHOS/RICH");
+ comment = comment.Append(" | Holes for PHOS/HMPID");
} else {
- comment = comment.Append(" | No holes for PHOS/RICH");
+ comment = comment.Append(" | No holes for PHOS/HMPID");
}
printf("\n \n Comment: %s \n \n", comment.Data());
-// Field (L3 0.4 T)
- AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., smag);
- field->SetL3ConstField(0); //Using const. field in the barrel
- rl->CdGAFile();
- gAlice->SetField(field);
+// Field (L3 +-0.2/0.5 T)
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",-1,-1, AliMagF::k5kG)); // +0.5/-0.7
+ /*
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1, 1, AliMagF::k5kG)); // -0.5/-0.7
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",-1,-1, AliMagF::k2kG)); // +0.2/-0.7
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1, 1, AliMagF::k2kG)); // -0.2/+0.7
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1, 0, AliMagF::k5kG)); // +-0.5/0
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1, 0, AliMagF::k2kG)); // +-0.2/0
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 0, 0, AliMagF::k5kGUniform)); // 0/0
+ */
+
//
Int_t iABSO = 0;
Int_t iDIPO = 0;
Int_t iPHOS = 0;
Int_t iPIPE = 1;
Int_t iPMD = 0;
- Int_t iRICH = 0;
+ Int_t iHMPID = 0;
Int_t iSHIL = 0;
- Int_t iSTART = 0;
+ Int_t iT0 = 0;
Int_t iTOF = 1;
Int_t iTPC = 1;
Int_t iTRD = 1;
Int_t iZDC = 0;
Int_t iEMCAL = 0;
Int_t iVZERO = 0;
- Int_t iCRT = 0;
+ Int_t iACORDE = 0;
//=================== Alice BODY parameters =============================
AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
}
- if (iRICH)
+ if (iHMPID)
{
- //=================== RICH parameters ===========================
- AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+ //=================== HMPID parameters ===========================
+ AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");
}
if (sgeo == kHoles) {
// With hole in front of PHOS
TRD->SetPHOShole();
- // With hole in front of RICH
- TRD->SetRICHhole();
+ // With hole in front of HMPID
+ TRD->SetHMPIDhole();
}
// Switch on TR
AliTRDsim *TRDsim = TRD->CreateTR();
AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
}
- if (iSTART)
+ if (iT0)
{
- //=================== START parameters ============================
- AliSTART *START = new AliSTARTv1("START", "START Detector");
+ //=================== T0 parameters ============================
+ AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
}
if (iEMCAL)
AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "G56_2_55_19");
}
- if (iCRT)
+ if (iACORDE)
{
- //=================== CRT parameters ============================
- AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+ //=================== ACORDE parameters ============================
+ AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
}
if (iVZERO)
{
- //=================== CRT parameters ============================
+ //=================== ACORDE parameters ============================
AliVZERO *VZERO = new AliVZEROv3("VZERO", "normal VZERO");
}