#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"
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)
{
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;