Additional protection, for the Savannah ticket #68145 (M. Ivanov)

[u/mrichter/AliRoot.git] / ITS / DoVerticesSPD.C

diff --git a/ITS/DoVerticesSPD.C b/ITS/DoVerticesSPD.C
index a1d4920d191e576f4ade43b8913eb64a436d469b..616e294109a7f7fc60b71d62ecd36a49244edf5d 100644 (file)
--- 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");