PndTrkLegendreTransform Class Reference

#include <PndTrkLegendreTransform.h>

Inheritance diagram for PndTrkLegendreTransform:

Public Member Functions
	PndTrkLegendreTransform ()
	~PndTrkLegendreTransform ()
void	SetUpLegendreHisto ()
void	SetUpLegendreHisto (double thetaNofBin, double thetaMin, double thetaMax, double rNofBin, double rMin, double rMax)
void	SetUpZoomHisto (double theta, double r, double deltatheta, double deltar)
void	SetUpZoomHisto ()
void	ResetLegendreHisto ()
void	ResetZoomHisto ()
void	FillHisto (TH2F *histo, double thetamin, double thetamax, double x, double y, double radius)
void	FillLegendreHisto (double x, double y, double radius)
void	FillZoomHisto (double x, double y, double radius)
void	ApplyThresholdLegendreHisto (double threshold)
int	ExtractLegendreMaximum (double &theta_max, double &r_max)
void	ExtractLegendreMaxima (int nmaxima, std::vector< double > &theta_max, std::vector< double > &r_max, std::vector< int > &content_max)
int	ExtractZoomMaximum (double &theta_max, double &r_max)
int	ExtractMaximumFromHisto (TH2F *histo, double &theta_max, double &r_max)
void	ExtractMaximaFromHisto (int nmaxima, TH2F *histo, std::vector< double > &theta_max, std::vector< double > &r_max, std::vector< int > &content_max)
void	ExtractLineParameters (double theta, double r, double &slope, double &intercept)
void	ExtractLegendreSingleLineParameters (double &slope, double &intercept)
void	ExtractZoomSingleLineParameters (double &slope, double &intercept)
void	DeleteZoneAroundXYLegendre (double x, double y)
void	DeleteZoneAroundXYZoom (double x, double y)
void	DeleteZoneAroundXY (TH2F *histo, double x, double y)
void	Draw ()
void	DrawZoom ()

Private Attributes
Double_t	fThetaNofBin
Double_t	fThetaMin
Double_t	fThetaMax
Double_t	fRNofBin
Double_t	fRMin
Double_t	fRMax
Double_t	fThetaMinZoom
Double_t	fThetaMaxZoom
Double_t	fRMinZoom
Double_t	fRMaxZoom
TH2F *	fhLegendre
TH2F *	fhLegendreZoom

Detailed Description

Definition at line 7 of file PndTrkLegendreTransform.h.

Constructor & Destructor Documentation

PndTrkLegendreTransform::PndTrkLegendreTransform

(

)

Definition at line 19 of file PndTrkLegendreTransform.cxx.

References fThetaMax, and fThetaMin.

                                                 {
  fThetaNofBin = 180; 
  fThetaMin    = 0;    
  fThetaMax    = 180;  
  fRNofBin     = 1000; 
  fRMin        = -0.07; 
  fRMax        = 0.07;  


  fThetaMinZoom    = 0;    
  fThetaMaxZoom    = 180;  
  fRMinZoom        = -0.07; 
  fRMaxZoom        = 0.07;  

  fhLegendre = NULL;
  fhLegendreZoom = NULL;
}

PndTrkLegendreTransform::~PndTrkLegendreTransform

(

)

Definition at line 40 of file PndTrkLegendreTransform.cxx.

                                                  { 

}

Member Function Documentation

void PndTrkLegendreTransform::ApplyThresholdLegendreHisto

(

double

threshold

)

void PndTrkLegendreTransform::DeleteZoneAroundXY	(	TH2F *	histo,
		double	x,
		double	y
	)

Definition at line 237 of file PndTrkLegendreTransform.cxx.

                                                                                {

  //  int bin = histo->GetMaximumBin();
  //  int bin1, bin2, bin3;
  //  histo->GetBinXYZ(bin, bin1, bin2, bin3);
  // cout << "MAX BIN BEFORE " << bin << " " << bin1 << " " << bin2 << " " << bin3 << endl;

  int binx = (int) ((x -  histo->GetXaxis()->GetXmin())/histo->GetXaxis()->GetBinWidth(1)) + 1;
  int biny = (int) ((y -  histo->GetYaxis()->GetXmin())/histo->GetYaxis()->GetBinWidth(1)) + 1;
  
  // delete the bin and a line and a row around it
  histo->SetBinContent(binx,     biny,     0);
  histo->SetBinContent(binx + 1, biny,     0);
  histo->SetBinContent(binx - 1, biny,     0);
  histo->SetBinContent(binx + 1, biny + 1, 0);
  histo->SetBinContent(binx,     biny + 1, 0);
  histo->SetBinContent(binx - 1, biny + 1, 0);
  histo->SetBinContent(binx + 1, biny - 1, 0);
  histo->SetBinContent(binx,     biny - 1, 0);
  histo->SetBinContent(binx - 1, biny - 1, 0);
   
  // cout << "deleting binx " << binx << " " << biny << endl;
  //  bin = histo->GetMaximumBin();
  //  histo->GetBinXYZ(bin, bin1, bin2, bin3);
  //  cout << "MAX BIN AFTER " << bin << " " << bin1 << " " << bin2 << " " << bin3 << endl;
}

