EdbDataStore Class Reference

#include <EdbDataStore.h>

Inheritance diagram for EdbDataStore:

Collaboration diagram for EdbDataStore:

Public Member Functions
void	AddPattern (EdbPattern *pat)
EdbSegP *	AddSegment (EdbSegP *seg, EdbSegmentCut *cut=0, int Plt0=0, int Plt1=100)
	add methods More...
void	AddTrack (EdbTrackP *tr)
void	AddVertex (EdbVertex *v)
void	Clear (bool hard=false)
	clear methods More...
void	ClearGeom ()
void	ClearRaw (bool hard=false)
void	ClearSeg (bool hard=false)
void	ClearTracks (bool hard=false)
void	ClearVTX ()
void	DoEfficiency (TF1 *eff_seg, TF1 *eff_mtk)
void	DoSmearing (EdbScanCond *cond_btk, EdbScanCond *cond_mtk=0)
	methods for simulation: More...
	EdbDataStore ()
EdbLayer *	FindLayer (int plate, int side=0)
EdbTrackP *	FindLongTrk (int nsmin=8)
EdbPattern *	FindPattern (int plate, int side=0)
EdbVertex *	FindPrimVtx ()
EdbTrackP *	FindTrack (int id)
	find methods More...
EdbVertex *	FindVertex (int id)
long	Gen_mtk_BG (long NBG, int Plate, int Side, TH2 *pdf_Ang, TH2 *pdf_WT=0)
EdbPattern *	GetPattern (int n, bool btk=true)
EdbPatternsVolume *	GetPV (bool btk=true)
EdbPattern *	GetRawPat (int n)
EdbPattern *	GetSegPat (int n)
EdbTrackP *	GetTrack (int n)
	get methods More...
EdbVertex *	GetVertex (int n)
void	LoadMCVertices (TObjArray *vtx)
	restore MC info methods More...
void	MakePattern (double z, int plate, int side)
int	Nplt ()
int	Nt ()
	count methods: More...
int	Nv ()
void	PrintBrief ()
	print methods More...
void	PrintPatterns ()
void	PrintTracks (int vlev=0)
void	Restore_PatFromGeom (int np0=0, int np1=1000)
void	Restore_PIDFromID ()
void	Restore_SegFromTrx (EdbSegmentCut *cut=0, int Plt0=0, int Plt1=1000)
void	Restore_TrxFromVtx ()
void	SavePlateToRaw (const char *fname, int PID, Option_t *option="RECREATE")
void	SaveToRaw (const char *dir="./", const EdbID &idset="0.0.0.0", Option_t *option="RECREATE", bool doaff=true)
	save methods: More...
void	SetOwnTracks (bool own=true)
	setown methods: More...
void	SetOwnTrkSegs ()
void	SetOwnVertices (bool own=true)
void	TransferGeometry (EdbDataStore *ds)
void	TransferTo (EdbDataStore *ds, char level, EdbSegmentCut *cut=0, int FromPlate=0, int ToPlate=57)
	transfer methods More...
	~EdbDataStore ()

Static Public Member Functions
static void	TransferSegs (EdbPatternsVolume *pv0, EdbPatternsVolume *pv1, EdbSegmentCut *cut=0, int FromPlate=0, int ToPlate=57)

Public Attributes
EdbBrickP *	eBrick
EdbPatternsVolume	eRawPV
	geometry More...
EdbPatternsVolume	eSegPV
TObjArray	eTracks
TObjArray	eVTX

Constructor & Destructor Documentation

EdbDataStore()

EdbDataStore::EdbDataStore

(

)

---

                          {
  eRawPV.Set0();
  eSegPV.Set0();
};

~EdbDataStore()

EdbDataStore::~EdbDataStore

(

)

---

{ Clear();};

Member Function Documentation

AddPattern()

void EdbDataStore::AddPattern

(

EdbPattern *

pat

)

