EdbDataPiece Class Reference

Edb raw data unit (scanned plate) associated with run file. More...

#include <EdbDataSet.h>

Inheritance diagram for EdbDataPiece:

Collaboration diagram for EdbDataPiece:

Public Member Functions
int	AcceptViewHeader (const EdbViewHeader *head)
void	AddCutCP (float var[6])
void	AddRCut (int layer, TCut &cut)
void	AddRunFile (const char *name)
void	AddSegmentCut (int layer, int xi, float min[5], float max[5])
void	AddSegmentCut (int layer, int xi, float var[10])
float	CalculateSegmentChi2 (EdbSegment *seg, float sx, float sy, float sz)
int	CheckCCD (int maxentr=2000)
void	CloseCPData ()
void	CloseRun ()
int	CorrectAngles ()
int	CorrectAngles (TTree *tree)
void	CorrectShrinkage (int layer, float shr)
	EdbDataPiece ()
	EdbDataPiece (int plate, int piece, char *file, int flag)
int	Flag () const
int	GetAreaData (EdbPVRec *ali, int const area, int const side)
EdbScanCond *	GetCond (int id)
int	GetCPData (EdbPattern *pat, EdbPattern *p1=0, EdbPattern *p2=0)
int	GetCPData_new (EdbPattern *pat, EdbPattern *p1=0, EdbPattern *p2=0, TIndex2 *trseg=0)
EdbSegmentCut *	GetCut (int layer, int i)
float	GetCutGR () const
EdbLayer *	GetLayer (int id)
int	GetLinkedSegEntr (int side, int aid, int vid, int sid, TArrayI &entr) const
EdbScanCond *	GetMakeCond (int id)
EdbLayer *	GetMakeLayer (int id)
const char *	GetNameCP () const
int	GetOUTPUT () const
int	GetRawData (EdbPVRec *ali)
float	GetRawSegmentPix (EdbSegment *seg)
TCut *	GetRCut (int layer)
const char *	GetRunFile (int i) const
int	InitCouplesInd ()
int	InitCouplesTree (const char *mode="READ")
int	MakeLinkListArea (int irun)
int	MakeLinkListCoord (int irun)
const char *	MakeName ()
const char *	MakeNameCP (const char *dir)
void	MakeNamePar (const char *dir)
int	NCuts (int layer)
int	Nruns () const
int	PassCutCP (float var[6])
int	PassCuts (int id, float var[5])
int	Plate () const
void	Print ()
int	ReadPiecePar (const char *file)
int	RemoveCCDPeak (TMatrix &matr)
void	Set0 ()
void	SetCouplesTree (TTree *tree)
void	SetCutGR (float chi)
void	SetOUTPUT (int out=1)
void	SetVolume0 (float x0, float y0, float z0, float tx=0, float ty=0)
void	SetVolumeA (float dx, float dy)
int	TakeCPSegment (EdbSegCouple &cp, EdbSegP &segP)
int	TakePiecePar ()
int	TakeRawSegment (EdbView *view, int id, EdbSegP &segP, int side)
int	UpdateAffPar (int layer, EdbAffine2D &aff)
int	UpdateAffTPar (int layer, EdbAffine2D &aff)
int	UpdateSegmentCut (EdbSegmentCut cut)
int	UpdateShrPar (int layer)
int	UpdateZPar (int layer, float z)
void	WriteCuts ()
virtual	~EdbDataPiece ()

Static Public Member Functions
static TTree *	InitCouplesTree (const char *file, const char *mode)

Public Attributes
Int_t	eAFID
	1-use fiducial marks transformations, 0 - do not More...
Int_t	eCLUST
	1-use clusters, 0 - do not More...
TIndexCell *	eCouplesInd
TTree *	eCouplesTree
EdbMask *	eEraseMask
	id's (entries) of segments to be erased when read couples tree More...
TString	eFileNameCP
	root-style text cuts More...
TString	eFileNamePar
	name of the parameters file More...
EdbRun *	eRun

Private Attributes
TIndexCell *	eAreas [3]
	base/up/down surface areas list More...
EdbScanCond *	eCond [3]
Float_t	eCutCP [6]
Float_t	eCutGR
	grain cut (chi) More...
TObjArray *	eCuts [3]
	array of cuts More...
Int_t	eFlag
	0-do nothing, 1-do something More...
EdbLayer *	eLayers [3]
	base(0),up(1),down(2) layers More...
Int_t	eOUTPUT
Int_t	ePiece
	piece id in this plate More...
Int_t	ePlate
	plate id More...
TCut *	eRCuts [3]
TObjArray	eRunFiles

Detailed Description

Edb raw data unit (scanned plate) associated with run file.

//////////////////////////////////////////////////////////////////////// // EdbDataSet // // OPERA data set definition&reconstruction // // ////////////////////////////////////////////////////////////////////////

Constructor & Destructor Documentation

EdbDataPiece() [1/2]

EdbDataPiece::EdbDataPiece

(

)

EdbDataPiece() [2/2]

EdbDataPiece::EdbDataPiece	(	int	plate,
		int	piece,
		char *	file,
		int	flag
	)

---

{
  for(int i=0; i<3; i++) eLayers[i]=0;
  Set0();
  ePlate=plate;
  ePiece=piece;
  AddRunFile(file);
  eFlag=flag;
}

~EdbDataPiece()

EdbDataPiece::~EdbDataPiece ( )

virtual

---

{
   int i;
   for(i=0; i<3; i++)  if(eLayers[i])   {delete eLayers[i]; eLayers[i]=0;}
   for(i=0; i<3; i++)  if(eAreas[i])    {delete eAreas[i];  eAreas[i] =0;}
   for(i=0; i<3; i++)  if(eCond[i])     {delete eCond[i];   eCond[i]  =0;}
   for(i=0; i<3; i++)  if(eCuts[i])     {delete eCuts[i];   eCuts[i]  =0;}
   for(i=0; i<3; i++)  if(eRCuts[i])    {delete eRCuts[i];  eRCuts[i] =0;}
   if( eCouplesInd) { delete eCouplesInd; eCouplesInd=0;}
   if( eEraseMask ) { delete eEraseMask;  eEraseMask=0; }
   eRunFiles.Delete();
   CloseRun();
   CloseCPData();
}

Member Function Documentation

AcceptViewHeader()

int EdbDataPiece::AcceptViewHeader

(

const EdbViewHeader *

head

)

---

{
  EdbSegP p;
  if(eAFID==0) {
    p.SetX( head->GetXview() );
    p.SetY( head->GetYview() );
  } else {
    p.SetX( 0. );
    p.SetY( 0. );
    p.Transform( head->GetAffine() );
  }
  p.SetX( p.X() - GetLayer(0)->X() );
  p.SetY( p.Y() - GetLayer(0)->Y() );
  p.Transform( GetLayer(0)->GetAffineXY() );
  if( p.X() < -GetLayer(0)->DX() )    return 0;
  if( p.X() >  GetLayer(0)->DX() )    return 0;
  if( p.Y() < -GetLayer(0)->DY() )    return 0;
  if( p.Y() >  GetLayer(0)->DY() )    return 0;
  return 1;
}

AddCutCP()

void EdbDataPiece::AddCutCP

(

float

var[6]

)

---

{
  for(int i=0; i<6; i++) eCutCP[i]=var[i];
}

AddRCut()

void EdbDataPiece::AddRCut	(	int	layer,
		TCut &	cut
	)

---

