LASGlobalData.h

Go to the documentation of this file.

#ifndef __LASGLOBALDATA_H
#define __LASGLOBALDATA_H

#include<vector>
#include<iostream>
#include "TNamed.h"


template <class T> 
class LASGlobalData : public TNamed
{

 public:
  // enums not in use so far...
  enum Subdetector { TECPLUS, TECMINUS, TIB, TOB };
  enum TecRing { RING4, RING6 };
  enum Beam { BEAM0, BEAM1, BEAM2, BEAM3, BEAM4, BEAM5, BEAM6, BEAM7 };
  enum TecDisk { DISK1, DISK2, DISK3, DISK4, DISK5, DISK6, DISK7, DISK8, DISK9 };
  enum TibTobPosition { MINUS3, MINUS2, MINUS1, PLUS1, PLUS2, PLUS3 };
  LASGlobalData();
  LASGlobalData(const T&);
  T& GetTECEntry( int subdetector, int tecRing, int beam, int tecDisk );
  T& GetTIBTOBEntry( int subdetector, int beam, int tibTobPosition );
  T& GetTEC2TECEntry( int subdetector, int beam, int tecDisk );
  void SetTECEntry( int subdetector, int tecRing, int beam, int tecDisk, T );
  void SetTIBTOBEntry( int subdetector, int beam, int tibTobPosition, T );
  void SetTEC2TECEntry( int subdetector, int beam, int tecDisk, T );
  //  LASGlobalData<T>& operator=( LASGlobalData<T>& );

 private:
  //void Init( void );
  void Init( const T& in=T());
  std::vector<std::vector<std::vector<T> > > tecPlusData; // ring<beam<disk<T>>>
  std::vector<std::vector<std::vector<T> > > tecMinusData; // ring<beam<disk<T>>>
  std::vector<std::vector<T> > tecPlusATData; // beam<disk<T>>
  std::vector<std::vector<T> > tecMinusATData; // beam<disk<T>>
  std::vector<std::vector<T> > tibData; // beam<pos<T>>
  std::vector<std::vector<T> > tobData; // beam<pos<T>>

  ClassDef( LASGlobalData, 2 );
};

// since this is a template
//#include "Alignment/LaserAlignment/src/LASGlobalData.cc"


template <class T>
LASGlobalData<T>::LASGlobalData() {

  Init();

}

template <class T>
LASGlobalData<T>::LASGlobalData(const T& in)
{
  Init(in);
}



template <class T>
T& LASGlobalData<T>::GetTECEntry( int theDetector, int theRing, int theBeam, int theDisk ) {
  
  // do a range check first
  if( !( ( theDetector == 0 || theDetector == 1 ) &&        // TEC+ or TEC-
     ( theRing == 0 || theRing == 1 )         &&        // ring4 or ring6
     ( theBeam >= 0 && theBeam < 8 )          &&        // eight beams in a TEC
     ( theDisk >= 0 && theDisk < 9 )             ) ) {  // disk1..disk9
    std::cerr << " [LASGlobalData::GetTECEntry] ** ERROR: illegal input coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", ring " << theRing << ", beam " << theBeam << ", disk " << theDisk << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 0 ) return( tecPlusData.at( theRing ).at( theBeam ).at( theDisk ) );
    else return( tecMinusData.at( theRing ).at( theBeam ).at( theDisk ) );
  }

}




template <class T>
T& LASGlobalData<T>::GetTIBTOBEntry( int theDetector, int theBeam, int thePosition ) {

  // do a range check first
  if( !( ( theDetector == 2 || theDetector == 3 ) &&        // TIB or TOB
     ( theBeam >= 0 && theBeam < 8 )          &&        // there are eight AT beams
     ( thePosition >= 0 && thePosition < 6 )     ) ) {  // z-pos -3 .. z-pos +3
    std::cerr << " [LASGlobalData::GetTIBTOBEntry] ** ERROR: illegal coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", beam " << theBeam << ", position " << thePosition << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 2 ) return( tibData.at( theBeam ).at( thePosition ) );
    else return( tobData.at( theBeam ).at( thePosition ) );
  }

}




template <class T>
T& LASGlobalData<T>::GetTEC2TECEntry( int theDetector, int theBeam, int theDisk ) {

  // do a range check first
  if( !( ( theDetector == 0 || theDetector == 1 ) &&        // TEC+ or TEC-
     ( theBeam >= 0 && theBeam < 8 )          &&        // eight AT beams in a TEC
     ( theDisk >= 0 && theDisk < 6 )     ) ) {          // disk1...disk5 are hit by AT
    std::cerr << " [LASGlobalData::GetTEC2TECEntry] ** ERROR: illegal coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", beam " << theBeam << ", disk " << theDisk << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 0 ) return( tecPlusATData.at( theBeam ).at( theDisk ) );
    else return( tecMinusATData.at( theBeam ).at( theDisk ) );
  }

}





template <class T>
void LASGlobalData<T>::SetTECEntry( int theDetector, int theRing, int theBeam, int theDisk, T theEntry ) {
  
  // do a range check first
  if( !( ( theDetector == 0 || theDetector == 1 ) &&        // TEC+ or TEC-
     ( theRing == 0 || theRing == 1 )         &&        // ring4 or ring6
     ( theBeam >= 0 && theBeam < 8 )          &&        // eight beams in a TEC
     ( theDisk >= 0 && theDisk < 9 )             ) ) {  // disk1..disk9
    std::cerr << " [LASGlobalData::SetTECEntry] ** ERROR: illegal coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", ring " << theRing << ", beam " << theBeam << ", disk " << theDisk << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 0 ) tecPlusData.at( theRing ).at( theBeam ).at( theDisk ) = theEntry;
    else tecMinusData.at( theRing ).at( theBeam ).at( theDisk ) = theEntry;
  }

}





