EdbScanSet Class Reference

#include <EdbScanSet.h>

Inheritance diagram for EdbScanSet:

Collaboration diagram for EdbScanSet:

Public Member Functions
Bool_t	AddID (EdbID *id, Int_t step)
Int_t	AssembleBrickFromPC ()
EdbBrickP &	Brick ()
const EdbBrickP &	Brick () const
void	Copy (EdbScanSet &sc)
	EdbScanSet (EdbID id)
	EdbScanSet (int brickid=0)
EdbID *	FindNextPlateID (Int_t p, Bool_t direction)
EdbID *	FindPlateID (Int_t p)
Bool_t	GetAffP2P (Int_t p1, Int_t p2, EdbAffine2D &aff)
Float_t	GetDZP2P (Int_t p1, Int_t p2)
EdbID *	GetID (Int_t i)
EdbLayer *	GetLayer (Int_t p, Int_t side)
EdbPlateP *	GetPlate (Int_t p)
void	MakeNominalSet (EdbID id, int from_plate=57, int to_plate=1, float z0=0, float dz=-1300, float shr=1, float dzbase=210., float dzem=45.)
Int_t	MakeParFiles (Int_t piece=0, const char *dir=".")
void	MakePIDList ()
void	Print ()
void	ReadGeom (const char *fname)
Int_t	ReadIDS (const char *file)
void	RemovePlate (int pid)
Bool_t	SetAsReferencePlate (Int_t pid)
void	SetID (EdbID id)
Int_t	ShiftBrickZ (Float_t z)
Int_t	TransformBrick (EdbAffine2D aff)
void	TransformBrick (EdbScanSet &ss)
Int_t	TransformSidesIntoBrickRS ()
void	UpdateBrickWithP2P (EdbLayer &la, int plate1, int plate2)
void	UpdateGap (float dz, int pid1, int pid2)
void	UpdateIDS (int brick, int ma, int mi)
bool	ValidSide (int side) const
void	WriteGeom (const char *fname)
Int_t	WriteIDS (const char *file=0)
float	Zlayer (int plate, int side)
virtual	~EdbScanSet ()

Public Attributes
EdbBrickP	eB
	all layers of the brick defined here More...
EdbID	eID
	id of the scanset itself More...
TList	eIDS
	list of the identifiers to be processed More...
TObjArray	ePC
	Plate Couples. Each couple represented as EdbPlateP where 2 sides in reality corresponds to 2 plates. More...
TIndexCell	ePID
	correspondance between index in eB and the plate id More...
Int_t	eReferencePlate

Constructor & Destructor Documentation

EdbScanSet() [1/2]

EdbScanSet::EdbScanSet

(

EdbID

id

)

EdbScanSet() [2/2]

EdbScanSet::EdbScanSet

(

int

brickid = 0

)

{
  eB.SetID(id);
  eID.eBrick=id;
  eReferencePlate=0;
}

~EdbScanSet()

virtual EdbScanSet::~EdbScanSet ( )

inlinevirtual

{}

Member Function Documentation

AddID()

Bool_t EdbScanSet::AddID	(	EdbID *	id,
		Int_t	step
	)

{
  if (FindPlateID(id->GetPlate())) return kFALSE;
  for (EdbID *it = (EdbID*)eIDS.First(); it ; it = (EdbID*)eIDS.After(it)) {
    if (step < 0 && id->GetPlate() > it->GetPlate()) {
      eIDS.AddBefore(it, id);
      return kTRUE;
    }
    if (step > 0 && id->GetPlate() < it->GetPlate()) {
      eIDS.AddBefore(it, id);
      return kTRUE;
    }
  }
  eIDS.Add(id);
  return kTRUE;
}

AssembleBrickFromPC()

Int_t EdbScanSet::AssembleBrickFromPC

(

)

input: ePC - couples of layers represented as "plates" output: brick with a first plate as a reference one

