EdbPatCell2 Class Reference

helper class to analyse the single pattern More...

#include <EdbPatCell2.h>

Inheritance diagram for EdbPatCell2:

Collaboration diagram for EdbPatCell2:

Public Member Functions
void	AcceptCorrections (EdbLayer &layer)
void	ApplyCorrections (EdbLayer &layer)
void	ApplyCorrections (EdbPattern &pat)
	EdbPatCell2 ()
int	FillCell (EdbPattern &pat)
	limits should be already defined More...
int	FillCell (TObjArray &pat)
	limits should be already defined More...
int	FillCombinations (EdbPatCell2 &cell, int ir2[2], TObjArray &arrC1, TObjArray &arrC2, bool doFill=1)
EdbH2 *	FillSelectedH2 ()
void	InitPat (EdbPattern &pat, int nx, int ny)
	define limits to fit the pattern More...
void	InitPatBin (EdbPattern &pat, float binx, float biny)
	define limits to fit the pattern More...
void	PrintCorr ()
void	ResetCorr ()
float	TXs (EdbSegP &s)
float	TYs (EdbSegP &s)
float	Xs (EdbSegP &s)
float	XsPos (EdbSegP &s)
float	XsRot (EdbSegP &s)
float	Ys (EdbSegP &s)
float	YsPos (EdbSegP &s)
float	YsRot (EdbSegP &s)
virtual	~EdbPatCell2 ()
Public Member Functions inherited from EdbCell2
bool	AddObject (float v[2], TObject *obj)
bool	AddObject (float x, float y, TObject *obj)
bool	AddObject (int ix, int iy, TObject *obj)
bool	AddObject (int j, TObject *obj)
int	CellLim ()
void	Copy (const EdbCell2 &cell)
void	Delete ()
	EdbCell2 ()
	EdbCell2 (const EdbCell2 &cell)
TObject *	GetObject (float x, float y, int ientr) const
TObject *	GetObject (int ix, int iy, int ientr) const
TObject *	GetObject (int j, int ientr) const
int	InitCell (EdbCell2 &c)
int	InitCell (int maxpercell, int n[2], float min[2], float max[2])
int	InitCell (int nx, float minx, float maxx, int ny, float miny, float maxy, int maxpercell)
void	InitEPC ()
void	PrintStat ()
void	Reset ()
int	SelectObjects (float min[2], float max[2], TObjArray &arr)
int	SelectObjects (int min[2], int max[2], TObjArray &arr)
int	SelectObjects (TObjArray &arr)
int	SelectObjectsC (float v[2], float r, TObjArray &arr)
int	SelectObjectsC (float v[2], int ir[2], TObjArray &arr)
int	SelectObjectsC (int iv[2], int ir[2], TObjArray &arr)
int	SelectObjectsCJ (int j, int ir, TObjArray &arr)
void	Set0 ()
	~EdbCell2 ()
Public Member Functions inherited from EdbH2
void	AddBin (int jcell, int n)
int	Bin (float x, float y) const
int	Bin (int iv[2]) const
int	Bin (int ix, int iy) const
int	Bin (int j) const
void	CleanCells ()
void	Copy (const EdbH2 &h)
void	Delete ()
int	DiscardHighCells (int nmax)
TH2F *	DrawH2 (const char *name="plot2d", const char *title="EdbH2plot2D")
TH1I *	DrawSpectrum (const char *name="plot1d", const char *title="EdbH2 DrawSpectrun")
	EdbH2 ()
	EdbH2 (const EdbH2 &h)
	EdbH2 (int nx, float minx, float maxx, int ny, float miny, float maxy)