---

                                            {
   assert(pat);
   assert(FindLayer(pat->Plate(),pat->Side()));
   if(pat->Side()==0)eSegPV.AddPattern(pat);
   else eRawPV.AddPattern(pat);
}

AddSegment()

EdbSegP * EdbDataStore::AddSegment	(	EdbSegP *	s,
		EdbSegmentCut *	cut = 0,
		int	Plt0 = 0,
		int	Plt1 = 100
	)

add methods

---

                                                                                  {
  if(s->Plate()<Plt0 || s->Plate()>Plt1)return 0;
      if(cut){
    float par[]={s->X(),s->Y(),s->TX(),s->TY(),s->W()};
    if(!cut->PassCut(par))return 0;
      }
      if(gEDBDEBUGLEVEL>5){
        s->PrintNice();
        printf("plt#%d side#%d\n",s->Plate(),s->Side());
        }
      EdbPattern* p=FindPattern(s->Plate(),s->Side());
      assert(p);
      EdbSegP* seg=p->AddSegment(*s);
      //p->SetID(s->Plate());
      seg->SetPID(p->ID());
      return seg;
}

AddTrack()

void EdbDataStore::AddTrack ( EdbTrackP *  tr )

inline

{assert(tr!=0); eTracks.Add(tr);}

AddVertex()

void EdbDataStore::AddVertex ( EdbVertex *  v )

inline

{assert(v!=0);  eVTX.Add(v);}

Clear()

void EdbDataStore::Clear ( bool  hard = false )

inline

clear methods

{ClearTracks(hard);ClearVTX();ClearRaw(hard);ClearSeg(hard);}

ClearGeom()

void EdbDataStore::ClearGeom

(

)

---

                            {
  if(eBrick==0)return;
  for(int np=0;np<eBrick->Npl();++np){
    eBrick->GetPlate(np)->Clear();
    delete eBrick->GetPlate(np);
  }
  eBrick->Clear();
}

ClearRaw()

void EdbDataStore::ClearRaw

(

bool

hard = false

)

---

                                    {
  EdbPatternsVolume* pv=&eRawPV;
 
  
  int N=pv->Npatterns();
  if(!N)return;
  TObjArray* Sgs=0;
  for(int k=0;k<N;k++) {
    Sgs=pv->GetPattern(k)->GetSegments();
    if(Sgs)Sgs->Delete();
//     delete pv->GetPattern(k);
  }
  if(hard){
    pv->ePatterns->SetOwner(1);
    pv->ePatterns->Clear();
    pv->Clear();
    pv->DropCell();
    pv->DropCouples();
  }
}

ClearSeg()

void EdbDataStore::ClearSeg

(

bool

hard = false

)

---

                                    {
  EdbPatternsVolume* pv=&eSegPV;
//   int N=pv->Npatterns();
//   if(!N)return;
//   TObjArray* Sgs=0;
//   for(int k=0;k<N;k++) {
//     Sgs=pv->GetPattern(k)->GetSegments();
//     if(Sgs)Sgs->Delete();
//   }
  int N=pv->Npatterns();
  if(!N)return;
  TObjArray* Sgs=0;
  for(int k=0;k<N;k++) {
    Sgs=pv->GetPattern(k)->GetSegments();
    if(Sgs)Sgs->Delete();
//     delete pv->GetPattern(k);
  }
  if(hard){
    pv->ePatterns->SetOwner(1);
    pv->DropCell();
    pv->DropCouples();
    pv->ePatterns->Clear();
    pv->Clear();
  }
}

ClearTracks()

void EdbDataStore::ClearTracks

(

bool

hard = false

)

---

                                       {
  eTracks.Clear();
}

ClearVTX()

void EdbDataStore::ClearVTX

(

)

---

                           {
  eVTX.Clear();
}

DoEfficiency()

void EdbDataStore::DoEfficiency	(	TF1 *	eff_seg,
		TF1 *	eff_mtk
	)

