if (!geom) return kFALSE;
printf("Geometry name: %s \n",(initgeom.GetGeometryName()).Data());
- AliITSDetTypeRec *fDetTypeRec = new AliITSDetTypeRec();
+ fDetTypeRec = new AliITSDetTypeRec();
fDetTypeRec->SetITSgeom(geom);
fDetTypeRec->SetDefaults();
fDetTypeRec->SetDefaultClusterFindersV2(kTRUE);
else {
// Run standard vertexer3d
AliITSVertexer3D *vertexer2 = new AliITSVertexer3D();
- AliESDVertex *vtx = vertexer2->FindVertexForCurrentEvent(clustersTree);
+ vtx = vertexer2->FindVertexForCurrentEvent(clustersTree);
AliMultiplicity *mult = vertexer2->GetMultiplicity();
delete vertexer2;
if(Filter(vtx,mult)) vtxOK = kTRUE;
#ifndef ALIITSMEANVERTEXER_H
#define ALIITSMEANVERTEXER_H
-#include <TString.h>
-
///////////////////////////////////////////////////////////////////////
// //
// Class to compute vertex position using SPD local reconstruction //
/* $Id$ */
-class TObject;
+#include <TObject.h>
+
class TH1F;
class TH2F;
class AliRawReader;
class AliMultiplicity;
class AliESDVertex;
+class AliITSDetTypeRec;
class AliITSMeanVertexer : public TObject {
Bool_t Reconstruct(AliRawReader *rawReader, Bool_t mode = kTRUE);
void WriteVertices(const char *filename);
+ const TH2F*GetVertexXY() const { return fVertexXY; }
+ const TH1F*GetVertexZ() const { return fVertexZ; }
private:
// copy constructor (NO copy allowed)
Int_t fFilterOnContributors; //! Numb. of contrib must be > fFilter...
Int_t fFilterOnTracklets; //! Numb. of tracklets must be > fFilterOnTr...
- ClassDef(AliITSMeanVertexer,0);
+ ClassDef(AliITSMeanVertexer,0)
};
#endif
#include <TPluginManager.h>
#include <TROOT.h>
+#include <TH1.h>
+#include <TH2.h>
#include "AliRawReaderDate.h"
#include "AliCDBManager.h"
+#include "AliITSMeanVertexer.h"
int main(int argc, char **argv) {
#ifdef ALI_AMORE
// send the histos to AMORE pool
- printf("AMORE send status: %d",vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY()));
- printf("AMORE send status: %d",vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ()));
+ printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY()));
+ printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ()));
#endif
delete mv;