#include <EdbPattern.h>

Inheritance diagram for EdbPattern:

[legend]

Collaboration diagram for EdbPattern:

[legend]

Public Member Functions
TIndexCell *	Cell () const

	EdbPattern ()

	EdbPattern (EdbPattern &p)

	EdbPattern (float x0, float y0, float z0, int n=0)

EdbPattern *	ExtractSubPattern (float min[5], float max[5], int MCevt=-1)

void	FillCell (float stepx, float stepy, float steptx, float stepty)

int	FindCompliments (EdbSegP &s, TObjArray &arr, float nsig, float nsigt)

EdbSegP *	FindSegment (int id)

int	ID () const

Int_t	NAff () const

int	PID () const

Int_t	Plate () const

void	Reset ()

EdbID	ScanID () const

void	Set0 ()

void	SetID (int id)

void	SetNAff (int n)

void	SetPID (int pid)

void	SetScanID (EdbID id)

void	SetSegmentsPID ()

void	SetSegmentsScanID (EdbID id)

void	SetSide (int side)

void	SetSigma (float sx, float sy, float stx, float sty)

void	SetStep (float stepx, float stepy, float steptx, float stepty)

Int_t	Side () const

float	StepTX () const

float	StepTY () const

float	StepX () const

float	StepY () const

float	SummaryPath ()

float	Xmean ()

float	Ymean ()

virtual	~EdbPattern ()

Public Member Functions inherited from EdbSegmentsBox
EdbSegP *	AddSegment (EdbSegP &s)

EdbSegP *	AddSegment (EdbSegP &s1, EdbSegP &s2)

EdbSegP *	AddSegment (int i, EdbSegP &s)

EdbSegP *	AddSegment (int id, float x, float y, float tx, float ty, float w=0, int flag=0)

EdbSegP *	AddSegmentNoDuplicate (EdbSegP &s)

EdbPoint *	At (int i) const

Int_t	CalculateAXAY (EdbSegmentsBox p, EdbAffine2D affA)

Int_t	CalculateXY (EdbSegmentsBox p, EdbAffine2D aff)

Float_t	Diff (EdbSegmentsBox &p)

Float_t	DiffAff (EdbAffine2D *aff)

Float_t	DZ () const

	EdbSegmentsBox (EdbSegmentsBox &box)

	EdbSegmentsBox (float x0, float y0, float z0, int nseg=0)

	EdbSegmentsBox (int nseg=0)

Int_t	GetN () const

EdbSegP *	GetSegment (int i) const

EdbSegP *	GetSegmentLast () const

TClonesArray *	GetSegments () const

void *	GetSegmentsAddr ()

Float_t	GetSize (Int_t XorY)

Float_t	GetSizeX ()

Float_t	GetSizeXY ()

Float_t	GetSizeY ()

Float_t	GetTrackDensity ()

Float_t	GetTrackDensitymm2 ()

Int_t	N () const

void	Print (Option_t *opt="") const
	other functions More...

void	ProjectTo (const float dz)

void	Reset ()

void	Set0 ()

void	SetSegmentsDZ (float dz)

void	SetSegmentsPlate (int plate)

void	SetSegmentsZ ()

void	SetX (float x)
	mandatory virtual functions: More...

void	SetY (float y)

void	SetZ (float z)

void	TransformA (const EdbAffine2D *affA)

void	TransformARot (const EdbAffine2D *affA)

void	TransformShr (const float shr)

Float_t	X () const

Float_t	Y () const

Float_t	Z () const

virtual	~EdbSegmentsBox ()

Public Member Functions inherited from EdbPointsBox2D
virtual EdbPoint *	At (int i) const =0

virtual Float_t	DeltaX ()

virtual Float_t	DeltaY ()

virtual void	DrawPoints (int style=23, int col=4, float size=1.)

	EdbPointsBox2D ()

	EdbPointsBox2D (const EdbPointsBox2D &pb)

virtual const EdbAffine2D *	GetKeep () const

virtual void	GetKeep (EdbAffine2D &aff)

virtual Int_t	N () const =0

virtual void	Print (Option_t *opt="") const

virtual void	Retransform ()

virtual void	Rotate (float angle)

virtual void	ScaleX (float scaleFactor)

virtual void	ScaleY (float scaleFactor)

virtual void	SetKeep (float a11, float a12, float a21, float a22, float b1, float b2)

virtual void	SetX (float x)