int	Fill (float x, float y)
int	Fill (float x, float y, int n)
int	InitH2 (const EdbH2 &h)
int	InitH2 (int n[2], float min[2], float max[2])
int	InitH2 (int nx, float minx, float maxx, int ny, float miny, float maxy)
Long_t	Integral ()
Long_t	Integral (int iv[2], int ir[2])
int	IX (float x) const
int	IX (int jcell) const
int	IY (float y) const
int	IY (int jcell) const
int	Jcell (float v[2]) const
int	Jcell (float x, float y) const
int	Jcell (int ix, int iy) const
int	MaxBin ()
Float_t	Mean ()
int	Ncell () const
int	NX () const
int	NY () const
void	PrintStat ()
EdbH1 *	ProjectionX ()
EdbH1 *	ProjectionY ()
void	Set0 ()
void	SetBin (int ix, int iy, int n)
void	SetBin (int j, int n)
float	X (int ix) const
float	Xbin () const
float	Xj (int j) const
float	Xmax () const
float	XmaxA (float level=0)
float	Xmin () const
float	XminA (float level=0)
float	Y (int iy) const
float	Ybin () const
float	Yj (int j) const
float	Ymax () const
float	YmaxA (float level=0)
float	Ymin () const
float	YminA (float level=0)
	~EdbH2 ()

Public Attributes
bool	eApplyCorr
TH1I *	eDoubletsRate
	to be filled in FillCombinations() More...
float	eDTX
float	eDTXlim
float	eDTY
	corrections to be applied if eApplyCorr==true More...
float	eDTYlim
float	eDX
float	eDXlim
float	eDY
float	eDYlim
	acceptance limits for the combinations selection More...
float	eDZ
	corrections to be applied if eApplyCorr==true More...
float	ePhi
	rotation angle (using dx,dy one can set-up the center of rotation) More...
float	eShr
	corrections to be applied if eApplyCorr==true More...
float	eXmarg
float	eYmarg
	margins for the cell definition More...

Additional Inherited Members
Protected Attributes inherited from EdbH2
Float_t	eBin [2]
	bin size More...
Float_t	eMax [2]
	max More...
Float_t	eMin [2]
	min More...
Int_t	eN [2]
	divisions More...
Int_t *	eNC
	[eNcell] number of objects/cell More...
Int_t	eNcell
	eNx*eNy More...

Detailed Description

helper class to analyse the single pattern

Constructor & Destructor Documentation

EdbPatCell2()

EdbPatCell2::EdbPatCell2

(

)

{
  eDXlim=3; eDYlim=3; eDTXlim=0.003; eDTYlim=0.003;
  eXmarg = eYmarg = 0.0001;        // margins for the cell definition
  eDoubletsRate=0;
  ResetCorr();
}

~EdbPatCell2()

virtual EdbPatCell2::~EdbPatCell2 ( )

inlinevirtual

{ResetCorr();}

Member Function Documentation

AcceptCorrections()

void EdbPatCell2::AcceptCorrections

(

EdbLayer &

layer

)

accept corrections from the layer

{
  eDX =  layer.GetAffineXY()->B1();
  eDY =  layer.GetAffineXY()->B1();
  eDZ =  layer.Zcorr();
  eDTX = layer.GetAffineTXTY()->B1();
  eDTY = layer.GetAffineTXTY()->B2();
  eShr = layer.Shr();
}

ApplyCorrections() [1/2]

void EdbPatCell2::ApplyCorrections

(

EdbLayer &

layer

)

propagate corrections to layer

{
  EdbAffine2D a2( 1,0,0,1, eDX,   eDY);
  EdbAffine2D at2(1,0,0,1, eDTX,  eDTY);
  layer.ApplyCorrections( eShr,  eDZ, a2, at2 );
  Log(3,"EdbPatCell2::ApplyCorrections","");
  if(gEDBDEBUGLEVEL>2)  PrintCorr();
 
}

ApplyCorrections() [2/2]

void EdbPatCell2::ApplyCorrections

(

EdbPattern &

pat

)

apply corrections to the segments of a given pat

{
  pat.TransformShr(eShr);
  EdbAffine2D aff(1,0,0,1, eDX, eDY);
  pat.Transform(&aff);
  EdbAffine2D affa(1,0,0,1, eDTX, eDTY);
  pat.TransformA(&affa);
}

FillCell() [1/2]

int EdbPatCell2::FillCell

(

EdbPattern &

pat

)

limits should be already defined

assume that the cell is already initialized

{
  EdbSegP *s=0;
  int n = pat.N();
  for(int i=0; i<n; i++) {
    s = pat.GetSegment(i);
    AddObject( Xs(*s), Ys(*s), (TObject*)s );
  }
  return eNcell;
}

