#include <EdbPattern.h>

Inheritance diagram for EdbPattern:

[legend]

Collaboration diagram for EdbPattern:

[legend]

Public Member Functions
Int_t	AddPattern (EdbPattern &p)

Int_t	Brick () const

TIndexCell *	Cell () const

	EdbPattern ()
	emulsion side 0/1/2, defined when possible More...

	EdbPattern (EdbPattern &p)

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

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

EdbPattern *	ExtractSubPatternXY (float xmin, float xmax, float ymin, float ymax)

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)

void	Expand (int n)

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

void	ProjectTo (const float dz)

void	RemoveAngMargins (float tx0, float ty0, float dtx, float dty)

void	RemovePosMargins (float x0, float y0, float dx, float dy)

void	Reset ()

void	Set0 ()

void	SetSegmentsDZ (float dz)

void	SetSegmentsFlag (int flag)

void	SetSegmentsPlate (int plate)

void	SetSegmentsW (float w)

void	SetSegmentsZ ()

void	SetX (float x)

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 Float_t	DeltaX ()

virtual Float_t	DeltaY ()

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

	EdbPointsBox2D ()

	EdbPointsBox2D (const EdbPointsBox2D &pb)

virtual const EdbAffine2D *	GetKeep () const

virtual void	GetKeep (EdbAffine2D &aff)

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	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 TH1F *	Xhist ()

virtual Float_t	Xmax () const

virtual Float_t	Xmin () const

virtual TH2F *	XYhist ()

virtual TH1F *	Yhist ()

virtual Float_t	Ymax () const

virtual Float_t	Ymin () const

virtual	~EdbPointsBox2D ()

Public Member Functions inherited from EdbPoint3D
virtual void	Substruct (EdbPoint *p)

virtual void	Test () const

virtual void	TestPoint3D () const

virtual void	Transform (const EdbAffine3D *a)

virtual	~EdbPoint3D ()

Public Member Functions inherited from EdbPoint2D
virtual void	TestPoint2D () const

virtual	~EdbPoint2D ()

Public Member Functions inherited from EdbPoint
virtual	~EdbPoint ()

Private Attributes
TIndexCell *	eCell

Int_t	eFlag

Int_t	eID

Int_t	eNAff

Int_t	ePID

EdbID	eScanID

Int_t	eSide
	main scanning ID for this pattern, defined when possible More...

Float_t	eSigma [4]

Float_t	eStepTX

Float_t	eStepTY

Float_t	eStepX
	associated with eSegments More...

Float_t	eStepY

Constructor & Destructor Documentation

◆ EdbPattern() [1/3]

EdbPattern::EdbPattern ( )

emulsion side 0/1/2, defined when possible

                        : 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

◆ AddPattern()

Int_t EdbPattern::AddPattern ( EdbPattern & p )

 {
   int n=p.N();
   Expand( N()+n );
   for(int i=0; i<n; i++)
     AddSegment( *(p.GetSegment(i)) );
   return n;
 }

◆ Brick()

Int_t EdbPattern::Brick ( ) const

inline

341 {return eScanID.eBrick;}

EdbID::eBrick

Int_t eBrick

Definition: EdbID.h:10

◆ Cell()

TIndexCell* EdbPattern::Cell ( ) const

inline

330 {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;
 }

◆ ExtractSubPatternXY()

EdbPattern * EdbPattern::ExtractSubPatternXY	(	float	xmin,
		float	xmax,
		float	ymin,
		float	ymax
	)

 {
   EdbPattern *pat = new EdbPattern( X(), Y(), Z() );
   EdbSegP *s;
   int nseg = N();
  
   for(int i=0; i<nseg; i++) {
     s = GetSegment(i);
     if(s->X() < xmin || s->X() > xmax || s->Y() < ymin || s->Y() > ymax) 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
	)

TODO: move this cycle into iterator

 {
   // return the array of segments compatible with the
   // prediction segment s with the accuracy of nsig (in sigmas)
  
   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

328 {return eID;}

◆ NAff()

Int_t EdbPattern::NAff ( ) const

inline

323 {return eNAff;}

◆ PID()

int EdbPattern::PID ( ) const

inline

329 {return ePID;}

◆ Plate()

Int_t EdbPattern::Plate ( ) const

inline

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

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

318 {eID=id;}

◆ SetNAff()

void EdbPattern::SetNAff ( int n )

inline

320 {eNAff=n;}

◆ SetPID()

void EdbPattern::SetPID ( int pid )

inline

319 {ePID=pid;}

pid

int pid[1000]

Definition: m2track.cpp:13

◆ SetScanID()

void EdbPattern::SetScanID ( EdbID id )

inline

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

321 {eSide=side;}

◆ SetSigma()

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

inline

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

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

◆ Side()

Int_t EdbPattern::Side ( ) const

inline

342 {return eSide;}

◆ StepTX()

float EdbPattern::StepTX ( ) const

inline

326 {return eStepTX;}

◆ StepTY()

float EdbPattern::StepTY ( ) const

inline

327 {return eStepTY;}

◆ StepX()

float EdbPattern::StepX ( ) const

inline

324 {return eStepX;}

◆ StepY()

float EdbPattern::StepY ( ) const

inline

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

◆ eFlag

Int_t EdbPattern::eFlag

private

◆ eID

Int_t EdbPattern::eID

private

◆ eNAff

Int_t EdbPattern::eNAff

private

◆ ePID

Int_t EdbPattern::ePID

private

◆ eScanID

EdbID EdbPattern::eScanID

private

◆ eSide

Int_t EdbPattern::eSide

private

main scanning ID for this pattern, defined when possible

◆ eSigma

Float_t EdbPattern::eSigma[4]

private

◆ eStepTX

Float_t EdbPattern::eStepTX

private

◆ eStepTY

Float_t EdbPattern::eStepTY

private

◆ eStepX

Float_t EdbPattern::eStepX

private

associated with eSegments

◆ eStepY

Float_t EdbPattern::eStepY

private

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

src/libEdr/EdbPattern.h
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

◆ AddPattern()

◆ Brick()

◆ Cell()

◆ ExtractSubPattern()

◆ ExtractSubPatternXY()

◆ 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