{
  if(layer<0) return;
  if(layer>2) return;
  
  if(!eRCuts[layer]) eRCuts[layer] = new TCut(cut);
  else (*(eRCuts[layer]))+=cut;
}

AddRunFile()

void EdbDataPiece::AddRunFile

(

const char *

name

)

---

{
  TObjString *str = new TObjString(name);
  eRunFiles.Add(str);
}

AddSegmentCut() [1/2]

void EdbDataPiece::AddSegmentCut	(	int	layer,
		int	xi,
		float	min[5],
		float	max[5]
	)

---

{
  if(!eCuts[layer])  eCuts[layer] = new TObjArray();
  EdbSegmentCut *cut=new EdbSegmentCut();
  cut->SetXI(xi);
  cut->SetMin(min);
  cut->SetMax(max);
  eCuts[layer]->Add( cut );
}

AddSegmentCut() [2/2]

void EdbDataPiece::AddSegmentCut	(	int	layer,
		int	xi,
		float	var[10]
	)

---

{
  if(!eCuts[layer])  eCuts[layer] = new TObjArray();
  eCuts[layer]->Add( new EdbSegmentCut(xi,var) );
}

CalculateSegmentChi2()

float EdbDataPiece::CalculateSegmentChi2	(	EdbSegment *	seg,
		float	sx,
		float	sy,
		float	sz
	)

---

TODO: remove this function from here (already in EdbViewRec)

{
 
    //assumed that clusters are attached to segments
  double chi2=0;
  EdbCluster *cl=0;
  TObjArray *clusters = seg->GetElements();
  if(!clusters) return 0;
  int ncl = clusters->GetLast()+1;
  if(ncl<=0)     return 0;
 
  float xyz1[3], xyz2[3];             // segment line parametrized as 2 points
  float xyz[3];
  bool inside=true;
 
  xyz1[0] = seg->GetX0() /sx;
  xyz1[1] = seg->GetY0() /sy;
  xyz1[2] = seg->GetZ0() /sz;
  xyz2[0] = (seg->GetX0() + seg->GetDz()*seg->GetTx()) /sx;
  xyz2[1] = (seg->GetY0() + seg->GetDz()*seg->GetTy()) /sy;
  xyz2[2] = (seg->GetZ0() + seg->GetDz())              /sz;
 
  double d;
  for(int i=0; i<ncl; i++ ) {
    cl = (EdbCluster*)clusters->At(i);
    xyz[0] = cl->GetX()/sx;
    xyz[1] = cl->GetY()/sy;
    xyz[2] = cl->GetZ()/sz;
    d = EdbMath::DistancePointLine3(xyz,xyz1,xyz2, inside);
    chi2 += d*d;
  }
 
  return TMath::Sqrt(chi2/ncl);
}

CheckCCD()

int EdbDataPiece::CheckCCD

(

int

maxentr = 2000

)

---

{
  if (eRun ) delete eRun;
  eRun =  new EdbRun( GetRunFile(0),"READ" );
  if(!eRun) { Log(1,"EdbDataPiece::CheckCCD","ERROR open file: %s",GetRunFile(0)); return -1; }
  EdbView    *view = eRun->GetView();
  EdbSegment *seg;
 
  int npeak=0;
  
  TMatrix matr(1000,1000);
  int ix,iy;
 
  int i,j;
  for(i=0; i<1000; i++) 
    for(j=0; j<1000; j++)
      matr[i][j]=0;
 
  int ncheck=0;
  int nentr = TMath::Min(maxentr,eRun->GetEntries());
  Log(2,"EdbDataPiece::CheckCCD","nentr=%d",nentr);
  for (i=0; i<nentr; i++) {
    view = eRun->GetEntry(i);
    int nseg=view->Nsegments();
    for (j=0; j<nseg; j++) {
      seg = view->GetSegment(j);
      if( seg->GetTx() >  .05 )   continue;
      if( seg->GetTx() < -.05 )   continue;
      if( seg->GetTy() >  .05 )   continue;
      if( seg->GetTy() < -.05 )   continue;
      ix = (Int_t)(seg->GetX0()+500.);
      iy = (Int_t)(seg->GetY0()+500.);
      (matr[ix][iy])++;
      ncheck++;
    }
  }
 
  Log(2,"EdbDataPiece::CheckCCD","ncheck=%d",ncheck);
  for(i=0; i<200; i++ ) {           //eliminate upto 200 CCD defects
    if(!RemoveCCDPeak(matr)) break;
    npeak++;
  }
  return npeak;
}

CloseCPData()

void EdbDataPiece::CloseCPData

(

)

---

{
  if (eCouplesTree) {
    TFile *f = eCouplesTree->GetDirectory()->GetFile();
    if (f) {
      if(f->IsWritable()) eCouplesTree->AutoSave();
      SafeDelete(eCouplesTree);
      SafeDelete(f);
    }
    eCouplesTree=NULL;
  }
}

CloseRun()

void EdbDataPiece::CloseRun

(

)

---

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

CorrectAngles() [1/2]

int EdbDataPiece::CorrectAngles

(

)

---

{
  int n=0;
  TTree *cptree=0;
  cptree=EdbDataPiece::InitCouplesTree(GetNameCP(),"READ");
  if( cptree )  n = CorrectAngles(cptree);
  Log(2,"CorrectAngles","in piece: %s  using %d basetracks",GetNameCP(),n);
  return n;
}

CorrectAngles() [2/2]

int EdbDataPiece::CorrectAngles

(

TTree *

tree

)

---

{
  EdbSegCouple    *cp = 0;
  EdbSegP         *s1 = 0;
  EdbSegP         *s2 = 0;
  EdbSegP         *s  = 0;
  tree->SetBranchAddress("cp"  , &cp  );
  tree->SetBranchAddress("s1." , &s1  );
  tree->SetBranchAddress("s2." , &s2  );
  tree->SetBranchAddress("s."  , &s   );
 
  EdbAffine2D *aff = new EdbAffine2D();
 
  int nentr = (int)(tree->GetEntries());
  TArrayF x(nentr);
  TArrayF y(nentr);
  TArrayF x1(nentr);
  TArrayF y1(nentr);
  TArrayF x2(nentr);
  TArrayF y2(nentr);
 
  int nseg = 0;
  Log(2,"EdbDataPiece::CorrectAngles","nentr = %d",nentr);
  for(int i=0; i<nentr; i++ ) {
    tree->GetEntry(i);
 
    if(cp->N1tot()>1)  continue;
    if(cp->N2tot()>1)  continue;
    if(cp->CHI2()>1.5) continue;
 
    x1[nseg] = s1->TX();
    y1[nseg] = s1->TY();
    x2[nseg] = s2->TX();
    y2[nseg] = s2->TY();
    x[nseg]  = (s2->X()-s1->X()) / (s2->Z()-s1->Z());
    y[nseg]  = (s2->Y()-s1->Y()) / (s2->Z()-s1->Z());
    nseg++;
  }
 
  aff->CalculateTurn( nseg,x1.fArray,y1.fArray,x.fArray,y.fArray );
  UpdateAffTPar(1,*aff);
  aff->CalculateTurn( nseg,x2.fArray,y2.fArray,x.fArray,y.fArray );
  UpdateAffTPar(2,*aff);
 
  delete aff;
 
  // Closing tree and file. A.C.
 
  TFile *f = tree->GetDirectory()->GetFile();
  SafeDelete(tree);
  if (f) SafeDelete(f);
 
  return nseg;
}

CorrectShrinkage()

