//_______________________________________________________________________
AliMagF::AliMagF(const char *name, const char* title, Double_t factorSol, Double_t factorDip,
- BMap_t maptype, BeamType_t btype, Double_t benergy, Int_t integ, Double_t fmax,
- const char* path):
+ BMap_t maptype, BeamType_t bt, Double_t be,Int_t integ, Double_t fmax, const char* path):
TVirtualMagField(name),
fMeasuredMap(0),
fMapType(maptype),
fSolenoid(0),
- fBeamType(btype),
- fBeamEnergy(benergy),
+ fBeamType(bt),
+ fBeamEnergy(be),
//
fInteg(integ),
fPrecInteg(1),
//
if (fBeamEnergy<=0 && fBeamType!=kNoBeamField) {
if (fBeamType == kBeamTypepp) fBeamEnergy = 7000.; // max proton energy
- else if (fBeamType == kBeamTypeAA) fBeamEnergy = 5500; // max PbPb energy
+ else if (fBeamType == kBeamTypeAA) fBeamEnergy = 2750; // max PbPb energy
AliInfo("Maximim possible beam energy for requested beam is assumed");
}
const char* parname = 0;
factorSol,(fMapType==k5kG||fMapType==k5kGUniform)?5.:2.,
fDipoleOFF ? "OFF":"ON",factorDip,fMapType==k5kGUniform?" |Constant Field!":""));
AliInfo(Form("Machine B fields for %s beam (%.0f GeV): QGrad: %.4f Dipole: %.4f",
- fBeamType==kBeamTypeAA ? "A-A":(fBeamType==kBeamTypepp ? "p-p":"OFF"),
- fBeamEnergy,fQuadGradient,fDipoleField));
+ bt==kBeamTypeAA ? "A-A":(bt==kBeamTypepp ? "p-p":"OFF"),be,fQuadGradient,fDipoleField));
}
//_______________________________________________________________________
//_____________________________________________________________________________
AliMagF* AliMagF::CreateFieldMap(Float_t l3Cur, Float_t diCur, Int_t convention, Bool_t uniform,
- Float_t beamenergy, const Char_t *beamtype, const Char_t *path)
+ Float_t sqrts, const Char_t *beamtype, const Char_t *path)
{
//------------------------------------------------
// The magnetic field map, defined externally...
// LHC and DCS08 conventions have opposite dipole polarities
if ( GetPolarityConvention() != convention) sclDip = -sclDip;
//
- return new AliMagF("MagneticFieldMap", ttl,sclL3,sclDip,map,btype,beamenergy,2,10.,path);
+ return new AliMagF("MagneticFieldMap", ttl,sclL3,sclDip,map,btype,sqrts/2,2,10.,path);
//
}
enum {kOverrideGRP=BIT(14)}; // don't recreate from GRP if set
//
AliMagF();
- AliMagF(const char *name, const char* title, Double_t factorSol=1., Double_t factorDip=1.,
- BMap_t maptype = k5kG, BeamType_t btype=kBeamTypepp,
- Double_t benergy=-1., Int_t integ=2, Double_t fmax=10,
- const char* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root");
+ AliMagF(const char *name, const char* title,Double_t factorSol=1., Double_t factorDip=1.,
+ BMap_t maptype = k5kG, BeamType_t btype=kBeamTypepp, Double_t benergy=-1,
+ Int_t integ=2, Double_t fmax=10,const char* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root");
AliMagF(const AliMagF& src);
AliMagF& operator=(const AliMagF& src);
virtual ~AliMagF();
static Int_t GetPolarityConvention() {return Int_t(fgkPolarityConvention);}
static AliMagF* CreateFieldMap(Float_t l3Current=-30000., Float_t diCurrent=-6000.,
Int_t convention=0, Bool_t uniform = kFALSE,
- Float_t benergy=7000., const Char_t* btype="pp",
+ Float_t sqrts=14000, const Char_t* btype="pp",
const Char_t* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root");
//
protected: