/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/Alignment/CommonAlignment/src/AlignableExtras.cc

Go to the documentation of this file.

#include "FWCore/Utilities/interface/Exception.h"

// Geometry
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"

// Alignment
#include "Alignment/CommonAlignment/interface/AlignableBeamSpot.h"

#include "CondFormats/Alignment/interface/Alignments.h"
#include "CondFormats/Alignment/interface/AlignmentErrors.h"
#include "CondFormats/Alignment/interface/AlignmentSorter.h"

#include "Alignment/CommonAlignment/interface/AlignableExtras.h"

//__________________________________________________________________________________________________
AlignableExtras::AlignableExtras()
{
  align::Alignables& alis = alignableLists_.get("BeamSpot");
  alis.push_back(new AlignableBeamSpot());
  components_.push_back(alis.back());
}

//__________________________________________________________________________________________________
void AlignableExtras::dump( void ) const
{
  Alignables comp = this->components();

  // Dump this
  edm::LogInfo("AlignableDump") 
    << " AlignableExtras knows " << comp.size() << " alignable(s)" << std::endl;

  // Dump components
  for ( Alignables::iterator i=comp.begin(); i!=comp.end(); i++ )
    (*i)->dump();
}

//__________________________________________________________________________________________________
Alignments* AlignableExtras::alignments( void ) const
{
  align::Alignables comp = this->components();
  Alignments* m_alignments = new Alignments();
  // Add components recursively
  for ( align::Alignables::iterator i=comp.begin(); i!=comp.end(); i++ )
    {
      Alignments* tmpAlignments = (*i)->alignments();
      std::copy( tmpAlignments->m_align.begin(), tmpAlignments->m_align.end(), 
                 std::back_inserter(m_alignments->m_align) );
          delete tmpAlignments;
    }

  std::sort( m_alignments->m_align.begin(), m_alignments->m_align.end(), 
             lessAlignmentDetId<AlignTransform>() );

  return m_alignments;
}

//__________________________________________________________________________________________________
AlignmentErrors* AlignableExtras::alignmentErrors( void ) const
{
  align::Alignables comp = this->components();
  AlignmentErrors* m_alignmentErrors = new AlignmentErrors();

  // Add components recursively
  for ( align::Alignables::iterator i=comp.begin(); i!=comp.end(); i++ )
    {
          AlignmentErrors* tmpAlignmentErrors = (*i)->alignmentErrors();
      std::copy( tmpAlignmentErrors->m_alignError.begin(), tmpAlignmentErrors->m_alignError.end(), 
                 std::back_inserter(m_alignmentErrors->m_alignError) );
          delete tmpAlignmentErrors;
    }
  
  std::sort( m_alignmentErrors->m_alignError.begin(), m_alignmentErrors->m_alignError.end(), 
             lessAlignmentDetId<AlignTransformError>() );

  return m_alignmentErrors;
}

void AlignableExtras::initializeBeamSpot(double x, double y, double z,
                                         double dxdz, double dydz)
{
  align::Alignables& alis = beamSpot();
  AlignableBeamSpot * aliBS = dynamic_cast<AlignableBeamSpot*>(alis.back());
  if (aliBS) {
    aliBS->initialize(x, y, z, dxdz, dydz);
  } else {
    edm::LogError("AlignableExtras") 
      << " AlignableBeamSpot not available. Cannot initialize!" << std::endl;
  }
}