EdbAlignment Class Reference

#include <EdbAlignment.h>

Inheritance diagram for EdbAlignment:

Collaboration diagram for EdbAlignment:

Public Member Functions
	EdbAlignment ()
	~EdbAlignment ()

Static Public Member Functions
static int	Calculate (const TObjArray &segarr1, const TObjArray &segarr2, EdbAffine2D &aff)
static int	Calculate (TArrayF &X1, TArrayF &Y1, TArrayF &X2, TArrayF &Y2, EdbAffine2D &aff)
static int	CalculateM (const TObjArray &segarr1, const TObjArray &segarr2, EdbAffine2D &aff, int flag=0)
static float	DTX (const TObjArray &segarr1, const TObjArray &segarr2)
static float	DTY (const TObjArray &segarr1, const TObjArray &segarr2)
static float	DVAR (const TObjArray &segarr1, const TObjArray &segarr2, int ivar)
static float	DX (const TObjArray &segarr1, const TObjArray &segarr2)
static float	DY (const TObjArray &segarr1, const TObjArray &segarr2)
static void	MakeAff (Double_t X0, Double_t Y0, Double_t dX, Double_t dY, Double_t dPhi, EdbAffine2D &aff)

Constructor & Destructor Documentation

EdbAlignment()

EdbAlignment::EdbAlignment ( )

inline

{}

~EdbAlignment()

EdbAlignment::~EdbAlignment ( )

inline

{}

Member Function Documentation

Calculate() [1/2]

int EdbAlignment::Calculate ( const TObjArray &  segarr1, const TObjArray &  segarr2, EdbAffine2D &  aff  )

static

aff applied to pattern 1 gives pattern 2

{
  int n = Min( segarr1.GetEntries(), segarr2.GetEntries() );
  if(n<1) return 0;
  TArrayF X1(n), Y1(n), X2(n), Y2(n);
  for(int i=0; i<n; i++) {
    X1[i] = ((EdbSegP *)segarr1.At(i))->X();
    Y1[i] = ((EdbSegP *)segarr1.At(i))->Y();
    X2[i] = ((EdbSegP *)segarr2.At(i))->X();
    Y2[i] = ((EdbSegP *)segarr2.At(i))->Y();
  }
  return Calculate( X1,Y1, X2,Y2, aff );
}

Calculate() [2/2]

int EdbAlignment::Calculate ( TArrayF &  X1, TArrayF &  Y1, TArrayF &  X2, TArrayF &  Y2, EdbAffine2D &  aff  )

static

this function do not use the matrix inversion for the precision reasons aff applied to pattern 1 gives pattern 2

{
 
  int n = Min( X1.GetSize(), X2.GetSize() );
  if(n<1) return 0;
 
  Double_t dX = 0, dY=0, dPhi=0;
  Double_t X0 = 0, Y0=0;
  
  for(int i=0; i<n; i++) {
    dX += X2[i] - X1[i];
    dY += Y2[i] - Y1[i];
    X0 += X1[i];
    Y0 += Y1[i];
  }
  dX /= n;  dY /= n;  X0 /= n;  Y0 /= n;
 
  if(n<2) return n;
  
  Double_t x1 = 0, y1=0, x2=0, y2=0;
  Double_t R=0, r=0;
  for(int i=0; i<n; i++) {
    x1 = X1[i] -X0;
    y1 = Y1[i] -Y0;
    x2 = X2[i] -X0 -dX;
    y2 = Y2[i] -Y0 -dY;
 
    r = Sqrt(x1*x1+y1*y1);
    dPhi += (ATan2(-y2,-x2) - ATan2(-y1,-x1)) * r;
    R += r;
  }
  dPhi /= R;
  
  MakeAff(X0,Y0, dX,dY, dPhi, aff);
 
  printf("X0,Y0 = %f %f  dX,dY= %f %f dPhi = %f\n",X0,Y0,dX,dY,dPhi);
  return n;
}

CalculateM()

int EdbAlignment::CalculateM ( const TObjArray &  segarr1, const TObjArray &  segarr2, EdbAffine2D &  aff, int  flag = 0  )

static

interface to the matrix procedure

{
  int n = Min( segarr1.GetEntries(), segarr2.GetEntries() );
  if(n<1) return 0;
  TArrayF X1(n), Y1(n), X2(n), Y2(n);
  for(int i=0; i<n; i++) {
    X1[i] = ((EdbSegP *)segarr1.At(i))->X();
    Y1[i] = ((EdbSegP *)segarr1.At(i))->Y();
    X2[i] = ((EdbSegP *)segarr2.At(i))->X();
    Y2[i] = ((EdbSegP *)segarr2.At(i))->Y();
  }
  return aff.Calculate( n, X1.fArray,Y1.fArray, X2.fArray,Y2.fArray , flag);
}

DTX()

static float EdbAlignment::DTX ( const TObjArray &  segarr1, const TObjArray &  segarr2  )

inlinestatic

{ return DVAR(segarr1,segarr2,3); }

DTY()

static float EdbAlignment::DTY ( const TObjArray &  segarr1, const TObjArray &  segarr2  )

inlinestatic

{ return DVAR(segarr1,segarr2,4); }

DVAR()

static float EdbAlignment::DVAR ( const TObjArray &  segarr1, const TObjArray &  segarr2, int  ivar  )

static

DX()

static float EdbAlignment::DX ( const TObjArray &  segarr1, const TObjArray &  segarr2  )

inlinestatic

{ return DVAR(segarr1,segarr2,1); }

DY()

static float EdbAlignment::DY ( const TObjArray &  segarr1, const TObjArray &  segarr2  )

inlinestatic

{ return DVAR(segarr1,segarr2,2); }

MakeAff()

void EdbAlignment::MakeAff ( Double_t  X0, Double_t  Y0, Double_t  dX, Double_t  dY, Double_t  dPhi, EdbAffine2D &  aff  )

static

{
  aff.Reset();
  aff.Rotate(dPhi);
  Double_t x1new = X0*Cos(dPhi)-Y0*Sin(dPhi);
  Double_t y1new = X0*Sin(dPhi)+Y0*Cos(dPhi);
  aff.ShiftX(X0-x1new+dX);
  aff.ShiftY(Y0-y1new+dY);
}

---

The documentation for this class was generated from the following files:

• /home/antonio/fedra_doxygen/src/libEdr/EdbAlignment.h
• /home/antonio/fedra_doxygen/src/libEdr/EdbAlignment.cxx