virtual void	SetY (float y)

virtual void	SetZ (float z)

virtual void	ShiftX (float offset)

virtual void	ShiftY (float offset)

virtual void	SmearXY (float sigmax, float sigmay)

virtual void	Substruct (EdbPointsBox2D *b)

virtual void	Transform (const EdbAffine2D *a)

virtual Float_t	X () const

virtual TH1F *	Xhist ()

virtual Float_t	Xmax () const

virtual Float_t	Xmin () const

virtual TH2F *	XYhist ()

virtual Float_t	Y () const

virtual TH1F *	Yhist ()

virtual Float_t	Ymax () const

virtual Float_t	Ymin () const

virtual Float_t	Z () const

virtual	~EdbPointsBox2D ()

Public Member Functions inherited from EdbPoint3D
virtual void	Print (Option_t *opt="") const

virtual void	SetZ (float z)=0

virtual void	Substruct (EdbPoint *p)

virtual void	Test () const

virtual void	TestPoint3D () const

virtual void	Transform (const EdbAffine3D *a)

virtual Float_t	Z () const =0

virtual	~EdbPoint3D ()

Public Member Functions inherited from EdbPoint2D
virtual void	Print (Option_t *opt="") const

virtual void	SetX (float x)=0

virtual void	SetY (float y)=0

virtual void	SetZ (float z)

virtual void	Substruct (EdbPoint *p)

virtual void	Test () const

virtual void	TestPoint2D () const

virtual void	Transform (const EdbAffine2D *a)

virtual Float_t	X () const =0

virtual Float_t	Y () const =0

virtual Float_t	Z () const

virtual	~EdbPoint2D ()

Public Member Functions inherited from EdbPoint
virtual void	SetX (float x)=0

virtual void	SetY (float y)=0

virtual void	SetZ (float z)=0

virtual void	Substruct (EdbPoint *p)=0

virtual void	Test () const

virtual void	Transform (const EdbAffine2D *a)

virtual void	Transform (const EdbAffine3D *a)

virtual Float_t	X () const =0

virtual Float_t	Y () const =0

virtual Float_t	Z () const =0

virtual	~EdbPoint ()

Private Attributes
TIndexCell *	eCell
	associated with eSegments More...

Int_t	eFlag
	pattern flag More...

Int_t	eID
	pattern id in the volume More...

Int_t	eNAff
	number of segments selected for affine calculation More...

Int_t	ePID
	correspond to the piece ID More...

EdbID	eScanID
	main scanning ID for this pattern, defined when possible More...

Int_t	eSide
	emulsion side 0/1/2, defined when possible More...

Float_t	eSigma [4]
	accuracy in comparison to neibour pattern More...

Float_t	eStepTX

Float_t	eStepTY
	bin size for the index cell More...

Float_t	eStepX

Float_t	eStepY
	bin size for the index cell More...

Constructor & Destructor Documentation

◆ EdbPattern() [1/3]

EdbPattern::EdbPattern ( )

                       : EdbSegmentsBox()
{
  eCell     = new TIndexCell();
  Set0();
}

◆ EdbPattern() [2/3]

EdbPattern::EdbPattern ( EdbPattern & p )

                                   :EdbSegmentsBox(p)
{
  eID = p.eID;
  ePID = p.ePID;
  eSide = p.eSide;
  eScanID = p.eScanID;
  eFlag = p.eFlag;
  eNAff = p.eNAff;
  eStepX = p.eStepX;  eStepY = p.eStepY;
  eStepTX = p.eStepTX;  eStepTY = p.eStepTY;
  for(int i=0; i<4; i++) eSigma[i] = p.eSigma[i];
  if(p.eCell) eCell = new TIndexCell(*(p.eCell));
}

◆ EdbPattern() [3/3]

EdbPattern::EdbPattern	(	float	x0,
		float	y0,
		float	z0,
		int	n = `0`
	)

                                                          : EdbSegmentsBox(x0,y0,z0,n) 
{
  Set0();
  eCell     = new TIndexCell();
}

◆ ~EdbPattern()

EdbPattern::~EdbPattern ( )

virtual

{
  if(eCell)  { delete eCell;  eCell=0; }
}

Member Function Documentation

◆ Cell()

TIndexCell * EdbPattern::Cell ( ) const

inline

321{return eCell;}

◆ ExtractSubPattern()