void PndTrkLegendreTransform::DeleteZoneAroundXYLegendre	(	double	x,
		double	y
	)

Definition at line 228 of file PndTrkLegendreTransform.cxx.

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

                                                                           {
  DeleteZoneAroundXY(fhLegendre, x, y);
}

void PndTrkLegendreTransform::DeleteZoneAroundXYZoom	(	double	x,
		double	y
	)

Definition at line 232 of file PndTrkLegendreTransform.cxx.

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

                                                                       {
 DeleteZoneAroundXY(fhLegendreZoom, x, y);
}

void PndTrkLegendreTransform::Draw

(

)

Definition at line 269 of file PndTrkLegendreTransform.cxx.

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

                                   {
  fhLegendre->Draw("colz");
}

void PndTrkLegendreTransform::DrawZoom

(

)

Definition at line 273 of file PndTrkLegendreTransform.cxx.

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

                                       {
  fhLegendreZoom->Draw("colz");
}

void PndTrkLegendreTransform::ExtractLegendreMaxima	(	int	nmaxima,
		std::vector< double > &	theta_max,
		std::vector< double > &	r_max,
		std::vector< int > &	content_max
	)

Definition at line 115 of file PndTrkLegendreTransform.cxx.

Referenced by PndTrkTrackFinder::Exec().

                                                                                                                                                           {
  ExtractMaximaFromHisto(nmaxima, fhLegendre, theta_max, r_max, content_max);
}

int PndTrkLegendreTransform::ExtractLegendreMaximum	(	double &	theta_max,
		double &	r_max
	)

Definition at line 111 of file PndTrkLegendreTransform.cxx.

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

                                                                                    {
  return ExtractMaximumFromHisto(fhLegendre, theta_max, r_max);
}

void PndTrkLegendreTransform::ExtractLegendreSingleLineParameters	(	double &	slope,
		double &	intercept
	)

Definition at line 120 of file PndTrkLegendreTransform.cxx.

References r, and theta.

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

                                                                                                  {
  double theta, r;
  ExtractLegendreMaximum(theta, r);
  ExtractLineParameters(theta, r, slope, intercept);
}

void PndTrkLegendreTransform::ExtractLineParameters	(	double	theta,
		double	r,
		double &	slope,
		double &	intercept
	)

Definition at line 216 of file PndTrkLegendreTransform.cxx.

References CAMath::Sin(), and CAMath::Tan().

Referenced by PndTrkTrackFinder::Exec().

                                                                                                            {
  // r = x cost + y sint --> y = r/sint - x/tant

  if(theta == 0.) theta = 1.e-6; // CHECK

  slope = - 1./TMath::Tan(TMath::DegToRad() * theta);
  intercept = r/TMath::Sin(TMath::DegToRad() * theta);

}

void PndTrkLegendreTransform::ExtractMaximaFromHisto	(	int	nmaxima,
		TH2F *	histo,
		std::vector< double > &	theta_max,
		std::vector< double > &	r_max,
		std::vector< int > &	content_max
	)

Definition at line 176 of file PndTrkLegendreTransform.cxx.

                                                                                                                                                                         {

  // PEAK Finder 2D
  int maxbins[nmaxima];
  for(int ibin = 0; ibin < nmaxima; ibin++) maxbins[ibin] = 0;
  
  for(int ibinx = 1; ibinx <= histo->GetNbinsX() ; ibinx++) {
    for(int ibiny = 1; ibiny <= histo->GetNbinsY() ; ibiny++) {
      
      int bincontent = histo->GetBinContent(ibinx, ibiny);
      if(bincontent == 0) continue; 
      
      int ibintotal = histo->GetBin(ibinx, ibiny);
      
      //      cout << "current bin: " << ibintotal << " contains " << bincontent << " entries" << endl;
      
      for(int ibinm = nmaxima - 1; ibinm >= 0; ibinm--) {
        // cout << "ibinm " << ibinm << " " << bincontent << " " << maxbins[ibinm] << " " << histo->GetBinContent(maxbins[ibinm]) << endl;
        if(bincontent >= histo->GetBinContent(maxbins[ibinm])) {
                  for(int jbinm = 0; jbinm < ibinm; jbinm++) maxbins[jbinm] = maxbins[jbinm + 1];
          maxbins[ibinm] = ibintotal;
          break;
        }
      }
    }
  }

  for(int ibin = 0; ibin < nmaxima; ibin++) {
    int bin = maxbins[ibin];
    int binx, biny, binz;
    histo->GetBinXYZ(bin, binx, biny, binz);
    theta_max.push_back(histo->GetXaxis()->GetBinCenter(binx));
    r_max.push_back(histo->GetYaxis()->GetBinCenter(biny));
    content_max.push_back(histo->GetBinContent(bin));

 }
   
}