---

                                                        {
  EdbSegP* s=0;
  double eff, ran, th;
  if(eff_seg){
    for(int np=0;np<eSegPV.Npatterns();++np){
      EdbPattern* pat=GetSegPat(np);
      for(int ns=0; ns< pat->N();++ns){
        s=pat->GetSegment(ns);
        th=s->Theta();
        eff=(th<=eff_seg->GetXmax())?eff_seg->Eval(th):0;
        ran=gRandom->Rndm();
        if(ran>eff)s->SetW(-s->W());
     } 
    }
  }  
  if(eff_mtk){
    for(int np=0;np<eRawPV.Npatterns();++np){
      int n1=0;
      EdbPattern* pat=GetRawPat(np);
      for(int ns=0; ns< pat->N();++ns){
        s=pat->GetSegment(ns);
        th=s->Theta();
        Log(5,"DoEfficiency",Form("mtk#%d/%d theta=%2.4f [%2.4f %2.4f]\n",ns,pat->N(),th,s->TX(),s->TY()));
        eff=(th<=eff_mtk->GetXmax())?eff_mtk->Eval(th):0;
        ran=gRandom->Rndm();
        Log(5,"DoEfficiency",Form("eff=%2.4f, ran=%2.4f => %s\n",eff,ran,eff>ran?"survive":"KILL"));
        if(ran>eff)s->SetW(-s->W());
        else ++n1;
     } 
     Log(2,"DoEfficiency",Form("Plate #%d side#%d - efficiency %d/%d\n",pat->Plate(),pat->Side(),n1,pat->N()));
    }
  }
}

DoSmearing()

void EdbDataStore::DoSmearing	(	EdbScanCond *	cond_btk,
		EdbScanCond *	cond_mtk = 0
	)

methods for simulation:

---

                                                                        {
  EdbPVGen pvg;
  if(cond_btk){
    pvg.SetVolume(&eSegPV);
    pvg.SetScanCond(cond_btk);
    pvg.SmearSegments();
  }
  if(cond_mtk){
    pvg.SetVolume(&eRawPV);
    pvg.SetScanCond(cond_mtk);
    pvg.SmearSegments();
  }
}

FindLayer()

EdbLayer * EdbDataStore::FindLayer	(	int	plate,
		int	side = 0
	)

---

                                                     {
  assert(plate>=0);
  assert(side>=0 && side<3);
  
  EdbPlateP* plt=0;
  for(int np=0;np<eBrick->Npl();++np){
    plt=eBrick->GetPlate(np);
    if(plt->ID()==plate){
      return plt->GetLayer(side);
    }
  }
  return 0;
}

FindLongTrk()

EdbTrackP * EdbDataStore::FindLongTrk

(

int

nsmin = 8

)

---

find longest track

                                             {
  int nsmax=0;
  int imax=0;
  EdbTrackP* trk=0;
  for(int nt=0;nt<Nt();++nt){
    trk=GetTrack(nt);
    if(trk->PDG()!=22 && trk->P()>1.&& trk->N()>nsmax){
      imax=nt;
      nsmax=trk->N();
    }
  }
  if(nsmax<nsmin)return 0;
  return GetTrack(imax);
}

FindPattern()

EdbPattern * EdbDataStore::FindPattern	(	int	plate,
		int	side = 0
	)

---

                                                        {
  Log(3,"EdbDataStore::FindPattern",Form("search plt#%d side=%d",plate,side));
  assert(plate>=0);
  assert(side>=0 && side<3);
  
  EdbPatternsVolume* pv=0;
  pv=(side==0)?&eSegPV:&eRawPV;
  EdbPattern* pat=0;
  for(int np=0;np<pv->Npatterns();++np){
    pat=pv->GetPattern(np);
    //Log(6,"pat",Form("#%d/%d plt#%d side#%d",np,pv->Npatterns(),pat->Plate(),pat->Side()));
    if(pat->Plate()==plate && pat->Side()==side)return pat;
  }
  return 0;
}

