}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps",2, -1.,-1., 10.,mag);
+ AliMagF* field = new AliMagF("Maps","Maps",-1.,-1., mag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iACORDE = 0;
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
rl->LoadgAlice();
gAlice = rl->GetAliRun();
- //AliMagF *magf = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG);
- AliMagF *magf = new AliMagF("Maps","Maps", 2, 1., 1.);
+ //AliMagF *magf = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG);
+ AliMagF *magf = new AliMagF("Maps","Maps", 1., 1.);
//
// Open ESD file and recoveries TTree of ESD objects.
//__________________________________________________________________
// Field (L3 0.4 T)
- AliMagF* field = new AliMagF("Maps","Maps",2, -1., -1., 10.,mag);
+ AliMagF* field = new AliMagF("Maps","Maps",-1., -1., mag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
//__________________________________________________________________
// Field (L3 0.4 T)
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,mag);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,mag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
// TODO: does the factor need adjustment?
const AliMagF* currentMap = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
if (!currentMap) {
- AliMagF* field = new AliMagF("MagneticFieldMap", "BMap", 2, 1., 1., 10., map);
+ AliMagF* field = new AliMagF("MagneticFieldMap", "BMap", 1., 1., map);
TGeoGlobalMagField::Instance()->SetField(field);
HLTInfo("Solenoid Field set to: %f map %d", SolenoidBz, map);
} else if (currentMap->GetMapType()!=map) {
// TODO: does the factor need adjustment?
const AliMagF* currentMap = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
if (!currentMap) {
- AliMagF* field = new AliMagF("MagneticFieldMap", "BMap", 2, 1., 1., 10., map);
+ AliMagF* field = new AliMagF("MagneticFieldMap", "BMap", 1., 1., map);
TGeoGlobalMagField::Instance()->SetField(field);
HLTInfo("Solenoid Field set to: %f map %d", SolenoidBz, map);
} else if (currentMap->GetMapType()!=map) {
AliMagF* field;
if (fMagneticField == 0){
// magnetic field OFF
- field = new AliMagF("Maps","Maps",2,0.,0., 10.,AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps",0.,0.,AliMagF::k5kGUniform);
TGeoGlobalMagField::Instance()->SetField(field);
HLTDebug("Magnetic field is OFF.");
}else{
// magnetic field ON
- field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
HLTDebug("Magnetic field is ON.");
if( fMagneticField < 0 )
Printf("------------------------------");
Printf(" Summary for B: %s",s.Data());
Printf("------------------------------");
- AliMagF* fld = new AliMagF("MagneticFieldMap", s.Data(), 2, fcL3, fcDip, 10., map, path,
- btype,beamenergy);
+ AliMagF* fld = new AliMagF("MagneticFieldMap", s.Data(), fcL3, fcDip, map, btype,beamenergy,path);
TGeoGlobalMagField::Instance()->SetField( fld );
TGeoGlobalMagField::Instance()->Lock();
//
WritePhys(pF); //physics processes
//Field
- if(fMagBG->GetButton(kFld0)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 0,1,10,AliMagF::k2kG));// NO field\n");
- else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k2kG));//0.2 Tesla field\n");
- else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k4kG));//0.4 Tesla field\n");
- else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k5kG));//0.5 Tesla field\n");
- else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k2kG));//-0.2 Tesla field\n");
- else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k4kG));//-0.4 Tesla field\n");
- else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k5kG));//-0.5 Tesla field\n");
+ if(fMagBG->GetButton(kFld0)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",0,0,AliMagF::k5kGUniform));// NO field\n");
+ else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k2kG));//0.2 Tesla field\n");
+ // else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k4kG));//0.4 Tesla field\n");
+ else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k5kG));//0.5 Tesla field\n");
+ else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k2kG));//-0.2 Tesla field\n");
+ // else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k4kG));//-0.4 Tesla field\n");
+ else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k5kG));//-0.5 Tesla field\n");
fprintf(pF," pAL->CdGAFile();\n\n"); //????
//Generator
f<<Form(" gMC->SetExternalDecayer(pDecayer);\n\n");
//Field
- if(fMagBG->GetButton(kFld0)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",0,1,1,10,0)); //no field\n\n");
- else if(fMagBG->GetButton(kFld2)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,0));//0.2 Tesla field\n\n");
- else if(fMagBG->GetButton(kFld4)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,1));//0.4 Tesla field\n\n");
- else if(fMagBG->GetButton(kFld5)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,2));//0.5 Tesla field\n\n");
- else if(fMagBG->GetButton(kFld_2)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,0));//-0.2 Tesla field\n\n");
- else if(fMagBG->GetButton(kFld_4)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,1));//-0.4 Tesla field\n\n");
- else if(fMagBG->GetButton(kFld_5)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,2));//-0.5 Tesla field\n\n");
+ if (fMagBG->GetButton(kFld0)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",0,0,AliMagF::k5kGUniform)); //no field\n\n");
+ else if(fMagBG->GetButton(kFld2)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k2kG));//0.2 Tesla field\n\n");
+// else if(fMagBG->GetButton(kFld4)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,1));//0.4 Tesla field\n\n");
+ else if(fMagBG->GetButton(kFld5)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k5kG));//0.5 Tesla field\n\n");
+ else if(fMagBG->GetButton(kFld_2)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",-1,-1,AliMagF::k2kG));//-0.2 Tesla field\n\n");
+ // else if(fMagBG->GetButton(kFld_4)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,1));//-0.4 Tesla field\n\n");
+ else if(fMagBG->GetButton(kFld_5)->GetState()) f<<Form(" gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",-1,-1,AliMagF::k5kG));//-0.5 Tesla field\n\n");
if(fProBG->GetButton(kDCAY)->GetState()) f<<Form(" gMC->SetProcess(\"DCAY\",1); ");else f<<Form(" gMC->SetProcess(\"DCAY\",0); ");
if(fProBG->GetButton(kPAIR)->GetState()) f<<Form(" gMC->SetProcess(\"PAIR\",1); ");else f<<Form(" gMC->SetProcess(\"PAIR\",0); ");
new TGeant3("C++ Interface to Geant3");
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
if ( !TGeoGlobalMagField::Instance()->GetField() ) {
AliInfo( "Loading field map..." );
- AliMagF* field = new AliMagF( "Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG );
+ AliMagF* field = new AliMagF( "Maps","Maps", 1., 1., AliMagF::k5kG );
TGeoGlobalMagField::Instance()->SetField( field );
TGeoGlobalMagField::Instance()->Lock();
}
}
//=============================================================
// Field (L3 0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1, 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1, AliMagF::k5kG));
//=============================================================
//=================== Alice BODY parameters =============================
AliBODY *BODY = new AliBODY("BODY","Alice envelop");
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
// set mag field
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1., AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
//Set mag field; waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1., AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
//Set mag field; waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1., AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
} printf("Loading field map...\n");
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
// waiting for mag field in CDB
if (!TGeoGlobalMagField::Instance()->GetField()) {
printf("Loading field map...\n");
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
TGeoGlobalMagField::Instance()->SetField(field);
}
// set the magnetic field for track extrapolations
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
// Field (L3 0.4 T)
//Zero magnetic field
- AliMagF* field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ AliMagF* field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
// AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., smag);
TGeoGlobalMagField::Instance()->SetField(field);
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 0;
Int_t iDIPO = 0;
//AliGeomManager::LoadGeometry("/lustre/alice/jacek/sim/HEADJB/flatPt_uniB/0/geometry.root");
// set magnetic field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
gROOT->LoadMacro("$ALICE_ROOT/PWG1/macros/RunPerformanceTaskMC.C");
RunPerformanceTaskMC(chain, kTRUE, kFALSE, kFALSE, 0);
// set magnetic field
// the best is to create macro MagField.C with the line:
- // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
gProof->Exec("gROOT->Macro(\"MagField.C\")");
gROOT->LoadMacro("$ALICE_ROOT/PWG1/macros/RunPerformanceTaskMC.C");
//AliGeomManager::LoadGeometry("/lustre/alice/local/TRDdata/SIM/P-Flat/TRUNK/test/RUN0/geometry.root");
// set magnetic field
- gProof->Exec("gROOT->ProcessLine(\"TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\", 2, 1., 1., 10., AliMagF::k5kG))\")");
+ gProof->Exec("gROOT->ProcessLine(\"TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\", 1., 1., AliMagF::k5kG))\")");
gROOT->LoadMacro("$ALICE_ROOT/PWG1/macros/RunPerformanceTaskMC.C");
RunPerformanceTaskMC(chain, kTRUE, kTRUE, kTRUE);
TGeoManager::Import("/u/miranov/proof/geometry.root");
AliGeomManager::LoadGeometry("/u/miranov/proof/geometry.root");
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
}
}
//_______________________________________________________________________
-AliMagF::AliMagF(const char *name, const char* title, Int_t integ,
- Double_t factorSol, Double_t factorDip,
- Double_t fmax, BMap_t maptype, const char* path,
- BeamType_t bt, Double_t be):
+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):
TVirtualMagField(name),
fMeasuredMap(0),
fMapType(maptype),
fSolenoid(0),
- fBeamType(bt),
- fBeamEnergy(be),
+ fBeamType(btype),
+ fBeamEnergy(benergy),
//
fInteg(integ),
fPrecInteg(1),
}
if (fInteg == 0) fPrecInteg = 0;
//
+ if (fBeamEnergy<=0 && fBeamType!=kNoBeamField) {
+ if (fBeamType == kBeamTypepp) fBeamEnergy = 7000.; // max proton energy
+ else if (fBeamType == kBeamTypeAA) fBeamEnergy = 5500; // max PbPb energy
+ AliInfo("Maximim possible beam energy for requested beam is assumed");
+ }
const char* parname = 0;
//
if (fMapType == k2kG) parname = fDipoleOFF ? "Sol12_Dip0_Hole":"Sol12_Dip6_Hole";
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",
- bt==kBeamTypeAA ? "A-A":(bt==kBeamTypepp ? "p-p":"OFF"),be,fQuadGradient,fDipoleField));
+ fBeamType==kBeamTypeAA ? "A-A":(fBeamType==kBeamTypepp ? "p-p":"OFF"),
+ fBeamEnergy,fQuadGradient,fDipoleField));
}
//_______________________________________________________________________
//_______________________________________________________________________
void AliMagF::InitMachineField(BeamType_t btype, Double_t benergy)
{
- if (btype==kNoBeamField || benergy<1.) {
+ if (btype==kNoBeamField) {
fQuadGradient = fDipoleField = fCCorrField = fACorr1Field = fACorr2Field = 0.;
return;
}
// LHC and DCS08 conventions have opposite dipole polarities
if ( GetPolarityConvention() != convention) sclDip = -sclDip;
//
- return new AliMagF("MagneticFieldMap", ttl, 2, sclL3, sclDip, 10., map, path,btype,beamenergy);
+ return new AliMagF("MagneticFieldMap", ttl,sclL3,sclDip,map,btype,beamenergy,2,10.,path);
//
}
enum {kOverrideGRP=BIT(14)}; // don't recreate from GRP if set
//
AliMagF();
- AliMagF(const char *name, const char* title, Int_t integ,
- Double_t factorSol=1., Double_t factorDip=1.,
- Double_t fmax=15, BMap_t maptype = k5kG,
- const char* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",
- BeamType_t btype=kBeamTypepp, Double_t benergy=7000.);
+ 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();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
// Field (L3 +-0.2/0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, -1,-1, 10,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",2, 1, 1, 10,AliMagF::k5kG)); // -0.5/-0.7
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, -1,-1, 10,AliMagF::k2kG)); // +0.2/-0.7
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, 1, 1, 10,AliMagF::k2kG)); // -0.2/+0.7
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, 1, 0, 10,AliMagF::k5kG)); // +-0.5/0
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, 1, 0, 10,AliMagF::k2kG)); // +-0.2/0
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",2, 0, 0, 10,AliMagF::k5kGUniform)); // 0/0
+ 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
*/
//
Float_t factor = bz/(5.); // default b filed in Cheb with minus sign
// was chenged in the Revision ???? (Ruben can you add here number)
- AliMagF* bmap = new AliMagF("MapsExB","MapsExB", 2,factor,1., 10.,AliMagF::k5kG,"$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root");
+ AliMagF* bmap = new AliMagF("MapsExB","MapsExB", factor,TMath::Sign(1.f,factor),AliMagF::k5kG);
AliTPCExBFirst *exb = new AliTPCExBFirst(bmap,0.88*2.6400e+04,50,50,50);
AliTPCExB::SetInstance(exb);
// currently scaling factor 0 used - no magnetic field
//
metaData = CreateMetaObject("AliTPCExB");
- AliMagF* field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kG));
+ AliMagF* field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kG));
AliTPCExBFirst *exbfirst = new AliTPCExBFirst(field,2.83000000000000000e+04,50,50,50);
StoreObject("TPC/Calib/ExB", exbfirst, metaData);
//
// Initialize
//
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG);
AliTPCExB::RegisterField(0,field);
- AliMagF* fieldC0 = new AliMagF("Maps","Maps", 2, 1, 1, 10., AliMagF::k5kG);
+ AliMagF* fieldC0 = new AliMagF("Maps","Maps", 1, 1, AliMagF::k5kG);
AliTPCExB::RegisterField(1,fieldC0);
- AliMagF* fieldC1 = new AliMagF("Maps","Maps", 2, 1, 1, 10., AliMagF::k5kG,kTRUE,"$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root");
+ AliMagF* fieldC1 = new AliMagF("Maps","Maps", 1, 1, AliMagF::k5kG);
AliTPCExB::RegisterField(2,fieldC1);
gSystem->Load("libSTAT.so");
gener->Init();
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
// detectors - only the TPC is relevant
bzfac=1;
bz=5;
}
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, bzfac, 1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", bzfac, 1., smag));
printf("\n\nSET EXB FIELD\t%f\n\n", -bz);
AliTPCcalibDB::Instance()->SetExBField(-bz);
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., smag));
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 0., 1., 10., AliMagF::k5kG));
+ // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., smag));
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
- //TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 0., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
+ //TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
rl->CdGAFile();
//
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
gAlice->SetField(field);
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., smag));
rl->CdGAFile();
//
// MAGNETIC FIELD IN THE BARREL
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., mag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., mag));
// TGeoGlobalMagField::Instance()->GetField()->SetL3ConstField(0); //Using const. field in the barrel
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
rootfile->cd();
}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., mag);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., mag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Magnetic field
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., smag);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., smag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
gener->Init(); // Initialization of the coctail generator
// Field (the last parameter is 1 => L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
// Make sure the current ROOT directory is in galice.root
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
// Field (L3 0.4 T)
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., smag);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., smag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
gGener = gener;
// MAGNETIC FIELD IN THE BARREL
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,mag);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1., mag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,mag);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,,mag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
}
//=============================================================
// Field (L3 0.5 T) outside dimuon spectrometer
- AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",1.,1.,AliMagF::k5kG);
field->SetL3ConstField(0); // Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
}
//=============================================================
// Field (L3 0.5 T) outside dimuon spectrometer
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 10., AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG);
field->SetL3ConstField(0); // Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., mag);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., mag);
field->SetL3ConstField(0); //Using const. field in the barrel
TGeoGlobalMagField::Instance()->SetField(field);
// =============================
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
comment = comment.Append(" | L3 field 0.5 T");
}
// OK
- AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., smag);
+ AliMagF* field = new AliMagF("Maps","Maps", -1., -1., smag);
TGeoGlobalMagField::Instance()->SetField(field);
printf("\n \n Comment: %s \n \n", comment.Data());
//
// create field map
- AliMagF* field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG);
// field-SetL3ConstField(1);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
// FIELD
//
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., mag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., mag));
printf("\n \n Comment: %s \n \n", comment.Data());
gener->Init();
// Field (L3 0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., mag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., mag));
Int_t iABSO = 1;
Int_t iDIPO = 1;
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;\r
if (mag == kNoField) {\r
comment = comment.Append(" | L3 field 0.0 T");\r
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);\r
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);\r
} else if (mag == k5kG) {\r
comment = comment.Append(" | L3 field 0.5 T");\r
// field = new AliMagWrapCheb("Maps","Maps", 2, 1., 10., AliMagWrapCheb::k5kG);\r
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
}\r
\r
printf("\n \n Comment: %s \n \n", comment.Data());\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
\r
//============================================================= \r
// Field (L3 0.5 T) outside dimuon spectrometer\r
- AliMagF *field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);\r
+ AliMagF *field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);\r
TGeoGlobalMagField::Instance()->SetField(field);\r
\r
Int_t iITS = 1;\r
AliMagF* field = 0x0;\r
if (mag == kNoField) {\r
comment = comment.Append(" | L3 field 0.0 T");\r
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k2kG,\r
- "$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",\r
- AliMagF::kBeamTypepp, energy/2.0);\r
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k2kG,AliMagF::kBeamTypepp, energy/2.0);\r
} else if (mag == k5kG) {\r
comment = comment.Append(" | L3 field 0.5 T");\r
- field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG,\r
- "$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",\r
- AliMagF::kBeamTypepp, energy/2.0);\r
+ field = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG,AliMagF::kBeamTypepp, energy/2.0);\r
}\r
printf("\n \n Comment: %s \n \n", comment.Data());\r
\r
static PprRun_t srun = kHijing_per2;
static PprRad_t srad = kGluonRadiation;
static AliMagF::BMap_t smag = AliMagF::k5kG;
+static AliMagF::BeamType_t beamType = AliMagF::kBeamTypeAA;
+static Double_t beamEnergy = 7000.*82./208;
static Int_t sseed = 12345; //Set 0 to use the current time
static PprTrigConf_t strig = kDefaultPbPbTrig; // default pp trigger configuration
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField
+ (new AliMagF("Maps","Maps", -1., -1., smag,beamType,beamEnergy));
rl->CdGAFile();
//
Int_t iABSO = 1;
}
printf("\n \n Comment: %s \n \n", comment.Data());
// to use constant field in the barrel use:
- // AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kGUniform);
- AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., mag);
+ // AliMagF* field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kGUniform);
+ AliMagF* field = new AliMagF("Maps","Maps", -1., -1., mag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
// AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform));
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
}
printf("\n \n Comment: %s \n \n", comment.Data());
// TGeoGlobalMagField::Instance()->SetField(field);
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform,
- "$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",
- AliMagF::kBeamTypepp, energy/2.0);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform,AliMagF::kBeamTypepp, energy/2.0);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG,
- "$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",
- AliMagF::kBeamTypepp, energy/2.0);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG, AliMagF::kBeamTypepp, energy/2.0);
}
printf("\n \n Comment: %s \n \n", comment.Data());
// Field (L3 0.4 T)
- AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., smag);
+ AliMagF* field = new AliMagF("Maps","Maps",-1., -1., smag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
//
//gener->SetVertexSmear(perTrack);
// Field (L3 0.5 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
Int_t iABSO = 1;
Int_t iDIPO = 1;
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., smag));
rl->CdGAFile();
//
Int_t iABSO = 1;
// Field (L3 0.4 T)
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps",-1., -1., smag));
rl->CdGAFile();
//
Int_t iABSO = 1;
AliMagF* field = 0x0;
if (mag == kNoField) {
comment = comment.Append(" | L3 field 0.0 T");
- field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k5kGUniform);
+ field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform);
} else if (mag == k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
- field = new AliMagF("Maps","Maps", 2, -1., -1., 10., AliMagF::k5kG);
+ field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
}
printf("\n \n Comment: %s \n \n", comment.Data());
TGeoGlobalMagField::Instance()->SetField(field);
// Field
- TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, -1., -1., 10., smag));
+ TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., smag));
rl->CdGAFile();
//
}
printf("\n \n Comment: %s \n \n", comment.Data());
- AliMagF* field = new AliMagF("Maps","Maps", 2, -1., -1., 10., mag);
+ AliMagF* field = new AliMagF("Maps","Maps", -1., -1., mag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();