PndTrkLegendreNew Class Reference

#include <PndTrkLegendreNew.h>

Inheritance diagram for PndTrkLegendreNew:

Public Member Functions
	PndTrkLegendreNew ()
	PndTrkLegendreNew (int verbose)
	~PndTrkLegendreNew ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
void	SetParContainers ()
void	Initialize ()
void	Reset ()
void	SwitchOnDisplay ()
void	DrawGeometry ()
void	DrawHits (PndTrkHitList *hitlist)
void	DrawLists ()
void	DrawNeighborings ()
void	DrawNeighboringsToHit (PndTrkHit *hit)
void	Refresh ()
void	RefreshConf ()
void	DrawGeometryConf (double x1, double x2, double y1, double y2)
void	DrawConfHit (double x, double y, double r, int marker=2)
void	LightCluster (PndTrkCluster *cluster)
void	DrawLegendreHisto ()
void	ComputePlaneExtremities (PndTrkCluster *cluster)
Int_t	FillConformalHitList (PndTrkCluster *cluster)
void	FillLegendreHisto (PndTrkCluster *cluster)
void	ComputeTraAndRot (PndTrkHit *hit, Double_t &delta, Double_t trasl[2])
PndTrkHit *	FindSttReferenceHit (int isec=-1)
PndTrkHit *	FindMvdPixelReferenceHit ()
PndTrkHit *	FindMvdStripReferenceHit ()
PndTrkHit *	FindMvdReferenceHit ()
PndTrkHit *	FindReferenceHit ()
PndTrkHit *	FindReferenceHit (PndTrkCluster *cluster)
void	RePrepareLegendre (PndTrkCluster *cluster)
Int_t	ApplyLegendre (PndTrkCluster *cluster, double &theta_max, double &r_max)
Int_t	ExtractLegendre (Int_t mode, double &theta_max, double &r_max)
void	FromConformalToRealTrack (double fitm, double fitp, double &x0, double &y0, double &R)
void	FromRealToConformalTrack (double x0, double y0, double R, double &fitm, double &fitp)
void	SearchSecondaryTracks ()
PndTrkClusterList	CreateFullClusterization ()
PndTrkClusterList	CreateFullClusterization2 ()
Int_t	CountTracksInCluster (PndTrkCluster *cluster)
Int_t	CountTracksInSkewSector (PndTrkCluster *cluster)
Int_t	CountTracksInCluster (PndTrkCluster *cluster, Int_t where)
Int_t	ClusterToConformal (PndTrkCluster *cluster)
PndTrkTrack *	LegendreFit (PndTrkCluster *cluster)
PndTrkCluster *	CreateClusterAroundTrack (PndTrkTrack *track)
void	AnalyticalFit (PndTrkCluster *cluster, double xc, double yc, double R, double &fitm, double &fitq)
void	AnalyticalFit2 (PndTrkCluster *cluster, double fitm, double fitp, double &fitm2, double &fip2)
void	IntersectionFinder (PndTrkHit *hit, double xc, double yc, double R)
void	IntersectionFinder (PndTrkConformalHit *chit, double fitm, double fitp)
void	FillHitMap ()

Private Member Functions
	ClassDef (PndTrkLegendreNew, 1)

Private Attributes
Int_t	fNofMvdPixHits
Int_t	fNofMvdStrHits
Int_t	fNofSttHits
Int_t	fNofTriplets
Int_t	fNofHits
TClonesArray *	fSttPointArray
TClonesArray *	fSttHitArray
TClonesArray *	fMvdPixelHitArray
TClonesArray *	fMvdStripHitArray
TClonesArray *	fTrackArray
TClonesArray *	fTrackCandArray
TClonesArray *	fTubeArray
PndGeoSttPar *	fSttParameters
char	fSttBranch [200]
char	fMvdPixelBranch [200]
char	fMvdStripBranch [200]
PndSttMapCreator *	fMapper
Int_t	fEventCounter
PndTrkSttHitList *	stthitlist
PndTrkSdsHitList *	mvdpixhitlist
PndTrkSdsHitList *	mvdstrhitlist
Double_t	fSttParalDistance
Double_t	fSttToMvdStripDistance
double	fDeltaThetaRad
PndTrkLegendreTransform *	legendre
Bool_t	fPersistence
Bool_t	fUseMVDPix
Bool_t	fUseMVDStr
Bool_t	fUseSTT
Bool_t	fSecondary
Bool_t	fInitDone
PndTrkConformalTransform *	conform
PndTrkConformalHitList *	conformalhitlist
PndTrkTools *	tools
PndTrkHit *	fRefHit
Double_t	fMvdPix_RealDistLimit
Double_t	fMvdStr_RealDistLimit
Double_t	fStt_RealDistLimit
Double_t	fMvdPix_ConfDistLimit
Double_t	fMvdStr_ConfDistLimit
Double_t	fStt_ConfDistLimit
std::vector< std::pair< double, double > >	fFoundPeaks
double	fTime
TStopwatch *	fTimer
PndTrkFitter *	fFitter
PndTrkNeighboringMap *	fHitMap
double	fUmin
double	fUmax
double	fVmin
double	fVmax
double	fRmin
double	fRmax
double	fThetamin
double	fThetamax
Bool_t	fDisplayOn
TH2F *	hxy
TH2F *	hxz
TH2F *	hzphi
TCanvas *	display
TH2F *	huv

Detailed Description

Definition at line 35 of file PndTrkLegendreNew.h.

Constructor & Destructor Documentation

PndTrkLegendreNew::PndTrkLegendreNew

(

)

Default constructor

Definition at line 56 of file PndTrkLegendreNew.cxx.

References fMvdPixelBranch, fMvdStripBranch, fSttBranch, and PndGeoHandling::Instance().

                                     : FairTask("secondary track finder", 0), fDisplayOn(kFALSE), fPersistence(kTRUE), fUseMVDPix(kTRUE), fUseMVDStr(kTRUE), fUseSTT(kTRUE), fSecondary(kFALSE), fMvdPix_RealDistLimit(1000), fMvdStr_RealDistLimit(1000), fStt_RealDistLimit(1000), fMvdPix_ConfDistLimit(1000), fMvdStr_ConfDistLimit(1000), fStt_ConfDistLimit(1000), fInitDone(kFALSE), fUmin(-0.07), fUmax(0.07), fVmin(-0.07), fVmax(0.07), fRmin(-1.5), fRmax(1.5), fThetamin(0), fThetamax(180) {
  sprintf(fSttBranch,"STTHit");
  sprintf(fMvdPixelBranch,"MVDHitsPixel");
  sprintf(fMvdStripBranch,"MVDHitsStrip");
  PndGeoHandling::Instance();
}

PndTrkLegendreNew::PndTrkLegendreNew

(

int

verbose

)

Definition at line 63 of file PndTrkLegendreNew.cxx.

References fMvdPixelBranch, fMvdStripBranch, fSttBranch, and PndGeoHandling::Instance().

                                                : FairTask("secondary track finder", verbose), fDisplayOn(kFALSE), fPersistence(kTRUE), fUseMVDPix(kTRUE), fUseMVDStr(kTRUE), fUseSTT(kTRUE), fSecondary(kFALSE), fMvdPix_RealDistLimit(1000), fMvdStr_RealDistLimit(1000), fStt_RealDistLimit(1000), fMvdPix_ConfDistLimit(1000), fMvdStr_ConfDistLimit(1000), fStt_ConfDistLimit(1000), fInitDone(kFALSE), fUmin(-0.07), fUmax(0.07), fVmin(-0.07), fVmax(0.07), fRmin(-1.5), fRmax(1.5), fThetamin(0), fThetamax(180) {
  sprintf(fSttBranch,"STTHit");
  sprintf(fMvdPixelBranch,"MVDHitsPixel");
  sprintf(fMvdStripBranch,"MVDHitsStrip");
  PndGeoHandling::Instance();
}

PndTrkLegendreNew::~PndTrkLegendreNew

(

)

Destructor

Definition at line 73 of file PndTrkLegendreNew.cxx.

                                      { 

}

Member Function Documentation

void PndTrkLegendreNew::AnalyticalFit	(	PndTrkCluster *	cluster,
		double	xc,
		double	yc,
		double	R,
		double &	fitm,
		double &	fitq
	)

Definition at line 1687 of file PndTrkLegendreNew.cxx.

References conform, display, fDisplayOn, fFitter, FromConformalToRealTrack(), PndTrkConformalTransform::GetConformalHit(), PndTrkCluster::GetHit(), PndTrkCluster::GetNofHits(), PndTrkConformalHit::GetPosition(), hit, IntersectionFinder(), PndTrkHit::IsSttParallel(), PndTrkHit::IsSttSkew(), Refresh(), PndTrkFitter::Reset(), PndTrkFitter::SetPointToFit(), sigma, and PndTrkFitter::StraightLineFit().

                                                                                                                       {
  
  
  // fit with analytical chi2 -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
  fFitter->Reset();
  if(fDisplayOn) {
    display->cd(1);
    Refresh();
  }
  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) 
    {
      PndTrkHit *hit = cluster->GetHit(ihit);
      if(hit->IsSttSkew()) continue;
      if(hit->IsSttParallel()) IntersectionFinder(hit, xc, yc, R);
      PndTrkConformalHit *chit = conform->GetConformalHit(hit);
      double sigma = 0.1; // CHECK
      fFitter->SetPointToFit(chit->GetPosition().X(), chit->GetPosition().Y(), sigma);
      //      cout << "set point to fit " << chit->GetPosition().X() << " " <<  chit->GetPosition().Y() << endl;
      //     if(fDisplayOn) {
      //        display->cd(1);
      //        TMarker *mrk = new TMarker(hit->GetPosition().X(), hit->GetPosition().Y(), 6);
      //        mrk->SetMarkerColor(kRed);
      //        mrk->Draw("SAME");

      //        display->cd(2);
      //        TMarker *mrk2 = new TMarker(chit->GetPosition().X(), chit->GetPosition().Y(), 6);
      //        mrk2->SetMarkerColor(kRed);
      //        mrk2->Draw("SAME");


      //        display->Update();
      //        display->Modified();
      //       } 
    }

  
  fFitter->StraightLineFit(fitm, fitq);
  cout << "previous " << xc << " " << yc << " " << R << endl;
  FromConformalToRealTrack(fitm, fitq, xc, yc, R);
  cout << "now " << xc << " " << yc << " " << R << endl;

  if(fDisplayOn) {
    display->cd(2);
    cout << "wanna see the line?" << endl;
    TLine *line = new TLine(-10.07, fitq + fitm * (-10.07), 10.07, fitq + fitm * (10.07));
    line->SetLineColor(2);
    line->Draw("SAME");
    char goOnChar;
    display->Update();
    display->Modified();
    cin >> goOnChar;
  }




}

void PndTrkLegendreNew::AnalyticalFit2	(	PndTrkCluster *	cluster,
		double	fitm,
		double	fitp,
		double &	fitm2,
		double &	fip2
	)

Definition at line 1746 of file PndTrkLegendreNew.cxx.

References conform, display, PndTrkConformalHit::Draw(), fDisplayOn, fFitter, FromConformalToRealTrack(), PndTrkConformalTransform::GetConformalSttHit(), PndTrkCluster::GetHit(), PndTrkConformalHit::GetIsochrone(), PndTrkCluster::GetNofHits(), PndTrkConformalHit::GetPosition(), PndTrkHit::GetPosition(), hit, IntersectionFinder(), PndTrkHit::IsSttSkew(), R, Refresh(), PndTrkFitter::Reset(), PndTrkFitter::SetPointToFit(), sigma, and PndTrkFitter::StraightLineFit().

Referenced by Exec().

                                                                                                                    {
  
  
  // fit with analytical chi2 -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
  fFitter->Reset();
  if(fDisplayOn) {
    display->cd(1);
    Refresh();
  }
  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) 
    {
      PndTrkHit *hit = cluster->GetHit(ihit);
      if(hit->IsSttSkew()) continue;
  
      PndTrkConformalHit *chit = conform->GetConformalSttHit(hit);
      IntersectionFinder(chit, fitm, fitp);
      
      double sigma = chit->GetIsochrone();
      fFitter->SetPointToFit(chit->GetPosition().X(), chit->GetPosition().Y(), sigma);
      //      cout << "set point to fit " << chit->GetPosition().X() << " " <<  chit->GetPosition().Y() << endl;
      if(fDisplayOn) {
        display->cd(1);
        TMarker *mrk = new TMarker(hit->GetPosition().X(), hit->GetPosition().Y(), 6);
        mrk->SetMarkerColor(kRed);
        mrk->Draw("SAME");

        display->cd(2);
        chit->Draw(1);

        display->cd(2);
        TMarker *mrk2 = new TMarker(chit->GetPosition().X(), chit->GetPosition().Y(), 6);
        mrk2->SetMarkerColor(kRed);
        mrk2->Draw("SAME");


        display->Update();
        display->Modified();
      } 
    }

  
  fFitter->StraightLineFit(fitm2, fitp2);

  double xc, yc, R;
  FromConformalToRealTrack(fitm, fitp, xc, yc, R);
  cout << "previous " << xc << " " << yc << " " << R << endl;
  FromConformalToRealTrack(fitm2, fitp2, xc, yc, R);
  cout << "now " << xc << " " << yc << " " << R << endl;

  if(fDisplayOn) {
    display->cd(2);
    cout << "wanna see the line?" << endl;
    TLine *line = new TLine(-10.07, fitp2 + fitm2 * (-10.07), 10.07, fitp2 + fitm2 * (10.07));
    line->SetLineColor(2);
    line->Draw("SAME");

    display->cd(1);
    TArc *aline = new TArc(xc, yc, R);
    aline->SetFillStyle(0);
    aline->SetLineColor(2);
    aline->Draw("SAME");


    char goOnChar;
    display->Update();
    display->Modified();
    cin >> goOnChar;
  }




}