void EdbDataPiece::CorrectShrinkage	(	int	layer,
		float	shr
	)

---

{
  GetLayer(layer)->SetShrinkage( shr*GetLayer(layer)->Shr() );
}

Flag()

int EdbDataPiece::Flag ( ) const

inline

{return eFlag;}

GetAreaData()

int EdbDataPiece::GetAreaData	(	EdbPVRec *	ali,
		int const	area,
		int const	side
	)

---

{
  TIndexCell *elist   = eAreas[side]->At(aid);
  if(!elist)  return 0;
  EdbPattern *pat = new EdbPattern( 0.,0., GetLayer(side)->Z() );
  pat->SetID(side);
 
 
  EdbSegP    segP;
  EdbView    *view = eRun->GetView();
  int   nseg=0, nrej=0;
  int   entry;
  int   nsegV;
 
  int niu=elist->N();
  for(int iu=0; iu<niu; iu++) {
    entry = elist->At(iu)->Value();
    if(eCLUST)       {
      view = eRun->GetEntry(entry,1,1,1);
      view->AttachClustersToSegments();
    }
    else             view = eRun->GetEntry(entry);
 
    nsegV = view->Nsegments();
 
    for(int j=0;j<nsegV;j++) {
      if(!TakeRawSegment(view,j,segP,side)) {
    nrej++;
    continue;
      }
      nseg++;
      segP.SetVid(entry,j);
      segP.SetAid(view->GetAreaID(),view->GetViewID(),side);
      pat->AddSegment( segP);
    }
  }
 
  Log(2,"EdbDataPiece::GetAreaData","Area: %d ( %d%%)  %d \t views: %d \t nseg: %d \t rejected: %d", 
       aid,100*aid/eAreas[side]->N(1),side,niu,nseg, nrej );
  pat->SetSegmentsZ();
  ali->AddPattern(pat);
  return nseg;
}

GetCond()

EdbScanCond * EdbDataPiece::GetCond ( int  id )

inline

    { if(eCond[id]) return (EdbScanCond *)eCond[id]; else return 0; }

GetCPData()

int EdbDataPiece::GetCPData	(	EdbPattern *	pat,
		EdbPattern *	p1 = 0,
		EdbPattern *	p2 = 0
	)

---

TO CHECK !!!

{
  Log(2,"EdbDataPiece::GetCPData","z = %f", GetLayer(0)->Z());
  pat->SetID(0);
  EdbSegP    segP;
 
  TTree *tree=eCouplesTree;
  EdbSegCouple    *cp = 0;
  EdbSegP         *s1 = 0;
  EdbSegP         *s2 = 0;
  EdbSegP         *s  = 0;
 
  TBranch *b_cp=0, *b_s=0, *b_s1=0, *b_s2=0;
  b_cp = tree->GetBranch("cp");
  b_s  = tree->GetBranch("s.");
  b_s1 = tree->GetBranch("s1.");
  b_s2 = tree->GetBranch("s2.");
 
  b_cp->SetAddress( &cp  );
  b_s->SetAddress(  &s   );
  b_s1->SetAddress( &s1  );
  b_s2->SetAddress( &s2  );
 
 
  int nseg = 0;
  int nentr = (int)(tree->GetEntries());
  for(int i=0; i<nentr; i++ ) {
    tree->GetEntry(i);
    b_cp->GetEntry(i);
    b_s->GetEntry(i);
    if( !TakeCPSegment(*cp,*s) )      continue;
    if(pat) {
      s->SetZ( s->Z() + pat->Z() );   
      //s->SetPID( ePlate*10 );       /// TO CHECK !!!
      EdbTraceBack::SetBaseTrackVid( *s, ePlate, ePiece, i );
      s->SetChi2(cp->CHI2P());
      pat->AddSegment( *s  );
      nseg++;
    }
    if(p1)  { 
      b_s1->GetEntry(i);                             // !!!
      s1->SetZ( s1->Z() + pat->Z() );
      //s1->SetPID( ePlate*10 +1 );    /// TO CHECK !!!
      p1->AddSegment(  *s1 ); 
      nseg++; 
    }
    if(p2)  { 
      b_s2->GetEntry(i);                             // !!!
      s2->SetZ( s2->Z() + pat->Z() );
      //s2->SetPID( ePlate*10 + 2 );   /// TO CHECK !!!
      p2->AddSegment(  *s2 ); 
      nseg++; 
    }
  }
 
  Log(2,"EdbDataPiece::GetCPData","%d (of %d) segments are readed", nseg,nentr );
 
  return nseg;
}

GetCPData_new()

int EdbDataPiece::GetCPData_new	(	EdbPattern *	pat,
		EdbPattern *	p1 = 0,
		EdbPattern *	p2 = 0,
		TIndex2 *	trseg = 0
	)

---

TO CHECK !!!

{
  Log(3,"EdbDataPiece::GetCPData_new","z = %f ", GetLayer(0)->Z());
  if(!eCouplesTree)  return  0;
  pat->SetID(0);
  EdbSegP    segP;
 
  TTree *tree=eCouplesTree;
  EdbSegCouple    *cp = 0;
  EdbSegP         *s1 = 0;
  EdbSegP         *s2 = 0;
  EdbSegP         *s  = 0;
 
  TBranch *b_cp=0, *b_s=0, *b_s1=0, *b_s2=0;                        // !!!
  b_cp = tree->GetBranch("cp");                   // !!!
  b_s  = tree->GetBranch("s.");                   // !!!
  b_s1 = tree->GetBranch("s1.");                  // !!!
  b_s2 = tree->GetBranch("s2.");                  // !!!
 
  b_cp->SetAddress( &cp  );
  b_s->SetAddress(  &s   );
  b_s1->SetAddress( &s1  );
  b_s2->SetAddress( &s2  );
 
  int nseg = 0;
  int nentr = (int)(tree->GetEntries());  if(nentr<1) return 0;
 
  TCut *cut = GetRCut(0);
  TEventList *lst =0;
  if(cut)       tree->Draw(">>lst", *cut );
  else          tree->Draw(">>lst", "" );
  lst = (TEventList*)(gDirectory->GetList()->FindObject("lst"));
 
  if(!lst) {Log(1,"EdbDataPiece::GetCPData_new","ERROR!: lst do not found! empty couples tree??"); return 0;}
 
  int nlst =lst->GetN();
 
  int entr=0;
  for(int i=0; i<nlst; i++ ) {
    entr = lst->GetEntry(i);
 
    if(trseg) {           //exclude segments participating in tracks
      if( (trseg->Find(ePlate*1000+ePiece,entr) >= 0) )  continue;
    }
 
    if(eEraseMask) if(eEraseMask->At(entr)) continue;
 
    b_cp->GetEntry(entr);                            // !!!
    b_s->GetEntry(entr);                             // !!!
    if( !TakeCPSegment(*cp,*s) )      continue;
    if(pat) {
      s->SetZ( s->Z() + pat->Z() );   
      //s->SetPID( ePlate*10 );       /// TO CHECK !!!
      s->SetVid(ePlate*1000+ePiece,entr);
      s->SetChi2(cp->CHI2P());
      pat->AddSegment( *s  );
      nseg++;
    }
    if(p1)  { 
      b_s1->GetEntry(entr);                             // !!!
      s1->SetZ( s1->Z() + pat->Z() );
      //s1->SetPID( ePlate*10 +1 );    /// TO CHECK !!!
      p1->AddSegment(  *s1 ); 
      nseg++; 
    }
    if(p2)  { 
      b_s2->GetEntry(entr);                             // !!!
      s2->SetZ( s2->Z() + pat->Z() );
      //s2->SetPID( ePlate*10 + 2 );   /// TO CHECK !!!
      p2->AddSegment(  *s2 ); 
      nseg++; 
    }
  }
 
  SafeDelete(lst);
 
  if(cut) Log(2,"EdbDataPiece::GetCPData_new","select %d of %d segments by cut %s",nlst, nentr, cut->GetTitle() );
  else    Log(2,"EdbDataPiece::GetCPData_new","%d (of %d) segments accepted", nseg,nentr );
 
  return nseg;
}

