PndTrkTrack Class Reference

#include <PndTrkTrack.h>

Inheritance diagram for PndTrkTrack:

Public Member Functions
	PndTrkTrack ()
	PndTrkTrack (PndTrkCluster *cluster)
	PndTrkTrack (PndTrkCluster *cluster, double x, double y, double radius)
	PndTrkTrack (PndTrkHit *hit, PndTrkCluster *cluster, double x, double y, double radius)
	PndTrkTrack (double x, double y, double radius)
	PndTrkTrack (PndTrack *trk)
	PndTrkTrack (const PndTrkTrack &track)
	~PndTrkTrack ()
PndTrkTrack &	operator= (const PndTrkTrack &track)
Bool_t	operator== (PndTrkTrack track)
void	SetRadius (double radius)
void	SetCenter (double x, double y)
void	SetCenter (TVector3 center)
void	SetTanL (double tanl)
void	SetZ0 (double z0)
void	SetCluster (PndTrkCluster *cluster)
void	SetRefHit (PndTrkHit *hit)
Double_t	GetRadius ()
TVector2	GetCenter ()
Double_t	GetTanL ()
Double_t	GetZ0 ()
PndTrkCluster	GetCluster ()
Double_t	ComputePhi (TVector3 hit)
Double_t	ComputePhiFrom (TVector3 hit, TVector3 from)
TVector3	ComputeMomentumAtPosition (TVector3 position, TVector3 &newposition)
void	ComputeCharge ()
PndTrackCand	ConvertToPndTrackCand ()
PndTrack	ConvertToPndTrack ()
Int_t	GetCharge ()
void	Draw (Color_t color=kBlack)
void	LightUp ()
void	Clear (Option_t *opt="")

Protected Member Functions
	ClassDef (PndTrkTrack, 1)

Protected Attributes
PndTrkCluster	fCluster
PndTrkHit *	fRefHit
double	fCenterX
double	fCenterY
double	fRadius
double	fTanL
double	fZ0
int	fCharge
double	fPhiMin
double	fPhiMax

Detailed Description

PndTrkTrack

Author

Lia Lavezzi

Definition at line 17 of file PndTrkTrack.h.

Constructor & Destructor Documentation

PndTrkTrack::PndTrkTrack

(

)

Definition at line 17 of file PndTrkTrack.cxx.

:  fCluster(PndTrkCluster()), fRefHit(NULL), fCenterX(0), fCenterY(0), fRadius(0), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360) {}

PndTrkTrack::PndTrkTrack

(

PndTrkCluster *

cluster

)

Definition at line 19 of file PndTrkTrack.cxx.

:  fCluster(*cluster), fRefHit(NULL), fCenterX(0), fCenterY(0), fRadius(0), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360){}

PndTrkTrack::PndTrkTrack	(	PndTrkCluster *	cluster,
		double	x,
		double	y,
		double	radius
	)

Definition at line 21 of file PndTrkTrack.cxx.

: fCluster(*cluster), fRefHit(NULL), fCenterX(x), fCenterY(y), fRadius(radius), fTanL(0), fZ0(0), fCharge(0) , fPhiMin(0), fPhiMax(360){}

PndTrkTrack::PndTrkTrack	(	PndTrkHit *	hit,
		PndTrkCluster *	cluster,
		double	x,
		double	y,
		double	radius
	)

Definition at line 23 of file PndTrkTrack.cxx.

: fCluster(*cluster), fRefHit(hit), fCenterX(x), fCenterY(y), fRadius(radius), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360) {}

PndTrkTrack::PndTrkTrack	(	double	x,
		double	y,
		double	radius
	)

Definition at line 25 of file PndTrkTrack.cxx.

: fCluster(PndTrkCluster()), fRefHit(NULL), fCenterX(x), fCenterY(y), fRadius(radius), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360) {}

PndTrkTrack::PndTrkTrack

(

PndTrack *

trk

)

Definition at line 28 of file PndTrkTrack.cxx.