FillCell() [2/2]

int EdbPatCell2::FillCell

(

TObjArray &

pat

)

limits should be already defined

assume that the cell is already initialized

{
  EdbSegP *s=0;
  int n = pat.GetEntriesFast();
  for(int i=0; i<n; i++) {
    s = (EdbSegP*)pat.UncheckedAt(i);
    AddObject( Xs(*s), Ys(*s), (TObject*)s );
  }
  return eNcell;
}

FillCombinations()

int EdbPatCell2::FillCombinations	(	EdbPatCell2 &	cell,
		int	ir2[2],
		TObjArray &	arrC1,
		TObjArray &	arrC2,
		bool	doFill = 1
	)

{
  int nout=0;
 
  TObjArray arr1(10000);
  TObjArray arr2(100);
  float v[2];
  EdbSegP *s1,*s2;
 
  int ncomb;                      // douplications rate counter
  int n1 = SelectObjects(arr1);
  for(int i=0; i<n1; i++) {
    ncomb=0;
    s1 = (EdbSegP*)arr1.UncheckedAt(i);
    if(s1->Flag()==-10)              continue;
    arr2.Clear();
    v[0] = Xs(*s1);
    v[1] = Ys(*s1);
    int n2 = cell.SelectObjectsC(v,ir2,arr2);
    if(n2<1)                         continue;
    float tx1 = TXs(*s1);
    float ty1 = TYs(*s1);
    for(int i2=0; i2<n2; i2++) {
      s2 = (EdbSegP*)arr2.UncheckedAt(i2);
      if(s2->Flag()==-10)                                    continue;
      if( s2==s1 )                                           continue;
      if( Abs(cell.Xs(*s2) - v[0]) > eDXlim )                continue;
      if( Abs(cell.Ys(*s2) - v[1]) > eDYlim )                continue;
      if( Abs(cell.TXs(*s2)- tx1)  > eDTXlim )               continue;
      if( Abs(cell.TYs(*s2)- ty1)  > eDTYlim )               continue;
 
      if(doFill) {
    arrC1.Add(s1);
    arrC2.Add(s2);
      }
      ncomb++;
    }
    nout+=ncomb;
    if(eDoubletsRate) eDoubletsRate->Fill(ncomb);
  }
  
  Log(3,"EdbPatCell2::FillCombinations","%7d couples found with the acceptance:(dx,dy): %6.2f %6.2f (dtx,dty): %6.4f %6.4f",
      nout,eDXlim,eDYlim,eDTXlim,eDTYlim );
  if(gEDBDEBUGLEVEL>2&&eApplyCorr)             PrintCorr();
  if(gEDBDEBUGLEVEL>2&&cell.eApplyCorr)   cell.PrintCorr();
  return nout;
}

FillSelectedH2()

EdbH2 * EdbPatCell2::FillSelectedH2

(

)

{
  EdbH2 *h = new EdbH2();
  h->InitH2(eN, eMin, eMax);
  EdbSegP *s;
  for(int i=0; i<eNcell; i++)
    for(int j=0; j<eNC[i]; j++) {
      s = (EdbSegP*)(GetObject(i,j));
      if(s->Flag()==-10) continue;
      h->AddBin(i,1);
    }
  return h;
}

InitPat()

void EdbPatCell2::InitPat	(	EdbPattern &	pat,
		int	nx,
		int	ny
	)

define limits to fit the pattern

{
  float min[2] = {pat.Xmin()-eXmarg, pat.Ymin()-eYmarg };
  float max[2] = {pat.Xmax()+eXmarg, pat.Ymax()+eYmarg };
  int   n[2]   = {nx,ny};
  int maxcell = pat.N()/n[0]/n[1]+10;
  maxcell += (int)(5*Sqrt(maxcell));
  Log(2, "EdbPatCell2::InitPat", "maxcell = %d\n",maxcell);
  InitCell(maxcell,n,min,max);
}

InitPatBin()

void EdbPatCell2::InitPatBin	(	EdbPattern &	pat,
		float	binx,
		float	biny
	)

define limits to fit the pattern

