PndSttGeometryMap Class Reference

#include <PndSttGeometryMap.h>

Inheritance diagram for PndSttGeometryMap:

Public Member Functions
	PndSttGeometryMap ()
	PndSttGeometryMap (TClonesArray *tubearray, Int_t geoType)
	~PndSttGeometryMap ()
	PndSttGeometryMap (const PndSttGeometryMap &)=delete
PndSttGeometryMap &	operator= (const PndSttGeometryMap &)=delete
void	SetGeneralParameters ()
void	GenerateStrawMap (Int_t map)
Bool_t	FillGeometryParameters ()
void	SetGeneralParametersGeoType1 ()
void	GenerateStrawMapGeoType1 (Int_t map)
Bool_t	FillGeometryParametersGeoType1 ()
void	GenerateStrawMapTubeIDGeoType1 ()
void	GenerateStrawMapAngleGeoType1 ()
int	GetRow (int strawindex) const
int	GetSector (int strawindex) const
const vector< int > &	GetStrawRow (int sector, int row) const
const vector< vector< int > > &	GetStrawSector (int sector) const
bool	IsEdgeStraw (int strawindex) const
int	IsSectorBorderStraw (int strawindex) const
bool	IsAxialStraw (int strawindex) const
bool	IsSkewedStraw (int strawindex) const
bool	IsAxialRow (int rowindex) const
bool	IsSkewedRow (int rowindex) const
TArrayI	FindNeighborings (PndSttTube *tube)
TArrayI	FindNeighborings (int tubeId)
TArrayI	GetNeighboringsByMap (int tubeId)
void	FillStrawNeighborsMap ()
Double_t	CalculateStrawPoca (PndSttHit *hit1, PndSttHit *hit2, TVector3 &poca)
double	GetAngleBetweenTubes (int tubeID1, int tubeID2) const
bool	InStraightLine (int tube1, int tube2, int tube3) const
void	SetVerbose (int ver)

Private Member Functions
void	GenerateAngles ()

Private Attributes
Int_t	fGeoType
Int_t	fVerbose
map< int, vector< vector< int > > >	fStrawIndex
map< int, TArrayI >	fStrawNeighbors
vector< int >	fSectorOfStraw
vector< int >	fLayerOfStraw
vector< bool >	fAxialStraw
vector< double >	fSectorStart
vector< double >	fSectorEnd
bool	fStrawMapInitialized
Int_t	fNLayers
Int_t	fNSectors
Int_t	fNTubes
Int_t	fNTubes_inner_parallel
Int_t	fNTubes_outer_parallel
Int_t	fNTubes_fillup_parallel
Int_t	fNTubes_skewed
Int_t	fNLayers_inner_parallel
Int_t	fNLayers_skew
Int_t	fNLayers_outer_parallel
Int_t	fNLayers_fillup_parallel
Int_t **	fStartTube
Int_t **	fEndTube
Int_t **	fShift
Int_t **	fShiftSkew
TClonesArray *	fTubeArray

Detailed Description

Definition at line 24 of file PndSttGeometryMap.h.

Constructor & Destructor Documentation

PndSttGeometryMap::PndSttGeometryMap

(

)

Definition at line 53 of file PndSttGeometryMap.cxx.

: fGeoType(-1), fVerbose(0), fTubeArray(0) {}

PndSttGeometryMap::PndSttGeometryMap	(	TClonesArray *	tubearray,
		Int_t	geoType
	)

Definition at line 56 of file PndSttGeometryMap.cxx.

References fGeoType, GenerateStrawMap(), and SetGeneralParameters().

                                                                           : fGeoType(geoType), fVerbose(0), fTubeArray(tubearray) {
  if(fGeoType != 1) cout << "ERROR PndSttGeometryMap: fGeoType " << fGeoType << " not supported" << endl;
  SetGeneralParameters();
  GenerateStrawMap(1);
}

PndSttGeometryMap::~PndSttGeometryMap

(

)

Definition at line 62 of file PndSttGeometryMap.cxx.

{}

PndSttGeometryMap::PndSttGeometryMap ( const PndSttGeometryMap &  )

delete

Member Function Documentation

Double_t PndSttGeometryMap::CalculateStrawPoca	(	PndSttHit *	hit1,
		PndSttHit *	hit2,
		TVector3 &	poca
	)

Definition at line 252 of file PndSttGeometryMap.cxx.

References CAMath::Abs(), Double_t, fTubeArray, fVerbose, PndSttTube::GetHalfLength(), PndSttHit::GetTubeID(), PndSttTube::GetWireDirection(), m, p1, p2, r, t1, and t2.

Referenced by PndSttSkewedCombineTask::Exec().