Int_t PndTrkLegendreNew::ApplyLegendre	(	PndTrkCluster *	cluster,
		double &	theta_max,
		double &	r_max
	)

Definition at line 905 of file PndTrkLegendreNew.cxx.

References ExtractLegendre(), and RePrepareLegendre().

                                                                                               {
  RePrepareLegendre(cluster);
  return ExtractLegendre(1, theta_max, r_max);
}

PndTrkLegendreNew::ClassDef ( PndTrkLegendreNew  , 1    )

private

Int_t PndTrkLegendreNew::ClusterToConformal

(

PndTrkCluster *

cluster

)

Definition at line 1466 of file PndTrkLegendreNew.cxx.

References ComputeTraAndRot(), conform, conformalhitlist, Double_t, FillConformalHitList(), FindReferenceHit(), fRefHit, fSecondary, and PndTrkConformalTransform::SetOrigin().

Referenced by Exec().

                                                                  {
  // ================ --> TO CONFORMAL PLANE
  conformalhitlist = new PndTrkConformalHitList(); // CHECK
  // translation and rotation
  Int_t nchits = 0;
  Double_t delta = 0, trasl[2] = {0., 0.};
  if(fSecondary) {
    // translation and rotation - CHECK
    //  cout << " REFERENCE HIT " << cluster->GetNofHits() << endl;
    fRefHit = FindReferenceHit(cluster);
    if(fRefHit == NULL)  {
      //          cout << "REFHIT " << fRefHit << endl;
      //        Reset();  
      return 0;
    }
    ComputeTraAndRot(fRefHit, delta, trasl);
  }
  //    
  cout << "DELTA " << delta << " TRASL " << trasl[0] << " " << trasl[1] << endl;
  conform->SetOrigin(trasl[0], trasl[1], delta);
  nchits = FillConformalHitList(cluster);
  
  return nchits;
}

void PndTrkLegendreNew::ComputePlaneExtremities

(

PndTrkCluster *

cluster

)

Definition at line 829 of file PndTrkLegendreNew.cxx.

References CAMath::ATan2(), conformalhitlist, CAMath::Cos(), fabs(), fRmax, fRmin, fUmax, fUmin, fVmax, fVmin, PndTrkConformalHitList::GetHit(), PndTrkConformalHit::GetIsochrone(), PndTrkConformalHitList::GetNofHits(), PndTrkConformalHit::GetU(), PndTrkConformalHit::GetV(), Pi, r1, r2, CAMath::Sin(), and v.

Referenced by Exec().

                                                                      {
  fUmin =  1000, fVmin =  1000, fRmin =  1000;
  fUmax = -1000, fVmax = -1000, fRmax = -1000;
  double rc_of_min, rc_of_max;
  
  
  for(int ihit = 0; ihit < conformalhitlist->GetNofHits(); ihit++) {
    PndTrkConformalHit *chit = conformalhitlist->GetHit(ihit);
    double u = chit->GetU();
    double v = chit->GetV();
    double rc = chit->GetIsochrone();
    u - rc < fUmin ? fUmin = u - rc : fUmin;
    v - rc < fVmin ? fVmin = v - rc : fVmin;
    u + rc > fUmax ? fUmax = u + rc : fUmax;
    v + rc > fVmax ? fVmax = v + rc : fVmax;
    
    double theta1 = TMath::ATan2(v, u);
    double theta2 = theta1 + TMath::Pi();
    
    double r1 = u * TMath::Cos(theta1) + v * TMath::Sin(theta1);
    double r2 = u * TMath::Cos(theta2) + v * TMath::Sin(theta2);
    
    double rimin, rimax;
    r1 < r2 ? (rimin = r1, rimax = r2) : (rimin = r2, rimax = r1);
    
    rimin < fRmin ? (rc_of_min = rc, fRmin = rimin) : fRmin;
    rimax > fRmax ? (rc_of_max = rc, fRmax = rimax) : fRmax;
  }
  
  fRmin -= rc_of_min;
  fRmax += rc_of_max;
  
  // to square the conformal plane
  double du = fUmax - fUmin;
  double dv = fVmax - fVmin;
  double delta = fabs(dv - du)/2.;
  du < dv ? (fUmin -= delta, fUmax += delta) : (fVmin -= delta, fVmax += delta);
  
  cout << "u_min " << fUmin << " u_max " << fUmax << endl;
  cout << "v_min " << fVmin << " v_max " << fVmax << endl;
  cout << "r_min " << fRmin << " r_max " << fRmax << endl;
  cout << "theta_min 0 theta_max 180" << endl;
  
  
}

void PndTrkLegendreNew::ComputeTraAndRot	(	PndTrkHit *	hit,
		Double_t &	delta,
		Double_t	trasl[2]
	)

Definition at line 820 of file PndTrkLegendreNew.cxx.

References PndTrkHit::GetPosition().

Referenced by ClusterToConformal().

                                                                                           {
  
  trasl[0] = hit->GetPosition().X();
  trasl[1] = hit->GetPosition().Y();

  delta = 0.; // TMath::ATan2(hit->GetPosition().Y() - 0., hit->GetPosition().X() - 0.); // CHECK 
  
}

Int_t PndTrkLegendreNew::CountTracksInCluster

(

PndTrkCluster *

cluster

)

Definition at line 1347 of file PndTrkLegendreNew.cxx.

Referenced by CountTracksInSkewSector().

                                                                    {
  return CountTracksInCluster(cluster, 0);
}

Int_t PndTrkLegendreNew::CountTracksInCluster	(	PndTrkCluster *	cluster,
		Int_t	where
	)

int noftracks = nofhitsinlay[tmplayid] - isneigh; if(tmplayid != -1) cout << "CLUSTER CONTAINS @ LAYER " << tmplayid << " ACTUALLY " << nofhitsinlay[tmplayid] << " - " << isneigh << " = " << noftracks << " TRACKS" << endl; if(noftracks > maxnoftracks) maxnoftracks = noftracks;

Definition at line 1356 of file PndTrkLegendreNew.cxx.

References counter, counter1, fTubeArray, PndTrkCluster::GetHit(), PndSttTube::GetLayerID(), PndTrkCluster::GetNofHits(), PndTrkHit::GetTubeID(), hit, PndSttTube::IsNeighboring(), PndTrkHit::IsSttParallel(), PndTrkHit::SetSortVariable(), and PndTrkCluster::Sort().

                                                                                 {
  // where means:
  // 0 all: parallel & skewed sectors
  // 1: only skewed

  // check how many neighboring tubes each skew 
  // tube on a layer has on that same layer
  // the total number of tubes on a layer minus the number 
  // of neighboring couples gives the number of tracks:
  // example 1 with 3 tracks:
  // OOO OO OOO are tube no.: 0 1 2   3 4  5 6 7
  // nof tubes on the layer = 8
  // calculation of neighborings:
  // 0 <--> 1
  // 1 <--> 2
  // 3 <--> 4
  // 5 <--> 6
  // 6 <--> 7
  // so, nof neigboging couples = 5
  // Then: noftubes (8) - nofcouples (5) = 3 tracks !!OK!!
  cout << "COUNT TRACKS IN SKEW SECTOR" << endl;
  int nofhitsinlay[30]; // CHECK initialize this
  for(int ilay = 0; ilay < 30; ilay++) nofhitsinlay[ilay] = 0;

  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
    PndTrkHit *hit = cluster->GetHit(ihit);
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());
    hit->SetSortVariable(tube->GetLayerID());
    nofhitsinlay[tube->GetLayerID()]++;
    //    cout << "hit " << ihit << " " << tube->GetLayerID() << " " << nofhitsinlay[tube->GetLayerID()] << endl;

  }
  cluster->Sort();

  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
    PndTrkHit *hit = cluster->GetHit(ihit);
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());
    //    cout << "SORTED " << ihit << " " << hit->GetHitID() << " " << tube->GetLayerID() << endl;
  }

  int maxnoftracks = 1;
  int tmplayid = -1;
  int counter = 0, counter1 = 0;;
  int isneigh = 0;
  // loop over cluster hits
  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
    PndTrkHit *hit = cluster->GetHit(ihit);
    counter++; 

    // which sector?
    if(where == 1 && hit->IsSttParallel() == kTRUE) continue;
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());

    int layid = tube->GetLayerID();
    if(nofhitsinlay[layid] <= 1) continue;

    // new layer?
    if(layid != tmplayid) {
      isneigh = 0;
      tmplayid = layid;
      counter1 = 0;
      //        continue; //    break;
    }
    //    cout << "hit " << ihit << " on layid " << layid << "/ " <<  nofhitsinlay[layid] << endl;

    // count processed hits 
    // in this same layer
    counter1++;

    // if it is the last hit ==> all its
    // neighborings have already been taken
    // into account
    if(counter1 == nofhitsinlay[layid]) continue;

    for(int jhit = counter; jhit < counter + nofhitsinlay[layid] -  counter1; jhit++) {
      PndTrkHit *hit2 = cluster->GetHit(jhit);
      int tubeid2 =   hit2->GetTubeID();
      PndSttTube *tube2 = (PndSttTube*) fTubeArray->At(tubeid2);



      if(tube->GetLayerID() != tube2->GetLayerID()) cout << "ERROR" << tube->GetLayerID()  << " " << tube2->GetLayerID() << endl;
      //      cout << "compare " << ihit << "(" << hit->GetHitID() << "- " << hit->GetTubeID() << ") with " << jhit << " (" << hit2->GetHitID() << "- " << tubeid2 << ") from " << counter << " to " << counter + nofhitsinlay[layid] - 1 << endl;
      if(tube->IsNeighboring(tubeid2) == kTRUE) {
        isneigh++;
        //      cout << "isneigh " << isneigh << endl;
        //      // break;
      }
    }

    // if all the hits in the layer have been processed
    if(counter1 == nofhitsinlay[layid] - 1) {
      int noftracks = nofhitsinlay[layid] - isneigh;
      cout << "CLUSTER CONTAINS @ LAYER " << layid << " ACTUALLY " << nofhitsinlay[layid] << " - " << isneigh << " = " << noftracks << " TRACKS" << endl;
      if(noftracks > maxnoftracks) maxnoftracks = noftracks;
    }

  }
    
  
  cout << "THIS CLUSTER HAS A TOTAL OF " << maxnoftracks << " TRACKS" << endl;
  return maxnoftracks;
}

Int_t PndTrkLegendreNew::CountTracksInSkewSector

(

PndTrkCluster *

cluster

)

Definition at line 1352 of file PndTrkLegendreNew.cxx.

References CountTracksInCluster().

Referenced by Exec().

                                                                       {
  return CountTracksInCluster(cluster, 1); 
}

PndTrkCluster * PndTrkLegendreNew::CreateClusterAroundTrack

(

PndTrkTrack *

track

)

Definition at line 1533 of file PndTrkLegendreNew.cxx.

References PndTrkCluster::AddHit(), conform, display, PndTrkCluster::DoesContain(), PndTrkTrack::Draw(), PndTrkCluster::Draw(), PndTrkHit::DrawTube(), fabs(), fDisplayOn, FromRealToConformalTrack(), fTubeArray, PndTrkTrack::GetCenter(), PndTrkTrack::GetCluster(), PndTrkConformalTransform::GetConformalHit(), PndTrkConformalTransform::GetConformalSttHit(), PndTrkHitList::GetHit(), PndTrkCluster::GetHit(), PndSttTube::GetLayerID(), PndTrkHitList::GetNofHits(), PndTrkCluster::GetNofHits(), PndTrkTrack::GetRadius(), PndSttTube::GetSectorID(), PndTrkHit::GetTubeID(), PndTrkConformalHit::GetU(), PndTrkConformalHit::GetV(), PndTrkHit::GetXYDistance(), hit, PndTrkHit::IsSttParallel(), R, CAMath::Sqrt(), and stthitlist.

