EdbVertexG Class Reference

Vertex. More...

#include <EdbDisplay.h>

Inheritance diagram for EdbVertexG:

Collaboration diagram for EdbVertexG:

Public Member Functions
virtual void	DeleteVertex ()
virtual void	DumpVertex ()
	EdbVertexG ()
	EdbVertexG (EdbDisplay *D)
virtual const char *	GetName () const
virtual char *	GetObjectInfo (int px, int py) const
virtual const char *	GetTitle () const
virtual void	InspectVertex ()
virtual void	RemoveKink ()
virtual void	SetAsWorking ()
void	SetVertex (EdbVertex *v)
virtual void	TestVertex ()
virtual	~EdbVertexG ()

Private Attributes
EdbDisplay *	eD
EdbVertex *	eV

Detailed Description

Vertex.

Constructor & Destructor Documentation

EdbVertexG() [1/2]

EdbVertexG::EdbVertexG ( )

inline

:TPolyMarker3D(1) {eV=0; eD=0;}

EdbVertexG() [2/2]

EdbVertexG::EdbVertexG ( EdbDisplay *  D )

inline

:TPolyMarker3D(1) {eV=0; eD=D;}

~EdbVertexG()

virtual EdbVertexG::~EdbVertexG ( )

inlinevirtual

{}

Member Function Documentation

DeleteVertex()

void EdbVertexG::DeleteVertex ( )

virtual

{
  EdbDisplay *eDs = 0;
  EdbVertex  *eVs = 0;
  eDs = eD;
  eVs = eV;
  if (eDs && eVs)
  {
    if (eDs->eWait_Answer) return;
    if (eDs->eWorking == eVs)
    {
    eDs->DialogModifiedVTX();
    return;
    }
    if (eDs->eVertex == eVs)
    {
        eDs->CancelModifiedVTX();
    }
    if (eDs->eArrV) 
    {
    eDs->eArrV->Remove((TObject *)eVs);
    eDs->eArrV->Compress();
    eDs->Draw();
    }
    else
    {
    if ((eDs->eArrV)->FindObject(eVs))
    {
        eDs->eArrV->Remove(eVs);
        eDs->eArrV->Compress();
        eDs->eVertex = 0;
        eDs->Draw();
    }
    }
    delete eVs;
    eDs->Draw();
  }
}

DumpVertex()

void EdbVertexG::DumpVertex ( )

virtual

{
  if (eV) eV->Print();
}

GetName()

const char * EdbVertexG::GetName ( ) const

virtual

{
    static char name[] = "Vertex";
    return name;
}

GetObjectInfo()

char * EdbVertexG::GetObjectInfo ( int  px, int  py  ) const

virtual

{
    static char coordinates[80];
    if (eV->V())
    sprintf(coordinates, "X = %.1f, Y = %.1f, Z = %.1f", eV->VX(), eV->VY(), eV->VZ());
    else
    sprintf(coordinates, "X = %.1f, Y = %.1f, Z = %.1f", eV->X(), eV->Y(), eV->Z());
    return coordinates;
}

GetTitle()

const char * EdbVertexG::GetTitle ( ) const

virtual

{
    static char title[80];
    if (eV->V())
    sprintf(title, "Vertex ID %d, Prongs %d, Prob %f", eV->ID(), eV->N(), eV->V()->prob());
    else
    sprintf(title, "Vertex ID %d, Prongs %d, Prob %f", eV->ID(), eV->N(), 1.);
    return title;
}

InspectVertex()

void EdbVertexG::InspectVertex ( )

virtual

{
  if (eV) eV->Inspect();
}

RemoveKink()

void EdbVertexG::RemoveKink ( )

virtual