References ComputePhi(), Double_t, fCenterX, fCenterY, fCharge, fPhiMax, fPhiMin, fRadius, fTanL, fZ0, PndTrack::GetParamFirst(), PndTrack::GetParamLast(), and pt().

                                       :  fCluster(PndTrkCluster()), fRefHit(NULL), fCenterX(0), fCenterY(0), fRadius(0), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360){

  TVector3 momentum = trk->GetParamFirst().GetMomentum();
  TVector3 position = trk->GetParamFirst().GetPosition();

  fCharge = TMath::Sign(1., trk->GetParamFirst().GetQp());

  Double_t pl = trk->GetParamFirst().GetMomentum().Z();
  Double_t pt = trk->GetParamFirst().GetMomentum().Perp();
  fRadius = pt/0.006;
  fTanL = pl/pt;

  TVector2 radius = momentum.XYvector();
  radius = radius.Unit();
  double rotx = radius.X(); 
  double roty = radius.Y(); 
  TVector2 myrad(fCharge * roty, - fCharge * rotx);
  myrad *= fRadius;

  TVector2 center = myrad + position.XYvector();
  fCenterX = center.X();
  fCenterY = center.Y();

  fZ0 = 0; // CHECK

  // ----------------------------------------------
  // cluster? // CHECK
 



  fPhiMin = ComputePhi(trk->GetParamFirst().GetPosition());
  if(fPhiMin > 180) fPhiMin -= 360;

  fPhiMax = ComputePhi(trk->GetParamLast().GetPosition());
  if(fPhiMax > 180) fPhiMax -= 360;

  if(fCharge > 0 && fPhiMin < fPhiMax) fPhiMin += 360;
  else if(fCharge < 0 && fPhiMin > fPhiMax) fPhiMax += 360;


}

PndTrkTrack::PndTrkTrack

(

const PndTrkTrack &

track

)

Definition at line 71 of file PndTrkTrack.cxx.

References track.

                                                 : TObject(track), fCluster(PndTrkCluster()), fRefHit(NULL), fCenterX(0), fCenterY(0), fRadius(0), fTanL(0), fZ0(0), fCharge(0), fPhiMin(0), fPhiMax(360) {
  *this = track;
}

PndTrkTrack::~PndTrkTrack

(

)

Definition at line 76 of file PndTrkTrack.cxx.

References fRefHit.

                          {
  delete fRefHit;

}

Member Function Documentation

PndTrkTrack::ClassDef ( PndTrkTrack  , 1    )

protected

void PndTrkTrack::Clear

(

Option_t *

opt = ""

)

Definition at line 102 of file PndTrkTrack.cxx.

References PndTrkCluster::Clear(), fCluster, and fRefHit.

                                     {
  delete fRefHit;
  fCluster.Clear(opt);
}

void PndTrkTrack::ComputeCharge

(

)

consider the centers Int_t nleft, nright = 0;

first PndTrkHit *hit1 = fCluster.GetHit(0); last PndTrkHit *hit2 = fCluster.GetHit(fCluster->GetNofHits() - 1);

vector from 1st to last hit TVector3 firsttolast = hit2->GetPosition() - hit1->GetPosition();

for(int ihit = 1; ihit < fCluster.GetNofHits() - 1; ihit++) { PndTrkHit *hit = fCluster.GetHit(ihit);

vector from 1st to this hit TVector3 position = hit->GetPosition() - hit1->GetPosition();

double crossz = position.Cross(firsttolast).Z();

if crossz > 0 hits stays 'on the right' (counterclockwise) –> negative otherwise it stays 'on the left' (clockwise) –> positive if (crossz > 0.) nright++; else nleft++; }

if(nright > nleft) fCharge = -1; else fCharge = 1;

first PndTrkHit hit1 = (PndTrkHit) fCluster.GetHit(0); TVector3 pos1 = hit1->GetPosition(); last PndTrkHit hit2 = (PndTrkHit) fCluster.GetHit(fCluster.GetNofHits() - 1); TVector3 pos2 = hit2->GetPosition();

double m = (pos2.Y() - pos1.Y())/(pos2.X() - pos1.X()); double q = pos1.Y() - m + pos1.X();

TVector3 posbar(0.5 * (pos1.X() + pos2.X()), 0.5 * (pos1.Y() + pos2.Y()));

double mortho = -1./m; double qortho = posbar.Y() - mortho * posbar.X();

solve the system (x - xc)**2 + (y - yc)**2 = R**2 y = mortho * x + qortho and find xI xII double alpha = 1 + mortho * mortho; double beta = mortho * (qortho - fCenterY) - fCenterX; double gamma = fCenterX * fCenterX + (qortho - fCenterY) * (qortho - fCenterY) - fRadius * fRadius;

double xI = (-beta + TMath::Sqrt(beta * beta - alpha * gamma))/alpha; double yI = xI * mortho + qortho; double xII = (-beta - TMath::Sqrt(beta * beta - alpha * gamma))/alpha; double yII = xII * mortho + qortho;

