EdbDSRec Class Reference

#include <EdbDataStore.h>

Inheritance diagram for EdbDSRec:

Collaboration diagram for EdbDSRec:

Public Member Functions
void	Clear (bool hard=false)
int	DoDecaySearch ()
int	DoFindBlkSeg (EdbVertex *v, int w0, double ImpMax=50., double RMax=3000, int Dpat=1)
	TODO. More...
int	DoMomEst ()
int	DoTracking (bool use_btk=true, int p0=0, int p1=100)
int	DoTracking0 (bool use_btk=true, int p0=0, int p1=100)
int	DoVertexing ()
	EdbDSRec ()
void	FillECovPV (EdbPatternsVolume *, EdbScanCond *cnd=0)
	prepare segments' cov matrix More...
void	FillECovSeg (EdbSegP *seg, EdbScanCond *cnd=0)
void	FillECovTrks ()
void	FillErrorsCOV ()
	~EdbDSRec ()
Public Member Functions inherited from EdbDataStore
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 ()

Public Attributes
EdbScanCond	eCond_b
EdbScanCond	eCond_m
EdbMomentumEstimator	eMomEst
EdbVertexRec	eVRec
Public Attributes inherited from EdbDataStore
EdbBrickP *	eBrick
EdbPatternsVolume	eRawPV
	geometry More...
EdbPatternsVolume	eSegPV
TObjArray	eTracks
TObjArray	eVTX

Additional Inherited Members
Static Public Member Functions inherited from EdbDataStore
static void	TransferSegs (EdbPatternsVolume *pv0, EdbPatternsVolume *pv1, EdbSegmentCut *cut=0, int FromPlate=0, int ToPlate=57)

Constructor & Destructor Documentation

EdbDSRec()

EdbDSRec::EdbDSRec

(

)

---

init Vertexing:

init momentum estimator:

                  {
  eVRec.eEdbTracks=&eTracks;
  eVRec.eDZmax=5000;
  eVRec.eProbMin=0.001;
  eVRec.eImpMax=50;
  eVRec.eUseMom=0;
  eVRec.eUseSegPar=1;
  eVRec.eQualityMode=0;
  eMomEst.SetParPMS_Mag();
  eMomEst.eAlg = 0;
  eMomEst.eMinEntr = 3;
}

~EdbDSRec()

EdbDSRec::~EdbDSRec ( )

inline

{};

Member Function Documentation

Clear()

void EdbDSRec::Clear

(

bool

hard = false

)

---

                             {
  EdbDataStore::Clear(hard);
  eVRec.Reset();
}

DoDecaySearch()

int EdbDSRec::DoDecaySearch

(

)

---

                           {
  return 0;
}

DoFindBlkSeg()

int EdbDSRec::DoFindBlkSeg	(	EdbVertex *	v,
		int	w0,
		double	ImpMax = 50.,
		double	RMax = 3000,
		int	Dpat = 1
	)

TODO.

---

segment assigned to track

not black enough

                                                                                   {
  int nn=0;
  v->ClearNeighborhood();
  EdbVTA* vta=0;
  printf("start\n");
  EdbPattern* pat=0;
  EdbSegP* seg=0;
  double dzmax=Dpat*1300;
  double dx,dy,dz;
  double imp=0;
  for(int np=0;np<eSegPV.Npatterns();++np){
    pat=eSegPV.GetPattern(np);
    dz=fabs(pat->Z()-v->Z());
    if(dz<0 || dz>dzmax)continue;
    for(int ns=0;ns<pat->N();++ns){
      seg=pat->GetSegment(ns);
      if(seg->Track()!=-1)continue; 
      if(seg->W()<w0)continue; 
      imp=v->DistSeg(seg);
      printf("imp =%2.1f\n",imp);
      if (imp > ImpMax) continue;
      dx=seg->X()-v->X();
      dy=seg->Y()-v->Y();
      vta = new EdbVTA((EdbTrackP *)seg, v);
      vta->SetZpos(1);
      vta->SetFlag(1);
      vta->SetImp(imp);
      vta->SetDist(sqrt(dx*dx+dy*dy+dz*dz));
      v->AddVTA(vta);
      ++nn;
    }
  }
  return nn;
}

DoMomEst()

int EdbDSRec::DoMomEst

(

)

---

                      {
  
  float p,p0,p1;
  float tl,tt;
  EdbTrackP* trk=0;
  EdbTrackP* clon=0;
  int Nest=0;
  for(int nt=0;nt<eTracks.GetEntries();++nt){
    trk=GetTrack(nt);
    clon=(EdbTrackP*)trk->Clone();
    Log(1,"EdbDSRec::DoMomEst",Form("Mom Est: track #%d (%dseg)",nt,trk->N()));
//     if(trk->N()<5){trk->SetP(0); printf("track #%d (%dseg) - skip!\n",nt,trk->N());continue;}
    p=eMomEst.PMSang(*clon);
//     if(eMomEst.eStatus==-1){printf("track #%d (%dseg) - skip!\n",nt,trk->N());continue;}
    printf("track #%d (%dseg) mom %f - status %d\n",nt,trk->N(),p,eMomEst.eStatus);
    p0=eMomEst.ePmin;
    p1=eMomEst.ePmax;
/*    tl=eMomEst.eGX->GetY()[0];
    tt=eMomEst.eGY->GetY()[0]; */
//     printf("track #%d (%dseg) mom %f [%2.2f <> %2.2f]\n",nt,trk->N(),p,p0,p1);
//     printf("dTL=%2.4g, dTT=%2.4g\n",tl,tt);
//     if(p<0 && (tl>0.015 || tt>0.015)){
//       printf("BAD!\n");
//       trk->SetP(-p);
//       continue;
//     }
//     if(p<0)p=0;
    trk->SetP(p);
    trk->SetPerr(p0,p1);
    delete clon;
  }
  return Nest;
}