Referenced by Exec().

                                                                              {

  double R = track->GetRadius();
  double xc = track->GetCenter().X();
  double yc = track->GetCenter().Y();
  double fitm, fitp;
  FromRealToConformalTrack(xc, yc, R, fitm, fitp);


 PndTrkCluster *cluster = track->GetCluster();

  // create cluster depending on fitting
  double rmin = R - R * 0.05; // CHECK 5%?
  double rmax = R + R * 0.05; // "      "
  
  
  if(fDisplayOn) {
    display->cd(1);
    track->Draw(kBlue);
    
    TArc *arcmin = new TArc(xc, yc, rmin);
    TArc *arcmax = new TArc(xc, yc, rmax);
      
    arcmin->SetFillStyle(0);
    arcmax->SetFillStyle(0);
    arcmin->SetLineColor(kGreen);
    arcmax->SetLineColor(kBlue);

 //    arcmin->Draw("SAME");
//     arcmax->Draw("SAME");

    display->Update();
    display->Modified();
    char goOnChar;
    cout << "want to go to new cluster?" << endl;
    cin >> goOnChar;
  }
    
  // create cluster depending on fitting
  PndTrkCluster *thiscluster = new PndTrkCluster();
  int startsecid = 1000, endsecid = -1, startlayid = 1000, endlayid = -1;
  double totaldistanceconf = 0,  chi2 = 0;
  // clean existing cluster
  for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
    PndTrkHit *hit = cluster->GetHit(ihit);
    PndTrkConformalHit *chit = NULL;
    if(hit->IsSttParallel()) chit = conform->GetConformalSttHit(hit);
    else chit = conform->GetConformalHit(hit); // CHECK if skew?
    double distanceconf = fabs((chit->GetV() - fitm * chit->GetU() - fitp)/ TMath::Sqrt(fitm * fitm + 1));
 

    double distance = hit->GetXYDistance(TVector3(xc, yc, 0.));
    //    cout << "distance " << distance << " " << rmin << " " << rmax << " " << distanceconf << endl;
    if(distance <= rmax && distance >= rmin) {
      thiscluster->AddHit(hit);
      if(fDisplayOn) {
        display->cd(1);
        hit->DrawTube(kGreen);
        display->Update();
        display->Modified();
        //        char goOnChar;
        //        cout << "want to go to next hitcluster2?" << endl;
        //        cin >> goOnChar;
      } 
      
      PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());
      
      if(tube->GetLayerID() < startlayid)  startlayid = tube->GetLayerID();
      if(tube->GetSectorID() < startsecid) startsecid = tube->GetSectorID();
      if(tube->GetLayerID() > endlayid)  endlayid = tube->GetLayerID();
      if(tube->GetSectorID() > endsecid) endsecid = tube->GetSectorID();
      
      //      cout << "hit " << ihit << " " << hit->GetHitID() << " " << endsecid << " " << endlayid << endl;
      totaldistanceconf += distanceconf;

    }
  }

  double meandistanceconf = totaldistanceconf/thiscluster->GetNofHits();

  cout << "START SECTOR " << startsecid << " END SECTOR " << endsecid << endl;
  cout << "START LAYER  " << startlayid << " END LAYER  " << endlayid << endl;

  if(fDisplayOn) {
    display->cd(1);
    thiscluster->Draw(kRed);
    display->Update();
    display->Modified();
    char goOnChar;
    cout << "want to go to next cluster1?" << endl;
    cin >> goOnChar;
  } 


  if(startlayid != 0 || endlayid != 23) {
    for(int ihit = 0; ihit < stthitlist->GetNofHits(); ihit++) {
      PndTrkHit *hit = stthitlist->GetHit(ihit);
      if(cluster->DoesContain(hit)) continue;     
      PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());
     
      double distance = hit->GetXYDistance(TVector3(xc, yc, 0.));
      if(distance <= rmax && distance >= rmin) {
        //      cout << endl;
        //      cout << "other sector " << tube->GetSectorID() << " " << tube->GetLayerID();

        if(tube->GetSectorID() == 0 || tube->GetSectorID() == 5) {
          if(startsecid != 5 && endsecid != 5 && startsecid != 0 && endsecid != 0) continue; 
        }
        else if(fabs(tube->GetSectorID() - startsecid) > 1 && fabs(tube->GetSectorID() - endsecid) > 1) continue;
          
        if(tube->GetLayerID() > startlayid && tube->GetLayerID() < endlayid) continue;
        PndTrkConformalHit *chit = NULL;
        if(hit->IsSttParallel()) chit = conform->GetConformalSttHit(hit);
        else chit = conform->GetConformalHit(hit); // CHECK if skew?
        double distanceconf = fabs((chit->GetV() - fitm * chit->GetU() - fitp)/ TMath::Sqrt(fitm * fitm + 1));
   
        //      cout << "->distance " << distance << " (" << meandistanceconf << ") " << rmin << " " << rmax << " " << distanceconf << endl;
 


        if(fDisplayOn) {
          display->cd(1);
          hit->DrawTube(kBlue);
          display->Update();
          display->Modified();
          char goOnChar;
          cout << "want to go to next?" << endl;
          // cin >> goOnChar;
        } 

        //      cout << "tubeid " << hit->GetTubeID() << " " << tube->GetLayerID() << endl;
        thiscluster->AddHit(hit);
        //      cout << " ***";


      }
    }
  }
  cout << endl;

  if(fDisplayOn) {
    display->cd(1);
    thiscluster->Draw(kRed);
    display->Update();
    display->Modified();
    char goOnChar;
    cout << "want to go to next cluster2?" << endl;
    cin >> goOnChar;
  } 
  // ---------------------------

  return thiscluster;
}

PndTrkClusterList PndTrkLegendreNew::CreateFullClusterization

(

)

Definition at line 1113 of file PndTrkLegendreNew.cxx.

References PndTrkClusterList::AddCluster(), PndTrkCluster::AddHit(), display, PndTrkCluster::DoesContain(), PndTrkHit::Draw(), fDisplayOn, fHitMap, fTubeArray, PndTrkNeighboringMap::GetCandseeds(), PndTrkCluster::GetHit(), PndSttTube::GetLayerID(), PndTrkNeighboringMap::GetNeighboringsToHit(), PndTrkHitList::GetNofHits(), PndTrkCluster::GetNofHits(), PndTrkNeighboringMap::GetSeeds(), PndTrkHit::GetTubeID(), PndTrkHit::IsUsed(), and stthitlist.

Referenced by Exec().

                                                              {
  PndTrkClusterList clusterlist;

  cout << "CLUSTERIZATION <---------------" << endl;

  // get seeds *********************************************8
  TObjArray seeds = fHitMap->GetSeeds();
  TObjArray *neighborings = NULL;
  int clusterizedhits = 0;
  // ----------------- loop over seeds
  for(int iseed = 0; iseed < seeds.GetEntriesFast(); iseed++) {
    PndTrkCluster *cluster = new PndTrkCluster();
    PndTrkHit *seedhit = (PndTrkHit*) seeds.At(iseed);
    
    // is it already used
    if(seedhit->IsUsed() == kTRUE) continue;
    
    int seedtubeID = seedhit->GetTubeID();
    PndSttTube *seedtube = (PndSttTube*) fTubeArray->At(seedtubeID);
    int seedlayerID = seedtube->GetLayerID();
    
    // add hit to cluster
    cluster->AddHit(seedhit);

    if(fDisplayOn) {
      display->cd(1);
      char goOnChar;
      cin >> goOnChar;
      cout << "SEED " << seedtubeID << endl; 
      seedhit->Draw(kRed);
      display->Update();
      display->Modified();
    }

    //    // add cluster to clusterlist
    //    clusterlist.AddCluster(cluster);

    int nlastadded = 1, addedcounter = 0;
    // cout << "nlastadded to " << seedhit->GetHitID() << "(" << seedtubeID << ")" << " " << nlastadded << endl;
    //    if(nlastadded == 0) continue;
    while(nlastadded > 0) {
      // loop on the last nlastadded hits to this cluster
      // example: add to a 5 hits cluster: 0 1 2 3 4
      // the hits no. 5, 6, 7
      // --> nlastadded = 3 & nof hits in cluster = 5 + 3 = 8
      // 7 6 5 = 8 - 3
      // here loop from hit 8 - 1 = 7 to hit 8 - 3 = 5

      addedcounter = 0;

      // cout << "@@@@@@@@@@@@@@@@@ loop on the last " << nlastadded << " hits of cluster" << endl;
      //   for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
      //        PndTrkHit *hit = cluster->GetHit(ihit);
      //        cout << " " << hit->GetHitID() ;
      //       }
      //       cout << endl;


      int nclusterhits = cluster->GetNofHits();
      for(int iadded = nclusterhits - 1; iadded >= (nclusterhits - nlastadded); iadded--) {
        PndTrkHit *addedhit = cluster->GetHit(iadded);
        neighborings = fHitMap->GetNeighboringsToHit(addedhit);
        if(neighborings->GetEntriesFast() == 0) continue;
        //      cout << "hit " << addedhit->GetHitID() << "(" << addedhit->GetTubeID() << ")" << " has " << neighborings->GetEntriesFast()  << " neighborigns: " << endl;

        // loop over the neighborings and add them all
        for(int ineigh = 0; ineigh < neighborings->GetEntriesFast(); ineigh++)
          {
            PndTrkHit *neighhit = (PndTrkHit*) neighborings->At(ineigh);
            //  cout << " " << neighhit->GetHitID() << "(" << neighhit->GetTubeID() << ")";
            if(cluster->DoesContain(neighhit) == kTRUE) {
              //  cout << "UN-ADDED, in cluster already" << endl;
              continue;
            }
            cluster->AddHit(neighhit);
            addedcounter++;
            //  cout << " - ADDED; ";

            if(fDisplayOn) {
              display->cd(1);
              char goOnChar;
              cin >> goOnChar;
              neighhit->Draw(kGreen);
                cout << "nof hits " << cluster->GetNofHits() << endl;
              // cluster->LightUp();
              display->Update();
              display->Modified();
              cin >> goOnChar;
            }
            
          }


      }
      //       cout << endl;
      nlastadded = addedcounter;
    }
    //     cout << "NEXT seed " << endl;
    // add cluster to clusterlist
    if(cluster->GetNofHits() > 3) {
      clusterlist.AddCluster(cluster); // CHECK
      clusterizedhits += cluster->GetNofHits();
    }
    
  }

  // -----------------------------------------
  cout << "NOF TOTAL HITS " << stthitlist->GetNofHits() << " NOF CLUSTERIZED HITS " <<  clusterizedhits << endl;
  if(stthitlist->GetNofHits() - clusterizedhits > 6) {
  // get candseeds *********************************************8
  TObjArray candseeds = fHitMap->GetCandseeds();
  // ----------------- loop over cand seeds

  cout << "we have " << candseeds.GetEntriesFast() << " canduidate seeds" << endl;
  for(int jseed = 0; jseed < candseeds.GetEntriesFast(); jseed++) {
    PndTrkCluster *cluster = new PndTrkCluster();
    PndTrkHit *cseedhit = (PndTrkHit*) candseeds.At(jseed);
    
    // is it already used

    if(fDisplayOn) {
      display->cd(1);
      char goOnChar;
      cin >> goOnChar;
      cout << "CSEED " << endl;
      cseedhit->Draw(kBlue);
      display->Update();
      display->Modified();
    }


    if(cseedhit->IsUsed() == kTRUE) { cout << "already" << endl ; continue; }
    
    int cseedtubeID = cseedhit->GetTubeID();
    PndSttTube *cseedtube = (PndSttTube*) fTubeArray->At(cseedtubeID);
    int cseedlayerID = cseedtube->GetLayerID();
    
    // add hit to cluster
    cluster->AddHit(cseedhit);

    if(fDisplayOn) {
      display->cd(1);
      char goOnChar;
      cin >> goOnChar;
      cout << "SEED " << cseedtubeID << endl; 
      cseedhit->Draw(kRed);
      display->Update();
      display->Modified();
    }

    //    // add cluster to clusterlist
    //    clusterlist.AddCluster(cluster);

    int nlastadded = 1, addedcounter = 0;
    // cout << "nlastadded to " << seedhit->GetHitID() << "(" << seedtubeID << ")" << " " << nlastadded << endl;
    //    if(nlastadded == 0) continue;
    while(nlastadded > 0) {
      // loop on the last nlastadded hits to this cluster
      // example: add to a 5 hits cluster: 0 1 2 3 4
      // the hits no. 5, 6, 7
      // --> nlastadded = 3 & nof hits in cluster = 5 + 3 = 8
      // 7 6 5 = 8 - 3
      // here loop from hit 8 - 1 = 7 to hit 8 - 3 = 5

      addedcounter = 0;

      // cout << "@@@@@@@@@@@@@@@@@ loop on the last " << nlastadded << " hits of cluster" << endl;
      //   for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
      //        PndTrkHit *hit = cluster->GetHit(ihit);
      //        cout << " " << hit->GetHitID() ;
      //       }
      //       cout << endl;


      int nclusterhits = cluster->GetNofHits();
      for(int iadded = nclusterhits - 1; iadded >= (nclusterhits - nlastadded); iadded--) {
        PndTrkHit *addedhit = cluster->GetHit(iadded);
        neighborings = fHitMap->GetNeighboringsToHit(addedhit);
        if(neighborings->GetEntriesFast() == 0) continue;
        //      cout << "hit " << addedhit->GetHitID() << "(" << addedhit->GetTubeID() << ")" << " has " << neighborings->GetEntriesFast()  << " neighborigns: " << endl;

        // loop over the neighborings and add them all
        for(int ineigh = 0; ineigh < neighborings->GetEntriesFast(); ineigh++)
          {
            PndTrkHit *neighhit = (PndTrkHit*) neighborings->At(ineigh);
            //  cout << " " << neighhit->GetHitID() << "(" << neighhit->GetTubeID() << ")";
            if(cluster->DoesContain(neighhit) == kTRUE) {
              //  cout << "UN-ADDED, in cluster already" << endl;
              continue;
            }
            cluster->AddHit(neighhit);
            addedcounter++;
            //  cout << " - ADDED; ";

            if(fDisplayOn) {
              display->cd(1);
              char goOnChar;
              cin >> goOnChar;
              neighhit->Draw(kGreen);
                cout << "nof hits " << cluster->GetNofHits() << endl;
              // cluster->LightUp();
              display->Update();
              display->Modified();
              cin >> goOnChar;
            }
            
          }


      }
      //       cout << endl;
      nlastadded = addedcounter;
    }   

    if(cluster->GetNofHits() > 3) {
      clusterlist.AddCluster(cluster); // CHECK
    }

  }
  }
  // -----------------------------------------
  return clusterlist;
}

PndTrkClusterList PndTrkLegendreNew::CreateFullClusterization2

(

)

 PndTrkClusterList clusterlist;

get seeds *********************************************8 TObjArray seeds = fHitMap->GetSeeds(); neighborings = NULL;

Definition at line 1337 of file PndTrkLegendreNew.cxx.

                                                               {
}

void PndTrkLegendreNew::DrawConfHit	(	double	x,
		double	y,
		double	r,
		int	marker = 2
	)

