[u/mrichter/AliRoot.git] / HLT / TPCLib / Ref / AliHLTTPCVertexFinder.cxx

// @(#) $Id$

// Author: Uli Frankenfeld <mailto:franken@fi.uib.no>
//*-- Copyright &copy ALICE HLT Group

#include "AliHLTTPCStandardIncludes.h"

#include "AliHLTTPCLogging.h"
#include "AliHLTTPCVertexArray.h"
#include "AliHLTTPCVertex.h"
#include "AliHLTTPCVertexFinder.h"
#include "AliHLTTPCSpacePointData.h"
#include "AliHLTTPCTransform.h"

/** \class AliHLTTPCVertexFinder
<pre>
//_____________________________________________________________
// AliHLTTPCVertexFinder
//
//   Implementation of AliHLTTPCArray 
//   usage:
// 
//   ResetSector();
//   for(n=0;n<NMEMSEC;n++)  
//     Read();
//   FindSectorVertex();
//   SetZ(GetZSector());
//   SetZErr(GetZErrSector());
// 
//
</pre>
*/

ClassImp(AliHLTTPCVertexFinder)
AliHLTTPCVertexFinder::AliHLTTPCVertexFinder(){
  //
  // default constructor for the AliHLTTPCVertexFinder class. 
  //

  //Set vertex to zero.
  SetX(0);
  SetY(0);
  SetZ(0);
  SetXErr(1);
  SetYErr(1);
  SetZErr(1);
  fR=0;
  fPhi=0;

  fMWxy = 1.;
    
}

AliHLTTPCVertexFinder::~AliHLTTPCVertexFinder(){
  //
  // destructor
  //
}

void AliHLTTPCVertexFinder::Reset(){
  ResetSector();
}


void AliHLTTPCVertexFinder::Read(Int_t ncluster, AliHLTTPCSpacePointData* hits ){
  //
  //  analyze sector 
  // 
  
  const Int_t seedRow = AliHLTTPCTransform::GetNRows() - 1; 
  const Int_t firstRow = seedRow-32;
  for(Int_t n=0;n<ncluster;n++){
    if(hits[n].fPadRow==seedRow)
      FillSectorSeed3D(hits[n].fX,hits[n].fY,hits[n].fZ);  //copy seeds in 3D
    if(hits[n].fPadRow<=seedRow && hits[n].fPadRow>=firstRow)
      FillSector3D(hits[n].fX,hits[n].fY,hits[n].fZ);  //copy data in 3D
  }
}

void AliHLTTPCVertexFinder::Analyze(){
  FindSectorVertex();
  SetZ(GetZSector());
  SetZErr(GetZSectorErr());
  LOG(AliHLTTPCLog::kInformational,"AliHLTTPCVertexFinder::Analyze","Result")
  <<AliHLTTPCLog::kDec<<"Vertex: "<<GetZ()<<"  RMS: "<<GetZErr()<<ENDLOG;
}

void AliHLTTPCVertexFinder::Write(AliHLTTPCVertex *vertex){
  vertex->SetX(GetX());
  vertex->SetY(GetZ());
  vertex->SetZ(GetZ());
  vertex->SetXErr(GetXErr());
  vertex->SetYErr(GetYErr());
  vertex->SetZErr(GetZErr());

  vertex->SetXYWeight(GetXYWeight());
}

void AliHLTTPCVertexFinder::Write(AliHLTTPCVertexData *vertex){
  vertex->fX=GetX();
  vertex->fY=GetZ();
  vertex->fZ=GetZ();
  vertex->fXErr=GetXErr();
  vertex->fYErr=GetYErr();
  vertex->fZErr=GetZErr();
}
Commit	Line	Data
78001a73	1	// @(#) $Id$
	2
	3	// Author: Uli Frankenfeld <mailto:franken@fi.uib.no>
	4	//*-- Copyright &copy ALICE HLT Group
	5
	6	#include "AliHLTTPCStandardIncludes.h"
	7
	8	#include "AliHLTTPCLogging.h"
	9	#include "AliHLTTPCVertexArray.h"
	10	#include "AliHLTTPCVertex.h"
	11	#include "AliHLTTPCVertexFinder.h"
	12	#include "AliHLTTPCSpacePointData.h"
	13	#include "AliHLTTPCTransform.h"
	14
	15	/** \class AliHLTTPCVertexFinder
	16	<pre>
	17	//_____________________________________________________________
	18	// AliHLTTPCVertexFinder
	19	//
	20	// Implementation of AliHLTTPCArray
	21	// usage:
	22	//
	23	// ResetSector();
	24	// for(n=0;n<NMEMSEC;n++)
	25	// Read();
	26	// FindSectorVertex();
	27	// SetZ(GetZSector());
	28	// SetZErr(GetZErrSector());
	29	//
	30	//
	31	</pre>
	32	*/
	33
	34	ClassImp(AliHLTTPCVertexFinder)
	35	AliHLTTPCVertexFinder::AliHLTTPCVertexFinder(){
	36	//
	37	// default constructor for the AliHLTTPCVertexFinder class.
	38	//
	39
	40	//Set vertex to zero.
	41	SetX(0);
	42	SetY(0);
	43	SetZ(0);
	44	SetXErr(1);
	45	SetYErr(1);
	46	SetZErr(1);
	47	fR=0;
	48	fPhi=0;
	49
	50	fMWxy = 1.;
	51
	52	}
	53
	54	AliHLTTPCVertexFinder::~AliHLTTPCVertexFinder(){
	55	//
	56	// destructor
	57	//
	58	}
	59
	60	void AliHLTTPCVertexFinder::Reset(){
	61	ResetSector();
	62	}
	63
	64
65	void AliHLTTPCVertexFinder::Read(Int_t ncluster, AliHLTTPCSpacePointData* hits ){
66	//
67	// analyze sector
68	//
69
70	const Int_t seedRow = AliHLTTPCTransform::GetNRows() - 1;
71	const Int_t firstRow = seedRow-32;
72	for(Int_t n=0;n<ncluster;n++){
73	if(hits[n].fPadRow==seedRow)
74	FillSectorSeed3D(hits[n].fX,hits[n].fY,hits[n].fZ); //copy seeds in 3D
75	if(hits[n].fPadRow<=seedRow && hits[n].fPadRow>=firstRow)
76	FillSector3D(hits[n].fX,hits[n].fY,hits[n].fZ); //copy data in 3D
77	}
78	}
79
80	void AliHLTTPCVertexFinder::Analyze(){
81	FindSectorVertex();
82	SetZ(GetZSector());
83	SetZErr(GetZSectorErr());
84	LOG(AliHLTTPCLog::kInformational,"AliHLTTPCVertexFinder::Analyze","Result")
85	<<AliHLTTPCLog::kDec<<"Vertex: "<<GetZ()<<" RMS: "<<GetZErr()<<ENDLOG;
86	}
87
88	void AliHLTTPCVertexFinder::Write(AliHLTTPCVertex *vertex){
89	vertex->SetX(GetX());
90	vertex->SetY(GetZ());
91	vertex->SetZ(GetZ());
92	vertex->SetXErr(GetXErr());
93	vertex->SetYErr(GetYErr());
94	vertex->SetZErr(GetZErr());
95
96	vertex->SetXYWeight(GetXYWeight());
97	}
98
99	void AliHLTTPCVertexFinder::Write(AliHLTTPCVertexData *vertex){
100	vertex->fX=GetX();
101	vertex->fY=GetZ();
102	vertex->fZ=GetZ();
103	vertex->fXErr=GetXErr();
104	vertex->fYErr=GetYErr();
105	vertex->fZErr=GetZErr();
106	}