#include "PYTHIA6/AliDecayerPythia.h"
#include "PYTHIA6/AliGenPythia.h"
#include "TDPMjet/AliGenDPMjet.h"
-#include "STEER/AliMagFCheb.h"
+#include "STEER/AliMagF.h"
#include "STRUCT/AliBODY.h"
#include "STRUCT/AliMAG.h"
#include "STRUCT/AliABSOv3.h"
"kPythia6", "kPhojet"
};
-//--- Magnetic Field ---
enum Mag_t
{
kNoField, k5kG, kFieldMax
// FIELD
//
- AliMagFCheb* field = 0x0;
+ AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagFCheb("Maps","Maps", 2, 0., 10., AliMagFCheb::k2kG);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagFCheb("Maps","Maps", 2, 1., 10., AliMagFCheb::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
+ TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
- gAlice->SetField(field);
-
-
Int_t iABSO = 1;
Int_t iACORDE= 0;
AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
AliTRDgeometry *geoTRD = TRD->GetGeometry();
- // Partial geometry: modules at 0,8,9,17
+ // Partial geometry: modules at 0,1,7,8,9,10,17
// starting at 3h in positive direction
- geoTRD->SetSMstatus(1,0);
geoTRD->SetSMstatus(2,0);
geoTRD->SetSMstatus(3,0);
geoTRD->SetSMstatus(4,0);
geoTRD->SetSMstatus(5,0);
geoTRD->SetSMstatus(6,0);
- geoTRD->SetSMstatus(7,0);
- geoTRD->SetSMstatus(10,0);
geoTRD->SetSMstatus(11,0);
geoTRD->SetSMstatus(12,0);
geoTRD->SetSMstatus(13,0);