Definition at line 642 of file PndTrkLegendreNew.cxx.

References display.

Referenced by FillLegendreHisto().

                                                                            {
  display->cd(2);
  if(r >= 0) {
    TArc *arc = new TArc(u, v, r);
    arc->SetFillStyle(0);
    arc->Draw("SAME");
  }
  else {
    TMarker *mrk = new TMarker(u, v, marker);
    mrk->Draw("SAME");
  }
  //  display->Update();
  //   display->Modified();
}

void PndTrkLegendreNew::DrawGeometry

(

)

Definition at line 498 of file PndTrkLegendreNew.cxx.

References display, fTubeArray, PndSttTube::GetPosition(), hxy, and PndSttTube::IsParallel().

Referenced by Refresh().

                                     {
  if(hxy == NULL)  hxy = new TH2F("hxy", "xy plane", 100, -43, 43, 100, -43, 43);
  else hxy->Reset();
  display->cd(1);
  hxy->SetStats(kFALSE);
  hxy->Draw();

  // draw all the tubes
  for(int itube = 1; itube < fTubeArray->GetEntriesFast(); itube++) {
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(itube);
    if(tube->IsParallel()) {
      TArc * arc = new TArc(tube->GetPosition().X(), tube->GetPosition().Y(), 0.5);
      arc->SetFillStyle(0);
      arc->SetLineColor(kCyan - 10);
      arc->Draw("SAME");
    }
    else {
      TMarker *mrk = new TMarker(tube->GetPosition().X(), tube->GetPosition().Y(), 6);
   
      mrk->SetMarkerColor(kCyan - 10);
      mrk->Draw("SAME");
    }
  }
  // ............................
 
  display->Update();
  display->Modified();  
 
}

void PndTrkLegendreNew::DrawGeometryConf	(	double	x1,
		double	x2,
		double	y1,
		double	y2
	)

Definition at line 528 of file PndTrkLegendreNew.cxx.

References display, and huv.

Referenced by Exec(), LegendreFit(), and RePrepareLegendre().

                                                                                   {
  // CHECK
  char goOnChar;
  //  cout << "DRAWING GEOMETRY CONF" << endl;
  // cin >> goOnChar;
  
  // CHECK previous calculations, now not used;
  if(huv == NULL)  huv = new TH2F("huv", "uv plane", 100, x1, x2, 100, y1, y2);
  else {
    huv->Reset(); 
    huv->GetXaxis()->SetLimits(x1, x2);
    huv->GetYaxis()->SetLimits(y1, y2);
   }
  display->cd(2);
  huv->Draw();
  display->Update();
  display->Modified();  
 
}

void PndTrkLegendreNew::DrawHits

(

PndTrkHitList *

hitlist

)

Definition at line 491 of file PndTrkLegendreNew.cxx.

References display, and PndTrkHitList::Draw().

Referenced by Refresh().

                                                       {
  display->cd(1);
  hitlist->Draw();
  display->Update();
  display->Modified();
}

void PndTrkLegendreNew::DrawLegendreHisto

(

)

Definition at line 476 of file PndTrkLegendreNew.cxx.

References display, PndTrkLegendreTransform::Draw(), and legendre.

Referenced by ExtractLegendre().

                                          {
  display->cd(3);
  legendre->Draw();
  display->Update();
  display->Modified();
}

void PndTrkLegendreNew::DrawLists

(

)

Definition at line 549 of file PndTrkLegendreNew.cxx.

References counter, display, PndTrkHit::DrawTube(), fHitMap, fTubeArray, PndTrkNeighboringMap::GetCandseeds(), PndTrkNeighboringMap::GetIndivisibles(), PndSttTube::GetLayerID(), PndTrkNeighboringMap::GetNeighboringsToHit(), PndTrkNeighboringMap::GetSeeds(), PndTrkNeighboringMap::GetStandalone(), PndTrkHit::GetTubeID(), i, and Refresh().

Referenced by FillHitMap().

                                  {
  char goOnChar;
  cout << "new hit?" << endl;
  cin >> goOnChar;
  Refresh(); 
  
  
  for(int i = 0; i < fHitMap->GetStandalone().GetEntriesFast(); i++) {
    PndTrkHit *hitA = (PndTrkHit*) fHitMap->GetStandalone().At(i);
    hitA->DrawTube(kGreen);
  }
  for(int i = 0; i < fHitMap->GetSeeds().GetEntriesFast(); i++) {
    PndTrkHit *hitA = (PndTrkHit*) fHitMap->GetSeeds().At(i);
    hitA->DrawTube(kRed);
  }
  for(int i = 0; i < fHitMap->GetCandseeds().GetEntriesFast(); i++) {
    PndTrkHit *hitA = (PndTrkHit*) fHitMap->GetCandseeds().At(i);
    hitA->DrawTube(kBlue);
  }
  for(int i = 0; i < fHitMap->GetIndivisibles().GetEntriesFast(); i++) {
    PndTrkHit *hitA = (PndTrkHit*) fHitMap->GetIndivisibles().At(i);
    hitA->DrawTube(kOrange);
    TObjArray *neighs = fHitMap->GetNeighboringsToHit(hitA);
    for(int j = 0; j < neighs->GetEntriesFast(); j++) {
      PndTrkHit *hitB = (PndTrkHit*) neighs->At(j);
      TObjArray *neighs2 = fHitMap->GetNeighboringsToHit(hitA);
      if(neighs2->GetEntriesFast() > 2) {
        int counter = 0;
        PndSttTube *tubeB = (PndSttTube* ) fTubeArray->At(hitB->GetTubeID());
        for(int k = 0; k < neighs2->GetEntriesFast(); k++) {
          PndTrkHit *hitC = (PndTrkHit*) neighs2->At(k);
          PndSttTube *tubeC = (PndSttTube* ) fTubeArray->At(hitC->GetTubeID());
          cout << "tubes " << tubeB << " "<< tubeC << endl;
          if(tubeB->GetLayerID() == tubeC->GetLayerID()) continue;
          counter++;
        }
        if(counter > 2) continue;
      }
      hitB->DrawTube(kYellow);
    }
        display->Update();
        display->Modified();
        cin >> goOnChar;

  }
  
  //     for(int i = 0; i < limits.GetEntriesFast(); i++) {
  //        PndTrkHit *hitA = (PndTrkHit*) limits.At(i);
  //        hitA->DrawTube(kYellow);
  //      }
  
  display->Update();
  display->Modified();
  cin >> goOnChar;
}

void PndTrkLegendreNew::DrawNeighborings

(

)

Definition at line 607 of file PndTrkLegendreNew.cxx.

References DrawNeighboringsToHit(), PndTrkHitList::GetHit(), PndTrkHitList::GetNofHits(), hit, and stthitlist.

Referenced by FillHitMap().

                                         {
  char goOnChar;
 
  for(int ihit = 0; ihit < stthitlist->GetNofHits(); ihit++) {
    cout << "new neigh hit?" << endl;
    cin >> goOnChar;
    PndTrkHit *hit = stthitlist->GetHit(ihit);
    DrawNeighboringsToHit(hit);
    cin >> goOnChar;
  }

}

void PndTrkLegendreNew::DrawNeighboringsToHit

(

PndTrkHit *

hit

)

Definition at line 620 of file PndTrkLegendreNew.cxx.

References display, PndTrkHit::DrawTube(), fHitMap, fTubeArray, PndTrkHit::GetHitID(), PndSttTube::GetLayerID(), PndTrkNeighboringMap::GetNeighboringsToHit(), PndTrkHit::GetTubeID(), i, and Refresh().

Referenced by DrawNeighborings().

                                                            {

 
  Refresh(); 
  hit->DrawTube(kYellow);
  PndSttTube *tube = (PndSttTube*) fTubeArray->At(hit->GetTubeID());
  TObjArray *neighs = fHitMap->GetNeighboringsToHit(hit);


  cout << "HIT " << hit->GetHitID() << "(" << hit->GetTubeID() << "/" << tube->GetLayerID() << ")" << " has " << neighs->GetEntriesFast() << " neighborings: ";
  for(int i = 0; i < neighs->GetEntriesFast(); i++) {
    PndTrkHit *hitA = (PndTrkHit*) neighs->At(i);
    hitA->DrawTube(kCyan);
    PndSttTube *tubeA = (PndSttTube*) fTubeArray->At(hitA->GetTubeID());
    cout << " " << hitA->GetHitID() << "(" << hitA->GetTubeID() << "/" << tubeA->GetLayerID() << ")";
  }
  cout << endl;
  display->Update();
  display->Modified();
  
}

void PndTrkLegendreNew::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 203 of file PndTrkLegendreNew.cxx.

References PndTrkCluster::AddHit(), PndTrkTrackList::AddTrack(), AnalyticalFit2(), PndTrkNeighboringMap::Clear(), ClusterToConformal(), ComputePlaneExtremities(), CountTracksInSkewSector(), CreateClusterAroundTrack(), CreateFullClusterization(), display, PndTrkCluster::DoesContain(), PndTrkTrack::Draw(), DrawGeometryConf(), fDisplayOn, fEventCounter, fHitMap, FillHitMap(), fMvdPixelHitArray, fMvdStripHitArray, fRmax, fRmin, FromConformalToRealTrack(), FromRealToConformalTrack(), fSecondary, fSttHitArray, fThetamax, fThetamin, fTrackArray, fTrackCandArray, fUmax, fUmin, fVerbose, fVmax, fVmin, PndTrkTrack::GetCenter(), PndTrkClusterList::GetCluster(), PndTrkTrack::GetCluster(), PndTrkCluster::GetHit(), PndTrkHit::GetHitID(), PndTrkClusterList::GetNofClusters(), PndTrkCluster::GetNofHits(), PndTrkTrackList::GetNofTracks(), PndTrkTrack::GetRadius(), PndTrkTrackList::GetTrack(), hit, Initialize(), legendre, LegendreFit(), PndTrkCluster::LightUp(), nhits, R, Refresh(), RefreshConf(), Reset(), PndTrkTrack::SetCenter(), PndTrkTrack::SetCluster(), PndTrkTrack::SetRadius(), PndTrkLegendreTransform::SetUpLegendreHisto(), and track.

                                       {
  // ############## I N I T I A L I Z A T I O N S ##############
  fTrackArray->Delete();
  fTrackCandArray->Delete();
  if(fVerbose > 0) cout << "*********************** " << fEventCounter << " ***********************" << endl;
 
  // initialize -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~
  Initialize();
  if(fVerbose > 1) {
    cout << "number of stt    hits " << fSttHitArray->GetEntriesFast() << endl;
    cout << "number of mvdpix hits " << fMvdPixelHitArray->GetEntriesFast() << endl;
    cout << "number of mvdstr hits " << fMvdStripHitArray->GetEntriesFast() << endl;
  }
  // initialize display -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-
  if(fDisplayOn)  {
    Refresh();
    char goOnChar;
    display->Update();
    display->Modified();
    cout << " STARTING" << endl;
    cin >> goOnChar;
    display->Update();
    display->Modified();
  }

  // initialize hit map -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-
  fHitMap->Clear();
  FillHitMap();
  // ##########################################################

  //  fDisplayOn = kFALSE;
  // clusterization -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----
  PndTrkClusterList clusterlist = CreateFullClusterization();
  int nofclusters = clusterlist.GetNofClusters(); 
  cout << "CLUSTERLIST " << nofclusters << endl;

  if(fDisplayOn)  {
    for(int iclus = 0; iclus < nofclusters; iclus++) {
      PndTrkCluster *cluster = clusterlist.GetCluster(iclus);
      cout << "CLUSTER " << iclus << ":";
      char goOnChar;
      cin >> goOnChar;
      Refresh();
      cluster->LightUp();
      display->Update();
      display->Modified();
    }
  }




  PndTrkTrackList tracklist;

  // loop on clusterlist -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~
  for(int iclus = 0; iclus < nofclusters; iclus++) {
    PndTrkCluster *cluster = clusterlist.GetCluster(iclus);
    // if(cluster->GetNofHits() < 3) continue;

    // print and display cluster -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
    cout << "CLUSTER " << iclus << ":";
    if(fDisplayOn)  {
      char goOnChar;
      cin >> goOnChar;
      Refresh();
      cluster->LightUp();
      display->Update();
      display->Modified();
    }
      
    for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
      PndTrkHit *hit = cluster->GetHit(ihit);
      cout << " " << hit->GetHitID();
    }
    
    cout << endl;

    // count tracks in skew sector -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~
    Int_t maxnoftracks =  CountTracksInSkewSector(cluster);

    //  if(maxnoftracks > 1) continue; // CHECK

    PndTrkClusterList partialcluslist;
    cout << "\033[1;36m ITERATIONS -----------------------> " << maxnoftracks << "\033[0m" << endl;
    for(int iter = 0; iter < maxnoftracks; iter++) {
      cout << "\033[1;36m ############### ITER "  << iter << "\033[0m" << endl;
      if(fDisplayOn)  {
        char goOnChar;
        cin >> goOnChar;
        Refresh();
        cluster->LightUp();
        display->Update();
        display->Modified();
      }

      // fitting procedure -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~
      Int_t nhits = ClusterToConformal(cluster);
      if(fSecondary) {
        cout << ">>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<" << endl;
        ComputePlaneExtremities(cluster);
        
        int nofsteps = (int) ((fRmax - fRmin)/0.04);
        cout << "nofsteps " << nofsteps << endl;
        legendre->SetUpLegendreHisto(1000, fThetamin, fThetamax, 1000, fRmin, fRmax);
      }
      PndTrkTrack *track = LegendreFit(cluster);
      if(track == NULL) continue;
      PndTrkCluster *thiscluster = CreateClusterAroundTrack(track);
      cout << "nof hits " << thiscluster->GetNofHits() << endl;
      track->SetCluster(thiscluster);
      tracklist.AddTrack(track);
      cout << "ADD CLUSTER TO TRACKLIST " << thiscluster->GetNofHits() << endl;
 
      PndTrkCluster *remainingcluster = new PndTrkCluster();
      for(int ihit = 0; ihit < cluster->GetNofHits(); ihit++) {
        PndTrkHit *hit = cluster->GetHit(ihit);
        if(thiscluster->DoesContain(hit)) continue;
        remainingcluster->AddHit(hit);
      }
      cluster = remainingcluster;
      cout << "partialclusterlist " << partialcluslist.GetNofClusters() << endl;
      cout << "remainning hits " << remainingcluster->GetNofHits() << endl;

      // if there are still a lot of hits not assigned, maybe 
      // there is another track...
      if(remainingcluster->GetNofHits() > 15 && iter == maxnoftracks - 1) {
        cout << "\033[1;31m LETS TRY ANOTHER ONE -------- BONUS ------------ \033[0m" << endl;
        
        maxnoftracks++;
      }
    }


    cout << "tracklist is " << tracklist.GetNofTracks() << endl;

  }
  // ---------------------------------------
  
  // loop over the tracklist -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~
  int noftrack = tracklist.GetNofTracks(); 
  cout << "\033[1;35m (AFTER LEGENDRE AND RECLUSTERING) NEW CLUSTERLIST " << noftrack << " \033[0m" << endl;

  // fDisplayOn = kTRUE;

  PndTrkTrackList tracklist2;
  // loop on clusterlist -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~
  for(int itrk = 0; itrk < noftrack; itrk++) {
    PndTrkTrack *track = tracklist.GetTrack(itrk);
    PndTrkCluster *cluster = track->GetCluster();
    // if(cluster->GetNofHits() < 3) continue;

    // print and display cluster -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
    cout << "CLUSTER " << itrk << ":";
    if(fDisplayOn)  {
      char goOnChar;
      cin >> goOnChar;
      Refresh();
      cluster->LightUp();
      display->Update();
      display->Modified();
    }
      
    // fit with legendre the new cluster -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
    Int_t nhits = ClusterToConformal(cluster);

    track = LegendreFit(cluster); 
    if(track == NULL) continue; 

    // create cluster around new fit
    PndTrkCluster *cluster2 = CreateClusterAroundTrack(track); 

    double fitm, fitp;
    FromRealToConformalTrack(track->GetCenter().X(), track->GetCenter().Y(), track->GetRadius(), fitm, fitp);
    if(fDisplayOn) {
      // draw old conformal line
      RefreshConf();
      DrawGeometryConf(fUmin, fUmax, fVmin, fVmax);
      display->cd(2);
      TLine *line = new TLine(-10.07, fitp + fitm * (-10.07), 10.07, fitp + fitm * (10.07));
      line->Draw("SAME");
      
    }

    // fit with analytical chi2 -----~~~~~-----~~~~~-----~~~~~-----~~~~~-----~~~~~--
    // analytical LS fit
    double fitm2, fitp2;
    AnalyticalFit2(cluster2, fitm, fitp, fitm2, fitp2);

    double xc, yc, R;
    FromConformalToRealTrack(fitm2, fitp2, xc, yc, R);
    track->SetCenter(xc, yc);
    track->SetRadius(R);

    // create final cluster on the refitted track
    PndTrkCluster *thiscluster = CreateClusterAroundTrack(track);
    track->SetCluster(thiscluster);
    
    tracklist2.AddTrack(track);
  }
  //  fDisplayOn = kTRUE;

  if(fDisplayOn) {
    
    Refresh();
    for(int itrk = 0; itrk < tracklist2.GetNofTracks(); itrk++) {
      PndTrkTrack *track = tracklist2.GetTrack(itrk);
      track->Draw(kGreen);
      display->Update();
      display->Modified();
    }
  }


  if(fDisplayOn) {
    char goOnChar;
    display->Update();
    display->Modified();
    cout << "Finish? ";
    cin >> goOnChar;
  }
  Reset();

}

