X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FDoVerticesSPD.C;h=616e294109a7f7fc60b71d62ecd36a49244edf5d;hb=7e1b73ddf64949f28d5bcfef55094b31ecba73b8;hp=a1d4920d191e576f4ade43b8913eb64a436d469b;hpb=3714ff25397611dfcc07cf4efffb7aa1bcd22188;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/DoVerticesSPD.C b/ITS/DoVerticesSPD.C index a1d4920d191..616e294109a 100644 --- a/ITS/DoVerticesSPD.C +++ b/ITS/DoVerticesSPD.C @@ -12,9 +12,11 @@ #include "AliTracker.h" #include "AliHeader.h" #include "AliITSLoader.h" +#include "AliITSsegmentationSPD.h" #include "AliVertexerTracks.h" #include "AliCDBManager.h" #include "AliGeomManager.h" +#include "AliGRPManager.h" #include "AliITSDetTypeRec.h" #include "AliITSVertexer3D.h" #include "AliITSVertexerZ.h" @@ -41,6 +43,7 @@ Bool_t DoVerticesSPD(Int_t pileupalgo=1, Int_t optdebug=0){ if (!man->IsDefaultStorageSet()) { printf("Setting a local default storage and run number 0\n"); man->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); + // man->SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd())); man->SetRun(0); } else { @@ -86,20 +89,15 @@ Bool_t DoVerticesSPD(Int_t pileupalgo=1, Int_t optdebug=0){ return kFALSE; } esd->ReadFromTree(tree); - tree->GetEvent(0); - Int_t fieldkG=(Int_t)(TMath::Abs(esd->GetMagneticField())+0.001); - - if(fieldkG==5){ - TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG)); - }else if(fieldkG==2){ - TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k2kG)); - }else{ - TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 0., 1., 10., AliMagF::k5kG)); - } - AliMagF* fld = (AliMagF*)TGeoGlobalMagField::Instance()->GetField(); - printf("Magnetic field set to %f\n",-fld->SolenoidField()); + AliGRPManager * grpMan=new AliGRPManager(); + grpMan->ReadGRPEntry(); + grpMan->SetMagField(); + printf("Magnetic field set to %f T\n",AliTracker::GetBz()/10.); AliITSDetTypeRec* detTypeRec = new AliITSDetTypeRec(); + AliITSsegmentation* seg = new AliITSsegmentationSPD(); + detTypeRec->SetSegmentationModel(0,seg); + Double_t xnom=0.,ynom=0.; AliITSVertexerZ *vertz = new AliITSVertexerZ(xnom,ynom); vertz->Init("default");