RegSelModule.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
 
/**************************************************************************
 **
 **   File:         RegSelModule.cxx        
 **
 **   Description:    
 **                   
 **                   
 ** 
 **   Author:       M.Sutton    
 **
 **   Created:      Tue Apr  3 11:16:12 BST 2007
 **   Modified:     
 **                   
 **                   
 **
 **************************************************************************/ 
 
 
 
#include "RegSelLUT/RegSelModule.h"
 
#include <iostream>
#include <cmath>
 
 
RegSelModule::RegSelModule() :
  ZRObject( 0, 0, 0, 0 ),
  m_phiMin(0), m_phiMax(0),
  m_layer(0), m_detector(0),
  m_robID(0), m_hashID(0) { }
 
 
 
RegSelModule::RegSelModule(double zmin,   double zmax, 
               double rmin,   double rmax, 
               double phimin, double phimax, 
               int layer, 
               int detector,
               uint32_t       robid, 
               IdentifierHash hashid) : 
  //  ZRObject( std::sqrt(rmin), std::sqrt(rmax), // zmin, zmax ),
  //        ( zmin<0 ? -1 : 1 )*std::sqrt(std::fabs(zmin)), 
  //        ( zmax<0 ? -1 : 1 )*std::sqrt(std::fabs(zmax)) ),
  //  ZRObject( std::log10(rmin)-1.6, std::log10(rmax)-1.6, // zmin, zmax ),
  //        ( zmin<0 ? -1 : 1 )*std::log10(std::fabs(zmin)), 
  //        ( zmax<0 ? -1 : 1 )*std::log10(std::fabs(zmax)) ),
  ZRObject( rmin, rmax, zmin, zmax ),
  m_phiMin(phimin), m_phiMax(phimax),
  m_layer(layer), m_detector(detector),
  m_robID(robid), m_hashID(hashid)
{
  while ( m_phiMax> M_PI ) m_phiMax-=2*M_PI; 
  while ( m_phiMax<-M_PI ) m_phiMax+=2*M_PI; 
  while ( m_phiMin> M_PI ) m_phiMin-=2*M_PI; 
  while ( m_phiMin<-M_PI ) m_phiMin+=2*M_PI; 
  //  std::cout << "RegSelModule() " << *this << std::endl;
} 
 
 
 
 
RegSelModule::RegSelModule(double zmin,   double zmax, 
               double z2min,  double z2max, 
               double rmin,   double rmax, 
               double phimin, double phimax, 
               int layer, 
               int detector,
               uint32_t       robid, 
               IdentifierHash hashid) : 
  //  ZRObject( std::sqrt(rmin), std::sqrt(rmax), // zmin, zmax ),
  //        ( zmin<0 ? -1 : 1 )*std::sqrt(std::fabs(zmin)), 
  //        ( zmax<0 ? -1 : 1 )*std::sqrt(std::fabs(zmax)) ),
  //  ZRObject( std::log10(rmin)-1.6, std::log10(rmax)-1.6, // zmin, zmax ),
  //        ( zmin<0 ? -1 : 1 )*std::log10(std::fabs(zmin)), 
  //        ( zmax<0 ? -1 : 1 )*std::log10(std::fabs(zmax)) ),
  ZRObject( rmin, rmax, zmin, zmax, z2min, z2max ),
  m_phiMin(phimin), m_phiMax(phimax),
  m_layer(layer), m_detector(detector),
  m_robID(robid), m_hashID(hashid)
{
  while ( m_phiMax> M_PI ) m_phiMax-=2*M_PI; 
  while ( m_phiMax<-M_PI ) m_phiMax+=2*M_PI; 
  while ( m_phiMin> M_PI ) m_phiMin-=2*M_PI; 
  while ( m_phiMin<-M_PI ) m_phiMin+=2*M_PI; 
  //  std::cout << "RegSelModule() " << *this << std::endl;
} 
 
 
 
std::ostream& operator<<(std::ostream& s, const RegSelModule& m)
{
  s << "[ lyr= "     << m.layer() 
    << " ,\tdet= "   << m.detector() 
    << " ,\tr= "     << m.rMin()            << " - " << m.rMax()
    << " ,\tphi= "   << m.phiMin()*180/M_PI << " - " << m.phiMax()*180/M_PI
    << " ,\tz= "     << m.zMin()            << " - " << m.zMax();
  if ( m.zMin()!=m.z2Min() || m.zMax()!=m.z2Max() )  s << " ,\tz2= " << m.z2Min() << " - " << m.z2Max();
  s << " ,\trob= 0x"  << std::hex << m.robID() 
    << " ,\thash= 0x" << std::hex << m.hashID() << std::dec << ( m.enabled() ? " ]" : " (disabled)]");
  return s;
}  
 
 
 
bool getModule(std::istream& ss, RegSelModule& m)
{
  char s[128], s1[128], s2[128]; 
 
  int layer;
  int detector;
 
  double rMin,   rMax;
  double zMin,   zMax;
  double z2Min,  z2Max;
  double phiMin, phiMax;
 
  uint32_t robid;
  //  IdentifierHash hash;
  unsigned int hashint;
 
  ss >> s >> s >> std::dec >> layer 
    >> s >> s >> std::dec >> detector 
    >> s >> s >> rMin    >> s >> rMax  
    >> s >> s >> phiMin  >> s >> phiMax  
    >> s >> s >> zMin    >> s >> zMax;
 
  z2Min = zMin;
  z2Max = zMax;
  
  ss >> s1 >> s2;

  if ( std::string(s2)=="z2=" ) { 
    ss >> z2Min    >> s >> z2Max;
    ss >> s >> s;
  }
 
  ss >> std::hex >> robid 
    >> s >> s >> std::hex >> hashint >> std::dec >> s;
 
  if ( ss.fail() ) return false;
  
  //  std::cout << "s.fail() " << s.fail() << std::endl; 
 
  phiMin *= M_PI/180;   
  phiMax *= M_PI/180;   
 
  RegSelModule tm( zMin,   zMax, 
           z2Min,  z2Max, 
           rMin,   rMax, 
           phiMin, phiMax, 
           layer, 
           detector, 
           robid, 
           IdentifierHash(hashint));
  
  m=tm;
 
  return true;
}