EdbFragmentAlignment Class Reference

#include <EdbFragmentAlignment.h>

Inheritance diagram for EdbFragmentAlignment:

Collaboration diagram for EdbFragmentAlignment:

Public Member Functions
void	AddPatternAt (EdbPattern *p, int i)
void	AlignAndShift (EdbMosaicPath &mp)
void	AlignFragment (EdbPattern &pf)
void	AlignFragment_dummy (EdbPattern &pf)
void	ApplyAff ()
void	CheckScale (EdbMosaicPath &mp, EdbAffine2D &aff)
float	CheckScaleX (float y0)
float	CheckScaleY (float x0)
	EdbFragmentAlignment ()
void	FillVC ()
void	FillVDT (EdbCouplesTree &vdt)
EdbViewHeader *	GetHeader (int i) const
TObjArray &	GetParr ()
EdbPattern *	GetPattern (int i) const
int	N () const
void	RealignAndShift (EdbMosaicPath &mp)
void	SetAlPar (const AlPar &ap, EdbPlateAlignment &al)
void	SetHarr (TObjArray &ha)
void	SetID (int id)
void	SetMinPeak (int mp)
void	SetPatternsOwner ()
void	SetSide (int s)
int	Side () const
int	ViewSideAl (EdbPattern &p1, EdbPattern &p2, EdbAffine2D &aff, bool do_shift)
virtual	~EdbFragmentAlignment ()

Public Attributes
AlPar	eAP
bool	eDoSaveAlignments
float	eR0

Private Attributes
EdbAffine2D	eAff
TObjArray	eHarr
int	eID
	centers after alignment More...
int	eMinPeak
int	eN
TObjArray	eParr
int	eSide
EdbPattern *	eVC
	original centers More...
EdbPattern *	eVC0

Constructor & Destructor Documentation

EdbFragmentAlignment()

EdbFragmentAlignment::EdbFragmentAlignment ( )

inline

{ eN=0; eVC0=0; eVC=0; }

~EdbFragmentAlignment()

virtual EdbFragmentAlignment::~EdbFragmentAlignment ( )

inlinevirtual

{ eParr.Delete(); SafeDelete(eVC); SafeDelete(eVC0);  }

Member Function Documentation

AddPatternAt()

void EdbFragmentAlignment::AddPatternAt	(	EdbPattern *	p,
		int	i
	)

{ 
  if(i>=0&&i<eParr.GetSize()) eParr.AddAt( p, i );
  else Log(1,"EdbFragmentAlignment::AddPatternAt","Error! i = %d is out of bounds(0:%d)",i,eParr.GetSize());
}

AlignAndShift()

void EdbFragmentAlignment::AlignAndShift

(

EdbMosaicPath &

mp

)

{
  Log(1,"EdbFragmentAlignment::AlignAndShift","");
  mp.SetOK( mp.I(0) );
  for(int i=1; i<eN; i++)
  {
    EdbPattern *p = GetPattern(mp.I(i));
    if(p->N()<eMinPeak)
    {
    Log(1,"EdbFragmentAlignment::AlignAndShift","Skip view with %d segments",p->N());
    continue;
    } 
    TArrayI narr(20);
    int nb = mp.GetAlignedNeighbours( mp.I(i), narr );
    if(nb>20) 
    { 
      Log(1,"EdbFragmentAlignment::AlignAndShift","Warning! too many neigbours: %d  in dr = %f  reset to 20",nb, mp.eR0);
      nb=20;
    }
    EdbPattern  alp;
    for(int ii=0; ii<nb; ii++) {
      alp.AddPattern( *GetPattern(narr[ii]) );
    }
    printf("align %d -> %d  at dist %.1f \n",p->N(), alp.N(), mp.Dist(mp.I(i)) );
    EdbAffine2D aff;
    int peak = ViewSideAl(*p, alp, aff, 1);
    EdbSegP *s = eVC->GetSegment( mp.I(i) );
    s->SetW(peak);
    if( peak > eMinPeak )
    {
      mp.SetOK( mp.I(i) );
      s->Transform(&aff);
    }
  }
}

