EdbPatCell2 Class Reference

#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)
int	FillCell (TObjArray &pat)
int	FillCombinations (EdbPatCell2 &cell, int ir2[2], TObjArray &arrC1, TObjArray &arrC2, bool doFill=1)
EdbH2 *	FillSelectedH2 ()
void	InitPat (EdbPattern &pat, int nx, int ny)
void	InitPatBin (EdbPattern &pat, float binx, float biny)
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	AddObjectOverlap (float x, float y, float dx, float dy, TObject *obj)
int	CellLim ()
void	Copy (const EdbCell2 &cell)
void	Delete ()
	EdbCell2 ()
	[eNcell] pointers to cells More...
	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 DrawSpectrum")
	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	Fill (int j, int n=1)
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
float	eDTX
float	eDTXlim
float	eDTY
float	eDTYlim
float	eDX
float	eDXlim
float	eDY
float	eDYlim
float	eDZ
float	ePhi
float	eShr
float	eXmarg
float	eYmarg

Additional Inherited Members
Protected Attributes inherited from EdbH2
Float_t	eBin [2]
Float_t	eMax [2]
Float_t	eMin [2]
Int_t	eN [2]
Int_t *	eNC
Int_t	eNcell

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

)

{
  // 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

)

{
  // 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
	)

{
  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
	)

{
  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

eDTX

float EdbPatCell2::eDTX

eDTXlim

float EdbPatCell2::eDTXlim

eDTY

float EdbPatCell2::eDTY

eDTYlim

float EdbPatCell2::eDTYlim

eDX

float EdbPatCell2::eDX

eDXlim

float EdbPatCell2::eDXlim

eDY

float EdbPatCell2::eDY

eDYlim

float EdbPatCell2::eDYlim

eDZ

float EdbPatCell2::eDZ

ePhi

float EdbPatCell2::ePhi

eShr

float EdbPatCell2::eShr

eXmarg

float EdbPatCell2::eXmarg

eYmarg

float EdbPatCell2::eYmarg

---

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

• src/libAlignment/EdbPatCell2.h
• src/libAlignment/EdbPatCell2.cxx