GetCut()

EdbSegmentCut * EdbDataPiece::GetCut ( int  layer, int  i  )

inline

    { return (EdbSegmentCut *)(eCuts[layer]->UncheckedAt(i)); }

GetCutGR()

float EdbDataPiece::GetCutGR ( ) const

inline

{return eCutGR;}

GetLayer()

EdbLayer * EdbDataPiece::GetLayer ( int  id )

inline

    { if(eLayers[id]) return (EdbLayer *)eLayers[id]; else return 0; }

GetLinkedSegEntr()

int EdbDataPiece::GetLinkedSegEntr	(	int	side,
		int	aid,
		int	vid,
		int	sid,
		TArrayI &	entr
	)		const

---

{
  if(!eCouplesInd) return 0;
  Long_t v[4];
  v[0]=side;  v[1]=aid;  v[2]=vid;  v[3]=sid;
  TIndexCell *c = eCouplesInd->Find(4,v);
  if(!c) return 0;
  int n = c->N();
  entr.Set(n);
  for(int i=0; i<n; i++) entr.AddAt( (int)(c->At(i)->Value()), i);
  return n;
}

GetMakeCond()

EdbScanCond * EdbDataPiece::GetMakeCond

(

int

id

)

---

{
  if(id<0) return 0;
  if(id>2) return 0;
  if(!GetCond(id))  eCond[id] = new EdbScanCond();
  return GetCond(id);
}

GetMakeLayer()

EdbLayer * EdbDataPiece::GetMakeLayer

(

int

id

)

---

{
  if(id<0) return 0;
  if(id>2) return 0;
  if(!GetLayer(id))  eLayers[id] = new EdbLayer();
  return GetLayer(id);
}

GetNameCP()

const char * EdbDataPiece::GetNameCP ( ) const

inline

{return eFileNameCP.Data();}

GetOUTPUT()

int EdbDataPiece::GetOUTPUT ( ) const

inline

{return eOUTPUT;}

GetRawData()

int EdbDataPiece::GetRawData

(

EdbPVRec *

ali

)

---

TODO: irun logic

{
 
  CloseRun();
  eRun =  new EdbRun( GetRunFile(0),"READ" );
  if(!eRun) { Log(1,"EdbDataPiece::GetRawData","ERROR open file: %s",GetRunFile(0)); return -1; }
 
  EdbPattern *pat1 = new EdbPattern( 0.,0., GetLayer(1)->Z() + GetLayer(0)->Z() );
  EdbPattern *pat2 = new EdbPattern( 0.,0., GetLayer(2)->Z() + GetLayer(0)->Z() );
 
  EdbViewHeader *head=0;
  EdbSegP        segP;
  EdbView       *view = eRun->GetView();
  int      nseg=0, nrej=0;
  int      nsegV=0;
  int      side=0;
 
  int nentr = eRun->GetEntries();
  for(int iv=0; iv<nentr; iv++ ) {
    head = eRun->GetEntryHeader(iv);
    if(      head->GetNframesTop()==0 ) side=2;
    else if( head->GetNframesBot()==0 ) side=1;
 
    if( !AcceptViewHeader(head) )  continue;
 
 
    if(eCLUST)       {
      view = eRun->GetEntry(iv,1,1,1);
      view->AttachClustersToSegments();
    }
    else             view = eRun->GetEntry(iv);
 
    nsegV = view->Nsegments();
    for(int j=0;j<nsegV;j++) {
      if(!TakeRawSegment(view,j,segP,side)) {
    nrej++;
    continue;
      }
      nseg++;
      segP.SetVid(iv,j);
      segP.SetAid(view->GetAreaID(),view->GetViewID(),side);
 
      if(side==1) pat1->AddSegment( segP );
      else if(side==2) pat2->AddSegment( segP );
    }
 
  }
  ali->AddPattern(pat1);
  ali->AddPattern(pat2);
  return nseg;
}

GetRawSegmentPix()

float EdbDataPiece::GetRawSegmentPix

(

EdbSegment *

seg

)

---

assumed that clusters are attached to segments

{
  float pix=0;
  EdbCluster *cl=0;
  TObjArray *clusters = seg->GetElements();
  if(!clusters) return 0;
  int ncl = clusters->GetLast()+1;
  for(int i=0; i<ncl; i++ ) {
    cl = (EdbCluster*)clusters->At(i);
    pix += cl->GetArea();
  }
  return pix;
}

GetRCut()

TCut * EdbDataPiece::GetRCut ( int  layer )

inline

{ return eRCuts[layer]; }

GetRunFile()

const char * EdbDataPiece::GetRunFile

(

int

i

)

const

---

{
  if(Nruns()<i+1) return 0;
  return ((TObjString *)eRunFiles.At(i))->GetName();
}

InitCouplesInd()

int EdbDataPiece::InitCouplesInd

(

)

---

{
  if(eCouplesInd) delete eCouplesInd;
  eCouplesInd = new TIndexCell();
 
  if(!eCouplesTree) if(!InitCouplesTree("READ")) return 0;
  EdbSegP         *s1  = 0;
  EdbSegP         *s2  = 0;
  eCouplesTree->SetBranchAddress("s1."  , &s1  );
  eCouplesTree->SetBranchAddress("s2."  , &s2  );
 
  Long_t v[5]; // side:aid:vid:sid:entry
 
  int nentr = (int)(eCouplesTree->GetEntries());
 
  for(int i=0; i<nentr; i++) {
    eCouplesTree->GetEntry(i);
    v[0] = 1;
    v[1] = s1->Aid(0);
    v[2] = s1->Aid(1);
    v[3] = s1->Vid(1);
    v[4] = i;
    eCouplesInd->Add(5,v);
    v[0] = 2;
    v[1] = s2->Aid(0);
    v[2] = s2->Aid(1);
    v[3] = s2->Vid(1);
    v[4] = i;
    eCouplesInd->Add(5,v);
  }
  eCouplesInd->Sort();
  return eCouplesInd->N();
}

InitCouplesTree() [1/2]

TTree * EdbDataPiece::InitCouplesTree ( const char *  file, const char *  mode  )

static