FindPrimVtx()

EdbVertex * EdbDataStore::FindPrimVtx

(

)

---

                                    {
  EdbTrackP* trk=FindLongTrk(8);
  return trk->VertexS();
}

FindTrack()

EdbTrackP * EdbDataStore::FindTrack

(

int

id

)

find methods

---

                                        {
  EdbTrackP* trk=0;
  for(int nt=0;nt<eTracks.GetEntries();++nt){
    trk=GetTrack(nt);
//     printf("id%d/%d=%d\n",nt,eTracks.GetEntries(),trk->ID());
    if(trk->ID()==id)return trk;
  }
  return 0;
}

FindVertex()

EdbVertex * EdbDataStore::FindVertex

(

int

id

)

---

                                         {
  EdbVertex* vtx=0;
  for(int nv=0;nv<eVTX.GetEntries();++nv){
    vtx=GetVertex(nv);
    if(vtx->ID()==id)return vtx;
  }
  return 0;
}

Gen_mtk_BG()

long EdbDataStore::Gen_mtk_BG	(	long	NBG,
		int	Plate,
		int	Side,
		TH2 *	pdf_Ang,
		TH2 *	pdf_WT = 0
	)

---

Generate background microtracks, using distributions: 1)from pdf_Ang(TY:TX) - should be normalized so that its maximum bin content = 1 2)from pdb_WT (W:tan(Theta)) - each row with fixed X(i.e. slope) should have maximum content = 1

Side==0 - this is basetrack. So we need to project each segment to RAW patterns

                                                                                     {
  EdbPattern* pat=FindPattern(Plate,Side);
  EdbPattern *p1=0,*p2=0;
  if(Side==0){
    p1=FindPattern(Plate,1);
    p2=FindPattern(Plate,2);
  }
  assert(pat!=0);
  EdbLayer  * lay=FindLayer  (Plate,Side);
  assert(lay!=0);
  float x,y,z,tx,ty,w=16,p; 
  z=lay->Z();
  int Ntot=0;
  EdbSegP *seg=0,*s1=0;
  for(long n=0; n<NBG; ++n){
    tx=gRandom->Uniform(pdf_Ang->GetXaxis()->GetXmin(),pdf_Ang->GetXaxis()->GetXmax());
    ty=gRandom->Uniform(pdf_Ang->GetYaxis()->GetXmin(),pdf_Ang->GetYaxis()->GetXmax());
    p=pdf_Ang->Interpolate(tx,ty);
    Log(5,"Gen_mtk_BG",Form("tt=[%4.2f %4.2f] p=%6.4f\n",tx,ty,p));
    if(gRandom->Uniform()>p){n--; continue;}
    Ntot++;
    x=gRandom->Uniform(lay->Xmin(),lay->Xmax());
    y=gRandom->Uniform(lay->Ymin(),lay->Ymax());
    seg=pat->AddSegment(-1,x,y,tx,ty,w,0);
    if(p1){
      s1=p1->AddSegment(*seg);
      s1->PropagateTo(p1->Z());
      }
    if(p2){
      s1=p2->AddSegment(*seg);
      s1->PropagateTo(p2->Z());
      }
  }
  Log(0,"Gen_mtk_BG",Form("Generated %ld/%ld\n",Ntot,NBG));
  return Ntot;
}

GetPattern()

EdbPattern* EdbDataStore::GetPattern ( int  n, bool  btk = true  )

inline

{return GetPV(btk)->GetPattern(n);}

GetPV()

EdbPatternsVolume* EdbDataStore::GetPV ( bool  btk = true )

inline

{return btk?(&eSegPV):(&eRawPV);}

GetRawPat()

EdbPattern* EdbDataStore::GetRawPat ( int  n )

inline

{return (n<eRawPV.Npatterns())?eRawPV.GetPattern(n):0;}

GetSegPat()

EdbPattern* EdbDataStore::GetSegPat ( int  n )