int PndTrkLegendreTransform::ExtractMaximumFromHisto	(	TH2F *	histo,
		double &	theta_max,
		double &	r_max
	)

Definition at line 166 of file PndTrkLegendreTransform.cxx.

                                                                                                  {
  int bin = histo->GetMaximumBin();
  int binx, biny, binz;
  histo->GetBinXYZ(bin, binx, biny, binz);
  theta_max = histo->GetXaxis()->GetBinCenter(binx);
  r_max = histo->GetYaxis()->GetBinCenter(biny);
  return histo->GetMaximum();
  //    cout << "== THETA, R max ======== " << theta_max << " " << r_max << endl;
}

int PndTrkLegendreTransform::ExtractZoomMaximum	(	double &	theta_max,
		double &	r_max
	)

Definition at line 155 of file PndTrkLegendreTransform.cxx.

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

                                                                                {
  return ExtractMaximumFromHisto(fhLegendreZoom, theta_max, r_max);
}

void PndTrkLegendreTransform::ExtractZoomSingleLineParameters	(	double &	slope,
		double &	intercept
	)

Definition at line 159 of file PndTrkLegendreTransform.cxx.

References r, and theta.

                                                                                              {
  double theta, r;
  ExtractZoomMaximum(theta, r);
  ExtractLineParameters(theta, r, slope, intercept);
}

void PndTrkLegendreTransform::FillHisto	(	TH2F *	histo,
		double	thetamin,
		double	thetamax,
		double	x,
		double	y,
		double	radius
	)

Definition at line 70 of file PndTrkLegendreTransform.cxx.

References CAMath::Cos(), Double_t, r, r1, r2, and CAMath::Sin().

                                                                                                                        {
  Double_t  thetaRad = thetamin * TMath::DegToRad();
  double deltaThetaRad = TMath::DegToRad() * (thetamax - thetamin)/fThetaNofBin;
 

  // computation
  //double deltar = deltaThetaRad * (y * TMath::Cos(thetaRad) - x * TMath::Sin(thetaRad)); //[R.K. 01/2017] unused variable
  //  cout << "cfr delta r " << x << " " << y << " " << radius << " " << deltar << " " << (fRMax - fRMin)/fRNofBin << " " << 2 * radius << endl; 


  while(thetaRad < thetamax * TMath::DegToRad()) {
    double thetaDeg = thetaRad * TMath::RadToDeg();
        
    // STT
    if(radius > 0.) {
      double r1 = x * TMath::Cos(thetaRad) + y * TMath::Sin(thetaRad) + radius;
      double r2 = x * TMath::Cos(thetaRad) + y * TMath::Sin(thetaRad) - radius;
      //   histo->Fill(thetaDeg, r1, radius);
      //   histo->Fill(thetaDeg, r2, radius);
      histo->Fill(thetaDeg, r1);
      histo->Fill(thetaDeg, r2);
      }
    else { // MVD, SCITIL, GEM
      double r = x * TMath::Cos(thetaRad) + y * TMath::Sin(thetaRad);
      histo->Fill(thetaDeg, r);
    } 
    
    thetaRad += deltaThetaRad;
    
  }
}

void PndTrkLegendreTransform::FillLegendreHisto	(	double	x,
		double	y,
		double	radius
	)

Definition at line 102 of file PndTrkLegendreTransform.cxx.

References fThetaMax, and fThetaMin.

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

                                                                                 {
  FillHisto(fhLegendre, fThetaMin, fThetaMax, x, y, radius);
}

void PndTrkLegendreTransform::FillZoomHisto	(	double	x,
		double	y,
		double	radius
	)

Definition at line 151 of file PndTrkLegendreTransform.cxx.

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

                                                                             {
  FillHisto(fhLegendreZoom, fThetaMinZoom, fThetaMaxZoom, x, y, radius);
}

void PndTrkLegendreTransform::ResetLegendreHisto

(

)

Definition at line 64 of file PndTrkLegendreTransform.cxx.