{
        TVector3 p1(hit1->GetX(), hit1->GetY(), hit1->GetZ());
        TVector3 p2(hit2->GetX(), hit2->GetY(), hit2->GetZ());


        PndSttTube* tube1 = (PndSttTube*)fTubeArray->At(hit1->GetTubeID());
        PndSttTube* tube2 = (PndSttTube*)fTubeArray->At(hit2->GetTubeID());

        TVector3 u1(tube1->GetWireDirection());
        TVector3 u2(tube2->GetWireDirection());

        if (fVerbose > 1) {
                std::cout << "PndSttGeometryMap::CalculateStrawPoca Hit1 center: (" << hit1->GetX() << "/" << hit1->GetY() << "/" << hit1->GetZ()
                                << ") Dir: (" << u1.x() << "/" << u1.y() << "/" << u1.z() << ") HL: " << tube1->GetHalfLength() << std::endl;

                std::cout << "PndSttGeometryMap::CalculateStrawPoca Hit2 center: (" << hit2->GetX() << "/" << hit2->GetY() << "/" << hit2->GetZ()
                                << ") Dir: (" << u2.x() << "/" << u2.y() << "/" << u2.z() << ") HL: " << tube2->GetHalfLength() << std::endl;
        }
        TVector3 p21 = p2 - p1;
            if (p21.Mag() < 0.000001) {
                                poca = p1;
                                if (fVerbose > 1) std::cout << "Poca: " << poca.x() << "/" << poca.y() << "/" << poca.z() << std::endl;
                                return 0;
                }
//      p21.Print();
                TVector3 m = u2.Cross(u1);
                if (m.Mag() < 0.000001) {
                                poca = p1 + 0.5*p21;
                                if (fVerbose > 1) std::cout << "Poca: " << poca.x() << "/" << poca.y() << "/" << poca.z() << std::endl;
                                return p21.Mag();
                }
//      m.Print();
                TVector3 mnorm = (1.0/m.Mag2())*m;
                TVector3 r = p21.Cross(mnorm);
                Double_t t1 = r.Dot(u2);
                Double_t t2 = r.Dot(u1);

                if (t1 > tube1->GetHalfLength()){
                        t1 = tube1->GetHalfLength();
                } else if (t1 < (-1)*tube1->GetHalfLength()) {
                        t1 = (-1)*tube1->GetHalfLength();
                }

                if (t2 > tube2->GetHalfLength()){
                        t2 = tube2->GetHalfLength();
                } else if (t2 < (-1)*tube2->GetHalfLength()) {
                        t2 = (-1)*tube1->GetHalfLength();
                }

                TVector3 q1 = p1+t1*u1;
                TVector3 q2 = p2+t2*u2;
                TVector3 q21 = q2 - q1;
                poca = q1 + 0.5*q21;

                if (fVerbose > 1) std::cout << "q1: " << q1.x() << "/" << q1.y() << "/" << q1.z() << " t1: " << t1 << std::endl;
                if (fVerbose > 1) std::cout << "q2: " << q2.x() << "/" << q2.y() << "/" << q2.z() << " t2: " << t2 << std::endl;

                if (fVerbose > 1) std::cout << "Poca: " << poca.x() << "/" << poca.y() << "/" << poca.z() << std::endl;
//      cout << "Crosspoints: " << endl;
//      q1.Print();
//      q2.Print();
//      cout << "Poca: " << endl;
//      poca.Print();
//              return (TMath::Abs(p21.Dot(m))/m.Mag());

           return (TMath::Abs(p21.Dot(m))/m.Mag());
}

Bool_t PndSttGeometryMap::FillGeometryParameters

(

)

Definition at line 74 of file PndSttGeometryMap.cxx.

References fGeoType, and FillGeometryParametersGeoType1().

Referenced by PndSttMapCreator::FillTubeArrayGeoType1(), and PndSttMapCreator::FillTubeArrayGeoType2().

                                                 {
  if(fGeoType == 1) return FillGeometryParametersGeoType1();
  else return kFALSE;
}

Bool_t PndSttGeometryMap::FillGeometryParametersGeoType1

(

)

Definition at line 106 of file PndSttGeometryMap.cxx.

References FindNeighborings(), fTubeArray, IsEdgeStraw(), IsSectorBorderStraw(), PndSttTube::SetLayerLimitFlag(), PndSttTube::SetNeighborings(), and PndSttTube::SetSectorLimitFlag().

Referenced by FillGeometryParameters().

                                                         {

  for (int itube = 1; itube < fTubeArray->GetEntriesFast(); itube++) {
    PndSttTube* tube = (PndSttTube*) fTubeArray->At(itube);
    tube->SetLayerLimitFlag(IsEdgeStraw(itube));
    tube->SetSectorLimitFlag(IsSectorBorderStraw(itube));

    bool found = true;
    tube->SetNeighborings(FindNeighborings(tube));
    if(!found) cout << "ERROR IN FINDIN NEIGHBORING " << endl;
  }
  return kTRUE;
}