AlignFragment()

void EdbFragmentAlignment::AlignFragment

(

EdbPattern &

pf

)

{
  EdbMosaicPath mp(eN);
  mp.eR0=eR0;
  EdbViewHeader *h = mp.FindNearest( eHarr, pf.X(), pf.Y() );
  if(h) {
    Log(1,"EdbFragmentAlignment::AlignFragment","with %d views at x0,y0 was %f %f:   %f %f",
    eN,h->GetXview(), h->GetYview(), pf.X(), pf.Y() );
    mp.InitArea( eHarr, h->GetXview(), h->GetYview(), eMinPeak );
  }
  else {
    Log(1,"EdbFragmentAlignment::AlignFragment","Warning! central view close to (%f %f) was not found! abandon this fragment",pf.X(), pf.Y());
    //mp.InitArea( eHarr, pf.X(), pf.Y(), eMinPeak );
    return;
  }
    
  AlignAndShift( mp );  
  RealignAndShift( mp );
  AlignAndShift( mp );
  RealignAndShift( mp );
  
  for(int i=0; i<eN; i++)
  {
    if( mp.OK(i) )   pf.AddPattern( *GetPattern(i) );
  }
  
  eAff.Reset();
  eAff.Calculate(eVC,eVC0);                     // from found to original
  Log(1,"EdbFragmentAlignment::AlignFragment","apply %s",eAff.AsString());
  eVC->Transform(&eAff);
  pf.Transform(&eAff);
}

AlignFragment_dummy()

void EdbFragmentAlignment::AlignFragment_dummy

(

EdbPattern &

pf

)

{
  EdbMosaicPath mp(eN);
  mp.eR0=eR0;
  EdbViewHeader *h = mp.FindNearest( eHarr, pf.X(), pf.Y() );
  if(h) {
    Log(1,"EdbFragmentAlignment::AlignFragment","with %d views at x0,y0 was %f %f:   %f %f",
    eN,h->GetXview(), h->GetYview(), pf.X(), pf.Y() );
    mp.InitArea( eHarr, h->GetXview(), h->GetYview(), eMinPeak );
  }
  else {
    Log(1,"EdbFragmentAlignment::AlignFragment","Warning! central view close to (%f %f) was not found! abandon this fragment",pf.X(), pf.Y());
    //mp.InitArea( eHarr, pf.X(), pf.Y(), eMinPeak );
    return;
  }
    
  //AlignAndShift( mp );  
  //RealignAndShift( mp );
  //AlignAndShift( mp );
  //RealignAndShift( mp );
  
  for(int i=0; i<eN; i++)
  {
    pf.AddPattern( *GetPattern(i) );
  }
  
  //eAff.Reset();
  //eAff.Calculate(eVC,eVC0);                     // from found to original
  Log(1,"EdbFragmentAlignment::AlignFragment_dummy","apply %s",eAff.AsString());
  //eVC->Transform(&eAff);
  //pf.Transform(&eAff);
}

ApplyAff()

void EdbFragmentAlignment::ApplyAff

(

)

{
  Log(1,"EdbFragmentAlignment::ApplyAff","%s",eAff.AsString());
  for( int i=0; i<eN; i++ ) GetPattern(i)->Transform(&eAff);
}

CheckScale()

void EdbFragmentAlignment::CheckScale	(	EdbMosaicPath &	mp,
		EdbAffine2D &	aff
	)

{
  int n = mp.N();
  mp.SetOK( (mp.I(0)) );
  
  for(int i=1; i<n; i++) 
  {
    EdbPattern *p = (EdbPattern *)(eParr.At(mp.I(i)));
    TArrayI narr(10);
    int nb = mp.GetAlignedNeighbours( mp.I(i), narr );
    if(nb>10)
    {
      Log(1,"EdbFragmentAlignment::CheckScale","Warning! too many neigbours: %d  in dr = %f  reset to 10",nb, mp.eR0);
      nb=10;
    }
    EdbPattern  alp;
    for(int ii=0; ii<nb; ii++) {
      alp.AddPattern( *(EdbPattern *)(eParr.At(narr[ii])) );
    }
    printf("align %d -> %d  at dist %.1f \n",p->N(), alp.N(), mp.Dist(mp.I(i)) );
    if( ViewSideAl(*p, alp, aff, 0) > eMinPeak ) mp.SetOK( mp.I(i) );
  }
}