{
 
  ePID.Delete();
  eB.Clear();
  EdbPlateP *plate=0;
 
  Float_t dz0=214,dz1=45,dz2=45;  //TODO!
  EdbAffine2D aff;
  Int_t npc = ePC.GetEntries();      if(npc<1) return 0;
  EdbPlateP *pc=0;                    //the couple of plates
  Float_t z;
  printf("\nAssembleBrickFromPC\n");
  for(Int_t i=-1; i<npc; i++) {
    plate = new EdbPlateP();
 
    if(i==-1) {
      pc = (EdbPlateP *)(ePC.At(0));
      plate->SetID(pc->GetLayer(1)->ID());
      z=0;
    }
    else {
      pc = (EdbPlateP *)(ePC.At(i));
      plate->SetID(pc->GetLayer(2)->ID());
      z = eB.GetPlate(eB.Npl()-1)->Z() - pc->GetLayer(1)->Z();
      aff.Reset();
      aff.Transform( pc->GetLayer(1)->GetAffineXY() );
      aff.Invert();
      aff.Transform( eB.GetPlate(eB.Npl()-1)->GetAffineXY() );
      plate->GetAffineXY()->Transform(&aff);
    }
    plate->SetZlayer(z, z - dz0/2 + dz1, z+dz0/2+dz2);                
    plate->GetLayer(0)->SetZlayer(0,-dz0/2,dz0/2);       // internal plate coord
    plate->GetLayer(2)->SetZlayer(-dz0/2,-dz0/2-dz2,-dz0/2);
    plate->GetLayer(1)->SetZlayer( dz0/2, dz0/2, dz0/2+dz1);
 
    eB.AddPlate(plate);    
  }
  MakePIDList();
 
  eReferencePlate = eB.GetPlate(0)->ID();
 
  return eB.Npl();
}

Brick() [1/2]

EdbBrickP & EdbScanSet::Brick ( )

inline

{return eB;}

Brick() [2/2]

const EdbBrickP & EdbScanSet::Brick ( ) const

inline

{return eB;}

Copy()

void EdbScanSet::Copy

(

EdbScanSet &

sc

)

copy from ss to this

{
  eB.Copy(sc.eB);
  eID = sc.eID;
 
  int npc = sc.ePC.GetEntries();
  for(int i=0; i<npc; i++) {
    EdbPlateP *p = new EdbPlateP();
    p->Copy( *((EdbPlateP *)sc.ePC.At(i)) );
    ePC.Add(p);
  }
 
  TIter next( &(sc.eIDS) );
  TObject *obj=0;
  while((obj = next()))
    eIDS.Add( new EdbID( *((EdbID*)obj) ) );
 
  MakePIDList();
  eReferencePlate = sc.eReferencePlate;
}

FindNextPlateID()

EdbID * EdbScanSet::FindNextPlateID	(	Int_t	p,
		Bool_t	direction
	)

< find next valid identifier in respect to plate in positive or negative direction (+-1)

{
  Log(2,"FindNextPlateID","plate = %d  direction = %d", plate, positive_direction);
  EdbID *id = 0;
  Int_t plmin=kMaxInt, plmax=kMinInt;
  for (Int_t i = 0; i < eIDS.GetEntries(); i++) {
    id = (EdbID *)eIDS.At(i);
    if (id->ePlate>plmax) plmax=id->ePlate;
    if (id->ePlate<plmin) plmin=id->ePlate;
  }
  int step = positive_direction? 1 : -1;
  int first = plate +step;
  int last = positive_direction? plmax:plmin;   last += step;
  printf("first=%d, last = %d step = %d\n",first,last, step);
  if( positive_direction &&  first>=last)  return 0;
  if( !positive_direction && first<=last)  return 0;
  for(int i=first; i!=last; i+=step) {
    id = FindPlateID(i);
    if(id) return id;
  }
  return 0;
}

FindPlateID()

EdbID * EdbScanSet::FindPlateID

(

Int_t

p

)

{
  EdbID *id;
  for (Int_t i = 0; i < eIDS.GetEntries(); i++) {
    id = (EdbID *)eIDS.At(i);
    if (id->ePlate == plate) return id;
  }
  return 0;
}

GetAffP2P()

Bool_t EdbScanSet::GetAffP2P	(	Int_t	p1,
		Int_t	p2,
		EdbAffine2D &	aff
	)

input: p1 - id of the "from-plate" p2 - id of the "to-plate" output: aff - the transformation - applied to "from" segment became "to" segment