Int_t PndTrkLegendreNew::ExtractLegendre	(	Int_t	mode,
		double &	theta_max,
		double &	r_max
	)

Definition at line 912 of file PndTrkLegendreNew.cxx.

References conformalhitlist, PndTrkLegendreTransform::DeleteZoneAroundXYLegendre(), PndTrkLegendreTransform::DeleteZoneAroundXYZoom(), display, DrawLegendreHisto(), PndTrkLegendreTransform::DrawZoom(), PndTrkLegendreTransform::ExtractLegendreMaximum(), PndTrkLegendreTransform::ExtractZoomMaximum(), fDisplayOn, fFoundPeaks, PndTrkLegendreTransform::FillZoomHisto(), fVerbose, PndTrkConformalHitList::GetHit(), PndTrkConformalHit::GetIsochrone(), PndTrkConformalHitList::GetNofHits(), PndTrkConformalHit::GetU(), PndTrkConformalHit::GetV(), legendre, PndTrkLegendreTransform::SetUpZoomHisto(), and tr.

Referenced by ApplyLegendre(), and LegendreFit().

                                                                                      {
  if(fDisplayOn) {
    char goOnChar;
    //      cin >> goOnChar;
    DrawLegendreHisto();
    //    cout << "LEGENDRE (nof conf hits = " <<  conformalhitlist->GetNofHits() << ")" << endl;
    display->cd();
    //      cin >> goOnChar;
    display->Update();
    display->Modified();
  }
  
  // FIND MAXIMUM IN LEGENDRE HISTO

  //  legendre->ApplyThresholdLegendreHisto(0.3);
  int  maxpeak = legendre->ExtractLegendreMaximum(theta_max, r_max);

  bool alreadythere = false;
  //  if(mode == 0) 
  {

    if(maxpeak <= 3) {
      //      if(fVerbose > 1)
      cout << "\033[1;31m MAXPEAK " << maxpeak <<  ", BREAK NOW! \033[0m" << endl;
      return maxpeak;
    }
  
    for(int ialready = 0; ialready < fFoundPeaks.size(); ialready++) {
      std::pair<double, double> foundthetar = fFoundPeaks.at(ialready);
      double foundtheta = foundthetar.first;
      double foundr = foundthetar.second;
      // IF THIS PEAK WAS ALREADY FOUND, DELETE THE PEAK AND GO ON (TO AVOID INFINITE LOOPS)
      if(theta_max == foundtheta && r_max == foundr) {
        legendre->DeleteZoneAroundXYLegendre(theta_max, r_max);
        maxpeak = legendre->ExtractLegendreMaximum(theta_max, r_max);
        alreadythere = true;
        if(fVerbose > 0) cout << "OH NO! THIS PEAK IS ALREADY THERE" << endl;
        return -1;
      }
    }
    
    if(alreadythere == false) {
      std::pair<double, double> tr(theta_max, r_max);
      fFoundPeaks.push_back(tr);
    }
  }
  // ZOOM LEGENDRE HISTO
  legendre->SetUpZoomHisto(theta_max, r_max, 3, 0.005);
  //   cout << "THETA/R " << theta_max << " " << r_max << " maxpeak " << maxpeak << endl;

  for(int ihit = 0; ihit < conformalhitlist->GetNofHits(); ihit++) {
    PndTrkConformalHit *chit = conformalhitlist->GetHit(ihit);
    legendre->FillZoomHisto(chit->GetU(), chit->GetV(), chit->GetIsochrone());
  }

  if(mode == 0 && alreadythere == true) {
     cout << "THIS PEAK IS ALREADY THERE" << endl;
    legendre->DeleteZoneAroundXYZoom(theta_max, r_max);
  }

 int maxpeakzoom = legendre->ExtractZoomMaximum(theta_max, r_max);
  //  cout << "THETA/R ZOOM " << theta_max << " " << r_max <<  " maxpeakzoom " << maxpeakzoom << endl;

  if(fDisplayOn) {
    char goOnChar;
    display->cd(3);
    TMarker *mrk = new TMarker(theta_max, r_max, 29);
    mrk->Draw("SAME");
    display->cd(4);
    legendre->DrawZoom();
    mrk->Draw("SAME");
    display->Update();
    display->Modified();
    //      cin >> goOnChar;
  }
  
  return maxpeak;
}

Int_t PndTrkLegendreNew::FillConformalHitList

(

PndTrkCluster *

cluster

)

Definition at line 657 of file PndTrkLegendreNew.cxx.

References PndTrkConformalHitList::AddHit(), conform, conformalhitlist, fRefHit, PndTrkConformalTransform::GetConformalSttHit(), PndTrkCluster::GetHit(), PndTrkConformalHitList::GetNofHits(), PndTrkCluster::GetNofHits(), hit, PndTrkHit::IsSttParallel(), and PndTrkConformalHitList::SetConformalTransform().

Referenced by ClusterToConformal().

                                                                    { 
  
  // set the conformal transformation, where the
  // translation and rotation must be already set
  conformalhitlist->SetConformalTransform(conform);
    
  // loop over the cluster hits 
  // and port them to conf plane
  for(int jhit = 0; jhit < cluster->GetNofHits(); jhit++) {
    PndTrkHit *hit = cluster->GetHit(jhit);
    if(hit == fRefHit) continue;
    PndTrkConformalHit * chit = NULL;
    //    cout << "HIT " << hit->GetHitID() << " " << hit->IsSttParallel() << " " << hit->IsSttSkew() << endl;
    if(hit->IsSttParallel() == kTRUE) chit = conform->GetConformalSttHit(hit);
    else continue; // CHECK  chit = conform->GetConformalHit(hit); // CHECK
    conformalhitlist->AddHit(chit);  
    // cout << hit->GetPosition().X() << " " << hit->GetPosition().Y() << " " << hit->GetIsochrone() << " " << hit->IsSttParallel() << " " << " to CONFORMAL " << chit->GetU() << " " << chit->GetV() << " " << chit->GetIsochrone() << endl;

  }
  return conformalhitlist->GetNofHits();
}

void PndTrkLegendreNew::FillHitMap

(

)

Definition at line 1030 of file PndTrkLegendreNew.cxx.

References PndTrkNeighboringMap::AddNeighboringsToHit(), display, PndTrkHit::Draw(), DrawLists(), DrawNeighborings(), fDisplayOn, fHitMap, fTubeArray, PndTrkHitList::GetHit(), PndTrkNeighboringMap::GetIndivisiblesToHit(), PndTrkHitList::GetNofHits(), PndTrkHit::GetTubeID(), hit, PndSttTube::IsNeighboring(), PndSttTube::IsSectorLimit(), PndTrkNeighboringMap::PrintIndivisibleMap(), Refresh(), PndTrkNeighboringMap::SetOwnerValue(), and stthitlist.

Referenced by Exec().

                                   {
  TObjArray limits;
  TObjArray sector[6];
  TObjArray *neighborings = NULL;


  //  PndTrkHit *hit, hit2;
  // Loop over all hits and look for:
  // 1 - tubes limiting sectors                               std::vector< int > limithit
  // 2 - tubes with no neighborings (will not use them)       std::vector< int > standalone
  // 3 - tubes with only 1 neighboring (will serve as seed)   std::vector< int > seeds
  // 4 - tubes with only 2 neighborings, one of which is on   std::vector< int > candseeds
  //     the same layer and has neighboring (will be candidate 
  //     to serve as seed, if needed)
  // 5 - Fill the sector std::vector according to sectors     std::vector< int > sector*

  fHitMap->SetOwnerValue(kTRUE); // CHECK
  for(int ihit = 0; ihit < stthitlist->GetNofHits(); ihit++) {

    PndTrkHit  *hit = stthitlist->GetHit(ihit);
    int tubeID = hit->GetTubeID();
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(tubeID);
    if(tube->IsSectorLimit() == kTRUE) limits.Add(hit);

    neighborings = new TObjArray();
   
    for(int jhit = 0; jhit < stthitlist->GetNofHits(); jhit++) {
      if(ihit == jhit) continue;
     
      PndTrkHit *hit2 = stthitlist->GetHit(jhit);
      int tubeID2 = hit2->GetTubeID();
     
      if(tube->IsNeighboring(tubeID2) == kTRUE) neighborings->Add(hit2);
     
    }
   
    //    cout << "HIT: " << hit->GetHitID() << " has " << neighborings->GetEntriesFast() << " hits" << endl;
//     if(fDisplayOn) { 
//       Refresh();
//     }
  
    fHitMap->AddNeighboringsToHit(hit, neighborings);

//     if(fDisplayOn) {
//       char goOnChar;
//       display->Update();
//       display->Modified();
//       cout << " go on?" << endl;
//       cin >> goOnChar;
//     }
  }
  //  neighborings = NULL;
  //   delete neighborings;

  if(fDisplayOn) {
    if(1 == 1)  DrawLists();
    if(1 == 2)  DrawNeighborings();
    // ================================================
    Refresh();
    char goOnChar;
    for(int ihit = 0; ihit < stthitlist->GetNofHits(); ihit++) {
      PndTrkHit *stthit = stthitlist->GetHit(ihit);
      TObjArray *indiv = fHitMap->GetIndivisiblesToHit(stthit);
      for(int jhit = 0; jhit < indiv->GetEntriesFast(); jhit++) {
        PndTrkHit *stthit2 = (PndTrkHit*) indiv->At(jhit);
        stthit2->Draw(kOrange);
        stthit->Draw(kOrange);
      }
    }

    display->Update();
    display->Modified();
    cout << " FILLHITM AP STARTING" << endl;
    cin >> goOnChar;
    display->Update();
    display->Modified();
    // ================================================
  }
  cout << "PRINT INDIVISIBILE MAP" << endl;
  fHitMap->PrintIndivisibleMap();
}