{
  EdbDisplay *eDs = 0;
  EdbVertex  *eVs = 0;
  eDs = eD;
  eVs = eV;
  if (eDs && eVs)
  {
    if (eDs->eWait_Answer) return;
    if (eVs->N() != 2)
    {
    printf("Wrong vertex type - not a kink-like!\n");
    fflush(stdout);
    return;
    }
    EdbTrackP *etr = 0;
    for (int it=0; it<eVs->N(); it++)
    {
    etr = eVs->GetTrack(it);
    if (etr->VTAS() && etr->VTAE())
    {
        if (etr->VertexS() == eVs && etr->VertexE()->Flag() >= 0)
        {
        printf("Vertex track (at zmax) belong to another vertex (ID=%d) too - impossible delete it!\n", etr->VertexE()->ID());
        fflush(stdout);
        return;
        }
        if (etr->VertexE() == eVs && etr->VertexS()->Flag() >= 0)
        {
        printf("Vertex track (at zmin) belong to another vertex (ID=%d) too - impossible delete it!\n", etr->VertexS()->ID());
        fflush(stdout);
        return;
        }
    }
    }
    if (eDs->eWorking == eVs)
    {
    eDs->DialogModifiedVTX();
    return;
    }
    if (eDs->eVertex == eVs)
    {
        eDs->CancelModifiedVTX();
    }
    if (eDs->eArrV) 
    {
    eDs->eArrV->Remove((TObject *)eVs);
    eDs->eArrV->Compress();
    eDs->Draw();
    }
    else
    {
    if ((eDs->eArrV)->FindObject(eVs))
    {
        eDs->eArrV->Remove(eVs);
        eDs->eArrV->Compress();
        eDs->eVertex = 0;
        eDs->Draw();
    }
    }
    EdbSegP *seg = 0;
    if (!eD->eTrack) eD->eTrack = new EdbTrackP();
    for (int it=0; it<eVs->N(); it++)
    {
    etr = eVs->GetTrack(it);
        for(int is=0; is<etr->N(); is++) {
        seg = etr->GetSegment(is);
        (eD->eTrack)->AddSegment((EdbSegP *)seg);
    }
    etr->SetFlag(-10);
    etr->SetSegmentsTrack(-2-(etr->ID()+1));
    if(eD->eArrTr) if((eD->eArrTr)->FindObject(etr))
    {
        eD->eArrTr->Remove(etr);
        eD->eArrTr->Compress();
        if (!(eD->eArrSegP)) eD->eArrSegP = new TObjArray();
            for(int is=0; is<etr->N(); is++) {
        seg = etr->GetSegment(is);
        if(eD->eArrSegP) if(!((eD->eArrSegP)->FindObject(seg))) eD->eArrSegP->Add(seg);
        }
    }
    if(eD->eArrTrSave) if((eD->eArrTrSave)->FindObject(etr))
    {
        eD->eArrTrSave->Remove(etr);
        eD->eArrTrSave->Compress();
        if (!(eD->eArrSegPSave)) eD->eArrSegPSave = new TObjArray();
            for(int is=0; is<etr->N(); is++) {
        seg = etr->GetSegment(is);
        if(eD->eArrSegPSave) if(!((eD->eArrSegPSave)->FindObject(seg))) eD->eArrSegPSave->Add(seg);
        }
    }
    }
    if(eD->eArrTr) eD->eArrTr->Compress();
    if(eD->eArrTrSave) eD->eArrTrSave->Compress();
    eVs->SetFlag(-11);
    eDs->Draw();
    float mass = etr->M();
    float momentum = etr->P();
    if (eD->eP > 0.) momentum = eD->eP;
    if (eD->eM > 0.) mass = eD->eM;
    (eD->eTrack)->SetM(mass);
    (eD->eTrack)->SetP(momentum);
    float X0 = 0.;
    if (eD->eVerRec) if ((eD->eVerRec)->ePVR) X0 = (((eD->eVerRec)->ePVR)->GetScanCond())->RadX0();
    (eD->eTrack)->FitTrackKFS(true, X0, 0);
    if (!eD->eArrTr) eD->eArrTr = new TObjArray();
    (eD->eArrTr)->Add(eD->eTrack);
    if (eD->eArrTrSave) (eD->eArrTrSave)->Add(eD->eTrack);
    eD->Draw();
  }
}

SetAsWorking()

void EdbVertexG::SetAsWorking ( )

virtual

{
  char text[512];
  EdbTrackP *tr = 0;
  EdbDisplay *eDs = 0;
  EdbVertex  *eVs = 0;
  eDs = eD;
  eVs = eV;
  if (eDs && eVs)
  {
    if (eDs->eVertex == eVs) return;
    if (eDs->eWait_Answer) return;
    if (eDs->eWorking)
    {
    eDs->DialogModifiedVTX();
    return;
    }
    if (eDs->eVertex)
    {
        eDs->CancelModifiedVTX();
    }
    if (eDs->eSegPM)
    { 
    eDs->ClearSegmentEnv();
    }
    SetMarkerColor(kGreen);
    if (!(eDs->eArrV)) 
    {
    eDs->eArrV = new TObjArray(20);
    eDs->eArrV->Add((TObject *)eVs);
    if (!(eDs->eArrTr))  eDs->eArrTr = new TObjArray(20);
    for (int i=0; i<eVs->N(); i++)
    {
        tr = eVs->GetTrack(i);
        if(!(eDs->eArrTr->FindObject(tr))) eDs->eArrTr->Add(tr);
    }
    eDs->Draw();
    }
    else
    {
    if (!((eDs->eArrV)->FindObject(eVs)))
    {
        eDs->eArrV->Add(eVs);
        if (!(eDs->eArrTr))  eDs->eArrTr = new TObjArray(20);
        for (int i=0; i<eVs->N(); i++)
        {
        tr = eVs->GetTrack(i);
        if(!(eDs->eArrTr->FindObject(tr))) eDs->eArrTr->Add(tr);
        }
        eDs->eVertex = 0;
        eDs->Draw();
    }
    }
    eDs->eSegment = 0;
    eDs->eVertex = eVs;
    eDs->ePrevious = 0;
    eDs->CreateCanvasVTX();
    eVs->V()->rmsDistAngle();
    sprintf(text,"Orig    %-4d  %-4d  %-8.1f   %-8.1f   %-8.1f   %-6.1f %-7.1f  %-7.5f",
    eVs->ID(), eVs->N(), eVs->VX(), eVs->VY(), eVs->VZ(), eVs->V()->dist(),
    eVs->V()->chi2()/eVs->V()->ndf(), eVs->V()->prob());
    eDs->DrawOldVTX(text);
    eDs->DrawVTXTracks("Original", eDs->eVertex);
    eDs->DrawEnv();
    if (eDs->eArrV) eDs->eIndVert = (eDs->eArrV)->IndexOf(eVs);
    else            eDs->eIndVert = -1;
    (eDs->eCreatedTracks).Clear();
    eDs->Draw();
  }
}

SetVertex()

void EdbVertexG::SetVertex ( EdbVertex *  v )

inline

{eV=v;}

TestVertex()

void EdbVertexG::TestVertex ( )

virtual

{
  printf("\n");
  eV->Print();
  EdbTrackP *t=0;
  EdbMomentumEstimator me;
  for(int i=0; i<eV->N(); i++) {
    t = eV->GetTrack(i);
    //t->PrintNice();
    printf("nseg = %d p = %f\n",t->N(), (float)(me.PMSang(*t)) );
  }
}

Member Data Documentation

eD

EdbDisplay* EdbVertexG::eD

private

eV

EdbVertex* EdbVertexG::eV

private

---

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

• /home/antonio/fedra_doxygen/src/libEdd/EdbDisplay.h
• /home/antonio/fedra_doxygen/src/libEdd/EdbDisplay.cxx