{
  static TString tree_name("couples");
  TTree *tree = NULL;
 
  // checking for the existing directory
 
  FileStat_t buf;
 
  if (gSystem->GetPathInfo(gSystem->DirName(file_name), buf)) {
    if (gEDBDEBUGLEVEL > 0) {
      cout << "ERROR! Directory " << gSystem->DirName(file_name) 
       << " does not exist.\n";
      cout << "Please run scanning procedure first for the plate " 
       << gSystem->BaseName(gSystem->DirName(file_name)) << endl;
    }
    return 0;
 
  }
 
 
  TFile *f = new TFile(file_name, mode);
 
  if (f->IsOpen()) tree = (TTree*)f->Get(tree_name);
  else return 0;
 
  if (!tree) {
    f->cd();
    tree = new TTree(tree_name, tree_name);
    tree->SetMaxTreeSize(15000000000LL);   //set 15 Gb file size limit)
    //tree->SetMaxVirtualSize( 512 * 1024 * 1024 ); // default is 64000000
 
    int pid1=0,pid2=0;
    float xv=0,yv=0;
    EdbSegCouple *cp=0;
    EdbSegP      *s1=0;
    EdbSegP      *s2=0;
    EdbSegP      *s=0;
      
    tree->Branch("pid1",&pid1,"pid1/I");
    tree->Branch("pid2",&pid2,"pid2/I");
    tree->Branch("xv",&xv,"xv/F");
    tree->Branch("yv",&yv,"yv/F");
    tree->Branch("cp","EdbSegCouple",&cp,32000,99);
    tree->Branch("s1.","EdbSegP",&s1,32000,99);
    tree->Branch("s2.","EdbSegP",&s2,32000,99);
    tree->Branch("s." ,"EdbSegP",&s,32000,99);
    tree->Write();
    // tree->SetAutoSave(2000000);
  }
 
  if(!tree) Log(1,"EdbDataPiece::InitCouplesTree","ERROR!!! InitCouplesTree: can't initialize tree at %s as %s",file_name,mode);
 
  return tree;
}

InitCouplesTree() [2/2]

int EdbDataPiece::InitCouplesTree

(

const char *

mode = "READ"

)

---

{
  if ((eCouplesTree=InitCouplesTree(eFileNameCP,mode))) return 1;
  return 0;
}

MakeLinkListArea()

int EdbDataPiece::MakeLinkListArea

(

int

irun

)

---

{
  if (eRun ) delete eRun;
  eRun =  new EdbRun( GetRunFile(irun),"READ" );
  if(!eRun) { Log(1," EdbDataPiece::MakeLinkListArea","ERROR open file: %s",GetRunFile(irun)); return -1; }
 
  for(int i=0; i<3; i++) {
    if(eAreas[i])    delete eAreas[i];
    eAreas[i]= new TIndexCell();
  }
 
  int nentr = eRun->GetEntries();
  Log(2,"EdbDataPiece::MakeLinkListArea","Make views entry map,  nentr = %d",nentr);
 
  Long_t v[2];   // areaID,entry
  EdbViewHeader *head=0;
 
  for(int iv=0; iv<nentr; iv++ ) {
    head = eRun->GetEntryHeader(iv);
    v[0]=head->GetAreaID();
    v[1]=iv;
    if(head->GetNframesTop()==0) {         // fill down views
      eAreas[2]->Add(2,v);
    }
    else if(head->GetNframesBot()==0) {    // fill up views
       eAreas[1]->Add(2,v);
    }
  }
 
  eAreas[1]->Sort();
  eAreas[2]->Sort();
  return eAreas[1]->N(1);
}

MakeLinkListCoord()

int EdbDataPiece::MakeLinkListCoord

(

int

irun

)

---

{
  if(eRun) delete eRun;
  eRun =  new EdbRun( GetRunFile(irun),"READ" );
  if(!eRun) { Log(1,"EdbDataPiece::MakeLinkListCoord","ERROR open file: %s",GetRunFile(0)); return -1; }
 
  
  for(int i=0; i<3; i++) {
    if(eAreas[i])    delete eAreas[i];
    eAreas[i]= new TIndexCell();
  }
 
  int nentr = eRun->GetEntries();
  Log(2,"EdbDataPiece::MakeLinkListCoord","Make views coordinate map,  nentr = %d",nentr);
 
  TIndexCell upc;
  TIndexCell downc;
  Long_t v[3];   // x,y,entry
  EdbViewHeader *head=0;
 
  Long_t xx=0, yy=0;
  float cx = 2000., cy = 2000.;  // 2x2 mm cells
  float dx = 400. , dy = 400.;   // 400 microns margins
  float xv=0,yv=0;
  int mx[9] = {0, 0, 0,-1,  1, -1,-1, 1, 1};
  int my[9] = {0,-1, 1, 0,  0, -1, 1,-1, 1};
 
  for(int iv=0; iv<nentr; iv++ ) {
 
    head = eRun->GetEntryHeader(iv);
    if( !AcceptViewHeader(head) )  continue;
 
    yv = head->GetXview();
    xv = head->GetYview();
    xx = (Long_t)(xv/cx);
    yy = (Long_t)(yv/cx);
    v[0] = xx;
    v[1] = yy;
    v[2] = iv;
 
    if(head->GetNframesBot()==0) {              // fill up views
      upc.Add(3,v);
    }
    else if(head->GetNframesTop()==0) {         // fill down views
      downc.Add(3,v);                           // add center
      int im;
      for( im=1; im<5; im++ ) {             // add sides margins
        v[0] = (Long_t)(( xv+dx*mx[im] ) / cx);
        v[1] = (Long_t)(( yv+dy*my[im] ) / cy);
        if( (v[0] != xx) || (v[1] != yy) )
            downc.Add(3,v);
      }
      for( im=5; im<9; im++ ) {             // add angles margins
        v[0] = (Long_t)(( xv+dx*mx[im] ) / cx);
        v[1] = (Long_t)(( yv+dy*my[im] ) / cy);
        if( (v[0] != xx) && (v[1] != yy) )
          //      if(!downc.Find(3,v))
            downc.Add(3,v);
      }
 
    }
 
  }
 
  upc.Sort();
  downc.Sort();
 
  TIndexCell *clx=0;
  TIndexCell *cly=0;
 
  int areac=0;
  int nix,niy,nie;
  nix=upc.N(1);
  for(int ix=0; ix<nix; ix++) {
    clx = upc.At(ix);
    xx = clx->Value();
    niy=clx->N(1);
    for(int iy=0; iy<niy; iy++) {
      cly  = clx->At(iy);
      yy = cly->Value();
      areac++;
 
      nie = cly->N(1);
      int ie;
      for(ie=0; ie<nie; ie++) {
        v[0]=areac;
        v[1]  = cly->At(ie)->Value();
        eAreas[1]->Add(2,v);
      }
 
      cly = downc.Find(xx)->Find(yy);
      if(!cly) continue;
      nie=cly->N(1);
      for(ie=0; ie<nie; ie++) {
        v[0] = areac;
        v[1] = cly->At(ie)->Value();
        eAreas[2]->Add(2,v);
      }
    }
  }
 
  eAreas[1]->Sort();
  eAreas[2]->Sort();
 
  return eAreas[1]->N(1);
}

MakeName()

const char * EdbDataPiece::MakeName

(

)

---

{
  char name[8];
  sprintf(name,"%2.2d_%3.3d", ePlate,ePiece);
  SetName(name);
  return GetName();
}

MakeNameCP()

const char * EdbDataPiece::MakeNameCP

(

const char *

dir

)

---

{
  eFileNameCP=dir;
  eFileNameCP+=MakeName();
  eFileNameCP+=".cp.root";
  return eFileNameCP.Data();
}

MakeNamePar()

void EdbDataPiece::MakeNamePar

(

const char *

dir

)

---

{
  eFileNamePar=dir;
  eFileNamePar += MakeName();
  eFileNamePar += ".par";
}

NCuts()

int EdbDataPiece::NCuts

(

int

layer

)

---

{
  if(!eCuts[layer])  return 0;
  return eCuts[layer]->GetEntriesFast();
}

Nruns()

int EdbDataPiece::Nruns ( ) const