void PndSttGeometryMap::FillStrawNeighborsMap

(

)

Definition at line 237 of file PndSttGeometryMap.cxx.

References FindNeighborings(), fNTubes, fStrawNeighbors, and i.

Referenced by GetNeighboringsByMap().

{
        for (int i = 1; i < fNTubes+1; i++){
                fStrawNeighbors[i] = FindNeighborings(i);
        }
}

TArrayI PndSttGeometryMap::FindNeighborings

(

PndSttTube *

tube

)

Definition at line 120 of file PndSttGeometryMap.cxx.

References fNLayers, fNLayers_inner_parallel, fNSectors, fTubeArray, PndSttTube::GetDistance(), PndSttTube::GetLayerID(), PndSttTube::GetSectorID(), GetStrawRow(), PndSttTube::GetTubeID(), PndSttTube::IsParallel(), and PndSttTube::IsSkew().

Referenced by FillGeometryParametersGeoType1(), FillStrawNeighborsMap(), and FindNeighborings().

                                                            {
        double tolerance0 = 1.2; // parallel/parallel && inner parallel/skew
        double tolerance1 = 1.3; // skew/skew
        double tolerance2 = 1.5;; // outer parallel/skew
        double tolerance = 1.5; // CHECK tolerance

        TArrayI neighboring;

        int isector = tube->GetSectorID();
        int ilayer = tube->GetLayerID();
        //  cout << "LAYER/SECTOR " << ilayer << " " << isector << endl;

        // possible lay/sec to check:
        // same layer/same sector, before, after
        // layer up/same sector, before, after
        // layer down/same sector, before, after
        int possible_lay_sec[9][2] = { { ilayer, isector }, { ilayer, isector + 1 }, { ilayer, isector - 1 }, { ilayer + 1, isector },
                                                                        { ilayer + 1, isector + 1 }, { ilayer + 1, isector - 1 }, { ilayer - 1, isector },
                                                                        { ilayer - 1, isector + 1 }, { ilayer - 1, isector - 1 } };

        std::vector<int> neigh_candidates;
        for (int itest = 0; itest < 9; itest++) {

                        PndSttTube *tube2 = NULL;
                        int tubeid = -1;
                        //    cout << "is " << possible_lay_sec[itest][0] << " " << possible_lay_sec[itest][1] << "good?" << endl;
                        // check it is an existing layer/sec
                        if (possible_lay_sec[itest][0] < 0 || possible_lay_sec[itest][0] >= fNLayers)
                                        continue;
                        if (possible_lay_sec[itest][1] < 0 || possible_lay_sec[itest][1] >= fNSectors)
                                        continue;
                        // skip che pipe
                        if ((isector == 0 && possible_lay_sec[itest][1] == 5)
                                                        || (isector == 5 && possible_lay_sec[itest][1] == 0)
                                                        || (isector == 2 && possible_lay_sec[itest][1] == 3)
                                                        || (isector == 3 && possible_lay_sec[itest][1] == 2))
                                        continue;
                        // if not @ limit of the sector
                        //                      if (tube->IsSectorLimit() == kFALSE && (possible_lay_sec[itest][1] == isector + 1 || possible_lay_sec[itest][1] == isector - 1))
                        //                                      continue;

                        //    cout << "CHECKING " << possible_lay_sec[itest][0] <<  " " << possible_lay_sec[itest][1] << endl;
                        neigh_candidates = GetStrawRow(possible_lay_sec[itest][1], possible_lay_sec[itest][0]);

                        for (size_t itube = 0; itube < neigh_candidates.size(); itube++) {
                                        tubeid = neigh_candidates.at(itube);
                                        if (tube->GetTubeID() == tubeid)
                                                        continue;
                                        if (tubeid != -1)
                                                        tube2 = (PndSttTube*) fTubeArray->At(tubeid);
                                        double distance = 1000;
                                        if (tube2) distance = tube->GetDistance(tube2);
                                        // pick the correct tolerance
                                        if (tube2->IsSkew() == kTRUE && tube->IsSkew() == kTRUE)
                                                tolerance = tolerance1;
                                        else if (tube2->IsParallel() == kTRUE && tube->IsParallel() == kTRUE)
                                                tolerance = tolerance0;
                                        else {
                                                if (ilayer < fNLayers_inner_parallel + 3)
                                                        tolerance = tolerance1;
                                                else
                                                        tolerance = tolerance2;
                                        }
                                        if (distance < tolerance) {
                                                        int size = neighboring.GetSize();
                                                        neighboring.Set(size + 1);
                                                        neighboring.AddAt(tubeid, size);
                                                        //      cout << "ADD " << tubeid << " " << distance << endl;
                                        }
                        }
        }
        return neighboring;
}