inline

{return (n<eSegPV.Npatterns())?eSegPV.GetPattern(n):0;}

GetTrack()

EdbTrackP* EdbDataStore::GetTrack ( int  n )

inline

get methods

{return (n<eTracks.GetEntries())?(EdbTrackP*)eTracks.At(n):0;}

GetVertex()

EdbVertex* EdbDataStore::GetVertex ( int  n )

inline

{return (n<eVTX.GetEntries())?(EdbVertex*)eVTX.At(n):0;}

LoadMCVertices()

void EdbDataStore::LoadMCVertices

(

TObjArray *

vtx

)

restore MC info methods

---

                                               {
  Log(2,"EdbDataStore::LoadMCVertices",Form("Reading %d vertices",vtx->GetEntries()));
  for(int nv=0; nv<vtx->GetEntries();++nv){
    AddVertex((EdbVertex*)vtx->At(nv));
  }
}

MakePattern()

void EdbDataStore::MakePattern	(	double	z,
		int	plate,
		int	side
	)

---

                                                         {
  Log(3,"EdbDataStore::MakePattern()",Form("plt#%d (side %d) z=%2.1f",plate,side,z));
  EdbPattern* p=new EdbPattern(0,0,z,0);
  p->SetID(plate);
  p->SetScanID(EdbID(0,plate,0,0));
  p->SetSide(side);
  p->SetPID(side?(plate*10+side):plate);
  if(side==0)eSegPV.AddPattern(p);
  else eRawPV.AddPattern(p);
   
}

Nplt()

int EdbDataStore::Nplt ( )

inline

{return eRawPV.Npatterns()/2;}

Nt()

int EdbDataStore::Nt ( )

inline

count methods:

{return eTracks.GetEntries();}

Nv()

int EdbDataStore::Nv ( )

inline

{return eVTX.GetEntries();}

PrintBrief()

void EdbDataStore::PrintBrief

(

)

print methods

---

                             {
  printf("EdbDataStore. We have here:\n");
  printf(" --%d vertcies\n --%d tracks\n --%d Plates (geometry)\n",eVTX.GetEntries(),eTracks.GetEntries(),(eBrick)?eBrick->Npl():0);
  printf(" --%d BT Patterns\n --%d MT Patterns\n",eSegPV.Npatterns(),eRawPV.Npatterns());
}

PrintPatterns()

void EdbDataStore::PrintPatterns

(

)

---

                                {
  printf("EdbDataStore. We have %d BT patterns:\n",eSegPV.Npatterns());
  EdbPattern* pat=0;
  for(int np=0;np<eSegPV.Npatterns();++np){
    pat=GetSegPat(np);
    printf(" (pat #%3d) [X,Y,Z]=[% 7.2f % 7.2f % 7.2f]\t Nseg=%d\n",
       pat->Plate(),pat->X(),pat->Y(),pat->Z(),pat->N());
  }
}

PrintTracks()

void EdbDataStore::PrintTracks

(

int

vlev = 0

)

---

                                      {
  printf("EdbDataStore. We have %d tracks\n",eTracks.GetEntries());
  if(vlev<1)return;
  EdbTrackP* trk=0;
  for(int nt=0;nt<eTracks.GetEntries();++nt){
    trk=GetTrack(nt);
    printf(" (trk #%3d) PDG=% 4d P=% 4.2f [X,Y,Z]=[% 7.2f % 7.2f % 7.2f]\t[TX,TY]=[% 2.2f % 2.2f] Nseg=%d\n",
    trk->ID(),trk->PDG(),trk->P(),trk->X(),trk->Y(),trk->Z(),trk->TX(),trk->TY(),trk->N());
    if(vlev<2)continue;
    for(int ns=0;ns<trk->N();++ns){
      EdbSegP* s=trk->GetSegment(ns);
      printf("    seg#%d  Plate=%d pid=%d flag=%d\n",s->ID(),s->Plate(),s->PID(),s->Flag());
    }
  }
}