EdbPattern * EdbPattern::ExtractSubPattern	(	float	min[5],
		float	max[5],
		int	MCevt = `-1`
	)

return the copy of selected segments

{
  EdbPattern *pat = new EdbPattern( X(), Y(), Z() );
  EdbSegP *s;
  int nseg = N();
 
  for(int i=0; i<nseg; i++) {
    s = GetSegment(i);
        
    if (s->MCEvt()!=MCEvt && s->MCEvt()>0 && MCEvt>=0) continue;
    // Do not continue not in case MCEvt was not specified at all.
    // This allows backward compability.
    
    if(s->X()  < min[0])   continue;
    if(s->Y()  < min[1])   continue;
    if(s->TX() < min[2])   continue;
    if(s->TY() < min[3])   continue;
    if(s->W()  < min[4])   continue;
 
    if(s->X()  > max[0])   continue;
    if(s->Y()  > max[1])   continue;
    if(s->TX() > max[2])   continue;
    if(s->TY() > max[3])   continue;
    if(s->W()  > max[4])   continue;
 
    pat->AddSegment(*s);
  }
  
  // Set ID() and PID() to have consistent values of this pattern:
  pat->SetID(this->ID());
  pat->SetPID(this->PID());
  
  return pat;
}

◆ FillCell()

void EdbPattern::FillCell	(	float	stepx,
		float	stepy,
		float	steptx,
		float	stepty
	)

fill cells with fixed size at z=zPat

{
 
  TIndexCell *cell = Cell();
  if(cell) cell->Drop();
  SetStep(stepx,stepy,steptx,stepty);
 
  float x,y,tx,ty,dz;
  Long_t  val[5];  // x,y,ax,ay,i
  EdbSegP *p;
  int npat = N();
  for(int i=0; i<npat; i++ ) {
    p = GetSegment(i);
    dz = Z() - p->Z();
    tx = p->TX();
    ty = p->TY();
    x  = p->X() + tx*dz;
    y  = p->Y() + ty*dz;
    val[0]= (Long_t)( x / stepx  );
    val[1]= (Long_t)( y / stepy  );
    val[2]= (Long_t)( tx/ steptx );
    val[3]= (Long_t)( ty/ stepty );
    val[4]= (Long_t)(i);
    cell->Add(5,val);
  }
  cell->Sort();
 
}

◆ FindCompliments()

int EdbPattern::FindCompliments	(	EdbSegP &	s,
		TObjArray &	arr,
		float	nsig,
		float	nsigt
	)

return the array of segments compatible with the prediction segment s with the accuracy of nsig (in sigmas)

TODO: move this cycle into iterator