TArrayI PndSttGeometryMap::FindNeighborings

(

int

tubeId

)

Definition at line 194 of file PndSttGeometryMap.cxx.

References FindNeighborings(), and fTubeArray.

                                                      {

        PndSttTube* tube = (PndSttTube*) fTubeArray->At(tubeId);
        return FindNeighborings(tube);
}

void PndSttGeometryMap::GenerateAngles ( )

private

Definition at line 332 of file PndSttGeometryMap.cxx.

References fSectorEnd, fSectorStart, and Pi.

Referenced by GenerateStrawMapAngleGeoType1().

{
  fSectorStart.push_back(1.57);
  fSectorEnd.push_back(2.62);
  fSectorStart.push_back(2.62);
  fSectorEnd.push_back(-2.62+2*TMath::Pi());
  fSectorStart.push_back(-2.62+2*TMath::Pi());
  fSectorEnd.push_back(-1.57+2*TMath::Pi());
  fSectorStart.push_back(-1.57+2*TMath::Pi());
  fSectorEnd.push_back(-0.52+2*TMath::Pi());
  fSectorStart.push_back(-0.52+2*TMath::Pi());
  fSectorEnd.push_back(0.52);
  fSectorStart.push_back(0.52);
  fSectorEnd.push_back(1.57);
}

void PndSttGeometryMap::GenerateStrawMap

(

Int_t

map

)

Definition at line 70 of file PndSttGeometryMap.cxx.

References fGeoType, and GenerateStrawMapGeoType1().

Referenced by PndSttGeometryMap().

                                                  {
  if(fGeoType == 1) GenerateStrawMapGeoType1(map);
}

void PndSttGeometryMap::GenerateStrawMapAngleGeoType1

(

)

Definition at line 348 of file PndSttGeometryMap.cxx.

References counter, fAxialStraw, fEndTube, fLayerOfStraw, fNLayers, fNSectors, fSectorEnd, fSectorOfStraw, fSectorStart, fShift, fShiftSkew, fStartTube, fStrawIndex, fStrawMapInitialized, fTubeArray, fVerbose, GenerateAngles(), PndSttTube::GetPosition(), GetStrawRow(), PndSttTube::GetWireDirection(), i, phi, Pi, row, PndSttTube::SetLayerID(), and PndSttTube::SetSectorID().

Referenced by GenerateStrawMapGeoType1().