void PndTrkLegendreNew::FillLegendreHisto

(

PndTrkCluster *

cluster

)

Definition at line 875 of file PndTrkLegendreNew.cxx.

References conformalhitlist, PndTrkCluster::DoesContain(), DrawConfHit(), fDisplayOn, PndTrkLegendreTransform::FillLegendreHisto(), PndTrkConformalHit::GetHit(), PndTrkConformalHitList::GetHit(), PndTrkConformalHit::GetIsochrone(), PndTrkConformalHitList::GetNofHits(), PndTrkConformalHit::GetU(), PndTrkConformalHit::GetV(), hit, and legendre.

Referenced by LegendreFit(), and RePrepareLegendre().

{
  // ---------------------------------------------------------------
  //  cout << "FILL LEGENDRE HISTO " << cluster->GetNofHits() << endl;

  for(int ihit = 0; ihit < conformalhitlist->GetNofHits(); ihit++) {
    PndTrkConformalHit *chit = conformalhitlist->GetHit(ihit);
    PndTrkHit *hit = chit->GetHit();
    if(cluster->DoesContain(hit) == kFALSE) continue;
    legendre->FillLegendreHisto(chit->GetU(), chit->GetV(), chit->GetIsochrone());
    if(fDisplayOn) {
      DrawConfHit(chit->GetU(), chit->GetV(), chit->GetIsochrone());
      //      cout << "conformal2: " << chit->GetU() << " " << chit->GetV() << " " << chit->GetIsochrone() << endl;   
    }
  }
}

PndTrkHit * PndTrkLegendreNew::FindMvdPixelReferenceHit

(

)

Definition at line 714 of file PndTrkLegendreNew.cxx.

References fVerbose, PndTrkHitList::GetHit(), PndTrkHit::GetHitID(), PndTrkHitList::GetNofHits(), hit, PndTrkHit::IsUsed(), and mvdpixhitlist.

Referenced by FindMvdReferenceHit().

  {
  if(mvdpixhitlist->GetNofHits() == 0) return NULL;
  // loop on mvd pix hits
  int tmphitid = -1;
  PndTrkHit *refhit = NULL;
  for(int jhit = 0; jhit < mvdpixhitlist->GetNofHits(); jhit++) {
    PndTrkHit *hit = mvdpixhitlist->GetHit(jhit);
    if(hit->IsUsed()) { 
      if(fVerbose > 1) cout << "already used V" << endl; 
      continue; 
    }
    tmphitid = jhit;
    break;
  }   
  if(tmphitid == -1)  return NULL;
  refhit = mvdpixhitlist->GetHit(tmphitid);
   if(fVerbose > 1) cout << "MVD PIXEL REFERENCE HIT " << refhit->GetHitID() << endl;
  return refhit;
}

PndTrkHit * PndTrkLegendreNew::FindMvdReferenceHit

(

)

Definition at line 756 of file PndTrkLegendreNew.cxx.

References FindMvdPixelReferenceHit(), and FindMvdStripReferenceHit().

Referenced by FindReferenceHit().

{
  PndTrkHit *refhit = NULL;
  refhit = FindMvdStripReferenceHit();
  // refhit = FindMvdPixelReferenceHit();
  if(refhit != NULL) return refhit;
  //  refhit = FindMvdStripReferenceHit();
  FindMvdPixelReferenceHit();
  return refhit;
}

PndTrkHit * PndTrkLegendreNew::FindMvdStripReferenceHit

(

)

Definition at line 735 of file PndTrkLegendreNew.cxx.

References fVerbose, PndTrkHitList::GetHit(), PndTrkHit::GetHitID(), PndTrkHitList::GetNofHits(), hit, PndTrkHit::IsUsed(), and mvdstrhitlist.

Referenced by FindMvdReferenceHit().

{
  if(mvdstrhitlist->GetNofHits() == 0) return NULL;
  // loop on mvd str hits
  int tmphitid = -1;
  PndTrkHit *refhit = NULL;
  for(int jhit = 0; jhit < mvdstrhitlist->GetNofHits(); jhit++) {
    PndTrkHit *hit = mvdstrhitlist->GetHit(jhit);
    if(hit->IsUsed()) { 
       if(fVerbose > 1) cout << "already used V" << endl; 
      continue; 
    }
    tmphitid = jhit;
    break;
  }   
  if(tmphitid == -1)  return NULL;
  refhit = mvdstrhitlist->GetHit(tmphitid);
   if(fVerbose > 1) cout << "MVD STRIP REFERENCE HIT " << refhit->GetHitID() << endl;
  return refhit;
}

PndTrkHit * PndTrkLegendreNew::FindReferenceHit

(

)

Definition at line 767 of file PndTrkLegendreNew.cxx.

References FindMvdReferenceHit(), and FindSttReferenceHit().

Referenced by ClusterToConformal().

{
  PndTrkHit *refhit = NULL;
  // refhit = FindMvdReferenceHit();
  refhit = FindSttReferenceHit();
  if(refhit != NULL) return refhit;
  // refhit = FindSttReferenceHit();
  FindMvdReferenceHit();

 return refhit;
}

PndTrkHit * PndTrkLegendreNew::FindReferenceHit

(

PndTrkCluster *

cluster

)

Definition at line 779 of file PndTrkLegendreNew.cxx.

References Double_t, fVerbose, PndTrkHit::GetDetectorID(), PndTrkCluster::GetHit(), PndTrkHit::GetHitID(), PndTrkHit::GetIsochrone(), PndTrkCluster::GetNofHits(), hit, PndTrkHit::IsStt(), and PndTrkHit::IsSttParallel().

                                                                     {
  int ntot = cluster->GetNofHits();
  //  cout << "FIND REFERENCE HIT " << ntot << endl;
  
  if(ntot == 0) return NULL;


  int tmphitid = -1;
  Double_t tmpiso = 1.;
  PndTrkHit *refhit = NULL;
  for(int jhit = 0; jhit < ntot; jhit++) {
    PndTrkHit *hit = cluster->GetHit(jhit);
    
    //   if(hit->IsUsed()) { 
    //       if(fVerbose > 1) cout << "STT hit " << jhit << "already used " << endl; 
    //       continue; }
    
    if(hit->IsStt()) {
      if(hit->IsSttParallel()) {
        if(hit->GetIsochrone() < tmpiso) {
          tmphitid = jhit;
          tmpiso = hit->GetIsochrone();
          refhit = hit;
        }
      }
    }
    else {
      tmphitid = jhit;
      break;
    }
  }

  if(tmphitid == -1)  return NULL;
  refhit = cluster->GetHit(tmphitid);
  if(fVerbose > 1) cout << "REFERENCE HIT " << refhit->GetHitID() << " " << refhit->GetDetectorID() << endl;
  return refhit;


}

PndTrkHit * PndTrkLegendreNew::FindSttReferenceHit

(

int

isec = -1

)

Definition at line 680 of file PndTrkLegendreNew.cxx.

References Double_t, fVerbose, PndTrkHitList::GetHit(), PndTrkSttHitList::GetHitFromSector(), PndTrkHit::GetIsochrone(), PndTrkHitList::GetNofHits(), PndTrkSttHitList::GetNofHitsInSector(), hit, PndTrkHit::IsSttSkew(), PndTrkHit::IsUsed(), and stthitlist.

Referenced by FindReferenceHit().

                                                          {
  int ntot = 0; 
  if(isec == -1) ntot = stthitlist->GetNofHits();
  else  ntot = stthitlist->GetNofHitsInSector(isec);
  
  if(ntot == 0) return NULL;


  int tmphitid = -1;
  Double_t tmpiso = 1.;
  PndTrkHit *refhit = NULL;
  for(int jhit = 0; jhit < ntot; jhit++) {
    PndTrkHit *hit = NULL;
    if(isec == -1) hit = stthitlist->GetHit(jhit);
    else  hit = stthitlist->GetHitFromSector(jhit, isec);
    
    if(hit->IsUsed()) { 
      if(fVerbose > 1) cout << "STT hit " << jhit << "already used " << endl; 
      continue; }
    if(hit->IsSttSkew()) continue;
    if(hit->GetIsochrone() < tmpiso) {
      tmphitid = jhit;
      tmpiso = hit->GetIsochrone();
      refhit = hit;
    }
  }   
  if(tmphitid == -1)  return NULL;
  
  // PndTrkHit *refhit =  &thishitlist[tmphitid];
  if(fVerbose > 1) cout << "STT REFERENCE HIT " <<  tmphitid << " " << refhit->GetIsochrone() << endl;
  return refhit;

}

void PndTrkLegendreNew::FromConformalToRealTrack	(	double	fitm,
		double	fitp,
		double &	x0,
		double &	y0,
		double &	R
	)

Definition at line 991 of file PndTrkLegendreNew.cxx.

References conformalhitlist, CAMath::Cos(), Double_t, PndTrkConformalHitList::GetConformalTransform(), PndTrkConformalTransform::GetRotation(), PndTrkConformalTransform::GetTranslation(), CAMath::Sin(), and sqrt().

Referenced by AnalyticalFit(), AnalyticalFit2(), Exec(), and LegendreFit().

                                                                                                            {
 // CHECK if this needs to be kept --> change xc0 to xc etc
  // center and radius
  Double_t xcrot0, ycrot0;
  ycrot0 = 1 / (2 * fitp);
  xcrot0 = - fitm * ycrot0;
  R =  sqrt(xcrot0 * xcrot0 + ycrot0 * ycrot0);
  // re-rotation and re-traslation of xc and yc
  // rotation    
  x0 = TMath::Cos(conformalhitlist->GetConformalTransform()->GetRotation())*xcrot0 - TMath::Sin(conformalhitlist->GetConformalTransform()->GetRotation())*ycrot0;
  y0 = TMath::Sin(conformalhitlist->GetConformalTransform()->GetRotation())*xcrot0 + TMath::Cos(conformalhitlist->GetConformalTransform()->GetRotation())*ycrot0;
  // traslation 
  x0 += conformalhitlist->GetConformalTransform()->GetTranslation().X();
  y0 += conformalhitlist->GetConformalTransform()->GetTranslation().Y();
}

void PndTrkLegendreNew::FromRealToConformalTrack	(	double	x0,
		double	y0,
		double	R,
		double &	fitm,
		double &	fitp
	)

|x| |c -s||xr| —> |xr| | c s||x| |y| |s c||yr| |yr| |-s c||y|

Definition at line 1007 of file PndTrkLegendreNew.cxx.

References conformalhitlist, CAMath::Cos(), Double_t, PndTrkConformalHitList::GetConformalTransform(), PndTrkConformalTransform::GetRotation(), PndTrkConformalTransform::GetTranslation(), and CAMath::Sin().

Referenced by CreateClusterAroundTrack(), and Exec().

                                                                                                           {
 // CHECK if this needs to be kept --> change xc0 to xc etc
  Double_t xcrot0, ycrot0;
 
  // traslation 
  x0 -= conformalhitlist->GetConformalTransform()->GetTranslation().X();
  y0 -= conformalhitlist->GetConformalTransform()->GetTranslation().Y();

  // rotation    
  xcrot0 = TMath::Cos(conformalhitlist->GetConformalTransform()->GetRotation())*x0 + TMath::Sin(conformalhitlist->GetConformalTransform()->GetRotation())*y0;
  ycrot0 = - TMath::Sin(conformalhitlist->GetConformalTransform()->GetRotation())*x0 + TMath::Cos(conformalhitlist->GetConformalTransform()->GetRotation())*y0;
  
  // yr = 1/ (2 p) --> p = 1 / (2 yr)
  // xr = - m / (2 p) --> m = - 2 p xr
  fitp = 1 / (2 * ycrot0);
  fitm = - 2 * fitp *  xcrot0;
}

InitStatus PndTrkLegendreNew::Init ( )

virtual

Virtual method Init

Definition at line 81 of file PndTrkLegendreNew.cxx.