{
  Int_t ip1,ip2;
  if(ePID.Find(p1)) ip1 = ePID.Find(p1)->At(0)->Value();  else return false;
  if(ePID.Find(p2)) ip2 = ePID.Find(p2)->At(0)->Value();  else return false;
  aff.Reset();
  aff.Transform( eB.GetPlate(ip2)->GetAffineXY() );
  aff.Invert();
  aff.Transform( eB.GetPlate(ip1)->GetAffineXY() );
  return true;
}

GetDZP2P()

Float_t EdbScanSet::GetDZP2P	(	Int_t	p1,
		Int_t	p2
	)

input: p1 - id of the "from-plate" p2 - id of the "to-plate" return: dz - segment of from-plate propagated to dz become in to-plate

{
  Int_t ip1,ip2;
  if(ePID.Find(p1)) ip1 = ePID.Find(p1)->At(0)->Value();  else return false;
  if(ePID.Find(p2)) ip2 = ePID.Find(p2)->At(0)->Value();  else return false;
  return  eB.GetPlate(ip2)->Z() - eB.GetPlate(ip1)->Z();
}

GetID()

EdbID * EdbScanSet::GetID ( Int_t  i )

inline

{return (EdbID *)(eIDS.At(i));}

GetLayer()

EdbLayer * EdbScanSet::GetLayer ( Int_t  p, Int_t  side  )

inline

{ if(GetPlate(p)) if(ValidSide(side)) return GetPlate(p)->GetLayer(side); return 0; }

GetPlate()

EdbPlateP * EdbScanSet::GetPlate ( Int_t  p )

inline

                                {
    if (ePID.Find(p)) 
      return eB.GetPlate(ePID.Find(p)->At(0)->Value()); 
    else return 0; 
  }

MakeNominalSet()

void EdbScanSet::MakeNominalSet	(	EdbID	id,
		int	from_plate = 57,
		int	to_plate = 1,
		float	z0 = 0,
		float	dz = -1300,
		float	shr = 1,
		float	dzbase = 210.,
		float	dzem = 45.
	)

{
  eID=id;
  Log(2,"EdbScanSet::MakeNominalSet","%s",eID.AsString());
  eReferencePlate = from_plate;
  eB.SetID(id.eBrick);
  Int_t step= (from_plate>to_plate)?-1:1;
  Int_t p=from_plate;
  Float_t z=z0;
  do {
    EdbPlateP *plate = new EdbPlateP();
    plate->SetID(p);
    plate->SetZ(z);
    plate->SetPlateLayout(dzbase,dzem,dzem);
 
//     plate->SetZlayer(z,z-dzbase/2-dzem,z+dzbase/2+dzem);
//     plate->GetLayer(0)->SetZlayer(0,-dzbase/2,dzbase/2);
//     plate->GetLayer(1)->SetZlayer(dzbase/2,dzbase/2,dzbase/2+dzem);
//     plate->GetLayer(2)->SetZlayer(-dzbase/2,-dzbase/2-dzem,-dzbase/2);
 
    plate->GetLayer(1)->SetShrinkage(shr);
    plate->GetLayer(2)->SetShrinkage(shr);
    eB.AddPlate(plate);
    eIDS.Add(new EdbID(id.eBrick,p,id.eMajor,id.eMinor));
    p += step;
    z += dz;
  } while( p!=(to_plate+step) );
  MakePIDList();
  if(gEDBDEBUGLEVEL>2) Print();
}

MakeParFiles()

Int_t EdbScanSet::MakeParFiles	(	Int_t	piece = 0,
		const char *	dir = "."
	)

prepare par files for the dataset for "recset-like" processing

{
 
  TString name;
  Log(2,"MakeParFiles","for brick %d with %d plates:\n", eB.ID(), eB.Npl());
  EdbPlateP *p=0;
  for(Int_t i=0; i<eB.Npl(); i++) {
    p = eB.GetPlate(i);
    name.Form("%s/%2.2d_%3.3d.par",dir,p->ID(),piece);
    Log(3,"MakeParFiles","%s\n",name.Data());
    FILE *f = fopen(name.Data(),"w");
    EdbAffine2D *a=p->GetAffineXY();
    fprintf(f,"INCLUDE default.par\n");
    fprintf(f,"ZLAYER 0  %15.3f 0. 0.\n", p->Z());
    fprintf(f,"AFFXY  0  %15.6f %9.6f %9.6f %9.6f %15.6f %15.6f\n",
        a->A11(),a->A12(),a->A21(),a->A22(),a->B1(),a->B2());
    fclose(f);
  }
  return eB.Npl();
}

