EdbMarksSet Class Reference

#include <EdbFiducial.h>

Inheritance diagram for EdbMarksSet:

Collaboration diagram for EdbMarksSet:

Public Member Functions
EdbAffine2D *	Abs2Stage () const
Int_t	ChangeMapStringSpacer (TString &str, char oldspacer, char newspacer)
void	DrawA (Option_t *opt="") const
void	DrawS (Option_t *opt="") const
	EdbMarksSet ()
	EdbMarksSet (const EdbMarksSet &ms)
EdbMarksBox *	GetAbsolute () const
Int_t	GetN () const
Int_t	GetNa () const
Int_t	GetNs () const
EdbMarksBox *	GetStage () const
void	Print (Option_t *opt="") const
Int_t	ReadDA (char *file)
Int_t	ReadDB (char *file)
Int_t	ReadDollar (char *file, EdbMarksBox *mbox)
Int_t	ReadMap (char *file, char spacer='_', char shape='S')
Int_t	ReadMap2 (char *file, char spacer='_', char shape='S')
EdbAffine2D *	Stage2Abs () const
Int_t	WriteMap (char *file, char spacer='_', char shape='S', int plate=1)
Int_t	WriteMap2 (char *file, char spacer='_', char shape='S', int plate=1)
virtual	~EdbMarksSet ()

Public Attributes
Long_t	eBrick
Float_t	eXmax
Float_t	eXmin
Float_t	eYmax
Float_t	eYmin

Private Attributes
EdbMarksBox *	eAbsolute
EdbMarksBox *	eStage

Constructor & Destructor Documentation

EdbMarksSet() [1/2]

EdbMarksSet::EdbMarksSet

(

)

{
  eXmin=eXmax=eYmin=eYmax=0;  eBrick=0;
  eAbsolute     = new EdbMarksBox();
  eStage        = new EdbMarksBox();
}

EdbMarksSet() [2/2]

EdbMarksSet::EdbMarksSet

(

const EdbMarksSet &

ms

)

{
  eXmin=ms.eXmin;
  eXmax=ms.eXmax;
  eYmin=ms.eYmin;
  eYmax=ms.eYmax;
  eBrick=ms.eBrick;
  if(ms.GetAbsolute()) eAbsolute = new EdbMarksBox( *(ms.GetAbsolute()) );
  else                 eAbsolute     = new EdbMarksBox();
  if(ms.GetStage())    eStage    = new EdbMarksBox( *(ms.GetStage()) );
  else                 eStage        = new EdbMarksBox();
}

~EdbMarksSet()

EdbMarksSet::~EdbMarksSet ( )

virtual

{
  if(eAbsolute)      delete  eAbsolute;
  if(eStage)         delete  eStage;
}

Member Function Documentation

Abs2Stage()

EdbAffine2D * EdbMarksSet::Abs2Stage

(

)

const

{ 
  EdbAffine2D *aff = new EdbAffine2D();
  aff->Calculate( (EdbPointsBox2D *)GetAbsolute(), (EdbPointsBox2D *)GetStage() );
  return aff;
}

ChangeMapStringSpacer()

Int_t EdbMarksSet::ChangeMapStringSpacer	(	TString &	str,
		char	oldspacer,
		char	newspacer
	)

{
  // Changes the spacer inside the map string with a new character
  for (int is=0;is<str.Length();is++) 
    if(str[is]==oldspacer) str[is]=newspacer;
 
  return(1);
}

DrawA()

void EdbMarksSet::DrawA

(

Option_t *

opt = ""

)

const

{ 
  if(eAbsolute) eAbsolute->Draw( 23, 4, 1. );
}

DrawS()

void EdbMarksSet::DrawS

(

Option_t *

opt = ""

)

const

{ 
  if(eStage) eStage->Draw( 22, 3, 1. );
}

GetAbsolute()

EdbMarksBox* EdbMarksSet::GetAbsolute ( ) const

inline

{ return eAbsolute; }

GetN()