{
  GenerateAngles();     

  vector<int> currentRow;
  int sector = 0;
  int lastsector = 0;
  int row = 0;
  fSectorOfStraw.push_back(-1);
  fLayerOfStraw.push_back(-1);
  fAxialStraw.push_back(false);
  if(fVerbose > 0)  cout << "Generating straw map for " << fTubeArray->GetEntriesFast() << " straws." << endl;
  for (int i = 1; i < fTubeArray->GetEntriesFast(); i++) {
    if(fVerbose > 0) cout << "Finding tube." << endl;
    PndSttTube* tube = (PndSttTube*) fTubeArray->At(i);
    if(fVerbose > 0) cout << "Tube address: " << tube << endl;
    bool isaxial = (tube->GetWireDirection().Theta() < 0.001);
    if(fVerbose > 0) {
      cout << "Axial Straw: " << isaxial << endl;
      cout << "Wire Direction: ";
      tube->GetWireDirection().Print();
    }
    double phi = tube->GetPosition().Phi();
    if(fVerbose > 0) cout << "Phi: " << phi << endl;
    if (phi < 0) phi += 2*TMath::Pi();
    if(fVerbose > 0) cout << "Checking sector." << endl;
    while ( !( (phi > fSectorStart.at(sector)) && (phi < fSectorEnd.at(sector)) )
            && !( (sector == 4) && ( (phi > fSectorStart.at(4)) || (phi < fSectorEnd.at(4)) ) ) ) {
      sector++;
      sector %= 6;
    }
    if (sector != lastsector) {
      fStrawIndex[lastsector].push_back(currentRow);
      currentRow.clear();
      if(fVerbose > 0) cout << "SECTOR COMPLETE: Row " << row << " added to sector " << lastsector << endl;
    }
    if (sector < lastsector) row++;
    lastsector = sector;
    currentRow.push_back(i);
    fSectorOfStraw.push_back(sector);
    fLayerOfStraw.push_back(row);
    fAxialStraw.push_back(isaxial);
    if(fVerbose > 0) cout << "Straw " << i << " added to " << sector << ", " << row << endl;

    tube->SetLayerID(row);
    tube->SetSectorID(sector); 
    

  }
  fStrawIndex[lastsector].push_back(currentRow);
  if(fVerbose > 0) cout << "STT COMPLETE. SECTOR COMPLETE: Row " << row << " added to sector " << lastsector << endl;

  fStartTube = (int**) malloc(sizeof(int*) * fNSectors);
  fEndTube   = (int**) malloc(sizeof(int*) * fNSectors);
  fShift     = (int**) malloc(sizeof(int*) * fNSectors);
  fShiftSkew = (int**) malloc(sizeof(int*) * fNSectors);
  for(int i = 0; i < fNSectors; i++) {
    fStartTube[i] = (int*) malloc(sizeof(int*) * fNLayers);
    fEndTube[i]   = (int*) malloc(sizeof(int*) * fNLayers);
    fShift[i]     = (int*) malloc(sizeof(int*) * fNLayers);
    fShiftSkew[i] = (int*) malloc(sizeof(int*) * fNLayers);
  }

  const int size = 26* 6; // fNLayers * fNSectors
  int shift_list[size] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 4,  7, 11, 16,
                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 5, 7, 10, 14, 19,
                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 5, 7, 10, 14, 19,
                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 4,  7, 11, 16,
                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 5, 7, 10, 14, 19,
                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 5, 7, 10, 14, 19};
  const int sizeskew = 8*6; //fNLayers_skew * fNSectors;
  int shift_list_skew[sizeskew] = {16, 19, 16, 16, 19, 16,
                                   17, 20, 17, 17, 20, 17,
                                   19, 21, 19, 19, 21, 19,
                                   20, 22, 20, 20, 22, 20,
                                   21, 24, 21, 21, 24, 21,
                                   22, 25, 22, 22, 25, 22,
                                   24, 26, 24, 24, 26, 24,
                                   25, 27, 25, 25, 27, 25};

  int counter = 0;
  for(int i = 0; i < fNSectors; i++) {
    for(int j = 0; j < fNLayers; j++) {
      fStartTube[i][j] = GetStrawRow(i,j).front();
      fEndTube[i][j]   = GetStrawRow(i,j).back();
      fShift[i][j]     = shift_list[counter];
      if (j<8 || j>15)
        fShiftSkew[i][j] = 0;
      else 
        fShiftSkew[i][j] = shift_list_skew[(j-8)*6+i];
      counter++;
    }
  }

  if(fVerbose > 0) {
    // =============== PRINT
    for(int irow = 0; irow < fNSectors; irow++) {
      for(int icol = 0; icol < fNLayers; icol++) {
        cout << " " << fStartTube[irow][icol];
      }
      cout << endl;
    }
    cout << endl;
    for(int irow = 0; irow < fNSectors; irow++) {
      for(int icol = 0; icol < fNLayers; icol++) {
        cout << " " << fEndTube[irow][icol];
      }
      cout << endl;
    }
    cout << endl;
    for(int irow = 0; irow < fNSectors; irow++) {
      for(int icol = 0; icol < fNLayers; icol++) {
        cout << " " << fShift[irow][icol];
      }
      cout << endl;
    }
    cout << endl;
  }

  fStrawMapInitialized = true;
}

void PndSttGeometryMap::GenerateStrawMapGeoType1

(

Int_t

map

)

Definition at line 325 of file PndSttGeometryMap.cxx.

References GenerateStrawMapAngleGeoType1(), and GenerateStrawMapTubeIDGeoType1().

Referenced by GenerateStrawMap().

                                                          {
  if(map == 0) GenerateStrawMapTubeIDGeoType1();
  else if(map == 1) GenerateStrawMapAngleGeoType1();
}

void PndSttGeometryMap::GenerateStrawMapTubeIDGeoType1

(

)

Definition at line 538 of file PndSttGeometryMap.cxx.

References fEndTube, fLayerOfStraw, fSectorOfStraw, fStartTube, fStrawIndex, fTubeArray, fVerbose, PndSttTube::GetWireDirection(), PndSttTube::SetLayerID(), and PndSttTube::SetSectorID().

Referenced by GenerateStrawMapGeoType1().