mean hit int mean = (int) (fCluster.GetNofHits() * 0.5); PndTrkHit hit_mean = (PndTrkHit) fCluster.GetHit(mean); double distI = hit_mean->GetXYDistance(TVector3(xI, yI, 0.)); double distII = hit_mean->GetXYDistance(TVector3(xII, yII, 0.));

double x, y; if(distI < distII) { x = xI; y = yI; } else { x = xII; y = yII; }

translation x -= pos1.X(); y -= pos1.Y();

double delta = (pos2 - pos1).Phi(); double xrot = x * TMath::Cos(delta) + y * TMath::Sin(delta); double yrot = -x * TMath::Sin(delta) + y * TMath::Cos(delta);

cout << "x, y " << xrot << " " << yrot << endl;

Definition at line 171 of file PndTrkTrack.cxx.

References PndTrkTools::ComputePocaToPointOnCircle3(), fCenterX, fCenterY, fCharge, fCluster, fRadius, PndTrkCluster::GetHit(), PndTrkCluster::GetNofHits(), PndTrkHit::GetPosition(), hit, phi, and Pi.

Referenced by ConvertToPndTrack(), PndTrkLegendreSecTask::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

                                { // CHECK!!!

   // consider the poca on track
   Int_t nmore = 0, nless = 0;
   
   // first
   PndTrkHit *hit1 = (PndTrkHit*) fCluster.GetHit(0);
   PndTrkTools tools;
   TVector3 position1 = tools.ComputePocaToPointOnCircle3(hit1->GetPosition().X(), hit1->GetPosition().Y(), fCenterX, fCenterY, fRadius);
   
   // phi of first hit with respect to center of curvature
   TVector3 direction1 = position1 - TVector3(fCenterX, fCenterY, 0);
   double tmpphi = direction1.Phi();
   //  
   if(direction1.Y()) tmpphi += 2 * TMath::Pi(); // CHECK
   
   for(int ihit = 1; ihit < fCluster.GetNofHits(); ihit++) {
   PndTrkHit *hit = fCluster.GetHit(ihit);
   TVector3 position = tools.ComputePocaToPointOnCircle3(hit->GetPosition().X(), hit->GetPosition().Y(), fCenterX, fCenterY, fRadius);
   
   // vector from 1st to this hit
   TVector3 direction = position - TVector3(fCenterX, fCenterY, 0);
   double phi = direction.Phi();
   if(ihit > 1) phi >= tmpphi ? nmore++ : nless++;
   tmpphi = phi;
   //    cout << "phi " << phi * TMath::RadToDeg() << " " << hit->GetHitID() << " " << hit->GetDetectorID() << endl;
   }
   
   if(nmore > nless) fCharge = -1;
   else fCharge = 1;
   
   //   cout << "fCharge " << fCharge << " " << nmore << " " << nless << endl;

}

TVector3 PndTrkTrack::ComputeMomentumAtPosition	(	TVector3	position,
		TVector3 &	newposition
	)

Definition at line 293 of file PndTrkTrack.cxx.

References CAMath::Abs(), PndTrkTools::ComputePocaToPointOnCircle3(), Double_t, fCenterX, fCenterY, fCharge, fRadius, fTanL, pos, and pt().

Referenced by ConvertToPndTrack().

                                                                                        {
  TVector3 momentum(-999, -999, -999);
  newposition = position;
  if(fCharge == 0) {
    // cout << "ComputeMomentumAtPos: CHARGE = 0!" << endl;
    return momentum;
  }

  TVector2 center(fCenterX, fCenterY);
  TVector2 pos = TVector2(position.X(), position.Y());
  TVector2 myrad = center - pos;
  double distance = TMath::Abs(myrad.Mod() - fRadius);
  if(distance > 0.5) {
    // cout << "ComputeMomentumAtPosition: POINT NOT ON THE TRACK " << distance << endl;
    // return momentum;  // CHECK
    // alternative
    // CHECK
    PndTrkTools tools;
    newposition = tools.ComputePocaToPointOnCircle3(pos.X(), pos.Y(), fCenterX, fCenterY, fRadius);
    myrad = center - pos;
    distance = TMath::Abs(myrad.Mod() - fRadius);
    //   if(distance > 0.5) {
    //       cout << "ComputeMomentumAtPosition: AGAIN POINT NOT ON THE TRACK " << distance << endl;
    //       return momentum;  
    //     }
  }


  Double_t rotx, roty;
  rotx = - fCharge * myrad.Y();
  roty = fCharge * myrad.X();


//   cout << "COMPUTE MOMENTUM " << fRadius << endl;
  Double_t pt = 0.006 * fRadius;

//   cout << "pt " << pt << endl;
  Double_t pl = -999;
  if(fTanL != -999) pl = pt * fTanL;
//   cout << "pl " << pl << " tanl " << fTanL << endl;

  //Double_t ptot = TMath::Sqrt(pt * pt + pl * pl); //[R.K. 01/2017] unused variable


//   cout << rotx << " " << roty << endl;
  momentum.SetX(rotx); // CHECK magnitude?
  momentum.SetY(roty); // CHECK magnitude?
  momentum.SetZ(0.); // CHECK magnitude?
  //   momentum.Print();
  momentum.SetMag(pt);
  //   momentum.Print();
  momentum.SetZ(pl);
  //  momentum.Print();
  return momentum;
}