inline

{ return eRunFiles.GetEntriesFast(); }

PassCutCP()

int EdbDataPiece::PassCutCP

(

float

var[6]

)

---

{
  if(eCutCP[0]<0) return 1;
  for(int i=0; i<6; i++)    if(var[i]>eCutCP[i])  return 0;
  return 1;
}

PassCuts()

int EdbDataPiece::PassCuts	(	int	id,
		float	var[5]
	)

---

{
  int nc = NCuts(id);
  for(int i=0; i<nc; i++)
    if( !(GetCut(id,i)->PassCut(var)) )  return 0;
  return 1;
}

Plate()

int EdbDataPiece::Plate ( ) const

inline

{return ePlate;}

Print()

void EdbDataPiece::Print

(

)

---

{
  printf("Piece: %s\n",GetName());
  printf("%d %d \n", ePlate,ePiece);
  int i;
  for(i=0; i<Nruns(); i++)
    printf("%s\n",GetRunFile(i));
  for(i=0; i<3; i++)  if(eLayers[i])  eLayers[i]->Print();
  for(i=0; i<3; i++)  if(eCond[i])    eCond[i]->Print();
  for(i=0; i<3; i++)
    if(eCuts[i])
      for(int j=0; j<NCuts(i); j++)  GetCut(i,j)->Print();
}

ReadPiecePar()

int EdbDataPiece::ReadPiecePar

(

const char *

file

)

---

read parameters from par-file
return: 0 if ok
-1 if file access failed

{
 
  char buf[512];
  char key[256];
  char name[512];
 
  FILE *fp = fopen(file,"r");
  if (!fp) {
    Log(2,"ReadPiecePar","ERROR open file: %s", file);
    return -1;
  }
  else Log(3,"ReadPiecePar","Read piece parameters from file: %s", file );
 
  int id,mode;
  float z,zmin,zmax,shr;
  float a11,a12,a21,a22,b1,b2;
  float x1,x2,x3,x4;
  float var[10];
 
  while (fgets(buf, sizeof(buf), fp)) {
    for (Int_t i = 0; i < (Int_t)strlen(buf); i++) 
      if (buf[i]=='#')  {
    buf[i]='\0';                       // cut out comments starting from #
    break;
      }
    
    if( sscanf(buf,"%s",key)!=1 )                             continue;
 
    if      ( !strcmp(key,"INCLUDE")   )
      {
    sscanf(buf+strlen(key),"%s",name);
    ReadPiecePar(name);
      }
    else if ( !strcmp(key,"ZLAYER")   )
      {
    sscanf(buf+strlen(key),"%d %f %f %f",&id,&z,&zmin,&zmax);
    GetMakeLayer(id)->SetZlayer(z,zmin,zmax);
      }
    else if ( !strcmp(key,"SHRINK")  )
      {
    sscanf(buf+strlen(key),"%d %f",&id,&shr);
    GetMakeLayer(id)->SetShrinkage(shr);
      }
    else if ( !strcmp(key,"AFFXY")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f %f %f",&id,&a11,&a12,&a21,&a22,&b1,&b2);
    GetMakeLayer(id)->SetAffXY(a11,a12,a21,a22,b1,b2);
      }
    else if ( !strcmp(key,"AFFTXTY")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f %f %f",&id,&a11,&a12,&a21,&a22,&b1,&b2);
    GetMakeLayer(id)->SetAffTXTY(a11,a12,a21,a22,b1,b2);
      }
    else if ( !strcmp(key,"SIGMA0")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f",&id,&x1,&x2,&x3,&x4);
    GetMakeCond(id)->SetSigma0(x1,x2,x3,x4);
      }
    else if ( !strcmp(key,"DEGRAD")  )
      {
    sscanf(buf+strlen(key),"%d %f",&id,&x1);
    GetMakeCond(id)->SetDegrad(x1);
      }
    else if ( !strcmp(key,"BINS")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f",&id,&x1,&x2,&x3,&x4);
    GetMakeCond(id)->SetBins(x1,x2,x3,x4);
      }
    else if ( !strcmp(key,"RAMP0")  )
      {
    sscanf(buf+strlen(key),"%d %f %f",&id,&x1,&x2);
    GetMakeCond(id)->SetPulsRamp0(x1,x2);
      }
    else if ( !strcmp(key,"RAMP04")  )
      {
    sscanf(buf+strlen(key),"%d %f %f",&id,&x1,&x2);
    GetMakeCond(id)->SetPulsRamp04(x1,x2);
      }
    else if ( !strcmp(key,"CHI2MAX")  )
      {
    sscanf(buf+strlen(key),"%d %f",&id,&x1);
    GetMakeCond(id)->SetChi2Max(x1);
      }
    else if ( !strcmp(key,"CHI2PMAX")  )
      {
    sscanf(buf+strlen(key),"%d %f",&id,&x1);
    GetMakeCond(id)->SetChi2PMax(x1);
      }
    else if ( !strcmp(key,"CHI2MODE")  )
      {
    sscanf(buf+strlen(key),"%d %d",&id,&mode);
    GetMakeCond(id)->SetChi2Mode(mode);
      }
    else if ( !strcmp(key,"OFFSET")  )
      {
    sscanf(buf+strlen(key),"%d %f %f",&id,&x1,&x2);
    GetMakeCond(id)->SetOffset(x1,x2);
      }
    else if ( !strcmp(key,"SIGMAGR")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f",&id,&x1,&x2,&x3,&x4);
    GetMakeCond(id)->SetSigmaGR(x1,x2,x3);
    SetCutGR(x4);
      }
    else if ( !strcmp(key,"RADX0")  )
      {
    sscanf(buf+strlen(key),"%d %f",&id,&x1);
    GetMakeCond(id)->SetRadX0(x1);
      }
    else if ( !strcmp(key,"XCUT")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f %f %f %f %f %f %f",&id,
           var,var+1,var+2,var+3,var+4,var+5,var+6,var+7,var+8,var+9);
    AddSegmentCut(id,0,var);
      }
    else if ( !strcmp(key,"ICUT")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f %f %f %f %f %f %f",&id,
           var,var+1,var+2,var+3,var+4,var+5,var+6,var+7,var+8,var+9);
    AddSegmentCut(id,1,var);
      }
    else if ( !strcmp(key,"RCUT")  )
      {
    char rcut[256];
    sscanf(buf+strlen(key),"%d %s",&id, rcut );
    TCut cut(rcut);
    AddRCut(id,cut);
      }
    else if ( !strcmp(key,"CUTCP")  )
      {
    sscanf(buf+strlen(key),"%d %f %f %f %f %f %f",&id,
           var,var+1,var+2,var+3,var+4,var+5);
    AddCutCP(var);
      }
    else if ( !strcmp(key,"AFID")  )
      {
    sscanf(buf+strlen(key),"%d",&id);
    eAFID=id;
      }
    else if ( !strcmp(key,"OUTPUT")  )
      {
    sscanf(buf+strlen(key),"%d",&id);
    eOUTPUT=id;
      }
    else if ( !strcmp(key,"VOLUME0")  )
      {
    float x0=0,y0=0,z0=0,tx=0,ty=0;
    sscanf(buf+strlen(key),"%f %f %f %f %f",&x0,&y0,&z0,&tx,&ty);
    SetVolume0(x0,y0,z0,tx,ty);
      }
    else if ( !strcmp(key,"VOLUMEA")  )
      {
    float dx=0,dy=0;
    sscanf(buf+strlen(key),"%f %f",&dx,&dy);
    SetVolumeA(dx,dy);
      }
    else if ( !strcmp(key,"CLUST")  )
      {
    int icl=0;
    sscanf(buf+strlen(key),"%d",&icl);
    eCLUST=icl;
      }
  }
  fclose(fp);
 
  return 0;
}