{
  vector<int> currentRow;
  int isector = 0, ilayer = 0;
  int lastsector = 0;
  fSectorOfStraw.push_back(-1);
  fLayerOfStraw.push_back(-1);

  if(fVerbose > 0) cout << "Generating straw map for " << fTubeArray->GetEntriesFast() << " straws." << endl;
  for (int itube = 1; itube < fTubeArray->GetEntriesFast(); itube++) {
    if(fVerbose > 0)  cout << "Finding tube." << endl;
    PndSttTube* tube = (PndSttTube*) fTubeArray->At(itube);
    bool isaxial = (tube->GetWireDirection().Theta() < 0.001);
    if(fVerbose > 0)  {
      cout << "Axial Straw: " << isaxial << endl;
      cout << "Wire Direction: ";
      tube->GetWireDirection().Print();
    }
   
    while(!(itube >= fStartTube[0][ilayer] && itube <= fEndTube[5][ilayer])) {
      ilayer++;
      isector = 0;
    }

    while(!(itube >= fStartTube[isector][ilayer] && itube <= fEndTube[isector][ilayer])) {
      isector++;
      isector %= 6;
    }
    
    if (isector != lastsector) {
      fStrawIndex[lastsector].push_back(currentRow);
      currentRow.clear();
      if(fVerbose > 0) cout << "SECTOR COMPLETE: Row " << ilayer << " added to sector " << lastsector << endl;
    }
    // if (isector < lastsector) ilayer++;
    lastsector = isector;
    currentRow.push_back(itube);
    fSectorOfStraw.push_back(isector);
    fLayerOfStraw.push_back(ilayer);

    tube->SetLayerID(ilayer); 
    tube->SetSectorID(isector); 
  
    if(fVerbose > 0) cout << "Straw " << itube << " added to " << isector << ", " << ilayer << endl;
  }
  fStrawIndex[lastsector].push_back(currentRow);
  if(fVerbose > 0) cout << "SECTOR COMPLETE: Row " << ilayer << " added to sector " << lastsector << endl;

}

double PndSttGeometryMap::GetAngleBetweenTubes	(	int	tubeID1,
		int	tubeID2
	)		const

Definition at line 200 of file PndSttGeometryMap.cxx.

References fTubeArray, PndSttTube::GetPosition(), phi, and Pi.

Referenced by InStraightLine().

                                                                             {

        PndSttTube *tube1, *tube2;
        tube1 = (PndSttTube*) fTubeArray->At(tubeID1);
        tube2 = (PndSttTube*) fTubeArray->At(tubeID2);

        TVector3 pos1, pos2, direction;
        pos1 = tube1->GetPosition();
        pos2 = tube2->GetPosition();

        direction = pos2 - pos1;

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

        return phi;
}

TArrayI PndSttGeometryMap::GetNeighboringsByMap

(

int

tubeId

)

Definition at line 244 of file PndSttGeometryMap.cxx.

References FillStrawNeighborsMap(), and fStrawNeighbors.

Referenced by PndSttSkewedCombineTask::Exec(), PndSttCellTrackFinderData::FindHitNeighborsEventBased(), PndSttCellTrackFinderData::FindHitNeighborsTimeBased(), PndSttCellTrackFinderTask::Init(), and PndSttCellTrackFinderData::PndSttCellTrackFinderData().

{
        if (fStrawNeighbors.size() == 0)
                FillStrawNeighborsMap();
        return fStrawNeighbors[tubeId];
}

int PndSttGeometryMap::GetRow ( int  strawindex ) const

inline

Definition at line 47 of file PndSttGeometryMap.h.

References fLayerOfStraw.

Referenced by IsEdgeStraw(), and IsSectorBorderStraw().

{return fLayerOfStraw.at(strawindex);}

int PndSttGeometryMap::GetSector ( int  strawindex ) const

inline

Definition at line 48 of file PndSttGeometryMap.h.

References fSectorOfStraw.

Referenced by IsEdgeStraw(), and IsSectorBorderStraw().

{return fSectorOfStraw.at(strawindex);}

const vector<int>& PndSttGeometryMap::GetStrawRow ( int  sector, int  row  ) const

inline

Definition at line 49 of file PndSttGeometryMap.h.

References fStrawIndex.

Referenced by FindNeighborings(), GenerateStrawMapAngleGeoType1(), IsAxialRow(), IsEdgeStraw(), IsSectorBorderStraw(), and IsSkewedRow().

{return (fStrawIndex.find(sector))->second.at(row);}

const vector<vector<int> >& PndSttGeometryMap::GetStrawSector ( int  sector ) const

inline

Definition at line 50 of file PndSttGeometryMap.h.

References fStrawIndex.

{return (fStrawIndex.find(sector))->second;}

bool PndSttGeometryMap::InStraightLine	(	int	tube1,
		int	tube2,
		int	tube3
	)		const

Definition at line 219 of file PndSttGeometryMap.cxx.

References GetAngleBetweenTubes(), and Pi.

                                                                            {

        double diff = abs(
                        GetAngleBetweenTubes(tube1, tube2)
                                        - GetAngleBetweenTubes(tube1, tube3));

        //accept a deviation of 1 percent
        //if ((TMath::Pi() * 0.99) < diff & diff < (TMath::Pi() * 1.01)) {
        if (((TMath::Pi() * 0.99) < diff) && (diff < (TMath::Pi() * 1.01))) {  //[R.K. 01/2017] more explicit parenthesis logic
                //tube1 is in the middle
                return true;
        } else if (diff < 0.01) {
                //tube1 is at the edge of the line
                return true;
        } else
                return false;
}

bool PndSttGeometryMap::IsAxialRow