CheckScaleX()

float EdbFragmentAlignment::CheckScaleX

(

float

y0

)

{
  EdbMosaicPath mp(eN);
  float length = mp.InitLineX(eHarr, y0, 100. );
  mp.eR0 = eR0;
  EdbAffine2D aff;
  CheckScale(mp,aff);
  return (length-aff.B1())/length;
}

CheckScaleY()

float EdbFragmentAlignment::CheckScaleY

(

float

x0

)

{
  EdbMosaicPath mp(eN);
  float length = mp.InitLineY(eHarr, x0, 100. );
  mp.eR0 = eR0;
  EdbAffine2D aff;
  CheckScale(mp,aff);
  return (length-aff.B2())/length;
}

FillVC()

void EdbFragmentAlignment::FillVC

(

)

{
  eVC0 = new EdbPattern(0,0,0,eN);
  eVC  = new EdbPattern(0,0,0,eN);
  for( int i=0; i<eN; i++ ) {
    EdbViewHeader *h = GetHeader(i);
    EdbSegP *s0 = eVC0->AddSegment(i, h->GetXview(), h->GetYview(), 0, 0, 0,0 );
    EdbSegP *s  = eVC->AddSegment( i, h->GetXview(), h->GetYview(), 0, 0, 0,0 );
//    int side=0;
//    if(     h->GetNframesTop()==0) side=2;         // 2- bottom
//    else if(h->GetNframesBot()==0) side=1;         // 1- top
//    s0->SetAid( h->GetAreaID(),h->GetViewID(), side);
//    s->SetAid( h->GetAreaID(),h->GetViewID(), side);
    s0->SetAid( h->GetAreaID(),h->GetViewID(), eSide);
    s->SetAid( h->GetAreaID(),h->GetViewID(), eSide);
  }
}

FillVDT()

void EdbFragmentAlignment::FillVDT

(

EdbCouplesTree &

vdt

)

{
  int n = eVC0->N();
  for(int i=0; i<n; i++)
  {
    vdt.Fill( eVC0->GetSegment(i), eVC->GetSegment(i), 0,0, 0,0, eID, eSide );
  }
}

GetHeader()

EdbViewHeader* EdbFragmentAlignment::GetHeader ( int  i ) const

inline

{ return (EdbViewHeader *)(eHarr.At(i)); }

GetParr()

TObjArray& EdbFragmentAlignment::GetParr ( )

inline

{ return eParr; }

GetPattern()

EdbPattern* EdbFragmentAlignment::GetPattern ( int  i ) const

inline

{ return    (EdbPattern *)(eParr.At(i)); }

N()

int EdbFragmentAlignment::N ( ) const

inline

{return eN;}

RealignAndShift()

void EdbFragmentAlignment::RealignAndShift

(

EdbMosaicPath &

mp

)

{
  //assume that most of views were already aligned
  Log(1,"EdbFragmentAlignment::RealignAndShift","");
  for(int i=0; i<eN; i++)
  {
    if(mp.OK(i)) continue;
    EdbPattern *p = GetPattern(i);
    if(p->N()<eMinPeak)
    {
    Log(1,"EdbFragmentAlignment::RealignAndShift","Skip view with %d segments",p->N());
    continue;
    } 
    TArrayI narr(20);
    int nb = mp.GetAlignedNeighbours( i, narr );
    if(nb>20) 
    { 
      Log(1,"EdbFragmentAlignment::AlignAndShift","Warning! too many neigbours: %d  in dr = %f  reset to 20",nb, mp.eR0);
      nb=20;
    }
    EdbPattern  alp;
    for(int ii=0; ii<nb; ii++) {
      alp.AddPattern( *GetPattern(narr[ii]) );
    }
    printf("align %d -> %d  at dist %.1f \n",p->N(), alp.N(), mp.Dist(i) );
    EdbAffine2D aff;
    int peak = ViewSideAl(*p, alp, aff, 1);
    EdbSegP *s = eVC->GetSegment( i );
    s->SetW(peak);
    if( peak > eMinPeak )
    {
      mp.SetOK( i );
      s->Transform(&aff);
    }
  }
}