{
  float min[2] = {pat.Xmin()-eXmarg, pat.Ymin()-eYmarg };
  float max[2] = {pat.Xmax()+eXmarg, pat.Ymax()+eYmarg };
  int   n[2]   = { (int)((max[0]-min[0])/binx+1), (int)((max[1]-min[1])/biny+1) };
  int maxcell = pat.N()/n[0]/n[1]+10;
  maxcell += (int)(5*Sqrt(maxcell)); 
  Log(2, "EdbPatCell2::InitPatBin", "maxcell = %d\n",maxcell);
  InitCell(maxcell,n,min,max);
}

PrintCorr()

void EdbPatCell2::PrintCorr

(

)

{
  printf("EdbPatCell2 corrections: (dx,dy,dz): %7.2f %7.2f %7.2f (shr): %6.4f (dtx,dty): %6.4f %6.4f\n",
     eDX,eDY,eDZ,eShr,eDTX,eDTY);
}

ResetCorr()

void EdbPatCell2::ResetCorr

(

)

{
  eApplyCorr = false;
  ePhi = eDX = eDY = eDZ = eDTX = eDTY = 0.;
  eShr = 1;
  SafeDelete(eDoubletsRate);
}

TXs()

float EdbPatCell2::TXs ( EdbSegP &  s )

inline

{ return eApplyCorr? s.eTX/eShr+eDTX : s.eTX; }

TYs()

float EdbPatCell2::TYs ( EdbSegP &  s )

inline

{ return eApplyCorr? s.eTY/eShr+eDTY : s.eTY; }

Xs()

float EdbPatCell2::Xs ( EdbSegP &  s )

inline

{ if(!eApplyCorr) return s.eX; if(ePhi<0.00001) return XsPos(s); else return XsRot(s); }

XsPos()

float EdbPatCell2::XsPos ( EdbSegP &  s )

inline

{ return eApplyCorr? eDX + s.eX + TXs(s)*eDZ : s.eX; }

XsRot()

float EdbPatCell2::XsRot ( EdbSegP &  s )

inline

{ return  XsPos(s)*TMath::Cos(ePhi)-YsPos(s)*TMath::Sin(ePhi); }

Ys()

float EdbPatCell2::Ys ( EdbSegP &  s )

inline

{ if(!eApplyCorr) return s.eY; if(ePhi<0.00001) return YsPos(s); else return YsRot(s); }

YsPos()

float EdbPatCell2::YsPos ( EdbSegP &  s )

inline

{ return eApplyCorr? eDY + s.eY + TYs(s)*eDZ : s.eY; }

YsRot()

float EdbPatCell2::YsRot ( EdbSegP &  s )

inline

{ return  XsPos(s)*TMath::Sin(ePhi)+YsPos(s)*TMath::Cos(ePhi); }

Member Data Documentation

eApplyCorr

bool EdbPatCell2::eApplyCorr

eDoubletsRate

TH1I* EdbPatCell2::eDoubletsRate

to be filled in FillCombinations()

eDTX

float EdbPatCell2::eDTX

eDTXlim

float EdbPatCell2::eDTXlim

eDTY

float EdbPatCell2::eDTY

corrections to be applied if eApplyCorr==true

eDTYlim

float EdbPatCell2::eDTYlim

eDX

float EdbPatCell2::eDX

eDXlim

float EdbPatCell2::eDXlim

eDY

float EdbPatCell2::eDY

eDYlim

float EdbPatCell2::eDYlim

acceptance limits for the combinations selection

eDZ

float EdbPatCell2::eDZ

corrections to be applied if eApplyCorr==true

ePhi

float EdbPatCell2::ePhi

rotation angle (using dx,dy one can set-up the center of rotation)

eShr

float EdbPatCell2::eShr

corrections to be applied if eApplyCorr==true

eXmarg

float EdbPatCell2::eXmarg

eYmarg

float EdbPatCell2::eYmarg

margins for the cell definition

---

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

• /home/antonio/fedra_doxygen/src/libAlignment/EdbPatCell2.h
• /home/antonio/fedra_doxygen/src/libAlignment/EdbPatCell2.cxx