(

int

rowindex

)

const

Definition at line 521 of file PndSttGeometryMap.cxx.

References GetStrawRow(), and IsAxialStraw().

{
  //only works if all straws in one row are of the same type
  //which is the case for the current geometry
  return IsAxialStraw(GetStrawRow(0, rowindex).at(0));
}

bool PndSttGeometryMap::IsAxialStraw ( int  strawindex ) const

inline

Definition at line 55 of file PndSttGeometryMap.h.

References fAxialStraw.

Referenced by IsAxialRow().

{return fAxialStraw.at(strawindex);}

bool PndSttGeometryMap::IsEdgeStraw

(

int

strawindex

)

const

Definition at line 470 of file PndSttGeometryMap.cxx.

References GetRow(), GetSector(), and GetStrawRow().

Referenced by FillGeometryParametersGeoType1().

{
  if (GetRow(strawindex) == 0) return true;
  int endstraw = 0; //there is no straw index 0
  switch (GetSector(strawindex)) {
  case 0: endstraw = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).front();
    break;
  case 2: endstraw = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).back();
    break;
  case 3: endstraw = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).front();
    break;
  case 5: endstraw = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).back();
    break;
  default: return false;
  };
  return (strawindex == endstraw);
}

int PndSttGeometryMap::IsSectorBorderStraw

(

int

strawindex

)

const

Definition at line 488 of file PndSttGeometryMap.cxx.

References fShiftSkew, fTubeArray, PndSttTube::GetHalfLength(), GetRow(), GetSector(), GetStrawRow(), and PndSttTube::IsParallel().

Referenced by FillGeometryParametersGeoType1().

{
  PndSttTube* tube = (PndSttTube*) fTubeArray->At(strawindex);
  if (tube->IsParallel()) {
    int endstrawcw  = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).front();
    int endstrawccw = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).back();
    if (strawindex == endstrawcw) return -1;
    if (strawindex == endstrawccw) return 1;
    return 0;
  }
  else {
    int endstrawcw  = GetStrawRow(GetSector(strawindex), GetRow(strawindex)).front();
    int endstrawscw = endstrawcw + fShiftSkew[GetSector(strawindex)][GetRow(strawindex)];
    int endstrawccw = endstrawscw -1;
    if (tube->GetHalfLength()==75) {
      if (strawindex == endstrawcw) return -1;
      if (strawindex == endstrawccw) return 1;
      return 0;
    }
    else {
      if (strawindex==endstrawccw+1 || strawindex==endstrawccw+3 || strawindex==endstrawccw+5 ||
          strawindex==endstrawccw+7 || strawindex==endstrawccw+9 || strawindex==endstrawccw+11 ||
          strawindex==endstrawccw+13 || strawindex==endstrawccw+15)
        return 1;
      else if (strawindex==endstrawscw +1 || strawindex==endstrawscw+3 || strawindex==endstrawscw+5 ||
          strawindex==endstrawscw+7 || strawindex==endstrawscw+9 || strawindex==endstrawscw+11 ||
          strawindex==endstrawscw+13 || strawindex==endstrawscw+15)
        return -1;
    }
  }
  return -1000;
}

bool PndSttGeometryMap::IsSkewedRow

(

int

rowindex

)

const

Definition at line 528 of file PndSttGeometryMap.cxx.

References GetStrawRow(), and IsSkewedStraw().

{
  //only works if all straws in one row are of the same type
  //which is the case for the current geometry
  return IsSkewedStraw(GetStrawRow(0, rowindex).at(0));
}

bool PndSttGeometryMap::IsSkewedStraw ( int  strawindex ) const

inline

Definition at line 56 of file PndSttGeometryMap.h.

References fAxialStraw.

Referenced by IsSkewedRow(), and PndSttCellTrackFinderData::PndSttCellTrackFinderData().

{return !(fAxialStraw.at(strawindex));}

PndSttGeometryMap& PndSttGeometryMap::operator= ( const PndSttGeometryMap &  )

delete

void PndSttGeometryMap::SetGeneralParameters

(

)

Definition at line 66 of file PndSttGeometryMap.cxx.

References fGeoType, and SetGeneralParametersGeoType1().

Referenced by PndSttGeometryMap().

                                             { 
  if(fGeoType == 1) SetGeneralParametersGeoType1();
}

void PndSttGeometryMap::SetGeneralParametersGeoType1

(

)

Definition at line 82 of file PndSttGeometryMap.cxx.

References fNLayers, fNLayers_fillup_parallel, fNLayers_inner_parallel, fNLayers_outer_parallel, fNLayers_skew, fNSectors, fNTubes, fNTubes_fillup_parallel, fNTubes_inner_parallel, fNTubes_outer_parallel, and fNTubes_skewed.