Double_t PndTrkTrack::ComputePhi

(

TVector3

hit

)

 TVector3 center(fCenterX, fCenterY, 0.);
 TVector3 fromcentertohit = hit - center;

I want the positive phi angle from x axis in range [0, 360[. I use TVector3::Phi() [fromcentertohit.Phi()]:

x y Phi use!

• + 0/90 phi
• + 90/180 phi
• - -180/-90 phi + 360
• - -90/0 phi + 360 cout << "phi " << hit.X() << " " << hit.Y() << endl;

   double phi = fromcentertohit.Phi();
   if(fromcentertohit.Y() < 0) phi += (2 * TMath::Pi());
  

  cout << "final phi in rad " << phi << endl; return phi * TMath::RadToDeg();

Definition at line 349 of file PndTrkTrack.cxx.

References alpha, CAMath::ATan2(), CAMath::Cos(), d, Double_t, fCenterX, fCenterY, fCharge, fRadius, p, phi, Pi, pi, CAMath::Sin(), CAMath::Sqrt(), v, x0, and y0.

Referenced by PndTrkLegendreSecTask::CreateSkewHitList(), PndTrkLegendreTask::CreateSkewHitList(), PndTrkLegendreSecTask2::CreateSkewHitList(), PndTrkTrackFinder::CreateSkewHitList(), Draw(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreSecTask::Exec(), PndTrkLegendreTask::Exec(), PndTrkTrackFinder::Exec(), and PndTrkTrack().

{
  // x0 y0
  Double_t d = TMath::Sqrt(fCenterX * fCenterX + fCenterY * fCenterY) - fRadius;
  Double_t phi =  TMath::ATan2(fCenterY, fCenterX);
  
  Double_t x0 = d * TMath::Cos(phi);
  Double_t y0 = d * TMath::Sin(phi);

  Double_t Phi0 = TMath::ATan2((y0 - fCenterY),(x0 - fCenterX));

  // CHECK :-)GOOD! ...
  TVector2 v(x0 - fCenterX, y0 - fCenterY); 
  double alpha = TMath::ATan2(hit.Y() - y0 + fRadius * TMath::Sin(Phi0), hit.X() - x0 + fRadius * TMath::Cos(Phi0));
  TVector2 p(hit.X() - fCenterX, hit.Y() - fCenterY);
 
  Double_t Fi = - fCharge *  TMath::ACos(v * p / (v.Mod() * p.Mod()));
  double pi = TMath::Pi();
  double pi2 = 2 * pi;
     
  // Fi = h * (pi2 - h * Fi) // should be correct
  if((fCharge > 0 && ((Phi0 > 0 && ((alpha > 0 && alpha > Phi0) ||
                                   (alpha < 0 && alpha < Phi0 - pi)))
                     ||
                     ((Phi0 < 0 && ((alpha > 0 && alpha < pi + Phi0) ||
                                    (alpha < 0 && alpha > Phi0)))) ))) Fi = - (pi2 + Fi)  ;
  else if((fCharge < 0 && ((Phi0 > 0 && ((alpha > 0 && alpha < Phi0) ||
                                        (alpha < 0 && alpha > Phi0 - pi)))
                          ||
                          ((Phi0 < 0 && ((alpha > 0 && alpha > pi + Phi0) ||
                                         (alpha < 0 && alpha < Phi0)))) ))) Fi = pi2 - Fi  ;

 //  cout << "PHI ----------------- " <<  Fi * TMath::RadToDeg() << endl; 
  
  return (Phi0 + Fi) * TMath::RadToDeg();

}

Double_t PndTrkTrack::ComputePhiFrom	(	TVector3	hit,
		TVector3	from
	)

Definition at line 408 of file PndTrkTrack.cxx.

References alpha, CAMath::ATan2(), CAMath::Cos(), Double_t, fCenterX, fCenterY, fCharge, fRadius, phi, phi0, Pi, pi, CAMath::Sin(), x0, and y0.

Referenced by PndTrkTrackFinder::Exec().

{
  //  cout << "phi " << hit.X() << " " << hit.Y() << endl;

  TVector3 center(fCenterX, fCenterY, 0.);
  TVector3 fromcentertofrom = from - center;
  double phi0 = fromcentertofrom.Phi();
  
  double xtr = hit.X() - fCenterX;
  double ytr = hit.Y() - fCenterY;
  double xp = TMath::Cos(phi0) * xtr + TMath::Sin(phi0) * ytr;
  double yp = -TMath::Sin(phi0) * xtr + TMath::Cos(phi0) * ytr;

  TVector3 trarothit = TVector3(xp, yp, 0.);

  // I want the positive phi angle from x axis
  // in range [0, 360[.
  // I use TVector3::Phi() [fromcentertohit.Phi()]:
  // x   y      Phi        use!
  // -------------------------
  // +   +      0/90       phi
  // -   +     90/180      phi
  // -   -   -180/-90    phi + 360
  // +   -    -90/0      phi + 360
  

  double phi = trarothit.Phi();
  //  if(trarothit.Y() < 0) phi += (2 * TMath::Pi());

  //  cout << "final phi in rad " << phi << endl;
  //  return phi * TMath::RadToDeg();

   Double_t Phi0 = phi0;

   Double_t x0 = from.X();
   Double_t y0 = from.Y();

//   Double_t Phi0 = TMath::ATan2((y0 - fCenterY),(x0 - fCenterX));

//   // CHECK :-)GOOD! ...
//   TVector2 v(x0 - fCenterX, y0 - fCenterY); 
   double alpha = TMath::ATan2(hit.Y() - y0 + fRadius * TMath::Sin(Phi0), hit.X() - x0 + fRadius * TMath::Cos(Phi0));
//   TVector2 p(hit.X() - fCenterX, hit.Y() - fCenterY);
 
//   Double_t Fi = - fCharge *  TMath::ACos(v * p / (v.Mod() * p.Mod()));
   double pi = TMath::Pi();
   double pi2 = 2 * pi;
     
//   // Fi = h * (pi2 - h * Fi) // should be correct

   Double_t Fi = phi;
   if((fCharge > 0 && ((Phi0 > 0 && ((alpha > 0 && alpha > Phi0) ||
                                    (alpha < 0 && alpha < Phi0 - pi)))
                      ||
                      ((Phi0 < 0 && ((alpha > 0 && alpha < pi + Phi0) ||
                                     (alpha < 0 && alpha > Phi0)))) ))) Fi = - (pi2 + Fi)  ;
  else if((fCharge < 0 && ((Phi0 > 0 && ((alpha > 0 && alpha < Phi0) ||
                                         (alpha < 0 && alpha > Phi0 - pi)))
                           ||
                           ((Phi0 < 0 && ((alpha > 0 && alpha > pi + Phi0) ||
                                          (alpha < 0 && alpha < Phi0)))) ))) Fi = pi2 - Fi  ;

  //  cout << "PHI ----------------- " <<  Fi * TMath::RadToDeg() << endl; 
 
  return Fi * TMath::RadToDeg();


}

PndTrack PndTrkTrack::ConvertToPndTrack

(

)

Definition at line 119 of file PndTrkTrack.cxx.

References ComputeCharge(), ComputeMomentumAtPosition(), ConvertToPndTrackCand(), fCharge, fCluster, PndTrkCluster::GetHit(), PndTrkCluster::GetNofHits(), PndTrkHit::GetPosition(), PndTrack::SetFlag(), and track.

Referenced by PndTrkTrackFinder::Exec(), PndTrkLegendreSecTask::RegisterTrack(), PndTrkLegendreTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

                                        {

  ComputeCharge();
  
  // first
  TVector3 pos1, mom1;
  PndTrkHit *hit1 = fCluster.GetHit(0);
  mom1 = ComputeMomentumAtPosition(hit1->GetPosition(), pos1);

  // last
  TVector3 pos2, mom2;
  PndTrkHit *hit2 = fCluster.GetHit(fCluster.GetNofHits() - 1);
  mom2 = ComputeMomentumAtPosition(hit2->GetPosition(), pos2);
    
  TVector3 dj(1, 0, 0), dk(0, 0, 1);   // CHECK
  TVector3 dpos1(2., 2., 10.);         // CHECK
  TVector3 dpos2(2., 2., 10.);         // CHECK
  TVector3 dmom1(0.2, 0.2, 0.5);       // CHECK
  TVector3 dmom2(0.2, 0.2, 0.5);       // CHECK

  FairTrackParP firstpar(pos1, mom1,
                         dpos1, dmom1, fCharge,
                         pos1, dj, dk);
    
  FairTrackParP lastpar(pos2, mom2,
                        dpos2, dmom2, fCharge,
                        pos2, dj, dk);

  PndTrackCand trkcand = ConvertToPndTrackCand();
  PndTrack track(firstpar, lastpar, trkcand);
  if(mom1.Z() == -999) track.SetFlag(-1);
  return track;
}

PndTrackCand PndTrkTrack::ConvertToPndTrackCand

(

)

Definition at line 108 of file PndTrkTrack.cxx.

References PndTrackCand::AddHit(), fCluster, PndTrkHit::GetDetectorID(), PndTrkCluster::GetHit(), PndTrkHit::GetHitID(), PndTrkCluster::GetNofHits(), and hit.

Referenced by ConvertToPndTrack().

                                                {
  PndTrackCand trkcand;
  for(int ihit = 0; ihit < fCluster.GetNofHits(); ihit++) {
    PndTrkHit *hit = fCluster.GetHit(ihit);
    int detid = hit->GetDetectorID();
    int hitid = hit->GetHitID();
    trkcand.AddHit(detid, hitid, ihit);
  }
  return trkcand;
}

void PndTrkTrack::Draw

(

Color_t

color = kBlack

)

Definition at line 477 of file PndTrkTrack.cxx.

References ComputePhi(), fCenterX, fCenterY, fCharge, fCluster, fPhiMax, fPhiMin, fRadius, PndTrkCluster::GetHit(), PndTrkCluster::GetNofHits(), PndTrkHit::GetPosition(), and track.

Referenced by PndTrkCombiLegendreTask::CleanTrack(), PndTrkLegendreNew::CreateClusterAroundTrack(), PndTrkTrackFinder::CreateClusterAroundTrack(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreNew::Exec(), PndTrkLegendreTask::Exec(), PndTrkLegendreSecTask::Exec(), PndTrkCombiLegendreTask::Exec(), PndTrkTrackFinder::Exec(), PndTrkCombiLegendreTask::LegendreFit(), LightUp(), PndTrkLegendreTask::RegisterTrack(), PndTrkLegendreSecTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

                                    {
  cout << "draw" << endl;

  if(fCluster.GetNofHits() > 0)    {
    PndTrkHit *hit0 = fCluster.GetHit(0);
    fPhiMin = ComputePhi(hit0->GetPosition());
  //   if(fPhiMin > 180) fPhiMin -= 360;
    hit0->GetPosition().Print();
    
    PndTrkHit *hitN = fCluster.GetHit(fCluster.GetNofHits() - 1);
    fPhiMax = ComputePhi(hitN->GetPosition());
 //    if(fPhiMax > 180) fPhiMax -= 360;
    hitN->GetPosition().Print();


 //    if(fCharge > 0 && fPhiMin < fPhiMax) fPhiMin += 360;
//     else if(fCharge < 0 && fPhiMin > fPhiMax) fPhiMax += 360;

  }
  TArc *track = new TArc(fCenterX, fCenterY, fRadius, fPhiMin, fPhiMax);




  cout << fCharge << " " << fCenterX << " " << fCenterY << " " << fRadius << " " << fPhiMin << " " << fPhiMax << endl;
  track->SetFillStyle(0);
  track->SetLineColor(color);
  track->Draw("only SAME");
}

TVector2 PndTrkTrack::GetCenter ( )

inline

Definition at line 44 of file PndTrkTrack.h.

References fCenterX, and fCenterY.

Referenced by PndTrkLegendreNew::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack(), PndTrkTrackFinder::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack2(), PndTrkCombiLegendreTask::CreateClusterAroundTrack3(), PndTrkLegendreTask::CreateSkewHitList(), PndTrkLegendreSecTask::CreateSkewHitList(), PndTrkLegendreSecTask2::CreateSkewHitList(), PndTrkTrackFinder::CreateSkewHitList(), PndTrkLegendreNew::Exec(), PndTrkCombiLegendreTask::Exec(), PndTrkTrackFinder::Exec(), operator==(), PndTrkLegendreTask::RegisterTrack(), PndTrkLegendreSecTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

{ return TVector2(fCenterX, fCenterY); }

Int_t PndTrkTrack::GetCharge ( )

inline

Definition at line 57 of file PndTrkTrack.h.

References fCharge.

Referenced by PndTrkLegendreSecTask::Exec(), PndTrkLegendreTask::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkTrackFinder::Exec(), PndTrkLegendreSecTask::RegisterTrack(), PndTrkLegendreTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

{ return fCharge; }

PndTrkCluster PndTrkTrack::GetCluster ( )

inline

Definition at line 48 of file PndTrkTrack.h.

References fCluster.

Referenced by PndTrkCombiLegendreTask::CleanTrack(), PndTrkLegendreNew::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack(), PndTrkTrackFinder::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack3(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreTask::Exec(), PndTrkLegendreNew::Exec(), PndTrkLegendreSecTask::Exec(), PndTrkCombiLegendreTask::Exec(), PndTrkTrackFinder::Exec(), and operator==().

{ return fCluster; }

Double_t PndTrkTrack::GetRadius ( )

inline

Definition at line 43 of file PndTrkTrack.h.

References fRadius.

Referenced by PndTrkLegendreNew::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack(), PndTrkTrackFinder::CreateClusterAroundTrack(), PndTrkCombiLegendreTask::CreateClusterAroundTrack2(), PndTrkCombiLegendreTask::CreateClusterAroundTrack3(), PndTrkLegendreTask::CreateSkewHitList(), PndTrkLegendreSecTask::CreateSkewHitList(), PndTrkLegendreSecTask2::CreateSkewHitList(), PndTrkTrackFinder::CreateSkewHitList(), PndTrkLegendreNew::Exec(), PndTrkCombiLegendreTask::Exec(), PndTrkTrackFinder::Exec(), operator==(), PndTrkLegendreTask::RegisterTrack(), PndTrkLegendreSecTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

{ return fRadius; }

Double_t PndTrkTrack::GetTanL ( )

inline

Definition at line 45 of file PndTrkTrack.h.

References fTanL.

Referenced by PndTrkTrackFinder::Exec(), PndTrkLegendreSecTask::RegisterTrack(), PndTrkLegendreTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

{ return fTanL; } 

Double_t PndTrkTrack::GetZ0 ( )

inline

Definition at line 46 of file PndTrkTrack.h.

References fZ0.

Referenced by PndTrkTrackFinder::Exec(), PndTrkLegendreSecTask::RegisterTrack(), PndTrkLegendreTask::RegisterTrack(), and PndTrkLegendreSecTask2::RegisterTrack().

{ return fZ0; }

void PndTrkTrack::LightUp

(

)

Definition at line 507 of file PndTrkTrack.cxx.

References Draw(), fCluster, and PndTrkCluster::LightUp().

                          {
  Draw(kYellow);
  fCluster.LightUp();
}

PndTrkTrack & PndTrkTrack::operator=

(

const PndTrkTrack &

track

)

Definition at line 83 of file PndTrkTrack.cxx.

References fCenterX, fCenterY, fCharge, fCluster, fPhiMax, fPhiMin, fRadius, fRefHit, fTanL, and fZ0.

                                                            {
  fRefHit = track.fRefHit;
  fCluster = track.fCluster;
  fCenterX = track.fCenterX;
  fCenterY = track.fCenterY;
  fRadius = track.fRadius;
  fTanL = track.fTanL;
  fZ0 = track.fZ0;
  fCharge = track.fCharge;
  fPhiMin = track.fPhiMin;
  fPhiMax = track.fPhiMax;
  return *this;
}

Bool_t PndTrkTrack::operator==

(

PndTrkTrack

track

)

Definition at line 97 of file PndTrkTrack.cxx.

References fCenterX, fCenterY, fRadius, GetCenter(), GetCluster(), PndTrkCluster::GetNofHits(), and GetRadius().

                                                 {
 int nofhits = GetCluster().GetNofHits();
  return nofhits == track.GetCluster().GetNofHits() && fRadius == track.GetRadius() && fCenterX == track.GetCenter().X() && fCenterY == track.GetCenter().Y(); // CHECK
}

void PndTrkTrack::SetCenter ( double  x, double  y  )

inline

Definition at line 37 of file PndTrkTrack.h.

References fCenterX, fCenterY, x, and y.

Referenced by PndTrkLegendreNew::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkCombiLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

{ fCenterX = x; fCenterY = y; }

void PndTrkTrack::SetCenter ( TVector3  center )

inline

Definition at line 38 of file PndTrkTrack.h.

References fCenterX, and fCenterY.

{ fCenterX = center.X(); fCenterY = center.Y(); }

void PndTrkTrack::SetCluster ( PndTrkCluster *  cluster )

inline

Definition at line 41 of file PndTrkTrack.h.

References fCluster.

Referenced by PndTrkLegendreNew::Exec(), PndTrkCombiLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

{ fCluster = *cluster; }

void PndTrkTrack::SetRadius ( double  radius )

inline

Definition at line 36 of file PndTrkTrack.h.

References fRadius.

Referenced by PndTrkLegendreNew::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkCombiLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

{ fRadius = radius; }

void PndTrkTrack::SetRefHit ( PndTrkHit *  hit )

inline

Definition at line 42 of file PndTrkTrack.h.

References fRefHit, and hit.

{ fRefHit = hit; }

void PndTrkTrack::SetTanL ( double  tanl )

inline

Definition at line 39 of file PndTrkTrack.h.

References fTanL, and tanl.

Referenced by PndTrkLegendreSecTask::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

{ fTanL = tanl; }

void PndTrkTrack::SetZ0 ( double  z0 )

inline

Definition at line 40 of file PndTrkTrack.h.

References fZ0, and z0.

Referenced by PndTrkLegendreSecTask::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreTask::Exec(), and PndTrkTrackFinder::Exec().

{ fZ0 = z0;}

Member Data Documentation

double PndTrkTrack::fCenterX

protected

Definition at line 68 of file PndTrkTrack.h.

Referenced by ComputeCharge(), ComputeMomentumAtPosition(), ComputePhi(), ComputePhiFrom(), Draw(), GetCenter(), operator=(), operator==(), PndTrkTrack(), and SetCenter().

double PndTrkTrack::fCenterY

protected

Definition at line 68 of file PndTrkTrack.h.

Referenced by ComputeCharge(), ComputeMomentumAtPosition(), ComputePhi(), ComputePhiFrom(), Draw(), GetCenter(), operator=(), operator==(), PndTrkTrack(), and SetCenter().

int PndTrkTrack::fCharge

protected

Definition at line 69 of file PndTrkTrack.h.

Referenced by ComputeCharge(), ComputeMomentumAtPosition(), ComputePhi(), ComputePhiFrom(), ConvertToPndTrack(), Draw(), GetCharge(), operator=(), and PndTrkTrack().

PndTrkCluster PndTrkTrack::fCluster

protected

Definition at line 66 of file PndTrkTrack.h.

Referenced by Clear(), ComputeCharge(), ConvertToPndTrack(), ConvertToPndTrackCand(), Draw(), GetCluster(), LightUp(), operator=(), and SetCluster().

double PndTrkTrack::fPhiMax

protected

Definition at line 71 of file PndTrkTrack.h.

Referenced by Draw(), operator=(), and PndTrkTrack().

double PndTrkTrack::fPhiMin

protected

Definition at line 70 of file PndTrkTrack.h.

Referenced by Draw(), operator=(), and PndTrkTrack().

double PndTrkTrack::fRadius

protected

Definition at line 68 of file PndTrkTrack.h.

Referenced by ComputeCharge(), ComputeMomentumAtPosition(), ComputePhi(), ComputePhiFrom(), Draw(), GetRadius(), operator=(), operator==(), PndTrkTrack(), and SetRadius().

PndTrkHit* PndTrkTrack::fRefHit

protected

Definition at line 67 of file PndTrkTrack.h.

Referenced by Clear(), operator=(), SetRefHit(), and ~PndTrkTrack().

double PndTrkTrack::fTanL

protected

Definition at line 68 of file PndTrkTrack.h.

Referenced by ComputeMomentumAtPosition(), GetTanL(), operator=(), PndTrkTrack(), and SetTanL().

double PndTrkTrack::fZ0

protected

Definition at line 68 of file PndTrkTrack.h.

Referenced by GetZ0(), operator=(), PndTrkTrack(), and SetZ0().

---

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

• PndTrkTrack.h
• PndTrkTrack.cxx