Referenced by PndTrkTrackFinder::Exec(), PndTrkLegendreNew::LegendreFit(), PndTrkTrackFinder::LegendreFit(), PndTrkLegendreSecTask::PrepareLegendre(), PndTrkLegendreTask::PrepareLegendre(), PndTrkLegendreSecTask2::PrepareLegendre(), PndTrkLegendreNew::RePrepareLegendre(), PndTrkTrackFinder::RePrepareLegendre(), PndTrkLegendreSecTask::RePrepareLegendre(), PndTrkLegendreTask::RePrepareLegendre(), and PndTrkLegendreSecTask2::RePrepareLegendre().

                                                 {
  fhLegendre->Reset();
  fhLegendre->SetStats(kFALSE);
}

void PndTrkLegendreTransform::ResetZoomHisto

(

)

Definition at line 146 of file PndTrkLegendreTransform.cxx.

                                             {
  fhLegendreZoom->Reset();
  fhLegendreZoom->SetStats(kFALSE);
}

void PndTrkLegendreTransform::SetUpLegendreHisto

(

)

Definition at line 46 of file PndTrkLegendreTransform.cxx.

References fThetaMax, and fThetaMin.

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

                                                 {
 fhLegendre = new TH2F("fhLegendre", "", fThetaNofBin, fThetaMin, fThetaMax, fRNofBin, fRMin, fRMax);
}

void PndTrkLegendreTransform::SetUpLegendreHisto	(	double	thetaNofBin,
		double	thetaMin,
		double	thetaMax,
		double	rNofBin,
		double	rMin,
		double	rMax
	)

Definition at line 50 of file PndTrkLegendreTransform.cxx.

References fThetaMax, and fThetaMin.

                                                                                                                                               {
  fThetaNofBin = thetaNofBin;
  fThetaMin    = thetaMin;
  fThetaMax    = thetaMax;
  fRNofBin     = rNofBin;
  fRMin        = rMin;
  fRMax        = rMax;
  if(fhLegendre) {
    fhLegendre->SetBins(fThetaNofBin, fThetaMin, fThetaMax, fRNofBin, fRMin, fRMax);
    ResetLegendreHisto();
  }
  else SetUpLegendreHisto();
}

void PndTrkLegendreTransform::SetUpZoomHisto	(	double	theta,
		double	r,
		double	deltatheta,
		double	deltar
	)

Definition at line 131 of file PndTrkLegendreTransform.cxx.

Referenced by PndTrkLegendreNew::ExtractLegendre(), PndTrkTrackFinder::ExtractLegendre(), PndTrkLegendreTask::ExtractLegendre(), PndTrkLegendreSecTask::ExtractLegendre(), PndTrkLegendreSecTask2::ExtractLegendre(), PndTrkLegendreSecTask2::Init(), PndTrkLegendreNew::Init(), PndTrkLegendreSecTask::Init(), and PndTrkLegendreTask::Init().

                                                                                                     {
  if(fhLegendreZoom) {
    ResetZoomHisto();
  }

  fThetaMinZoom = theta - deltatheta;
  fThetaMaxZoom = theta + deltatheta;
  fRMinZoom = r - deltar;
  fRMaxZoom = r + deltar;
  
  fhLegendreZoom->GetXaxis()->SetLimits(fThetaMinZoom, fThetaMaxZoom);
  fhLegendreZoom->GetYaxis()->SetLimits(fRMinZoom, fRMaxZoom);

}

void PndTrkLegendreTransform::SetUpZoomHisto

(

)

Definition at line 127 of file PndTrkLegendreTransform.cxx.

                                             {
  fhLegendreZoom = new TH2F("fhLegendreZoom", "", fThetaNofBin, fThetaMinZoom, fThetaMaxZoom, fRNofBin, fRMin, fRMax);
}

Member Data Documentation

TH2F* PndTrkLegendreTransform::fhLegendre

private

Definition at line 51 of file PndTrkLegendreTransform.h.

TH2F * PndTrkLegendreTransform::fhLegendreZoom

private

Definition at line 51 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fRMax

private

Definition at line 49 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fRMaxZoom

private

Definition at line 50 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fRMin

private

Definition at line 49 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fRMinZoom

private

Definition at line 50 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fRNofBin

private

Definition at line 49 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fThetaMax

private

Definition at line 49 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fThetaMaxZoom

private

Definition at line 50 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fThetaMin

private

Definition at line 49 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fThetaMinZoom

private

Definition at line 50 of file PndTrkLegendreTransform.h.

Double_t PndTrkLegendreTransform::fThetaNofBin

private

Definition at line 49 of file PndTrkLegendreTransform.h.

---

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

• PndTrkLegendreTransform.h
• PndTrkLegendreTransform.cxx