References conform, display, fDisplayOn, fEventCounter, fFitter, fHitMap, PndSttMapCreator::FillTubeArray(), fMapper, fMvdPix_ConfDistLimit, fMvdPix_RealDistLimit, fMvdPixelBranch, fMvdPixelHitArray, fMvdStr_ConfDistLimit, fMvdStr_RealDistLimit, fMvdStripBranch, fMvdStripHitArray, fPersistence, fSecondary, fStt_ConfDistLimit, fStt_RealDistLimit, fSttBranch, fSttHitArray, fSttParameters, fTimer, fTrackArray, fTrackCandArray, fTubeArray, fVerbose, legendre, PndTrkLegendreTransform::SetUpLegendreHisto(), PndTrkLegendreTransform::SetUpZoomHisto(), and tools.

                                   {
  
  fEventCounter = 0;

  fMvdPix_RealDistLimit = 0.5; // CHECK limits
  fMvdStr_RealDistLimit = 0.5; // CHECK limits
  fStt_RealDistLimit =  1.5 * 0.5; // CHECK limits
  fMvdPix_ConfDistLimit = 0.003; // CHECK limits
  fMvdStr_ConfDistLimit = 0.003; // CHECK limits
  fStt_ConfDistLimit =  0.001; // CHECK limits
  if(fSecondary) {
    fMvdPix_RealDistLimit = 1.; // CHECK limits
    fMvdStr_RealDistLimit = 1.; // CHECK limits
    fMvdPix_ConfDistLimit = 0.007; // CHECK limits
    fMvdStr_ConfDistLimit = 0.007; // CHECK limits
   }
   

  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- PndTrkLegendreNew::Init: "
         << "RootManager not instantiated, return!" << endl;
    return kFATAL;
  }

  // -- HITS -------------------------------------------------
  //
  // STT
  fSttHitArray = (TClonesArray*) ioman->GetObject(fSttBranch);
  if ( ! fSttHitArray ) {
    cout << "-W- PndTrkLegendreNew::Init: "
         << "No STTHit array, return!" << endl;
    return kERROR;
  }
  //
  // MVD PIXEL
  fMvdPixelHitArray = (TClonesArray*) ioman->GetObject(fMvdPixelBranch);
  if ( !fMvdPixelHitArray){
    std::cout << "-W- PndTrkLegendreNew::Init: " << "No MVD Pixel hitArray, return!" << std::endl;
    return kERROR;
  }
  //
  // MVD STRIP
  fMvdStripHitArray = (TClonesArray*) ioman->GetObject(fMvdStripBranch);
  if ( !fMvdStripHitArray){
    std::cout << "-W- PndTrkLegendreNew::Init: " << "No MVD Strip hitArray, return!" << std::endl;
    return kERROR;
  }

  fTrackArray = new TClonesArray("PndTrack");
  fTrackCandArray = new TClonesArray("PndTrackCand");
  ioman->Register("Track", "pr", fTrackArray, fPersistence); // CHECK
  ioman->Register("TrackCand",  "pr", fTrackCandArray, fPersistence); // CHECK

  
  // ----------------------------------------   maps of STT tubes
  fMapper = new PndSttMapCreator(fSttParameters);
  fTubeArray = fMapper->FillTubeArray();
  // ----------------------------------------------------  end map
 
  if(fDisplayOn) {
    display = new TCanvas("display", "display", 0, 0, 800, 800); // CHECK
    display->Divide(2, 2);
  }

  legendre = new PndTrkLegendreTransform();

  if(fSecondary == kFALSE) legendre->SetUpLegendreHisto();
  // else it will be set up cluster by cluster
  legendre->SetUpZoomHisto();

  conform = new PndTrkConformalTransform();

  tools = new PndTrkTools();
  fFitter = new PndTrkFitter(fVerbose);
  fHitMap = new PndTrkNeighboringMap(fTubeArray);
  fTimer = new TStopwatch();

  return kSUCCESS;

}

void PndTrkLegendreNew::Initialize

(

)

Definition at line 174 of file PndTrkLegendreNew.cxx.

References PndTrkHitList::AddTCA(), fFoundPeaks, fInitDone, fMvdPixelBranch, fMvdPixelHitArray, fMvdStripBranch, fMvdStripHitArray, fSttBranch, fSttHitArray, fTubeArray, fUseMVDPix, fUseMVDStr, fUseSTT, PndTrkSttHitList::Instanciate(), PndTrkSdsHitList::InstanciatePixel(), PndTrkSdsHitList::InstanciateStrip(), MVDPIXEL, mvdpixhitlist, mvdstrhitlist, MVDSTRIP, and stthitlist.

Referenced by Exec().

                                   {
  
  stthitlist = new PndTrkSttHitList(fTubeArray);
  mvdpixhitlist = new PndTrkSdsHitList(MVDPIXEL);
  mvdstrhitlist = new PndTrkSdsHitList(MVDSTRIP);

  if(fUseSTT) { 
    stthitlist->AddTCA(FairRootManager::Instance()->GetBranchId(fSttBranch), fSttHitArray);
    stthitlist->Instanciate();
  }

  if(fUseMVDPix) {
    mvdpixhitlist->AddTCA(FairRootManager::Instance()->GetBranchId(fMvdPixelBranch), fMvdPixelHitArray);
    mvdpixhitlist->InstanciatePixel();
  }

  if(fUseMVDStr) {
    mvdstrhitlist->AddTCA(FairRootManager::Instance()->GetBranchId(fMvdStripBranch), fMvdStripHitArray);
    mvdstrhitlist->InstanciateStrip();
  }

  //   conformalhitlist = new PndTrkConformalHitList();
  fFoundPeaks.clear();

  fInitDone = kTRUE;
  //  stthitlist->PrintSectors();
}

void PndTrkLegendreNew::IntersectionFinder	(	PndTrkHit *	hit,
		double	xc,
		double	yc,
		double	R
	)

Definition at line 1835 of file PndTrkLegendreNew.cxx.

References Double_t, fabs(), fTubeArray, PndTrkHit::GetIsochrone(), PndSttTube::GetPosition(), PndTrkHit::GetTubeID(), m, point, R, PndTrkHit::SetPosition(), sqrt(), and vec.

Referenced by AnalyticalFit(), and AnalyticalFit2().

                                                                                         {

  TVector2 vec(xc, yc);
  
  // tubeID  CHECK added
    Int_t tubeID = hit->GetTubeID();
    PndSttTube *tube = (PndSttTube*) fTubeArray->At(tubeID);

    // [xp, yp] point = coordinates xy of the centre of the firing tube
    TVector2 point(tube->GetPosition().X(), tube->GetPosition().Y());
    double radius = hit->GetIsochrone();

    // the coordinates of the point are taken from the intersection
    // between the circumference from the drift time and the R radius of
    // curvature. -------------------------------------------------------
    // 2. find the intersection between the little circle and the line // R
    // 2.a
    // find the line passing throught [xc, yc] (centre of curvature) and [xp, yp] (first wire)
    // y = mx + q
    Double_t m = (point.Y() - vec.Y())/(point.X() - vec.X());
    Double_t q = point.Y() - m*point.X();

    Double_t x1 = 0, y1 = 0,
      x2 = 0, y2 = 0,
      xb1 = 0, yb1 = 0,
      xb2 = 0, yb2 = 0;

    // CHECK the vertical track
    if(fabs(point.X() - vec.X()) < 1e-6) {
      
      // 2.b
      // intersection little circle and line --> [x1, y1]
      // + and - refer to the 2 possible intersections
      // +
      x1 = point.X();
      y1 = point.Y() + sqrt(radius * radius - (x1 - point.X()) * (x1 - point.X()));
      // - 
      x2 = x1;
      y2 = point.Y() - sqrt(radius * radius - (x2 - point.X()) * (x2 - point.X()));
      
      // 2.c intersection between line and circle
      // +
      xb1 = vec.X();
      yb1 = vec.Y() + sqrt(R * R - (xb1 - vec.X()) * (xb1 - vec.X()));
      // -
      xb2 = xb1;
      yb2 = vec.Y() - sqrt(R * R - (xb2 - vec.X()) * (xb2 - vec.X()));

    }    // END CHECK
    else {
      
      // 2.b
      // intersection little circle and line --> [x1, y1]
      // + and - refer to the 2 possible intersections
      // +
      x1 = (-(m*(q - point.Y()) - point.X()) + sqrt((m*(q - point.Y()) - point.X())*(m*(q - point.Y()) - point.X()) - (m*m + 1)*((q - point.Y())*(q - point.Y()) + point.X()*point.X() - radius*radius))) / (m*m + 1);
      y1 = m*x1 + q;
      // - 
      x2 = (-(m*(q - point.Y()) - point.X()) - sqrt((m*(q - point.Y()) - point.X())*(m*(q - point.Y()) - point.X()) - (m*m + 1)*((q - point.Y())*(q - point.Y()) + point.X()*point.X() - radius*radius))) / (m*m + 1);
      y2 = m*x2 + q;
      
      // 2.c intersection between line and circle
      // +
      xb1 = (-(m*(q - vec.Y()) - vec.X()) + sqrt((m*(q - vec.Y()) - vec.X())*(m*(q - vec.Y()) - vec.X()) - (m*m + 1)*((q - vec.Y())*(q - vec.Y()) + vec.X()*vec.X() - R * R))) / (m*m + 1);
      yb1 = m*xb1 + q;
      // -
      xb2 = (-(m*(q - vec.Y()) - vec.X()) - sqrt((m*(q - vec.Y()) - vec.X())*(m*(q - vec.Y()) - vec.X()) - (m*m + 1)*((q - vec.Y())*(q - vec.Y()) + vec.X()*vec.X() - R * R))) / (m*m + 1);
      yb2 = m*xb2 + q;
    }
    
    // calculation of the distance between [xb, yb] and [xp, yp]
    Double_t distb1 = sqrt((yb1 - point.Y())*(yb1 - point.Y()) + (xb1 - point.X())*(xb1 - point.X()));
    Double_t distb2 = sqrt((yb2 - point.Y())*(yb2 - point.Y()) + (xb2 - point.X())*(xb2 - point.X()));
    
    // choice of [xb, yb]
    TVector2 xyb;
    if(distb1 > distb2) xyb.Set(xb2, yb2); 
    else xyb.Set(xb1, yb1); 

    // calculation of the distance between [x, y] and [xb. yb]
    Double_t dist1 = sqrt((xyb.Y() - y1)*(xyb.Y() - y1) + (xyb.X() - x1)*(xyb.X() - x1));
    Double_t dist2 = sqrt((xyb.Y() - y2)*(xyb.Y() - y2) + (xyb.X() - x2)*(xyb.X() - x2));

    // choice of [x, y]
    TVector2 *xy;
    if(dist1 > dist2) xy = new TVector2(x2, y2);
    else xy = new TVector2(x1, y1);   // <========= THIS IS THE NEW POINT to be used for the fit
    
    hit->SetPosition(TVector3(xy->X(), xy->Y(), 0.0));
    
    delete xy;
}

void PndTrkLegendreNew::IntersectionFinder	(	PndTrkConformalHit *	chit,
		double	fitm,
		double	fitp
	)

Definition at line 1819 of file PndTrkLegendreNew.cxx.

References fabs(), PndTrkConformalHit::GetIsochrone(), PndTrkConformalHit::GetU(), PndTrkConformalHit::GetV(), PndTrkConformalHit::SetPosition(), and CAMath::Sqrt().

                                                                                             {
  
  double xi1 = chit->GetU() + fitm * chit->GetIsochrone()/ TMath::Sqrt(fitm * fitm + 1);
  double yi1 = chit->GetV() - chit->GetIsochrone() / TMath::Sqrt(fitm * fitm + 1);
  
  double xi2 = chit->GetU() - fitm * chit->GetIsochrone() / TMath::Sqrt(fitm * fitm + 1);
  double yi2 = chit->GetV() + chit->GetIsochrone()/ TMath::Sqrt(fitm * fitm + 1);
  
  double xi = 0, yi = 0;
  
  fabs(yi1 - (fitm * xi1 + fitp)) < fabs(yi2 - (fitm * xi2 + fitp)) ? (yi = yi1, xi = xi1) : (yi = yi2, xi = xi2);
  
  chit->SetPosition(xi, yi);

}

PndTrkTrack * PndTrkLegendreNew::LegendreFit

(

PndTrkCluster *

cluster

)

Definition at line 1491 of file PndTrkLegendreNew.cxx.

References display, Double_t, DrawGeometryConf(), ExtractLegendre(), PndTrkLegendreTransform::ExtractLegendreSingleLineParameters(), fDisplayOn, FillLegendreHisto(), FromConformalToRealTrack(), fUmax, fUmin, fVmax, fVmin, PndTrkCluster::GetNofHits(), legendre, R, RefreshConf(), PndTrkLegendreTransform::ResetLegendreHisto(), and track.

Referenced by Exec().

                                                                   {

  cout << "APPLY LEGENDRE =======================" << endl;
  cout << "nof hits " << cluster->GetNofHits() << endl;

  // reset the legendre histo for a new legendre fit
  legendre->ResetLegendreHisto();
  
  if(fDisplayOn) {
    RefreshConf();
    DrawGeometryConf(fUmin, fUmax, fVmin, fVmax);
  }
  
  // fill legendre histo with the cluster hits
  FillLegendreHisto(cluster); 
  double theta_max, r_max;
  // get the peak
  int maxpeak = ExtractLegendre(1, theta_max, r_max); // CHECK mode??
          
  if(maxpeak < 4) return NULL; // CHECK
  // from theta/r to line parameters in CONFORMAL plane 
  double fitm, fitq;
  legendre->ExtractLegendreSingleLineParameters(fitm, fitq);
  if(fDisplayOn) {
    display->cd(2);
    TLine *line = new TLine(-10.07, fitq + fitm * (-10.07), 10.07, fitq + fitm * (10.07));
    line->Draw("SAME");
  }
    
  // from line parameters to center/radius in REAL plane
  Double_t xc, yc, R;
  FromConformalToRealTrack(fitm, fitq, xc, yc, R);
  cout << "\033[1;33m MAXPEAK " << maxpeak << " XR, YC, R: " << xc << " " << yc << " " << R << "\033[0m" << endl;
        //      cout << "start hit " << ihit << " " << hit->GetHitID() << " " << endsecid << " " << endlayid << endl;

  // create a track from the cluster
  PndTrkTrack *track = new PndTrkTrack(cluster, xc, yc, R);

  return track;
}

void PndTrkLegendreNew::LightCluster

(

PndTrkCluster *

cluster

)

void PndTrkLegendreNew::Refresh

(

)

Definition at line 459 of file PndTrkLegendreNew.cxx.

References DrawGeometry(), DrawHits(), fVerbose, mvdpixhitlist, mvdstrhitlist, and stthitlist.

Referenced by AnalyticalFit(), AnalyticalFit2(), DrawLists(), DrawNeighboringsToHit(), Exec(), and FillHitMap().

                               {
  // CHECK
  char goOnChar;
  //  cout << "Refresh?" << endl;
  //  cin >> goOnChar;
  //  cout << "REFRESHING" << endl;
  DrawGeometry();
  if(fVerbose)  cout << "Refresh stt" << endl;
  DrawHits(stthitlist);
  if(fVerbose)  cout << "Refresh pixel" << endl;
  DrawHits(mvdpixhitlist);
  if(fVerbose)  cout << "Refresh strip" << endl;
  DrawHits(mvdstrhitlist);
  if(fVerbose)  cout << "Refresh stop" << endl;

}