Int_t EdbMarksSet::GetN

(

)

const

{ 
  return TMath::Min( GetNa(), GetNs() );
}

GetNa()

Int_t EdbMarksSet::GetNa

(

)

const

{ 
  if(eAbsolute) return eAbsolute->GetN(); else return 0;
}

GetNs()

Int_t EdbMarksSet::GetNs

(

)

const

{ 
  if(eStage)    return eStage->GetN();    else return 0;
}

GetStage()

EdbMarksBox* EdbMarksSet::GetStage ( ) const

inline

{ return eStage;    }

Print()

void EdbMarksSet::Print

(

Option_t *

opt = ""

)

const

{
  printf("EdbMarksSet: eBrick = %ld\n", eBrick);
  printf("EdbMarksSet: eXmin, eYmin, eXmax, eYmax = %f %f %f %f\n", eXmin, eYmin, eXmax, eYmax);
  printf("EdbMarksSet: Absolute mark set:\n");
  if(eAbsolute)      eAbsolute->Print();
  printf("EdbMarksSet: Stage mark set:\n");
  if(eStage)         eStage->Print();
}

ReadDA()

Int_t EdbMarksSet::ReadDA

(

char *

file

)

{ 
  int     mod,  pln,   prg,  plt,  pos;
  float   sftx, sfty,  gap,  rot;
  float   pnm,  eml,   bas,  shr0,  shr1;
  int     n;
  float   x0,y0,z0, x1,y1,z1;
 
  int nmarks=0;
  int ncols=0;
 
  FILE *fp = fopen( file,"r");
  char            buffer[256]="";
 
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%d %d %d %d %d",  &mod, &pln, &prg, &plt, &pos   );
  if( ncols < 4 ) 
    { printf("ERROR: EdbMarksSet::ReadDA: wrong format\n");    return 0; }
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%f %f %f %f",     &sftx, &sfty,  &gap,  &rot     );
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%f %f %f %f %f",  &pnm, &eml, &bas, &shr0, &shr1 );
  if( ncols < 4 ) 
    { printf("ERROR: EdbMarksSet::ReadDA: wrong format\n");    return 0; }
 
  for(int i=0; i<10000; i++) {
    if (fgets (buffer, 256, fp) == NULL) break;
    ncols = sscanf( buffer,"%d %f %f %f %f %f %f",  
             &n, &x0, &y0, &z0, &x1, &y1, &z1 );
    if( ncols < 7 )                                        break;
    eAbsolute->AddMark(n,x0,y0);
    eStage->AddMark(n,x1,y1);
    nmarks++;
  }
 
  printf("%d marks are read from file %s\n", nmarks, file);
  return nmarks;
}

ReadDB()

Int_t EdbMarksSet::ReadDB

(

char *

file

)

{ 
  return ReadDollar( file, eStage );
}

ReadDollar()

Int_t EdbMarksSet::ReadDollar	(	char *	file,
		EdbMarksBox *	mbox
	)