RemoveCCDPeak()

int EdbDataPiece::RemoveCCDPeak

(

TMatrix &

matr

)

---

{
  double mean=0;
  int filled=0;
  float max=0;
  int ix=0, iy=0;
  int nc = matr.GetNcols();
  int nr = matr.GetNrows();
  for(int i=0; i<nc; i++) 
    for(int j=0; j<nr; j++) {
      if(matr[i][j]<=0) continue;
      filled++;
      mean+=matr[i][j];
      if(max<matr[i][j]) {
    max = matr[i][j];
    ix = i;
    iy = j;
      }
    }
  mean/=filled;
  Log(2,"EdbDataPiece::RemoveCCDPeak","mean = %f \t max[%d,%d]=%d",mean,ix,iy,(int)max);
 
  float vmin[5],vmax[5];
  EdbSegmentCut cut;
  if( max > 10.*mean ) {
    matr[ix][iy]=0;
    vmin[0] = ix-500.;    vmax[0] = vmin[0]+1;
    vmin[1] = iy-500.;    vmax[1] = vmin[1]+1;
    vmin[2] = -.6;        vmax[2] = .6;
    vmin[3] = -.6;        vmax[3] = .6;
    vmin[4] =   0;        vmax[4] = 17;
    cut.SetXI(0);
    cut.SetMin(vmin);
    cut.SetMax(vmax);
    UpdateSegmentCut(cut);
    return 1;
  }
  return 0;
}

Set0()

void EdbDataPiece::Set0

(

)

---

{
  eOUTPUT=0;
  ePlate=0;
  ePiece=0;
  eFlag=0;
  eAFID=0;
  eCLUST=0;
  eCutCP[0]=-1;
  int i;
  for(i=0; i<3; i++) { 
    if(eLayers[i]) delete eLayers[i];
    eLayers[i]=new EdbLayer();
  }
  for(i=0; i<3; i++) eCond[i]=0;
  for(i=0; i<3; i++) eAreas[i]=0;
  for(i=0; i<3; i++) eCuts[i]=0;
  for(i=0; i<3; i++) eRCuts[i]=0;
  eRun    = 0;
  eCouplesTree=0;
  eCouplesInd=0;
  eEraseMask=0;
}

SetCouplesTree()

void EdbDataPiece::SetCouplesTree ( TTree *  tree )

inline

{eCouplesTree=tree;}

SetCutGR()

void EdbDataPiece::SetCutGR ( float  chi )

inline

{eCutGR=chi;}

SetOUTPUT()

void EdbDataPiece::SetOUTPUT ( int  out = 1 )

inline

{eOUTPUT=out;}

SetVolume0()

void EdbDataPiece::SetVolume0	(	float	x0,
		float	y0,
		float	z0,
		float	tx = 0,
		float	ty = 0
	)

---

{
  float z = GetLayer(0)->Z();
  GetLayer(0)->SetXY( x0+(z-z0)*tx, y0+(z-z0)*ty );
  GetLayer(0)->SetTXTY(tx,ty);
}

SetVolumeA()

void EdbDataPiece::SetVolumeA ( float  dx, float  dy  )

inline

{ GetLayer(0)->SetDXDY( dx, dy); }

TakeCPSegment()

int EdbDataPiece::TakeCPSegment	(	EdbSegCouple &	cp,
		EdbSegP &	seg
	)

---

{
  float var[6];
 
  var[0] = cp.N1();
  var[1] = cp.N1tot();
  var[2] = cp.N2();
  var[3] = cp.N2tot();
  var[4] = cp.CHI2();
  var[5] = cp.CHI2P();
 
  if( !PassCutCP(var) )     return 0;
 
  var[0] = seg.X();
  var[1] = seg.Y();
  var[2] = seg.TX();
  var[3] = seg.TY();
  var[4] = seg.W();
 
  if( !PassCuts(0,var) )     return 0;
 
  return 1;
}

TakePiecePar()

int EdbDataPiece::TakePiecePar

(

)

---

{
  return ReadPiecePar( eFileNamePar.Data() );
}

TakeRawSegment()

int EdbDataPiece::TakeRawSegment	(	EdbView *	view,
		int	id,
		EdbSegP &	segP,
		int	side
	)

---

{
  EdbSegment *seg = view->GetSegment(id);
 
  float var[5];
  var[0] = seg->GetX0();
  var[1] = seg->GetY0();
  var[2] = seg->GetTx();
  var[3] = seg->GetTy();
  var[4] = seg->GetPuls();
 
  if( !PassCuts(side,var) )     return 0;
 
  float pix, chi2;
  if(eCLUST) {
    pix = GetRawSegmentPix(seg);
    segP.SetVolume( pix );
    chi2 = CalculateSegmentChi2( seg,
                 GetCond(1)->SigmaXgr(),  //TODO: side logic
                 GetCond(1)->SigmaYgr(), 
                 GetCond(1)->SigmaZgr());
 
    if(chi2>GetCutGR())  return 0;
    segP.SetChi2( chi2 );    
  }
 
  EdbLayer  *layer=GetLayer(side);
  seg->Transform(layer->GetAffineXY());                          //internal view transformation
  if(eAFID) seg->Transform( view->GetHeader()->GetAffine() );
 
  float x,y,z,tx,ty,puls;
  tx   = seg->GetTx()/layer->Shr();
  ty   = seg->GetTy()/layer->Shr();
  x    = seg->GetX0() + layer->Zmin()*tx;
  y    = seg->GetY0() + layer->Zmin()*ty;
  z    = layer->Z() + layer->Zmin();
  if(eAFID==0) {
    x+=view->GetXview();
    y+=view->GetYview();
  }
  puls = seg->GetPuls();
 
  EdbAffine2D *aff = layer->GetAffineTXTY();
  float txx = aff->A11()*tx+aff->A12()*ty+aff->B1();
  float tyy = aff->A21()*tx+aff->A22()*ty+aff->B2();
  segP.Set( seg->GetID(),x,y,txx,tyy,puls,0);
  segP.SetZ( z );
  segP.SetDZ( seg->GetDz()*layer->Shr() );
  segP.SetW( puls );
  segP.SetVolume( seg->GetVolume() );
  segP.SetChi2( seg->GetSigmaY() );   // make sence in case of fedra tracking 
  segP.SetMC( view->GetHeader()->GetEvent(),view->GetHeader()->GetTrack() );
  return 1;
}

UpdateAffPar()

int EdbDataPiece::UpdateAffPar	(	int	layer,
		EdbAffine2D &	aff
	)

---

{
  const char *file=eFileNamePar.Data();
 
  FILE *fp=fopen(file,"a");
  if (!fp) {
    Log(1,"EdbDataPiece::UpdateAffPar","ERROR open file: %s", file);
    return -1;
  } else Log(2,"EdbDataPiece::UpdateAffPar","Update parameters file with AFFXY: %s", file );
 
  char str[124];
  sprintf(str,"AFFXY \t %d \t %f %f %f %f %f %f\n",layer,
     aff.A11(),aff.A12(),aff.A21(),aff.A22(),aff.B1(),aff.B2() );
  fprintf(fp,"\n%s",str);
 
  fclose(fp);
  return 1;
}

