2 #if !defined(__CINT__) || defined(__MAKECINT__)
3 #include "AliCDBManager.h"
4 #include "AliCDBStorage.h"
6 #include "AliCDBMetaData.h"
7 #include "AliESDVertex.h"
14 void UpdateCDBVertexDiamond(Double_t xmed = 0., Double_t ymed = 0., Double_t sigx = 0.005, Double_t sigy = 0.005, Double_t sigz = 5.3) {
15 // produce the mean vertex with the current AliRoot and store it in the
18 AliCDBManager* man = AliCDBManager::Instance();
19 man->SetDefaultStorage("local://$ALICE_ROOT");
21 AliCDBId id("GRP/Calib/MeanVertex",0,AliCDBRunRange::Infinity());
22 AliCDBId idTPC("GRP/Calib/MeanVertexTPC",0,AliCDBRunRange::Infinity());
23 AliCDBMetaData *metadata= new AliCDBMetaData();
25 // Get root and AliRoot versions
26 const char* rootv = gROOT->GetVersion();
27 TString av(ALIROOT_SVN_BRANCH);
28 Int_t revnum = ALIROOT_SVN_REVISION;
30 metadata->SetResponsible("prino@to.infn.it");
31 metadata->SetComment("Default mean vertex position");
32 metadata->SetAliRootVersion(av.Data());
33 metadata->SetComment(Form("Default mean vertex produced with root version %s and AliRoot %s, revision number %d",rootv,av.Data(),revnum));
36 Printf(Form("Storing in CDB the default mean vertex produced with root version %s and"
37 "AliRoot version %s, revision number %d", rootv, av.Data(), revnum));
39 Double_t resolx=35./10000.;
40 Double_t resoly=35./10000.;
41 Double_t sigma[3],position[3];
45 sigma[0]=TMath::Sqrt(sigx*sigx+resolx*resolx);
46 sigma[1]=TMath::Sqrt(sigy*sigy+resoly*resoly);
49 AliESDVertex *vertex = new AliESDVertex(position,sigma,"vtxmean");
50 vertex->PrintStatus();
52 man->Put(vertex,id,metadata);
57 sigma[0]=TMath::Sqrt(sigx*sigx+resolx*resolx);
58 sigma[1]=TMath::Sqrt(sigy*sigy+resoly*resoly);
61 AliESDVertex *vertexTPC = new AliESDVertex(position,sigma,"vtxmean");
62 vertexTPC->PrintStatus();
64 man->Put(vertexTPC,idTPC,metadata);