{ 
  //$B - Fiducial marks (Chorus coordinates): -> NFID
  //
  //    mod pln prd plt   pos
  //     2   5   11  -1                   <- plate line
  //    sftx  sfty  gap  rot
  //     0.    0.    0.   0.              <- position correction line
  //    pnm  eml   bas  shr0  shr1
  //     0.  100.  850.  2.               <- emulsion thikness line
  //     n     x         y     z           
  //     10 -339667. -209560.  0.         <- X-ray marks line
  //     12 -39704.  -209630.  0.
  //     13 -339646. -59477.   0.
  //     15 -39746.  -59533.   0.
 
  int     mod,  pln,   prg,  plt,  pos;
  float   sftx, sfty,  gap,  rot;
  float   pnm,  eml,   bas,  shr0,  shr1;
  int     n;
  float   x,y,z;
 
  int nmarks=0;
  int ncols=0;
 
  FILE *fp = fopen( file,"r");
  char            buffer[256]="";
 
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%d %d %d %d %d",  &mod, &pln, &prg, &plt, &pos   );
  if( ncols < 4 ) 
    { printf("ERROR: EdbMarksSet::ReadDA: wrong format\n");    return 0; }
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%f %f %f %f",     &sftx, &sfty,  &gap,  &rot     );
  if (fgets (buffer, 256, fp) == NULL) return 0;
  ncols = sscanf( buffer,"%f %f %f %f %f",  &pnm, &eml, &bas, &shr0, &shr1 );
  if( ncols < 4 ) 
    { printf("ERROR: EdbMarksSet::ReadDA: wrong format\n");    return 0; }
 
  for(int i=0; i<10000; i++) {
    if (fgets (buffer, 256, fp) == NULL)                  break;
    ncols = sscanf( buffer,"%d %f %f %f",  &n, &x, &y, &z );
    if( ncols < 3 )                                       break;
    mbox->AddMark(n,x,y);
    nmarks++;
  }
 
  printf("%d marks are read from file %s\n", nmarks, file);
  return nmarks;
}

ReadMap()

Int_t EdbMarksSet::ReadMap	(	char *	file,
		char	spacer = '_',
		char	shape = 'S'
	)

{
  // Reads map file and adds informations to eAbsolute
  // If shape is 'S', then the "Front X-Ray Mark Map or Optical Mark Map" format is assumed
  // If shape is 'L', then the "Lateral X-Ray Mark Map" format is assumed
  Int_t nmarks = 0;
  Int_t mark, flag;
  Float_t x, y;
 
  FILE *fp = fopen( file,"r");
 
  TString str_header;
  TString str_mark;
 
  if (shape=='S'||shape=='X') {
    str_header = "mapext:_%ld_%*d_%*d_%*d;_%d_%f_%f_%f_%f";
    str_mark = ";_%d_%f_%f_%*f_%*f_%*d_%*d_%d";
  } else if (shape=='L') {
    str_header = "mapx:_%ld_%*d_%*d_%*d;_%d_%f_%f_%f_%f";
    str_mark = ";_%d_%f_%f_%*f_%*f_%*d_%*d_%d";
  } else {
    Log(1,"EdbMarksSet::ReadMap","ERROR: shape '%c' is not implemented\n", shape); return 0;
  }
 
  ChangeMapStringSpacer(str_header,'_',spacer);
  ChangeMapStringSpacer(str_mark,'_',spacer);
 
  if(!fp)
    { Log(1,"EdbMarksSet::ReadMap","ERROR: can not open file %s\n", file); return 0; }
 
  if ( fscanf(fp, str_header.Data(), &eBrick,&nmarks,&eXmin,&eYmin,&eXmax,&eYmax) != 6 )
    { Log(1,"EdbMarksSet::ReadMap","ERROR: file %s is empty or wrong format\n", file); return 0; }
 
  GetAbsolute()->GetMarks()->Clear("C");
 
  for (int imarks=0;imarks<nmarks;imarks++)
    {
      if (fscanf(fp, str_mark.Data(), &mark, &x, &y, &flag) != 4)
    { Log(1,"EdbMarksSet::ReadMap","ERROR: file %s is empty or wrong format\n", file); return 0; }
      GetAbsolute()->AddMark(mark,x,y,flag);
    }
 
  fclose(fp);
 
  Log(3,"EdbMarksSet::ReadMap","%d marks are read from file %s\n", nmarks, file);
  return nmarks;
}

ReadMap2()

Int_t EdbMarksSet::ReadMap2	(	char *	file,
		char	spacer = '_',
		char	shape = 'S'
	)