Referenced by SetGeneralParameters().

                                                     {   //  CHECK whether it depends on geometry or not
 
  //   fNLayers = 18;
  //   fNSectors = 6;
  //   const int size = 18* 6; // CHECK
 
  fNLayers = 26;
  fNSectors = 6;
  //const int size = 26* 6; // CHECK //[R.K. 01/2017] unused variable?
 
 
  fNTubes_inner_parallel = 1000;
  fNTubes_outer_parallel = 884;
  fNTubes_fillup_parallel = 914;
  fNTubes_skewed = 1744;
  fNTubes = fNTubes_inner_parallel + fNTubes_outer_parallel + fNTubes_fillup_parallel + fNTubes_skewed;

  fNLayers_inner_parallel = 8;
  fNLayers_skew = 8;
  fNLayers_outer_parallel = 4;
  fNLayers_fillup_parallel = 6;

}

void PndSttGeometryMap::SetVerbose ( int  ver )

inline

Definition at line 76 of file PndSttGeometryMap.h.

References fVerbose.

{ fVerbose = ver; }

Member Data Documentation

vector<bool> PndSttGeometryMap::fAxialStraw

private

Definition at line 88 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), IsAxialStraw(), and IsSkewedStraw().

Int_t ** PndSttGeometryMap::fEndTube

private

Definition at line 96 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), and GenerateStrawMapTubeIDGeoType1().

Int_t PndSttGeometryMap::fGeoType

private

Definition at line 80 of file PndSttGeometryMap.h.

Referenced by FillGeometryParameters(), GenerateStrawMap(), PndSttGeometryMap(), and SetGeneralParameters().

vector<int> PndSttGeometryMap::fLayerOfStraw

private

Definition at line 87 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), GenerateStrawMapTubeIDGeoType1(), and GetRow().

Int_t PndSttGeometryMap::fNLayers

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by FindNeighborings(), GenerateStrawMapAngleGeoType1(), and SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNLayers_fillup_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNLayers_inner_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by FindNeighborings(), and SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNLayers_outer_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNLayers_skew

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNSectors

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by FindNeighborings(), GenerateStrawMapAngleGeoType1(), and SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNTubes

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by FillStrawNeighborsMap(), and SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNTubes_fillup_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNTubes_inner_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNTubes_outer_parallel

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

Int_t PndSttGeometryMap::fNTubes_skewed

private

Definition at line 95 of file PndSttGeometryMap.h.

Referenced by SetGeneralParametersGeoType1().

vector<double> PndSttGeometryMap::fSectorEnd

private

Definition at line 90 of file PndSttGeometryMap.h.

Referenced by GenerateAngles(), and GenerateStrawMapAngleGeoType1().

vector<int> PndSttGeometryMap::fSectorOfStraw

private

Definition at line 86 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), GenerateStrawMapTubeIDGeoType1(), and GetSector().

vector<double> PndSttGeometryMap::fSectorStart

private

Definition at line 89 of file PndSttGeometryMap.h.

Referenced by GenerateAngles(), and GenerateStrawMapAngleGeoType1().

Int_t ** PndSttGeometryMap::fShift

private

Definition at line 96 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1().

Int_t ** PndSttGeometryMap::fShiftSkew

private

Definition at line 96 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), and IsSectorBorderStraw().

Int_t** PndSttGeometryMap::fStartTube

private

Definition at line 96 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), and GenerateStrawMapTubeIDGeoType1().

map<int, vector< vector<int> > > PndSttGeometryMap::fStrawIndex

private

Definition at line 84 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1(), GenerateStrawMapTubeIDGeoType1(), GetStrawRow(), and GetStrawSector().

bool PndSttGeometryMap::fStrawMapInitialized

private

Definition at line 91 of file PndSttGeometryMap.h.

Referenced by GenerateStrawMapAngleGeoType1().

map<int, TArrayI> PndSttGeometryMap::fStrawNeighbors

private

Definition at line 85 of file PndSttGeometryMap.h.

Referenced by FillStrawNeighborsMap(), and GetNeighboringsByMap().

TClonesArray* PndSttGeometryMap::fTubeArray

private

Definition at line 99 of file PndSttGeometryMap.h.

Referenced by CalculateStrawPoca(), FillGeometryParametersGeoType1(), FindNeighborings(), GenerateStrawMapAngleGeoType1(), GenerateStrawMapTubeIDGeoType1(), GetAngleBetweenTubes(), and IsSectorBorderStraw().

Int_t PndSttGeometryMap::fVerbose

private

Definition at line 80 of file PndSttGeometryMap.h.

Referenced by CalculateStrawPoca(), GenerateStrawMapAngleGeoType1(), GenerateStrawMapTubeIDGeoType1(), and SetVerbose().

---

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

• PndSttGeometryMap.h
• PndSttGeometryMap.cxx