{
 
  float dz = Z()-s.Z();
  float x = s.X() + s.TX()*dz;
  float y = s.Y() + s.TY()*dz;
  float sx = TMath::Sqrt( s.SX() + s.STX()*dz*dz );
  float sy = TMath::Sqrt( s.SY() + s.STY()*dz*dz );
 
  float stx = s.STX();
  if (stx <= 0.) stx = 1.;
  float sty = s.STY();
  if (sty <= 0.) sty = 1.;
 
//   printf("Step: %f %f %f %f\n",StepX(),StepY(),StepTX(),StepTY());
//   printf("Sigm: %f %f %f %f\n",sx,sy,TMath::Sqrt(stx),TMath::Sqrt(sty));
 
 
  long vcent[4] = { (long)(x/StepX()),
            (long)(y/StepY()),
            (long)(s.TX()/StepTX()),
            (long)(s.TY()/StepTY())  };
  long vdiff[4] = { (long)(sx*nsigx/StepX()+1),
            (long)(sy*nsigx/StepY()+1),
            (long)(TMath::Sqrt(stx)*nsigt/StepTX()+1),
            (long)(TMath::Sqrt(sty)*nsigt/StepTY()+1) };
 
  sy*=sy;
  sx*=sx;
 
  float dx,dy,dtx,dty;
 
  int nseg=0;
  EdbSegP *seg=0;
 
 
  //printf("find compliments in pattern %d\n",ID());
  //s.Print();
 
  long vmin[4],vmax[4];
  for(int i=0; i<4; i++) {
    vmin[i] = vcent[i]-vdiff[i];
    vmax[i] = vcent[i]+vdiff[i];
    //printf("%f +- %f \t",(float)vcent[i],(float)vdiff[i]);
  }
  //printf("\n");
 
 
  TIndexCell *c1=0,*c2=0,*c3=0,*c4=0;
  for(vcent[0]=vmin[0]; vcent[0]<=vmax[0]; vcent[0]++) {
    c1 = eCell->Find(vcent[0]);
    if(!c1) continue;
    for(vcent[1]=vmin[1]; vcent[1]<=vmax[1]; vcent[1]++) {
      c2 = c1->Find(vcent[1]);
      if(!c2) continue;
      for(vcent[2]=vmin[2]; vcent[2]<=vmax[2]; vcent[2]++) {
    c3 = c2->Find(vcent[2]);
    if(!c3) continue;
    for(vcent[3]=vmin[3]; vcent[3]<=vmax[3]; vcent[3]++) {
      c4 = c3->Find(vcent[3]);
      if(!c4) continue;
 
      for(int i=0; i<c4->N(); i++) {
        seg = GetSegment(c4->At(i)->Value());
        dtx=s.TX()-seg->TX();
        if( dtx*dtx > stx*nsigt*nsigt )    continue;
        dty=s.TY()-seg->TY();
        if( dty*dty > sty*nsigt*nsigt )    continue;
        dz = seg->Z()-s.Z();
        dx=s.X()+s.TX()*dz-seg->X();
        if( dx*dx > sx*nsigx*nsigx )           continue;
        dy=s.Y()+s.TY()*dz-seg->Y();
        if( dy*dy > sy*nsigx*nsigx )           continue;
        //if(!s.IsCompatible(*seg,nsigx,nsigt)) continue;
        arr.Add(seg);
        nseg++;
      }
 
    }
      }
    }
  }
 
 
//    TIndexCellIterV itr( eCell, 4, vcent, vdiff );
//    const TIndexCell *c=0;
//    while( (c=itr.Next()) ) 
//      for(int i=0; i<c->N(); i++) {
//        seg = GetSegment(c->At(i)->Value());
//        if(!s.IsCompatible(*seg,nsigx,nsigt)) continue;
//        arr.Add(seg);
//        nseg++;
//      }
 
  return nseg;
}

◆ FindSegment()

EdbSegP * EdbPattern::FindSegment ( int id )

{
  for(int i=0; i<N(); i++) if(GetSegment(i)->ID()==id) return GetSegment(i);
  return 0;
}

◆ ID()

int EdbPattern::ID ( ) const

inline

319{return eID;}

◆ NAff()

Int_t EdbPattern::NAff ( ) const

inline

314{return eNAff;}

◆ PID()

int EdbPattern::PID ( ) const

inline

320{return ePID;}

◆ Plate()

Int_t EdbPattern::Plate ( ) const

inline

327{return eScanID.ePlate;}

EdbID::ePlate

Int_t ePlate

Definition: EdbID.h:11

◆ Reset()

void EdbPattern::Reset ( )

{
  ((EdbSegmentsBox *)this)->Reset();
  Set0();
  if(eCell)     eCell->Drop();
}

◆ ScanID()

EdbID EdbPattern::ScanID ( ) const

inline

326{return eScanID;}

◆ Set0()

void EdbPattern::Set0 ( )

{
  eID   = 0;
  ePID  = 0;
  eFlag = 0;
  eNAff = 0;
  eStepX=eStepY=eStepTX=eStepTY=0;
  for(int i=0; i<4; i++) eSigma[i]=0;
  eSide=0;
}

◆ SetID()

void EdbPattern::SetID ( int id )

inline

309{eID=id;}

◆ SetNAff()

void EdbPattern::SetNAff ( int n )

inline

311{eNAff=n;}

◆ SetPID()

void EdbPattern::SetPID ( int pid )

inline

310{ePID=pid;}

pid

int pid[1000]

Definition: m2track.cpp:13

◆ SetScanID()

void EdbPattern::SetScanID ( EdbID id )

inline

296{eScanID=id;}

◆ SetSegmentsPID()

void EdbPattern::SetSegmentsPID ( )

{
  int nseg = N();
  for(int i=0; i<nseg; i++) {
    GetSegment(i)->SetPID(ID());  //PID of the segment must be ID of the pattern!
  }
}

◆ SetSegmentsScanID()

void EdbPattern::SetSegmentsScanID ( EdbID id )

{
  for(int i=0; i<N(); i++)     GetSegment(i)->SetScanID(id);
}

◆ SetSide()

void EdbPattern::SetSide ( int side )

inline

312{eSide=side;}

◆ SetSigma()