{
  // Reads map file and adds informations to eAbsolute and eStage
  // If shape is 'S', then the "Front X-Ray Mark Map or Optical Mark Map" format is assumed
  // If shape is 'L', then the "Lateral X-Ray Mark Map" format is assumed
  Int_t nmarks = 0;
  Int_t mark, flag;
  Float_t xabs, yabs, xst, yst;
 
  FILE *fp = fopen( file,"r");
 
  TString str_header;
  TString str_mark;
 
  if (shape=='S'||shape=='X') {
    str_header = "mapext:_%ld_%*d_%*d_%*d;_%d_%f_%f_%f_%f";
    str_mark = ";_%d_%f_%f_%f_%f_%*d_%*d_%d";
  } else if (shape=='L') {
    str_header = "mapx:_%ld_%*d_%*d_%*d;_%d_%f_%f_%f_%f";
    str_mark = ";_%d_%f_%f_%f_%f_%*d_%*d_%d";
  } else {
    Log(1,"EdbMarksSet::ReadMap2","ERROR: shape '%c' is not implemented\n", shape); return 0;
  }
 
  ChangeMapStringSpacer(str_header,'_',spacer);
  ChangeMapStringSpacer(str_mark,'_',spacer);
 
  if(!fp)
  { Log(1,"EdbMarksSet::ReadMap2","ERROR: can not open file %s\n", file); return 0; }
 
  if ( fscanf(fp, str_header.Data(), &eBrick,&nmarks,&eXmin,&eYmin,&eXmax,&eYmax) != 6 )
  { Log(1,"EdbMarksSet::ReadMap2","ERROR: file %s is empty or wrong format\n", file); return 0; }
 
  GetAbsolute()->GetMarks()->Clear("C");
  GetStage()->GetMarks()->Clear("C");
 
  for (int imarks=0;imarks<nmarks;imarks++)
  {
    if (fscanf(fp, str_mark.Data(), &mark, &xabs, &yabs, &xst, &yst, &flag) != 6)
    { Log(1,"EdbMarksSet::ReadMap2","ERROR: file %s is empty or wrong format\n", file); return 0; }
    GetAbsolute()->AddMark(mark,xabs,yabs,flag);
    GetStage()->AddMark(mark,xst,yst,flag);
  }
 
  fclose(fp);
 
  Log(3,"EdbMarksSet::ReadMap","%d marks are read from file %s\n", nmarks, file);
  return nmarks;
}

Stage2Abs()

EdbAffine2D * EdbMarksSet::Stage2Abs

(

)

const

{
  EdbAffine2D *aff = new EdbAffine2D();
  aff->Calculate( (EdbPointsBox2D *)GetStage(), (EdbPointsBox2D *)GetAbsolute() );
  return aff;
}

WriteMap()

Int_t EdbMarksSet::WriteMap	(	char *	file,
		char	spacer = '_',
		char	shape = 'S',
		int	plate = 1
	)

{
  // Create a map string starting from the informations stored in this class
  // and in the EdbMarksBox object eAbsolute
  // If shape is 'S', then the "Front X-Ray Mark Map or Optical Mark Map" format is assumed
  // If shape is 'L', then the "Lateral X-Ray Mark Map" format is assumed
 
  FILE *fp = fopen( file,"w");
 
  TString str_header;
  TString str_mark;
 
  if (shape=='S'||shape=='X') {
    str_header = "mapext:_%ld_%02d_0_0;_%d_%.4f_%.4f_%.4f_%.4f";
    str_mark = ";_%d_%0.0f_%0.0f_%0.0f_%0.0f_1_1_%d";
  } else if (shape=='L') {
    str_header = "mapx:_%ld_%02d_0_0;_%d_%.4f_%.4f_%.4f_%.4f";
    str_mark = ";_%d_%0.0f_%0.0f_%0.0f_%0.0f_1_1_%d";
  } else {
    Log(1,"EdbMarksSet::ReadMap","ERROR: shape '%c' is not implemented\n", shape); return 0;
  }
 
  ChangeMapStringSpacer(str_header,'_',spacer);
  ChangeMapStringSpacer(str_mark,'_',spacer);
 
  if(!fp)
    { Log(1,"EdbMarksSet::ReadMap","ERROR: can not open file %s\n", file); return 0; }
  
  int nmarks = GetAbsolute()->GetN();
  
  fprintf(fp, str_header.Data(), eBrick,plate,nmarks,eXmin,eYmin,eXmax,eYmax);
  
  for (int imarks=0;imarks<nmarks;imarks++)
    {
      EdbMark *mark = GetAbsolute()->GetMark(imarks);
      fprintf(fp, str_mark.Data(), mark->GetID(),mark->GetX(),mark->GetY(),mark->GetX(),mark->GetY(),mark->Flag());
    }
  
  fclose(fp);
 
  Log(3,"EdbMarksSet::WriteMap","%d marks are written as string map in the file %s\n", nmarks, file);
  return nmarks;
}