MakePIDList()

void EdbScanSet::MakePIDList

(

)

{
  ePID.Delete();
  Long_t v[2];
  for(Int_t i=0; i<eB.Npl(); i++) {
    v[0]= eB.GetPlate(i)->ID();
    v[1]= i;
    ePID.Add(2,v);
  }
  ePID.Sort();
}

Print()

void EdbScanSet::Print

(

)

{
  EdbPlateP *p=0;
  printf("EdbScanSet: %s\n", eID.AsString() );
//   int npc = ePC.GetEntries();
//   printf("%d couples\n",npc);
//   for(Int_t i=0; i<npc; i++) {
//     p = (EdbPlateP *)(ePC.At(i));
//     p->GetLayer(1)->Print();
//     p->GetLayer(2)->Print();
//  }
  printf("Brick %d with %d plates:\n", eB.ID(), eB.Npl());
  for(Int_t i=0; i<eB.Npl(); i++) { 
    p = eB.GetPlate(i);
    if(!p) continue;
    EdbAffine2D *a=p->GetAffineXY();
    if(!a) continue;
    printf("%3d  %12.2f       %9.6f %9.6f %9.6f %9.6f %15.6f %15.6f\n",
       p->ID(), p->Z(), a->A11(),a->A12(),a->A21(),a->A22(),a->B1(),a->B2());
  }
  printf("for this brick %d identifiers are defined:\n", eIDS.GetEntries());
  for(Int_t i=0; i<eIDS.GetEntries(); i++)  
  {
    EdbID *id = (EdbID*)eIDS.At(i);
    TString str = id->AsString();
    EdbPlateP *p = GetPlate(id->ePlate);
    if(p) str += Form("  shr1/base/shr2 = %.3f / %.1f / %.3f", p->GetLayer(1)->Shr(),p->GetLayer(0)->DZ(), p->GetLayer(2)->Shr());
    printf("%s\n",str.Data() );
  }
}

ReadGeom()

void EdbScanSet::ReadGeom

(

const char *

fname

)

todo

{
 
  char line[256];
  FILE *f = fopen (file, "r"); 
  if(!f) {Log(1,"EdbScanSet::ReadGeom","Can not open file %s", file); return; }
 
  int b,p,mi,ma;
  int n=0;
 
  if(fgets(line, 256, f) != NULL) if( sscanf(line, "%d", &n) != 1 ) return;
  printf("%d\n",n);
  if(!n) return;
 
  ePID.Delete();
  eB.Clear();
 
  EdbID *id1=0, *id2=0;
  float dz=0;
  for(int i=0; i<n; i++) {
    if(i==0) {
      if( fgets(line, 256, f) == NULL )  return;
      if( sscanf(line, "%d.%d.%d.%d", &b, &p, &mi, &ma) != 4 )  break;
      id1 = new EdbID(b,p,mi,ma);
      if(n==1) eIDS.Add(id1);
    }
    if( n>1 && i<n-1) {
      if( fgets(line, 256, f) == NULL )  return;
      if( sscanf(line, "%f", &dz) != 1 ) break;
      if( fgets(line, 256, f) == NULL )  return;
      if( sscanf(line, "%d.%d.%d.%d", &b, &p, &mi, &ma) != 4 )  break;
      id2 = new EdbID(b,p,mi,ma);
      eIDS.Add(id2);
    }
  }
  fclose(f);
}

ReadIDS()

Int_t EdbScanSet::ReadIDS

(

const char *

file

)

{
  char line[256];
  FILE *f = fopen (file, "rt");
  Int_t b,p,mi,ma;
  Int_t cnt=0;
  while(fgets(line, 256, f) != NULL)
    {
      if( sscanf(line, "%d, %d, %d, %d", &b, &p, &mi, &ma) != 4 )  break;
      eIDS.Add(new EdbID(b,p,mi,ma));
      cnt++;
    }
  fclose(f);
  return cnt;
}

RemovePlate()

void EdbScanSet::RemovePlate