DoTracking()

int EdbDSRec::DoTracking	(	bool	use_btk = true,
		int	p0 = 0,
		int	p1 = 100
	)

---

init tracker:

doing tracking using microtracks or basetracks

                                                    {
  EdbTrackAssembler tracker;
  tracker.eCond=use_btk?eCond_b:eCond_m;
//   .SetSigma0(3,3,0.005,0.005);
  tracker.eDTmax=0.07;
  tracker.eDRmax=45.;
//   tracker.eDZGapMax=10000;
  tracker.InitTrZMap(  2400, eBrick->Xmin(), eBrick->Xmax(),   2000, eBrick->Ymin(), eBrick->Ymax(),   30 );
  int nsegmin=use_btk?2:3;
  SetOwnTracks(1);
  EdbPatternsVolume* pv=use_btk?(&eSegPV):(&eRawPV);
  int npl=pv->Npatterns();
  EdbPattern *pat=0;
  EdbLayer *plate=0;
  for(int ipass=0; ipass<2; ipass++) {
    printf("\n\n*************** ipass=%d ************\n",ipass);
    tracker.eCollisionsRate=0;
    for(int i=0; i<npl; i++) {
      if(i<p0 || i>p1)continue;
      pat=pv->GetPattern(i);
      if(pat->N()==0)continue;
      Log(3,Form("EdbDSRec::DoTracking(%d,%d,%d)",use_btk,p0,p1),Form("plate=%d at Z=%4.1f side=%d [%d segm]\n",pat->Plate(),pat->Z(),pat->Side(),pat->N()));
      plate = eBrick->GetPlate(pat->Plate());
//       printf("pattern z = (%2.4f vs %2.4f)\n", pat->Z(),plate->Z());
      pat->SetSegmentsZ();
      pat->SetSegmentsPID();
      pat->Transform(    plate->GetAffineXY()   );
      pat->TransformShr( plate->Shr() );
      pat->TransformA(   plate->GetAffineTXTY() );
      pat->SetSegmentsPlate(pat->Plate());
//       printf("flag0=%d\n",pat->GetSegment(0)->Flag());
      
      //       pat=ds_rec->eRawPV.GetPattern(i);
      
//       printf("pat #%d has %d segments\n",i,pat->N());
      if(i>0) tracker.ExtrapolateTracksToZ(pat->Z());
      tracker.FillTrZMap();
      tracker.AddPattern(*pat);
//       printf("ok\n");
    }
  }
 
  int Ntr=0;
  int ntr = tracker.Tracks().GetEntries();
  for( int i=0; i<ntr; i++ )     {
    EdbTrackP *t = (EdbTrackP*)(tracker.Tracks().At(i));
    int nseg=t->N();
//     printf("trk %d/%d: nseg=%d (min %d)\n",i,ntr,nseg);
    if(nseg>=nsegmin) {
      for(int j=0; j<nseg; j++) {
    EdbSegP *s = t->GetSegment(j);
    s->SetErrors();
    tracker.eCond.FillErrorsCov(s->TX(),s->TY(),s->COV());
      }
      t->SetP(1.);
      t->FitTrackKFS(0);
      //fit.FitTrackLine(*t);
      tracker.RecalculateSegmentsProb(*t);
      t->SetCounters();
      AddTrack(new EdbTrackP(*t));
      
      Ntr++;
    }
//     else
    //if(t->N()>2)   t->PrintNice();
  }
  return Ntr;
}

DoTracking0()

int EdbDSRec::DoTracking0	(	bool	use_btk = true,
		int	p0 = 0,
		int	p1 = 100
	)

---

transfer patterns

propagate