UpdateAffTPar()

int EdbDataPiece::UpdateAffTPar	(	int	layer,
		EdbAffine2D &	aff
	)

---

{
 
  EdbAffine2D *a = GetLayer(layer)->GetAffineTXTY();
  a->Transform(&aff);
 
  const char *file=eFileNamePar.Data();
 
  FILE *fp=fopen(file,"a");
  if (!fp) {
    Log(1,"EdbDataPiece::UpdateAffTPar","ERROR open file: %s", file);
    return -1;
  }
  else Log(2,"EdbDataPiece::UpdateAffTPar","\nUpdate parameters file with AFFTXTY: %s\n\n", file );
 
  char str[124];
  sprintf(str,"AFFTXTY \t %d \t %f %f %f %f %f %f\n",layer,
     a->A11(),a->A12(),a->A21(),a->A22(),a->B1(),a->B2() );
  fprintf(fp,"\n%s",str);
 
  fclose(fp);
  return 1;
}

UpdateSegmentCut()

int EdbDataPiece::UpdateSegmentCut

(

EdbSegmentCut

cut

)

---

{
  const char *file=eFileNamePar.Data();
 
  FILE *fp=fopen(file,"a");
  if (!fp) {
    Log(1,"EdbDataPiece::UpdateSegmentCut","ERROR open file: %s", file);
    return -1;
  }
 
  char str[124];
  if(cut.XI()==0) {
    Log(2,"EdbDataPiece::UpdateSegmentCut","Update parameters file with XCUT: %s", file );
    sprintf(str,"XCUT \t %d \t %f %f %f %f %f %f %f %f %f %f\n",1,
        cut.Min(0),cut.Max(0), 
        cut.Min(1),cut.Max(1), 
        cut.Min(2),cut.Max(2), 
        cut.Min(3),cut.Max(3), 
        cut.Min(4),cut.Max(4)
        );
    fprintf(fp,"\n%s",str);
    sprintf(str,"XCUT \t %d \t %f %f %f %f %f %f %f %f %f %f\n",2,
        cut.Min(0),cut.Max(0), 
        cut.Min(1),cut.Max(1), 
        cut.Min(2),cut.Max(2), 
        cut.Min(3),cut.Max(3), 
        cut.Min(4),cut.Max(4)
        );
    fprintf(fp,"\n%s",str);
  } else  if(cut.XI()==1) {
    Log(2,"EdbDataPiece::UpdateSegmentCut","Update parameters file with ICUT: %s", file );
    sprintf(str,"ICUT \t %d \t %f %f %f %f %f %f %f %f %f %f\n",1,
        cut.Min(0),cut.Max(0), 
        cut.Min(1),cut.Max(1), 
        cut.Min(2),cut.Max(2), 
        cut.Min(3),cut.Max(3), 
        cut.Min(4),cut.Max(4)
        );
    fprintf(fp,"\n%s",str);
    sprintf(str,"XCUT \t %d \t %f %f %f %f %f %f %f %f %f %f\n",2,
        cut.Min(0),cut.Max(0), 
        cut.Min(1),cut.Max(1), 
        cut.Min(2),cut.Max(2), 
        cut.Min(3),cut.Max(3), 
        cut.Min(4),cut.Max(4)
        );
    fprintf(fp,"\n%s",str);
  }
 
  fclose(fp);
  return 1;
}

UpdateShrPar()

int EdbDataPiece::UpdateShrPar

(

int

layer

)

---

{
  const char *file=eFileNamePar.Data();
 
  FILE *fp=fopen(file,"a");
  if (fp==NULL)   {
    Log(1,"EdbDataPiece::UpdateShrPar","ERROR open file: %s", file);
    return -1;
  }else Log(2,"EdbDataPiece::UpdateShrPar","Update parameters file with SHRINK %d: %s", layer, file );
 
  char str[64];
  sprintf(str,"SHRINK \t %d \t %f \n",layer, GetLayer(layer)->Shr() );
  fprintf(fp,"\n%s",str);
 
  fclose(fp);
  return 1;
}

UpdateZPar()

int EdbDataPiece::UpdateZPar	(	int	layer,
		float	z
	)

---

{
  const char *file=eFileNamePar.Data();
 
  FILE *fp=fopen(file,"a");
  if (!fp) {
    Log(1,"EdbDataPiece::UpdateZPar", "ERROR open file: %s", file);
    return -1;
  }
  else Log(2,"EdbDataPiece::UpdateZPar", "Update parameters file with ZLAYER: %s", file );
 
  char str[124];
  sprintf(str,"ZLAYER \t %d \t %f %f %f\n",layer,
     z,0.,0. );
  fprintf(fp,"\n%s",str);
 
  fclose(fp);
  return 1;
}

WriteCuts()

void EdbDataPiece::WriteCuts

(

)

---

{
  TString file = eFileNamePar+".C";
 
  FILE *fp=fopen(file,"w");
  if (!fp) {
    Log(1," EdbDataPiece::WriteCuts","ERROR open file: %s", file.Data());
    return;
  }
  else Log(2,"EdbDataPiece::WriteCuts","Put Cuts to file: %s", file.Data() );
 
  fprintf(fp,"{\n");
 
  char str[256];
  for(int i=0; i<3; i++)
    if(eCuts[i])
      for(int j=0; j<NCuts(i); j++)  {
    GetCut(i,j)->CutLine(str,i,j);
    fprintf(fp,"%s",str);
      }
 
  fprintf(fp,"}\n");
  fclose(fp);
}

Member Data Documentation

eAFID

Int_t EdbDataPiece::eAFID

1-use fiducial marks transformations, 0 - do not

eAreas

TIndexCell* EdbDataPiece::eAreas[3]

private

base/up/down surface areas list

eCLUST

Int_t EdbDataPiece::eCLUST

1-use clusters, 0 - do not

eCond

EdbScanCond* EdbDataPiece::eCond[3]

private

eCouplesInd

TIndexCell* EdbDataPiece::eCouplesInd

eCouplesTree

TTree* EdbDataPiece::eCouplesTree

eCutCP

Float_t EdbDataPiece::eCutCP[6]

private

eCutGR

Float_t EdbDataPiece::eCutGR

private

grain cut (chi)

eCuts

TObjArray* EdbDataPiece::eCuts[3]

private

array of cuts

eEraseMask

EdbMask* EdbDataPiece::eEraseMask

id's (entries) of segments to be erased when read couples tree

eFileNameCP

TString EdbDataPiece::eFileNameCP

root-style text cuts

name of the couples data file

eFileNamePar

TString EdbDataPiece::eFileNamePar

name of the parameters file

eFlag

Int_t EdbDataPiece::eFlag

private

0-do nothing, 1-do something

eLayers

EdbLayer* EdbDataPiece::eLayers[3]

private

base(0),up(1),down(2) layers

eOUTPUT

Int_t EdbDataPiece::eOUTPUT

private

ePiece

Int_t EdbDataPiece::ePiece

private

piece id in this plate

ePlate

Int_t EdbDataPiece::ePlate

private

plate id

eRCuts

TCut* EdbDataPiece::eRCuts[3]

private

eRun

EdbRun* EdbDataPiece::eRun

eRunFiles

TObjArray EdbDataPiece::eRunFiles

private

---

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

• /home/antonio/fedra_doxygen/src/libEIO/EdbDataSet.h
• /home/antonio/fedra_doxygen/src/libEIO/EdbDataSet.cxx