(

int

pid

)

Remove both from eB and from eIDS all data corresponding to plate id note that gaps and affine transformations may become inconsistent and had to be updated after this operation

{
  eB.RemovePlate(pid);
  EdbID   *id = FindPlateID(pid);
  if(id) {
    eIDS.Remove(id);
    MakePIDList();
  }
}

SetAsReferencePlate()

Bool_t EdbScanSet::SetAsReferencePlate

(

Int_t

pid

)

{
  EdbPlateP *p = GetPlate(pid);
  if(!p) return kFALSE;
  eReferencePlate=pid;
  Float_t  z = p->Z();
  EdbAffine2D aff( *(p->GetAffineXY()) );
  aff.Invert();
  TransformBrick(aff);
  ShiftBrickZ(-z);
  return kTRUE;
}

SetID()

void EdbScanSet::SetID ( EdbID  id )

inline

{eID=id;}

ShiftBrickZ()

Int_t EdbScanSet::ShiftBrickZ

(

Float_t

z

)

{
  Int_t npl=eB.Npl();
  if(npl>0)  for(Int_t i=0; i<npl; i++) eB.GetPlate(i)->ShiftZ(z);
  return npl;
}

TransformBrick() [1/2]

Int_t EdbScanSet::TransformBrick

(

EdbAffine2D

aff

)

{
  Int_t npl=eB.Npl();
  if(npl>0)
    for(Int_t i=0; i<npl; i++)  eB.GetPlate(i)->GetAffineXY()->Transform(aff);
  return npl;
}

TransformBrick() [2/2]

void EdbScanSet::TransformBrick

(

EdbScanSet &

ss

)

{
  int n = eIDS.GetEntries();
  for(int i=0; i<n; i++) {
    EdbID *id   = GetID(i);                  if(!id) continue;
    EdbPlateP *p  = GetPlate(id->ePlate);    if(!p)  continue;
    EdbPlateP *p1 = ss.GetPlate(id->ePlate); if(!p1) continue;
    p->GetAffineXY()->Transform( p1->GetAffineXY() );
    p->SetZlayer(p1->Z(),p1->Zmin(),p1->Zmax());
    // ToDo: add also AFFTXTY?
  }
}

TransformSidesIntoBrickRS()

Int_t EdbScanSet::TransformSidesIntoBrickRS

(

)

to put layers 1&2 of each plate into brick reference system - to have absolute transformations for microtracks

{
  Int_t npl=eB.Npl();
  if(npl>0)
    for(Int_t i=0; i<npl; i++) 
       eB.GetPlate(i)->TransformSidesIntoPlateRS();
  return npl;
}

UpdateBrickWithP2P()

void EdbScanSet::UpdateBrickWithP2P	(	EdbLayer &	la,
		int	plate1,
		int	plate2
	)

update the brick geometry with aff & dz defined in layer (plate1->plate2) in this version assumed that plate1 and plate2 are consequtive in the brick Note eReferencePlate should be correctly defined

{
 
  int pfrom = TMath::Abs(plate1-eReferencePlate)<TMath::Abs(plate2-eReferencePlate)? plate2 : plate1;
  int step  = ((pfrom-eReferencePlate)>0)? 1: -1;  //the corrections to be propagated in the opposite direction from ref
  
  Log(3,"UpdateBrickWithP2P","from,step: %d %d   ref: %d  Z(%d)=%f",pfrom,step,eReferencePlate,pfrom,GetPlate(pfrom)->Z());
 
  float dz           = (GetDZP2P(plate1,plate2) - la.Zcorr());
 
  EdbAffine2D *affxy   = la.GetAffineXY();
  if(pfrom==plate2) { 
    dz *= -1.;
    affxy->Invert();
  }
  //Log(1,"UpdateBrickWithP2P","%d->%d  dz= %f",plate1,plate2, dz);
  int ipl=pfrom;
  for(int i=0; i<eB.Npl(); i++) {
    EdbPlateP *p = GetPlate(ipl);    if(!p) continue;
    p->SetZlayer( p->Z() + dz, p->Zmin(), p->Zmax() );
    p->GetAffineXY()->Transform(affxy);
    ipl+=step;
  }
  
  EdbPlateP *plate = GetPlate(plate1);
  plate->SetShrinkage( GetPlate(plate1)->Shr() * la.Shr() );
  //plate->GetAffineTXTY()->Print();
  plate->GetAffineTXTY()->Transform( la.GetAffineTXTY() );
  //plate->GetAffineTXTY()->Print();
  
  plate->ApplyCorrectionsLocal(la.Map());
  
  //Log( 1, "UpdateBrickWithP2P","after: Z(%d)=%f", pfrom, GetPlate(pfrom)->Z() );
 
}