SetAlPar()

void EdbFragmentAlignment::SetAlPar	(	const AlPar &	ap,
		EdbPlateAlignment &	al
	)

{
  av.eNoScaleRot  = ap.NoScaleRot;         // calculate shift only
  av.eOffsetMax   = ap.OffsetMax;
  av.eDZ          = ap.DZ;
  av.eDPHI        = ap.DPHI;
  av.eDoFine      = ap.DoFine;
  av.eSaveCouples = ap.DoSaveCouples;
  av.SetSigma( ap.SigmaR, ap.SigmaT );
  for(int i=0; i<4; i++) av.eDoublets[i] = ap.Doublets[i];
}

SetHarr()

void EdbFragmentAlignment::SetHarr

(

TObjArray &

ha

)

{
  eN = ha.GetEntries();
  eHarr.Expand(eN);
  for( int i=0; i<eN; i++ ) eHarr.AddAt( ha.At(i),i );
  eParr.Expand(eN);
  printf("*** eN=%d\n",eN);
  FillVC();
  printf("*** eN=%d\n",eN);
}

SetID()

void EdbFragmentAlignment::SetID ( int  id )

inline

{eID=id;}

SetMinPeak()

void EdbFragmentAlignment::SetMinPeak ( int  mp )

inline

{eMinPeak=mp;}

SetPatternsOwner()

void EdbFragmentAlignment::SetPatternsOwner ( )

inline

{ eParr.SetOwner(); }

SetSide()

void EdbFragmentAlignment::SetSide ( int  s )

inline

{eSide=s;}

Side()

int EdbFragmentAlignment::Side ( ) const

inline

{return eSide;}

ViewSideAl()

int EdbFragmentAlignment::ViewSideAl	(	EdbPattern &	p1,
		EdbPattern &	p2,
		EdbAffine2D &	aff,
		bool	do_shift
	)

{
  // transform aff instead of pattern
  EdbPlateAlignment av;
  SetAlPar( eAP, av );
  av.eSaveCouples=0;
  if(eDoSaveAlignments) av.InitOutputFile( Form( "p%.3d/%d_%d.al.vsa.root", p1.ScanID().ePlate, p1.ID(), p2.ID() ) ); 
  av.Align( p1, p2, 0);
  EdbAffine2D *affXY = av.eCorrL[0].GetAffineXY();
  EdbAffine2D *affTXTY = av.eCorrL[0].GetAffineTXTY();
  if(av.eNcoins > eMinPeak ) 
  {
    if(do_shift) p1.Transform( affXY );
    aff.Transform( affXY );
  }
  if(eDoSaveAlignments) av.CloseOutputFile();
  return av.eNcoins;  
}

Member Data Documentation

eAff

EdbAffine2D EdbFragmentAlignment::eAff

private

eAP

AlPar EdbFragmentAlignment::eAP

eDoSaveAlignments

bool EdbFragmentAlignment::eDoSaveAlignments

eHarr

TObjArray EdbFragmentAlignment::eHarr

private

eID

int EdbFragmentAlignment::eID

private

centers after alignment

eMinPeak

int EdbFragmentAlignment::eMinPeak

private

eN

int EdbFragmentAlignment::eN

private

eParr

TObjArray EdbFragmentAlignment::eParr

private

eR0

float EdbFragmentAlignment::eR0

eSide

int EdbFragmentAlignment::eSide

private

eVC

EdbPattern* EdbFragmentAlignment::eVC

private

original centers

eVC0

EdbPattern* EdbFragmentAlignment::eVC0

private

---

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

• src/libMosaic/EdbFragmentAlignment.h
• src/libMosaic/EdbFragmentAlignment.cxx