void EdbPattern::SetSigma	(	float	sx,
		float	sy,
		float	stx,
		float	sty
	)

inline

298 {eSigma[0]=sx; eSigma[1]=sy; eSigma[2]=stx; eSigma[3]=sty;}

◆ SetStep()

void EdbPattern::SetStep	(	float	stepx,
		float	stepy,
		float	steptx,
		float	stepty
	)

inline

300 {eStepX=stepx; eStepY=stepy; eStepTX=steptx; eStepTY=stepty;}

◆ Side()

Int_t EdbPattern::Side ( ) const

inline

328{return eSide;}

◆ StepTX()

float EdbPattern::StepTX ( ) const

inline

317{return eStepTX;}

◆ StepTY()

float EdbPattern::StepTY ( ) const

inline

318{return eStepTY;}

◆ StepX()

float EdbPattern::StepX ( ) const

inline

315{return eStepX;}

◆ StepY()

float EdbPattern::StepY ( ) const

inline

316{return eStepY;}

◆ SummaryPath()

float EdbPattern::SummaryPath ( )

calculate the microscope path of the eventual prediction scan of this pattern

{
  float sum=0,dx,dy;
  if(N()<2) return sum;
  for(int i=0; i<N()-1; i++ ) {
    dx = GetSegment(i+1)->X()-GetSegment(i)->X();
    dy = GetSegment(i+1)->Y()-GetSegment(i)->Y();
    sum += Sqrt(dx*dx+dy*dy);
  }
  return sum;
}

◆ Xmean()

float EdbPattern::Xmean ( )

{
  Double_t s=0;
  int n=N(); if(n<1) return 0;
  for(int i=0; i<n; i++) s += GetSegment(i)->eX;
  s /= n;
  return (float)s;
}

◆ Ymean()

float EdbPattern::Ymean ( )

{
  Double_t s=0;
  int n=N(); if(n<1) return 0;
  for(int i=0; i<n; i++) s += GetSegment(i)->eY;
  s /= n;
  return (float)s;
}

Member Data Documentation

◆ eCell

TIndexCell* EdbPattern::eCell

private

associated with eSegments

◆ eFlag

Int_t EdbPattern::eFlag

private

pattern flag

◆ eID

Int_t EdbPattern::eID

private

pattern id in the volume

◆ eNAff

Int_t EdbPattern::eNAff

private

number of segments selected for affine calculation

◆ ePID

Int_t EdbPattern::ePID

private

correspond to the piece ID

◆ eScanID

EdbID EdbPattern::eScanID

private

main scanning ID for this pattern, defined when possible

◆ eSide

Int_t EdbPattern::eSide

private

emulsion side 0/1/2, defined when possible

◆ eSigma

Float_t EdbPattern::eSigma[4]

private

accuracy in comparison to neibour pattern

◆ eStepTX

Float_t EdbPattern::eStepTX

private

◆ eStepTY

Float_t EdbPattern::eStepTY

private

bin size for the index cell

◆ eStepX

Float_t EdbPattern::eStepX

private

◆ eStepY

Float_t EdbPattern::eStepY

private

bin size for the index cell

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

/home/antonio/fedra_doxygen/src/libEdr/EdbPattern.h
/home/antonio/fedra_doxygen/src/libEdr/EdbPattern.cxx

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ EdbPattern() [1/3]

◆ EdbPattern() [2/3]

◆ EdbPattern() [3/3]

◆ ~EdbPattern()

Member Function Documentation

◆ Cell()

◆ ExtractSubPattern()

◆ FillCell()

◆ FindCompliments()

◆ FindSegment()

◆ ID()

◆ NAff()

◆ PID()

◆ Plate()

◆ Reset()

◆ ScanID()

◆ Set0()

◆ SetID()

◆ SetNAff()

◆ SetPID()

◆ SetScanID()

◆ SetSegmentsPID()

◆ SetSegmentsScanID()

◆ SetSide()

◆ SetSigma()

◆ SetStep()

◆ Side()

◆ StepTX()

◆ StepTY()

◆ StepX()

◆ StepY()

◆ SummaryPath()

◆ Xmean()

◆ Ymean()

Member Data Documentation

◆ eCell

◆ eFlag

◆ eID

◆ eNAff

◆ ePID

◆ eScanID

◆ eSide

◆ eSigma

◆ eStepTX

◆ eStepTY

◆ eStepX

◆ eStepY