UpdateGap()

void EdbScanSet::UpdateGap	(	float	dz,
		int	pid1,
		int	pid2
	)

update gap between plates pid1 and pid2. All updates will be done in a way to have the reference plate unchanged

{
  EdbLayer la; la.SetZcorr(dz);
  UpdateBrickWithP2P(la,pid1,pid2);
}

UpdateIDS()

void EdbScanSet::UpdateIDS	(	int	brick,
		int	ma,
		int	mi
	)

update all identifiers in dataset with (brik,ma,mi) keeping plateids unchanged

{
  
  eID.Set(brick, eID.ePlate, ma, mi);
  TIter next( &(eIDS) );
  TObject *obj=0;
  while((obj = next())) {
    EdbID *id=(EdbID*)obj;
    id->Set(brick, id->ePlate, ma, mi);
  }
}

ValidSide()

bool EdbScanSet::ValidSide ( int  side ) const

inline

{ if(side>-1&&side<3) return 1; return 0;}

WriteGeom()

void EdbScanSet::WriteGeom

(

const char *

fname

)

{
  FILE *f = 0;
  if(file) f = fopen(file, "w");
  if(!f) {Log(1,"EdbScanSet::WriteGeom","Can not open file %s", file); return; }
  EdbID *id1=0, *id2=0;
  int n = eIDS.GetEntries();
  fprintf(f,"%d\n", n);
  for( Int_t i=0; i<eIDS.GetEntries(); i++ ) {
    id1 = (EdbID *)eIDS.At(i);
    if(i==0) fprintf(f,"%d.%d.%d.%d\n", id1->eBrick, id1->ePlate, id1->eMajor, id1->eMinor );
    if( n>1 && i<n-1) id2 = (EdbID *)eIDS.At(i+1);    else continue;
    fprintf(f,"%f\n", GetDZP2P(id1->ePlate, id2->ePlate));
    fprintf(f,"%d.%d.%d.%d\n", id2->eBrick, id2->ePlate, id2->eMajor, id2->eMinor );
  }
 
  if(file) fclose(f);
}

WriteIDS()

Int_t EdbScanSet::WriteIDS

(

const char *

file = 0

)

{
  FILE *f = 0;
  if(file) f = fopen(file, "w");
  EdbID *id;
  for( Int_t i=0; i<eIDS.GetEntries(); i++ ) {
    id = (EdbID *)eIDS.At(i);
    if(file) fprintf(f,"%d.%d.%d.%d\n", id->eBrick, id->ePlate, id->eMajor, id->eMinor );
    else      printf("%d.%d.%d.%d\n", id->eBrick, id->ePlate, id->eMajor, id->eMinor );
  }
 
  if(file) fclose(f);
  return eIDS.GetEntries();
}

Zlayer()

float EdbScanSet::Zlayer ( int  plate, int  side  )

inline

{ if(GetLayer(plate,side)) return GetLayer(plate,side)->Z() + GetPlate(plate)->Z(); else return 0; }

Member Data Documentation

eB

EdbBrickP EdbScanSet::eB

all layers of the brick defined here

eID

EdbID EdbScanSet::eID

id of the scanset itself

eIDS

TList EdbScanSet::eIDS

list of the identifiers to be processed

ePC

TObjArray EdbScanSet::ePC

Plate Couples. Each couple represented as EdbPlateP where 2 sides in reality corresponds to 2 plates.

ePID

TIndexCell EdbScanSet::ePID

correspondance between index in eB and the plate id

eReferencePlate

Int_t EdbScanSet::eReferencePlate

---

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

• /home/antonio/fedra_doxygen/src/libEbase/EdbScanSet.h
• /home/antonio/fedra_doxygen/src/libEbase/EdbScanSet.cxx