geo::CRU Namespace Reference

Functions
static bool	sortAuxDet (const AuxDetGeo &ad1, const AuxDetGeo &ad2)
static bool	sortAuxDetSensitive (const AuxDetSensitiveGeo &ad1, const AuxDetSensitiveGeo &ad2)
static bool	sortCryo (const CryostatGeo &c1, const CryostatGeo &c2)
static bool	sortTPC (const TPCGeo &t1, const TPCGeo &t2)
static bool	sortPlane (const PlaneGeo &p1, const PlaneGeo &p2)
bool	sortWire (WireGeo const &w1, WireGeo const &w2)

Function Documentation

static bool geo::CRU::sortAuxDet ( const AuxDetGeo &  ad1, const AuxDetGeo &  ad2  )

static

Definition at line 32 of file GeoObjectSorterCRU.cxx.

    {

      // sort based off of GDML name, assuming ordering is encoded
      std::string ad1name = (ad1.TotalVolume())->GetName();
      std::string ad2name = (ad2.TotalVolume())->GetName();

      // assume volume name is "volAuxDet##"
      int ad1Num = atoi( ad1name.substr( 9, ad1name.size()).c_str() );
      int ad2Num = atoi( ad2name.substr( 9, ad2name.size()).c_str() );

      return ad1Num < ad2Num;

    }

static bool geo::CRU::sortAuxDetSensitive ( const AuxDetSensitiveGeo &  ad1, const AuxDetSensitiveGeo &  ad2  )

static

Definition at line 49 of file GeoObjectSorterCRU.cxx.

    {

      // sort based off of GDML name, assuming ordering is encoded
      std::string ad1name = (ad1.TotalVolume())->GetName();
      std::string ad2name = (ad2.TotalVolume())->GetName();

      // assume volume name is "volAuxDetSensitive##"
      int ad1Num = atoi( ad1name.substr( 9, ad1name.size()).c_str() );
      int ad2Num = atoi( ad2name.substr( 9, ad2name.size()).c_str() );

      return ad1Num < ad2Num;

    }

static bool geo::CRU::sortCryo ( const CryostatGeo &  c1, const CryostatGeo &  c2  )

static

Definition at line 66 of file GeoObjectSorterCRU.cxx.

    {
      double xyz1[3] = {0.}, xyz2[3] = {0.};
      double local[3] = {0.};
      c1.LocalToWorld(local, xyz1);
      c2.LocalToWorld(local, xyz2);

      return xyz1[0] < xyz2[0];
    }

static bool geo::CRU::sortPlane ( const PlaneGeo &  p1, const PlaneGeo &  p2  )

static

Definition at line 99 of file GeoObjectSorterCRU.cxx.

    {
      double xyz1[3] = {0.};
      double xyz2[3] = {0.};
      double local[3] = {0.};
      p1.LocalToWorld(local, xyz1);
      p2.LocalToWorld(local, xyz2);

      // drift direction is negative, plane number increases in drift direction
      return xyz1[0] > xyz2[0];
    }

static bool geo::CRU::sortTPC ( const TPCGeo &  t1, const TPCGeo &  t2  )

static

Definition at line 79 of file GeoObjectSorterCRU.cxx.

    {
      double xyz1[3] = {0.};
      double xyz2[3] = {0.};
      double local[3] = {0.};
      t1.LocalToWorld(local, xyz1);
      t2.LocalToWorld(local, xyz2);

      // First sort all TPCs into same-z groups
      if(xyz1[2]<xyz2[2]) return true;

      // Within a same-z group, sort TPCs into same-y groups
      if(xyz1[2] == xyz2[2] && xyz1[1] < xyz2[1]) return true;

      return false;
    }

bool geo::CRU::sortWire	(	WireGeo const &	w1,
		WireGeo const &	w2
	)

Definition at line 113 of file GeoObjectSorterCRU.cxx.

                                                       {
      // wire sorting algorithm
      // z_low -> z_high
      // for same-z group 
      //  the sorting either from bottom (y) left (z) corner up for strip angles > pi/2
      //  and horizontal wires
      //  or from top corner to bottom otherwise
      //  we assume all wires in the plane are parallel

      //  w1 geo info
      auto center1 = w1.GetCenter<geo::Point_t>();
      auto start1  = w1.GetStart<geo::Point_t>();
      auto end1    = w1.GetEnd<geo::Point_t>();
      auto Delta   = end1-start1;
      //double dx1   = Delta.X();
      double dy1   = Delta.Y();
      double dz1   = Delta.Z();
      
      // w2 geo info
      auto center2 = w2.GetCenter<geo::Point_t>();
      
      auto CheckTol = [](double val, double tol = 1.E-4){ 
        return (std::abs( val ) < tol); 
      };
      
      if( CheckTol(dz1) ){ // wires perpendicular to z axis
        return (center1.Z() < center2.Z());
      }
      
      if( CheckTol(dy1) ){ // wires perpendicular to y axis
        return (center1.Y() < center2.Y());
      }
      
      // wires at angle with same z center
      if( CheckTol( center2.Z() - center1.Z() ) ){
        if( dz1 < 0 ){ dy1 = -dy1; } // always compare here upstream - downstream ends
        if( dy1 < 0 ){ return (center1.Y() < center2.Y()); }
        if( dy1 > 0 ){ return (center1.Y() > center2.Y()); }
      }

      // otherwise sorted in z
      return ( center1.Z() < center2.Z() );
    }