template <class T>
void LASGlobalData<T>::SetTIBTOBEntry( int theDetector, int theBeam, int thePosition, T theEntry ) {

  // do a range check first
  if( !( ( theDetector == 2 || theDetector == 3 ) &&        // TIB or TOB
     ( theBeam >= 0 && theBeam < 8 )          &&        // there are eight AT beams
     ( thePosition >= 0 && thePosition < 6 )     ) ) {  // pos-3..pos+3
    std::cerr << " [LASGlobalData::SetTIBTOBEntry] ** ERROR: illegal coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", beam " << theBeam << ", position " << thePosition << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 2 ) tibData.at( theBeam ).at( thePosition ) = theEntry;
    else tobData.at( theBeam ).at( thePosition ) = theEntry;
  }

}





template <class T>
void LASGlobalData<T>::SetTEC2TECEntry( int theDetector, int theBeam, int theDisk, T theEntry ) {

  // do a range check first
  if( !( ( theDetector == 0 || theDetector == 1 ) &&        // TEC+ or TEC-
     ( theBeam >= 0 && theBeam < 8 )          &&        // eight beams in a TEC
     ( theDisk >= 0 && theDisk < 6 )             ) ) {  // disk1..disk5 for TEC AT
    std::cerr << " [LASGlobalData::SetTEC2TECEntry] ** ERROR: illegal coordinates:" << std::endl;
    std::cerr << "   detector " << theDetector << ", beam " << theBeam << ", disk " << theDisk << "." << std::endl;
    throw   "   Bailing out."; // @@@ REPLACE THIS BY cms::Exception (<FWCore/Utilities/interface/Exception.h> in 1_3_6)
  }
  else {
    if( theDetector == 0 ) tecPlusATData.at( theBeam ).at( theDisk ) = theEntry;
    else tecMinusATData.at( theBeam ).at( theDisk ) = theEntry;
  }

}





// ///
// /// element wise assignment operator
// ///
// template <class T>
// LASGlobalData<T>& LASGlobalData<T>::operator=( LASGlobalData<T>& anotherGlobalData ) {

//   // TEC copy
//   for( int det = 0; det < 2; ++det ) {
//     for( int ring = 0; ring < 2; ++ring ) {
//       for( int beam = 0; beam < 8; ++ beam ) {
//  for( int disk = 0; disk < 9; ++ disk ) {
//    this->SetTECEntry( det, ring, beam, disk, anotherGlobalData.GetTECEntry( det, ring, beam, disk ) );
//  }
//       }
//     }
//   }

//   // TIBTOB copy
//   for( int det = 2; det < 4; ++det ) {
//     for( int beam = 0; beam < 8; ++ beam ) {
//       for( int pos = 0; pos < 6; ++ pos ) {
//  this->SetTIBTOBEntry( det, beam, pos, anotherGlobalData.GetTIBTOBEntry( det, beam, pos ) );
//       }
//     }
//   }
  
//   // TEC2TEC copy
//   for( int det = 2; det < 4; ++det ) {
//     for( int beam = 0; beam < 8; ++ beam ) {
//       for( int disk = 0; disk < 9; ++ disk ) {
//  this->SetTEC2TECEntry( det, beam, disk, anotherGlobalData.GetTEC2TECEntry( det, beam, disk ) );
//       }
//     }
//   }

// }



template <class T>
void LASGlobalData<T>::Init( const T& in ) {

  // create TEC+ subdetector "multi"-vector of T
  tecPlusData.resize( 2 ); // create ring4 and ring6
  for( unsigned int ring = 0; ring < tecPlusData.size(); ++ring ) {
    tecPlusData.at( ring ).resize( 8 ); // create 8 beams for each ring
    for( unsigned int beam = 0; beam < tecPlusData.at( ring ).size(); ++beam ) {
      tecPlusData.at( ring ).at( beam ).resize( 9 , in); // create 9 disks for each beam
    }
  }

  // same for TEC-
  tecMinusData.resize( 2 ); // create ring4 and ring6
  for( unsigned int ring = 0; ring < tecMinusData.size(); ++ring ) {
    tecMinusData.at( ring ).resize( 8 ); // create 8 beams for each ring
    for( unsigned int beam = 0; beam < tecMinusData.at( ring ).size(); ++beam ) {
      tecMinusData.at( ring ).at( beam ).resize( 9, in ); // create 9 disks for each beam
    }
  }
  
  // same for TEC+ AT
  tecPlusATData.resize( 8 ); // create 8 beams
  for( unsigned int beam = 0; beam < tecPlusATData.size(); ++beam ) {
    tecPlusATData.at( beam ).resize( 5, in ); // five TEC disks hit by each AT beam
  }

  // same for TEC- AT
  tecMinusATData.resize( 8 ); // create 8 beams
  for( unsigned int beam = 0; beam < tecMinusATData.size(); ++beam ) {
    tecMinusATData.at( beam ).resize( 5, in ); // five TEC disks hit by each AT beam
  }

  // same for TIB..
  tibData.resize( 8 ); // create 8 beams
  for( unsigned int beam = 0; beam < tibData.size(); ++ beam ) {
    tibData.at( beam ).resize( 6, in ); // six TIB modules hit by each beam
  }

  // ..and for TOB
  tobData.resize( 8 ); // create 8 beams
  for( unsigned int beam = 0; beam < tobData.size(); ++ beam ) {
    tobData.at( beam ).resize( 6, in ); // six TOB modules hit by each beam
  }
}


#endif