WriteMap2()

Int_t EdbMarksSet::WriteMap2	(	char *	file,
		char	spacer = '_',
		char	shape = 'S',
		int	plate = 1
	)

{
  // Create a map string starting from the informations stored in this class: eAbsolute and eStage
  // If shape is 'S', then the "Front X-Ray Mark Map or Optical Mark Map" format is assumed
  // If shape is 'L', then the "Lateral X-Ray Mark Map" format is assumed
 
  FILE *fp = fopen( file,"w");
 
  TString str_header;
  TString str_mark;
 
  if (shape=='S'||shape=='X') {
    str_header = "mapext:_%ld_%02d_0_0;_%d_%.4f_%.4f_%.4f_%.4f";
    str_mark = ";_%d_%0.0f_%0.0f_%0.0f_%0.0f_1_1_%d";
  } else if (shape=='L') {
    str_header = "mapx:_%ld_%02d_0_0;_%d_%.4f_%.4f_%.4f_%.4f";
    str_mark = ";_%d_%0.0f_%0.0f_%0.0f_%0.0f_1_1_%d";
  } else {
    Log(1,"EdbMarksSet::WriteMap2","ERROR: shape '%c' is not implemented\n", shape); return 0;
  }
 
  ChangeMapStringSpacer(str_header,'_',spacer);
  ChangeMapStringSpacer(str_mark,'_',spacer);
 
  if(!fp)
  { Log(1,"EdbMarksSet::WriteMap2","ERROR: can not open file %s\n", file); return 0; }
  
  int nmarksAbs = GetAbsolute()->GetN();
  int nmarksSt  = GetStage()->GetN();
  if(nmarksAbs!=nmarksSt) 
  {
    Log(1,"EdbMarksSet::WriteMap2","ERROR: nmarksAbs != nmarksSt: %d %d use only Abs!\n", nmarksAbs, nmarksSt ); 
    return 0;
  };
  
  fprintf(fp, str_header.Data(), eBrick,plate,nmarksAbs,eXmin,eYmin,eXmax,eYmax);
  
  for (int imarks=0;imarks<nmarksAbs;imarks++)
  {
    EdbMark *markAbs = GetAbsolute()->GetMark(imarks);
    EdbMark *markSt  = 0;
    if(nmarksAbs!=nmarksSt) markSt = markAbs;
    else markSt  = GetStage()->GetMark(imarks);
    fprintf(fp, str_mark.Data(), markAbs->GetID(),markAbs->GetX(),markAbs->GetY(),markSt->GetX(),markSt->GetY(),markAbs->Flag());
  }
  
  fclose(fp);
 
  Log(3,"EdbMarksSet::WriteMap2","%d marks are written as string map in the file %s\n", nmarksAbs, file);
  return nmarksAbs;
}

Member Data Documentation

eAbsolute

EdbMarksBox* EdbMarksSet::eAbsolute

private

eBrick

Long_t EdbMarksSet::eBrick

eStage

EdbMarksBox* EdbMarksSet::eStage

private

eXmax

Float_t EdbMarksSet::eXmax

eXmin

Float_t EdbMarksSet::eXmin

eYmax

Float_t EdbMarksSet::eYmax

eYmin

Float_t EdbMarksSet::eYmin

---

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

• src/libEdb/EdbFiducial.h
• src/libEdb/EdbFiducial.cxx