Restore_PatFromGeom()

void EdbDataStore::Restore_PatFromGeom	(	int	np0 = 0,
		int	np1 = 1000
	)

---

fill patterns array from eBrick plates

                                                      {
  assert(eBrick);
  EdbPlateP* plt=0;
  for(int np=0;np<eBrick->Npl();++np){
    if(np<np0 || np>np1)continue;
    plt=eBrick->GetPlate(np);
    assert(plt);
    for(int side=0;side<3;++side)
      MakePattern(plt->Z()+plt->GetLayer(side)->Z(),np,side);
  }
}

Restore_PIDFromID()

void EdbDataStore::Restore_PIDFromID

(

)

---

restore patterns' "plate" and "side" from their ID

                                    {
  EdbPattern* p;
  int plate, side;
  for(int np=0;np<eSegPV.Npatterns();np++){
    p=eSegPV.GetPattern(np);
    plate=p->ID();
    p->SetPID(p->ID());
    p->SetScanID(EdbID(0,plate,0,0));
    p->SetSide(0);
  }
  for(int np=0;np<eRawPV.Npatterns();np++){
    p=eRawPV.GetPattern(np);
    plate=p->ID()/2;
    p->SetPID(p->ID());
    side=1+p->ID()%2;
    p->SetScanID(EdbID(0,plate,0,0));
    p->SetSide(side);
  }
};

Restore_SegFromTrx()

void EdbDataStore::Restore_SegFromTrx	(	EdbSegmentCut *	cut = 0,
		int	Plt0 = 0,
		int	Plt1 = 1000
	)

---

restore segments from tracks

                                                                          {
  EdbTrackP* t=0;
   EdbSegP* s=0;
  for(int nt=0;nt<eTracks.GetEntries();++nt){
    t=GetTrack(nt);
    if(t->N()==0)continue;
    for(int ns=0;ns<t->N();++ns){
      s=t->GetSegment(ns);
      s->SetMC(t->MCEvt(),t->MCTrack());
//       s->SetFlag(1);
      AddSegment(s,cut,Plt0,Plt1);
      }
  }
}

Restore_TrxFromVtx()

void EdbDataStore::Restore_TrxFromVtx

(

)

---

restore tracks from vertex

                                     {
  ClearTracks();
  EdbVertex* v;
  for(int nv=0;nv<eVTX.GetEntries();nv++){
    v=GetVertex(nv);
    for(int nt=0;nt<v->N();nt++){
      EdbTrackP* trk=v->GetTrack(nt);
      trk->AddVTA(v->GetVTa(nt));
      if(v->GetVTa(nt)->Zpos()==1)AddTrack(trk);
    }
  }
}

SavePlateToRaw()

void EdbDataStore::SavePlateToRaw	(	const char *	fname,
		int	PID,
		Option_t *	option = "RECREATE"
	)

---

Convert this event to RAW data format!

Get plate geometry:

---

