AlignmentParametersIO.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "Alignment/CommonAlignment/interface/Alignable.h"
#include "Alignment/CommonAlignmentParametrization/interface/RigidBodyAlignmentParameters.h"

// this class's header
#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentParametersIO.h"

//--------------------------------------------------------------------------------------------------
// write one set of original parameters
int AlignmentParametersIO::writeOneOrigRigidBody(Alignable *ali)
{

  AlignmentParameters *par = ali->alignmentParameters();
  AlignmentParameters *parBack = (par ? par->clone(par->parameters(), par->covariance()) : 0);

  ali->setAlignmentParameters(new RigidBodyAlignmentParameters(ali, true));
  int iret = this->writeOne(ali);

  ali->setAlignmentParameters(parBack); // deletes the above created RigidBodyAlignmentParameters

  return iret;
}


//-----------------------------------------------------------------------------
// write many parameters
int 
AlignmentParametersIO::write(const align::Alignables& alivec, 
                             bool validCheck) 
{
  int icount=0;
  for(align::Alignables::const_iterator it=alivec.begin();
      it!=alivec.end(); ++it) {
    if ((*it)->alignmentParameters()->isValid() || !(validCheck)) {
      icount++;
      int iret=writeOne(*it);
      if (iret!=0) return iret;
    }
  }
  edm::LogInfo("Alignment") << "@SUB=AlignmentParametersIO::write"
                            << "Wrote " << icount << " out of " << alivec.size() << " parameters";
  return 0;
}


//-----------------------------------------------------------------------------
// write many original parameters
int 
AlignmentParametersIO::writeOrigRigidBody(const align::Alignables& alivec, bool validCheck)
{
  int icount = 0;
  for(align::Alignables::const_iterator it = alivec.begin(); it != alivec.end(); ++it) {
    if (!validCheck || (*it)->alignmentParameters()->isValid()) {
      ++icount;
      int iret = this->writeOneOrigRigidBody(*it);
      if (iret != 0) return iret;
    }
  }
  edm::LogInfo("Alignment") << "@SUB=AlignmentParametersIO::writeOrigRigidBody"
                            << "Wrote " << icount << " out of " << alivec.size()
                            << " original parameters.";
  return 0;
}

//-----------------------------------------------------------------------------
// read many parameters

align::Parameters 
AlignmentParametersIO::read(const align::Alignables& alivec, int& ierr) 
{
  align::Parameters retvec;
  int ierr2;
  int icount=0;
  for(align::Alignables::const_iterator it=alivec.begin();
    it!=alivec.end(); ++it) {
    AlignmentParameters* ad=readOne(*it, ierr2);
    if (ad!=0 && ierr2==0) { retvec.push_back(ad); icount++; }
  }
  edm::LogInfo("Alignment") << "@SUB-AlignmentParametersIO::write"
                            << "Read " << icount << " out of " << alivec.size() << " parameters";
  return retvec;
}