add track here

                                                    {
  float  momentum = 0.3;
  float  mass     = 0.139;  // particle mass
  float  ProbMinP = 0.001;   // minimal probability to accept segment on propagation
  int    nsegmin  = 1;      // minimal number of segments to propagate this track
  int    ngapmax  = 3;      // maximal gap for propagation
  
  EdbPVRec gAli;
  EdbPatternsVolume* pv=use_btk?&eSegPV:&eRawPV;
  
  EdbPattern* pat,*pat1;
  for(int np=0; np<pv->Npatterns(); ++np){
    pat=pv->GetPattern(np);
    if(pat->Plate()<p0 || pat->Plate()>p1)continue;
    pat1=(EdbPattern*)pat->Clone();
    pat1->SetSegmentsZ();
    pat1->SetID(np);
    pat1->SetSegmentsPID();
    Log(1,"EdbDSRec",Form("add pattern %d with %d segments",pat1->PID(),pat1->N()));
    gAli.AddPattern(pat1);
  }
  
  gAli.SetScanCond(use_btk?&eCond_b:&eCond_m);
  gAli.SetCouplesPeriodic(0,1);
  int ntr = EdbDataProc::LinkTracksWithFlag( &gAli, momentum, ProbMinP, nsegmin, ngapmax, 0 );
  gAli.FitTracks( momentum, mass );
  Log(2,"EdbDSRec::DoTracking0",Form("ntr = %d\n",ntr));
  
  EdbTrackFitter tf;
  
  EdbTrackP *tr=0;
  ntr = gAli.eTracks->GetEntries();
  
  for(int i=0; i<ntr; i++) {
    tr = (EdbTrackP*)(gAli.eTracks->At(i));
    tr->SetID(i);
    tr->SetSegmentsTrack();
    tr->SetFlag(0);
  }
  
  int nadd = 0;
  int nseg=0;
  int Ntr=0;
//   nadd=gAli.PropagateTracks(0,57,0.001,3,0);
  for(int i=0; i<ntr; i++) {
    tr = (EdbTrackP*)(gAli.eTracks->At(i));
 
    float p=momentum;
    
    tr->SetErrorP(0.2*0.2*p*p);
    
    nseg = tr->N();
    tr->SetP(p);
    printf("pid=%d\n",tr->GetSegmentFirst()->PID());
    if(tr->Flag()<0) continue;
    //     tr->PrintNice();
    nadd += gAli.PropagateTrack( *tr, true, 0.001, 3, 0 );
    if(tr->Flag()<0) printf("%d flag=%d\n",i,tr->Flag());
    //if(tr->N() != nseg) printf("%d nseg=%d (%d) \t p = %f\n",i,tr->N(),nseg,tr->P());
    AddTrack((EdbTrackP*)tr->Clone());
    Ntr++;
  }
  printf("nadd = %d\n",nadd);
      
  return Ntr;
}

DoVertexing()

int EdbDSRec::DoVertexing

(

)

---

                         {
  SetOwnVertices(1);
  // performing vertexing
  if(gEDBDEBUGLEVEL>=1) printf("%d tracks for vertexing\n",  eVRec.eEdbTracks->GetEntries() );
  int nvtx = eVRec.FindVertex();
  if(gEDBDEBUGLEVEL>=1) printf("%d 2-track vertexes was found\n",nvtx);
  if(nvtx == 0) return 0;
  eVRec.ProbVertexN();
  for(int nv=0;nv<eVRec.Nvtx();nv++){
    EdbVertex* v=eVRec.GetVertex(nv);
    if(v->Flag()<0)delete v;
    else eVTX.Add(v);
  }
  return eVRec.Nvtx();
}

FillECovPV()

void EdbDSRec::FillECovPV	(	EdbPatternsVolume *	pv,
		EdbScanCond *	cnd = 0
	)

prepare segments' cov matrix

                                                               {
  assert(pv);
  EdbPattern* pat=0;
  EdbSegP* seg;
  for(int i=0;i<pv->Npatterns();i++){
    pat=pv->GetPattern(i);
    pat->SetSegmentsPID();
    for(int j=0;j<pat->GetN();j++){
      seg=pat->GetSegment(j);
      seg->SetID(j);
      FillECovSeg(seg);
    }
  }
}

FillECovSeg()

void EdbDSRec::FillECovSeg	(	EdbSegP *	seg,
		EdbScanCond *	cnd = 0
	)

---

if cnd=0 - decide which to use

                                                       {
  assert(seg);
  seg->SetErrors();
  seg->SetErrorP(1.);
  if(cnd==0) 
    cnd=(seg->DZ()>100)?&eCond_b:&eCond_m;
  cnd->FillErrorsCov(seg->TX(),seg->TY(),seg->COV());
}

FillECovTrks()

void EdbDSRec::FillECovTrks

(

)

---

                           {
    EdbTrackP* trk=0;
    TObjArray trx1;
    for(int nt=0;nt<Nt();nt++){
      trk=GetTrack(nt);
      for(int ns=0;ns<trk->N();ns++) FillECovSeg(trk->GetSegment(ns));
      trk->FitTrackKFS(0);
    }
}

FillErrorsCOV()

void EdbDSRec::FillErrorsCOV

(

)

---

                            {
  FillECovPV(&eRawPV,&eCond_m);
  FillECovPV(&eSegPV,&eCond_b);
  FillECovTrks();
}

Member Data Documentation

eCond_b

EdbScanCond EdbDSRec::eCond_b

eCond_m

EdbScanCond EdbDSRec::eCond_m

eMomEst

EdbMomentumEstimator EdbDSRec::eMomEst

eVRec

EdbVertexRec EdbDSRec::eVRec

---

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

• /home/antonio/fedra_doxygen/src/libEMC/EdbDataStore.h
• /home/antonio/fedra_doxygen/src/libEMC/EdbDataStore.cxx