maximum segments per view. When this is reached - create a new view.

                                                                           {
    EdbSegP* seg=0;
    EdbSegment* seg1=0;
    EdbPattern* ptop=FindPattern(PID,2);
    EdbPattern* pbot=FindPattern(PID,1);
    if(ptop->N()==0 && pbot->N()==0)return;
    
    EdbPlateP* p=eBrick->GetPlate(PID);
    if(strcmp(option,"UPDATE")==0){
      //test if file exists
      FILE *file= fopen(fname, "r");
      if (file) fclose(file);
      else{ 
        option="RECREATE";}
    }
    EdbRun run(fname,option);
    Log(1,"EdbDataStore::SavePlateToRaw",Form("Open file \"%s\" opt=\"%s\"",fname,option));
    Log(2,"EdbDataStore::SavePlateToRaw","created run");
    Log(2,"EdbDataStore::SavePlateToRaw",Form("cycle through %d TOP segs",ptop->N()));
    int Vid=0;
    const int NsegPerView=1000; 
    int Ns=0;
    run.GetView()->SetNframes(1,0);
    run.GetView()->GetHeader()->SetViewID(++Vid);
    for(int ns=0;ns<ptop->N();++ns){
      if(Ns==NsegPerView){
    Ns=0;
    run.AddView();
    run.GetView()->Clear();
    run.GetView()->GetHeader()->SetViewID(++Vid);
      }
      seg=ptop->GetSegment(ns);
      if(seg->W()<0)continue;
      seg1=run.GetView()->AddSegment(p->Xp(*seg),p->Yp(*seg),seg->Z(),p->TXp(*seg),p->TYp(*seg),seg->DZ(),0,(int)(seg->W()),seg->ID());
      seg1->SetSigma(seg->P(),1);
      ++Ns;
    }
    run.AddView();
    run.GetView()->Clear();
 
    Log(2,"EdbDataStore::SavePlateToRaw",Form("cycle through %d BOT segs",pbot->N()));
    run.GetView()->SetNframes(0,1);
    run.GetView()->GetHeader()->SetViewID(++Vid);
    Ns=0;
    for(int ns=0;ns<pbot->N();++ns){
      if(Ns==NsegPerView){
    Ns=0;
    run.AddView();
    run.GetView()->Clear();
    run.GetView()->GetHeader()->SetViewID(++Vid);
      }
      seg=pbot->GetSegment(ns);
      if(seg->W()<0)continue;
      seg1=run.GetView()->AddSegment(p->Xp(*seg),p->Yp(*seg),seg->Z(),p->TXp(*seg),p->TYp(*seg),seg->DZ(),1,(int)(seg->W()),seg->ID());      
      seg1->SetSigma(seg->P(),1);
      ++Ns;
    }
    run.AddView();
    run.GetView()->Clear();
    Log(2,"EdbDataStore::SavePlateToRaw","Close plate");
    run.Save();
    run.Close();    
};

SaveToRaw()

void EdbDataStore::SaveToRaw	(	const char *	dir = "./",
		const EdbID &	idset = "0.0.0.0",
		Option_t *	option = "RECREATE",
		bool	doaff = true
	)

save methods:

---

Save the raw data to FEDRA brick structure

generate plate filename

save plate.

save affine tranformations

save affine parameters

                                                                                           {
  EdbScanProc sp;
 
  sp.eProcDirClient=dir;
  sp.CheckBrickDir(idset);
  sp.CheckAFFDir(idset.eBrick);
 
  EdbID PL=idset;
  EdbID PL0=idset;
  EdbPlateP* p=0;
 
  float z,z0;
  for(int i=0;i<eBrick->Npl();i++){
    p=eBrick->GetPlate(i);
    z=p->Z();
    PL.ePlate=i+1;
    sp.CheckPlateDir(PL);
    TString fnm;
    sp.MakeFileName(fnm,PL,"raw.root",true); 
    SavePlateToRaw(fnm.Data(),i,option); 
    if(doaff && i>0){
      sp.MakeAffName(fnm,PL,PL0);
      FILE* afF=fopen(fnm.Data(),"w"); 
      fprintf(afF,"ZLAYER 0   %4f 0 0\n",z0-z);
      fclose(afF);
    }
    PL0=PL;
    z0=z;
  }
}

SetOwnTracks()

void EdbDataStore::SetOwnTracks ( bool  own = true )

inline

setown methods:

{eTracks.SetOwner(own);}

SetOwnTrkSegs()

void EdbDataStore::SetOwnTrkSegs ( )

inline

{for(int nt=0;nt<Nt();++nt)GetTrack(nt)->SetOwner();}

SetOwnVertices()

void EdbDataStore::SetOwnVertices ( bool  own = true )

inline

{eVTX.SetOwner(own);}

TransferGeometry()

void EdbDataStore::TransferGeometry

