/data/refman/pasoursint/CMSSW_6_1_2_SLHC4_patch1/src/Alignment/LaserAlignment/src/LASConstants.cc

Go to the documentation of this file.

#include "Alignment/LaserAlignment/interface/LASConstants.h"
#include "FWCore/Utilities/interface/Exception.h"

LASConstants::LASConstants() :
  atRadius(0.), tecBsZPosition(0.), atZPosition(0.)
{
}





LASConstants::LASConstants( std::vector<edm::ParameterSet> const& theConstConf ) {

  InitContainers();

  for( std::vector<edm::ParameterSet>::const_iterator iter = theConstConf.begin(); iter < theConstConf.end(); ++iter ) {

    if( iter->getParameter<std::string>( "PSetName" ) == "BeamsplitterKinks" ) FillBsKinks( *iter );
    else if( iter->getParameter<std::string>( "PSetName" ) == "Radii" ) FillRadii( *iter );
    else if( iter->getParameter<std::string>( "PSetName" ) == "ZPositions" ) FillZPositions( *iter );
    else {
      std::cerr << " [] ** WARNING: Cannot process unknown parameter set named: " << iter->getParameter<std::string>( "PSetName" ) << "." << std::endl;
    }

  }

}





LASConstants::~LASConstants() {
}





double LASConstants::GetEndcapBsKink( unsigned int det, unsigned int ring , unsigned int beam ) const {
  
  if( ! ( ( det == 0 || det == 1 ) && ( ring == 0 || ring == 1 ) && ( beam < 8U ) ) ) { // beam >= 0, since beam is unsigned
    throw cms::Exception( " [LASConstants::GetEndcapBsKink]" ) << " ** ERROR: no such element: det " << det << ", ring " << ring << ", beam " << beam << "." << std::endl;
  }

  return endcapBsKinks.at( det ).at( ring ).at( beam );
  
}





double LASConstants::GetAlignmentTubeBsKink( unsigned int beam ) const {
  
  if( beam >= 8U ) { // beam >= 0, since beam is unsigned
    throw cms::Exception( " [LASConstants::GetAlignmentTubeBsKink]" ) << " ** ERROR: no such beam: " << beam << "." << std::endl;
  }

  return alignmentTubeBsKinks.at( beam );

}





double LASConstants::GetTecRadius( unsigned int ring ) const {

  if( ring > 1U ) { // ring >= 0, since ring is unsigned
    throw cms::Exception( " [LASConstants::GetTecRadius]" ) << " ** ERROR: no such ring: " << ring << "." << std::endl;
  }

  return tecRadii.at( ring );

}





double LASConstants::GetAtRadius( void ) const {

  return atRadius;

}





double LASConstants::GetTecZPosition( unsigned int det, unsigned int disk ) const {

  if( ( det > 1 ) || ( disk > 8 ) ) {
    throw cms::Exception( " [LASConstants::GetTecZPosition]" ) << " ** ERROR: no such element: det " << det << ", disk " << disk << "." << std::endl;
  }

  if( det == 0 ) return tecZPositions.at( disk ); // tec+
  else return -1. * tecZPositions.at( disk ); // tec-

}





double LASConstants::GetTibZPosition( unsigned int pos ) const {

  if( pos > 5 ) {
    throw cms::Exception( " [LASConstants::GetTibZPosition]" ) << " ** ERROR: no such position: " << pos << "." << std::endl;
  }

  return tibZPositions.at( pos );

}





double LASConstants::GetTobZPosition( unsigned int pos ) const {

  if( pos > 5 ) {
    throw cms::Exception( " [LASConstants::GetTobZPosition]" ) << " ** ERROR: no such position: " << pos << "." << std::endl;
  }

  return tobZPositions.at( pos );

}





double LASConstants::GetTecBsZPosition( unsigned int det ) const {

  return tecBsZPosition;

}





double LASConstants::GetAtBsZPosition( void ) const {

  return atZPosition;

}





void LASConstants::InitContainers( void ) {

  // beam splitter kinks

  endcapBsKinks.resize( 2 ); // create two dets
  for( int det = 0; det < 2; ++det ) {
    endcapBsKinks.at( det ).resize( 2 ); // create two rings per det
    for( int ring = 0; ring < 2; ++ring ) {
      endcapBsKinks.at( det ).at( ring ).resize( 8 ); // 8 beams per ring
    }
  }

  alignmentTubeBsKinks.resize( 8 ); // 8 beams


  // radii
  tecRadii.resize( 2 );

  // z positions
  tecZPositions.resize( 9 );
  tibZPositions.resize( 6 );
  tobZPositions.resize( 6 );

}





void LASConstants::FillBsKinks( edm::ParameterSet const&  theBsKinkConf ) {

  // tec+
  endcapBsKinks.at( 0 ).at( 0 ) = theBsKinkConf.getParameter<std::vector<double> >( "LASTecPlusRing4BsKinks" );
  endcapBsKinks.at( 0 ).at( 1 ) = theBsKinkConf.getParameter<std::vector<double> >( "LASTecPlusRing6BsKinks" );

  // apply global offsets
  for( unsigned int ring = 0; ring < 2; ++ring ) {
    for( unsigned int beam = 0; beam < 8; ++beam ) {
      endcapBsKinks.at( 0 ).at( ring ).at( beam ) += theBsKinkConf.getParameter<double>( "TecPlusGlobalOffset" );
    }
  }

  // tec-
  endcapBsKinks.at( 1 ).at( 0 ) = theBsKinkConf.getParameter<std::vector<double> >( "LASTecMinusRing4BsKinks" );
  endcapBsKinks.at( 1 ).at( 1 ) = theBsKinkConf.getParameter<std::vector<double> >( "LASTecMinusRing6BsKinks" );

  // apply global offsets
  for( unsigned int ring = 0; ring < 2; ++ring ) {
    for( unsigned int beam = 0; beam < 8; ++beam ) {
      endcapBsKinks.at( 1 ).at( ring ).at( beam ) += theBsKinkConf.getParameter<double>( "TecMinusGlobalOffset" );
    }
  }

  // at
  alignmentTubeBsKinks = theBsKinkConf.getParameter<std::vector<double> >( "LASAlignmentTubeBsKinks" );

}





void LASConstants::FillRadii( edm::ParameterSet const&  theRadiiConf ) {
  
  tecRadii = theRadiiConf.getParameter<std::vector<double> >( "LASTecRadius" );
  atRadius = theRadiiConf.getParameter<double>( "LASAtRadius" );

}





void LASConstants::FillZPositions( edm::ParameterSet const& theZPosConf ) {

  tecZPositions  = theZPosConf.getParameter<std::vector<double> >( "LASTecZPositions" );
  tibZPositions  = theZPosConf.getParameter<std::vector<double> >( "LASTibZPositions" );
  tobZPositions  = theZPosConf.getParameter<std::vector<double> >( "LASTobZPositions" );
  tecBsZPosition = theZPosConf.getParameter<double>( "LASTecBeamSplitterZPosition" );
  atZPosition    = theZPosConf.getParameter<double>( "LASAtBeamsplitterZPosition" );

}