void PndTrkLegendreNew::RefreshConf

(

)

Definition at line 483 of file PndTrkLegendreNew.cxx.

Referenced by Exec(), LegendreFit(), and RePrepareLegendre().

                                    { // CHECK delete
  // CHECK
  char goOnChar;
  //  cout << "RefreshConf?" << endl;
  //  cin >> goOnChar;
  //  cout << "REFRESHING CONF" << endl;
}

void PndTrkLegendreNew::RePrepareLegendre

(

PndTrkCluster *

cluster

)

Definition at line 892 of file PndTrkLegendreNew.cxx.

References DrawGeometryConf(), fDisplayOn, FillLegendreHisto(), fUmax, fUmin, fVmax, fVmin, legendre, RefreshConf(), and PndTrkLegendreTransform::ResetLegendreHisto().

Referenced by ApplyLegendre().

                                                                {

  //    cout << "RESETTING LEGENDRE HISTO" << endl;
  legendre->ResetLegendreHisto();
 
  if(fDisplayOn) {
    RefreshConf();
    DrawGeometryConf(fUmin, fUmax, fVmin, fVmax);
  }
  // cout << "%%%%%%%%%%%%%%%%%%%% XY FINDE %%%%%%%%%%%%%%%%%%%%%%%%%%" << endl;
  FillLegendreHisto(cluster);
}

void PndTrkLegendreNew::Reset

(

)

Definition at line 427 of file PndTrkLegendreNew.cxx.

References display, fDisplayOn, fEventCounter, fInitDone, fTime, fTimer, fVerbose, mvdpixhitlist, mvdstrhitlist, and stthitlist.

Referenced by Exec().

{
  
  if(fDisplayOn) {
    char goOnChar;
    display->Update();
    display->Modified();
    cout << "Finish? ";
    cin >> goOnChar;
    cout << "FINISH" << endl;
  }
  
  if(fInitDone) {
    if(stthitlist) delete stthitlist;
    if(mvdpixhitlist)  delete mvdpixhitlist; 
    if(mvdstrhitlist)  delete mvdstrhitlist;
    if(fTimer) {
      fTimer->Stop();
      fTime += fTimer->RealTime();
      
      if(fVerbose > 0) cerr << fEventCounter << " Real time " << fTime << " s" << endl;
    }
  }

  fInitDone = kFALSE;
}

void PndTrkLegendreNew::SearchSecondaryTracks ( )

inline

Definition at line 95 of file PndTrkLegendreNew.h.

References fSecondary.

{ fSecondary = kTRUE; }

void PndTrkLegendreNew::SetParContainers

(

)

Definition at line 166 of file PndTrkLegendreNew.cxx.

References fSttParameters, and rtdb.

                                         {
  FairRuntimeDb* rtdb = FairRunAna::Instance()->GetRuntimeDb();
  fSttParameters = (PndGeoSttPar*) rtdb->getContainer("PndGeoSttPar");
}

void PndTrkLegendreNew::SwitchOnDisplay ( )

inline

Definition at line 61 of file PndTrkLegendreNew.h.

References fDisplayOn.

{ fDisplayOn = kTRUE; }

Member Data Documentation

PndTrkConformalTransform* PndTrkLegendreNew::conform

private

Definition at line 163 of file PndTrkLegendreNew.h.

Referenced by AnalyticalFit(), AnalyticalFit2(), ClusterToConformal(), CreateClusterAroundTrack(), FillConformalHitList(), and Init().

PndTrkConformalHitList* PndTrkLegendreNew::conformalhitlist

private

Definition at line 164 of file PndTrkLegendreNew.h.

Referenced by ClusterToConformal(), ComputePlaneExtremities(), ExtractLegendre(), FillConformalHitList(), FillLegendreHisto(), FromConformalToRealTrack(), and FromRealToConformalTrack().

TCanvas* PndTrkLegendreNew::display

private

Definition at line 182 of file PndTrkLegendreNew.h.

Referenced by AnalyticalFit(), AnalyticalFit2(), CreateClusterAroundTrack(), CreateFullClusterization(), DrawConfHit(), DrawGeometry(), DrawGeometryConf(), DrawHits(), DrawLegendreHisto(), DrawLists(), DrawNeighboringsToHit(), Exec(), ExtractLegendre(), FillHitMap(), Init(), LegendreFit(), and Reset().

double PndTrkLegendreNew::fDeltaThetaRad

private

Definition at line 156 of file PndTrkLegendreNew.h.

Bool_t PndTrkLegendreNew::fDisplayOn

private

Definition at line 180 of file PndTrkLegendreNew.h.

Referenced by AnalyticalFit(), AnalyticalFit2(), CreateClusterAroundTrack(), CreateFullClusterization(), Exec(), ExtractLegendre(), FillHitMap(), FillLegendreHisto(), Init(), LegendreFit(), RePrepareLegendre(), Reset(), and SwitchOnDisplay().

Int_t PndTrkLegendreNew::fEventCounter

private

Definition at line 147 of file PndTrkLegendreNew.h.

Referenced by Exec(), Init(), and Reset().

PndTrkFitter* PndTrkLegendreNew::fFitter

private

Definition at line 174 of file PndTrkLegendreNew.h.

Referenced by AnalyticalFit(), AnalyticalFit2(), and Init().

std::vector< std::pair<double, double> > PndTrkLegendreNew::fFoundPeaks

private

Definition at line 171 of file PndTrkLegendreNew.h.

Referenced by ExtractLegendre(), and Initialize().

PndTrkNeighboringMap* PndTrkLegendreNew::fHitMap

private

Definition at line 176 of file PndTrkLegendreNew.h.

Referenced by CreateFullClusterization(), DrawLists(), DrawNeighboringsToHit(), Exec(), FillHitMap(), and Init().

Bool_t PndTrkLegendreNew::fInitDone

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by Initialize(), and Reset().

PndSttMapCreator* PndTrkLegendreNew::fMapper

private

Definition at line 145 of file PndTrkLegendreNew.h.

Referenced by Init().

Double_t PndTrkLegendreNew::fMvdPix_ConfDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

Double_t PndTrkLegendreNew::fMvdPix_RealDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

char PndTrkLegendreNew::fMvdPixelBranch[200]

private

Definition at line 141 of file PndTrkLegendreNew.h.

Referenced by Init(), Initialize(), and PndTrkLegendreNew().

TClonesArray* PndTrkLegendreNew::fMvdPixelHitArray

private

Input array of MvdPixelHitArray

Definition at line 131 of file PndTrkLegendreNew.h.

Referenced by Exec(), Init(), and Initialize().

Double_t PndTrkLegendreNew::fMvdStr_ConfDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

Double_t PndTrkLegendreNew::fMvdStr_RealDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

char PndTrkLegendreNew::fMvdStripBranch[200]

private

Definition at line 141 of file PndTrkLegendreNew.h.

Referenced by Init(), Initialize(), and PndTrkLegendreNew().

TClonesArray* PndTrkLegendreNew::fMvdStripHitArray

private

Input array of MvdStripHitArray

Definition at line 133 of file PndTrkLegendreNew.h.

Referenced by Exec(), Init(), and Initialize().

Int_t PndTrkLegendreNew::fNofHits

private

Definition at line 121 of file PndTrkLegendreNew.h.

Int_t PndTrkLegendreNew::fNofMvdPixHits

private

Definition at line 121 of file PndTrkLegendreNew.h.

Int_t PndTrkLegendreNew::fNofMvdStrHits

private

Definition at line 121 of file PndTrkLegendreNew.h.

Int_t PndTrkLegendreNew::fNofSttHits

private

Definition at line 121 of file PndTrkLegendreNew.h.

Int_t PndTrkLegendreNew::fNofTriplets

private

Definition at line 121 of file PndTrkLegendreNew.h.

Bool_t PndTrkLegendreNew::fPersistence

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by Init().

PndTrkHit* PndTrkLegendreNew::fRefHit

private

Definition at line 167 of file PndTrkLegendreNew.h.

Referenced by ClusterToConformal(), and FillConformalHitList().

double PndTrkLegendreNew::fRmax

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), and Exec().

double PndTrkLegendreNew::fRmin

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), and Exec().

Bool_t PndTrkLegendreNew::fSecondary

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by ClusterToConformal(), Exec(), Init(), and SearchSecondaryTracks().

Double_t PndTrkLegendreNew::fStt_ConfDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

Double_t PndTrkLegendreNew::fStt_RealDistLimit

private

Definition at line 169 of file PndTrkLegendreNew.h.

Referenced by Init().

char PndTrkLegendreNew::fSttBranch[200]

private

Definition at line 141 of file PndTrkLegendreNew.h.

Referenced by Init(), Initialize(), and PndTrkLegendreNew().

TClonesArray* PndTrkLegendreNew::fSttHitArray

private

Input array of PndSttHit

Definition at line 127 of file PndTrkLegendreNew.h.

Referenced by Exec(), Init(), and Initialize().

Double_t PndTrkLegendreNew::fSttParalDistance

private

Definition at line 153 of file PndTrkLegendreNew.h.

PndGeoSttPar* PndTrkLegendreNew::fSttParameters

private

Definition at line 140 of file PndTrkLegendreNew.h.

Referenced by Init(), and SetParContainers().

TClonesArray* PndTrkLegendreNew::fSttPointArray

private

Input array of PndSttPoints

Definition at line 125 of file PndTrkLegendreNew.h.

Double_t PndTrkLegendreNew::fSttToMvdStripDistance

private

Definition at line 153 of file PndTrkLegendreNew.h.

double PndTrkLegendreNew::fThetamax

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by Exec().

double PndTrkLegendreNew::fThetamin

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by Exec().

double PndTrkLegendreNew::fTime

private

Definition at line 172 of file PndTrkLegendreNew.h.

Referenced by Reset().

TStopwatch* PndTrkLegendreNew::fTimer

private

Definition at line 173 of file PndTrkLegendreNew.h.

Referenced by Init(), and Reset().

TClonesArray* PndTrkLegendreNew::fTrackArray

private

Definition at line 135 of file PndTrkLegendreNew.h.

Referenced by Exec(), and Init().

TClonesArray * PndTrkLegendreNew::fTrackCandArray

private

Definition at line 135 of file PndTrkLegendreNew.h.

Referenced by Exec(), and Init().

TClonesArray* PndTrkLegendreNew::fTubeArray

private

Definition at line 137 of file PndTrkLegendreNew.h.

Referenced by CountTracksInCluster(), CreateClusterAroundTrack(), CreateFullClusterization(), DrawGeometry(), DrawLists(), DrawNeighboringsToHit(), FillHitMap(), Init(), Initialize(), and IntersectionFinder().

double PndTrkLegendreNew::fUmax

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), Exec(), LegendreFit(), and RePrepareLegendre().

double PndTrkLegendreNew::fUmin

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), Exec(), LegendreFit(), and RePrepareLegendre().

Bool_t PndTrkLegendreNew::fUseMVDPix

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by Initialize().

Bool_t PndTrkLegendreNew::fUseMVDStr

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by Initialize().

Bool_t PndTrkLegendreNew::fUseSTT

private

Definition at line 159 of file PndTrkLegendreNew.h.

Referenced by Initialize().

double PndTrkLegendreNew::fVmax

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), Exec(), LegendreFit(), and RePrepareLegendre().

double PndTrkLegendreNew::fVmin

private

Definition at line 177 of file PndTrkLegendreNew.h.

Referenced by ComputePlaneExtremities(), Exec(), LegendreFit(), and RePrepareLegendre().

TH2F* PndTrkLegendreNew::huv

private

Definition at line 183 of file PndTrkLegendreNew.h.

Referenced by DrawGeometryConf().

TH2F* PndTrkLegendreNew::hxy

private

Definition at line 181 of file PndTrkLegendreNew.h.

Referenced by DrawGeometry().

TH2F * PndTrkLegendreNew::hxz

private

Definition at line 181 of file PndTrkLegendreNew.h.

TH2F * PndTrkLegendreNew::hzphi

private

Definition at line 181 of file PndTrkLegendreNew.h.

PndTrkLegendreTransform* PndTrkLegendreNew::legendre

private

Definition at line 158 of file PndTrkLegendreNew.h.

Referenced by DrawLegendreHisto(), Exec(), ExtractLegendre(), FillLegendreHisto(), Init(), LegendreFit(), and RePrepareLegendre().

PndTrkSdsHitList* PndTrkLegendreNew::mvdpixhitlist

private

Definition at line 151 of file PndTrkLegendreNew.h.

Referenced by FindMvdPixelReferenceHit(), Initialize(), Refresh(), and Reset().

PndTrkSdsHitList* PndTrkLegendreNew::mvdstrhitlist

private

Definition at line 152 of file PndTrkLegendreNew.h.

Referenced by FindMvdStripReferenceHit(), Initialize(), Refresh(), and Reset().

PndTrkSttHitList* PndTrkLegendreNew::stthitlist

private

Definition at line 150 of file PndTrkLegendreNew.h.

Referenced by CreateClusterAroundTrack(), CreateFullClusterization(), DrawNeighborings(), FillHitMap(), FindSttReferenceHit(), Initialize(), Refresh(), and Reset().

PndTrkTools* PndTrkLegendreNew::tools

private

Definition at line 165 of file PndTrkLegendreNew.h.

Referenced by Init().

---

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

• PndTrkLegendreNew.h
• PndTrkLegendreNew.cxx