(

EdbDataStore *

ds

)

---

                                                   {
  Log(2,"transfer geometry","");
  ds->eBrick->Clear();
  ds->eBrick->Copy(*eBrick);
}

TransferSegs()

void EdbDataStore::TransferSegs ( EdbPatternsVolume *  pv0, EdbPatternsVolume *  pv1, EdbSegmentCut *  cut = 0, int  FromPlate = 0, int  ToPlate = 57  )

static

---

find corresponding DEST pattern

                                                                                                                          {
  Log(1,"EdbDataStore::TransferSegs()","Transfer segments");
  assert(pv0);
  assert(pv1);
  EdbPattern* pat0=0;
  EdbPattern* pat1=0;
  EdbSegP* seg=0;
  int npat=pv0->Npatterns();
 
  if(pv1->Npatterns()==0)pv0->PassProperties(*pv1);
 
  int nseg=0;
  float p[5];
  for(int np=0;np<npat;++np){
 
    pat0=pv0->GetPattern(np);
    if(pat0->Plate()<FromPlate)continue;
    if(pat0->Plate()>ToPlate)continue;
    pat1=pv1->GetPattern(np); 
    assert(pat1);
    assert(pat1->Z()==pat0->Z());
 
    nseg=pat0->N();
    for(int ns=0;ns<nseg;++ns){
      seg=pat0->GetSegment(ns);
      p[0]=seg->X();
      p[1]=seg->Y();
      p[2]=seg->TX();
      p[3]=seg->TY();
      p[4]=seg->W();
      if(cut && !cut->PassCut(p))continue;
      seg=pat1->AddSegment(*seg);
      seg->SetTrack(0);
    };
//     printf("dest.pat#%d has #%d segments\n",pat1->ID(),pat1->N());
    pat1->SetID(np);
    pat1->SetSegmentsPID();
    
  }
  
}

TransferTo()

void EdbDataStore::TransferTo	(	EdbDataStore *	ds,
		char	level,
		EdbSegmentCut *	cut = 0,
		int	FromPlate = 0,
		int	ToPlate = 57
	)

transfer methods

---

                                                                                                       {
  assert(level>=0 && level<16);
  
  if(level&0x1){
    Log(2,"EdbDataStore::TransferTo","transfer MTK segments\n");
    TransferSegs(&eRawPV,&(ds->eRawPV),cut,FromPlate,ToPlate);
    Log(2,"EdbDataStore::TransferTo","done");
  }
  
  if(level&0x2){
    Log(2,"EdbDataStore::TransferTo","transfer BTK segments\n");
    TransferSegs(&eSegPV,&(ds->eSegPV),cut,FromPlate,ToPlate);
    Log(2,"EdbDataStore::TransferTo","done");
  }
  
  if(level&0x4){
    Log(2,"EdbDataStore::TransferTo","transfer Tracks\n");
    for(int nt=0;nt<eTracks.GetEntries();++nt){
      if(GetTrack(nt)->N())ds->AddTrack(GetTrack(nt));
    }
    ds->SetOwnTracks(0);
    Log(2,"EdbDataStore::TransferTo","done");
  }
  if(level&0x8){
    Log(2,"EdbDataStore::TransferTo","transfer Vertices\n");
    for(int nv=0;nv<eVTX.GetEntries();++nv)ds->AddVertex(GetVertex(nv));
    ds->SetOwnVertices(0);
    Log(2,"EdbDataStore::TransferTo","done");
  }
}

Member Data Documentation

eBrick

EdbBrickP* EdbDataStore::eBrick

eRawPV

EdbPatternsVolume EdbDataStore::eRawPV

geometry

eSegPV

EdbPatternsVolume EdbDataStore::eSegPV

eTracks

TObjArray EdbDataStore::eTracks

eVTX

TObjArray EdbDataStore::eVTX

---

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

• src/libEMC/EdbDataStore.h
• src/libEMC/EdbDataStore.cxx