AlignmentParameterBuilder.cc

Go to the documentation of this file.

// This class's header should be first:
#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentParameterBuilder.h"

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

#include "Alignment/CommonAlignment/interface/Alignable.h"
#include "Alignment/CommonAlignment/interface/AlignmentParameters.h"

#include "Alignment/CommonAlignmentParametrization/interface/AlignmentParametersFactory.h"

#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentParameterSelector.h"
#include "Alignment/CommonAlignmentAlgorithm/interface/SelectionUserVariables.h"

#include <algorithm>

using namespace AlignmentParametersFactory;

//__________________________________________________________________________________________________
AlignmentParameterBuilder::AlignmentParameterBuilder(AlignableTracker *alignableTracker,
                                                     AlignableExtras *alignableExtras)
    : theAlignables(),
      theAlignableTracker(alignableTracker),
      theAlignableMuon(nullptr),
      theAlignableExtras(alignableExtras) {}

//__________________________________________________________________________________________________
AlignmentParameterBuilder::AlignmentParameterBuilder(AlignableTracker *alignableTracker,
                                                     AlignableMuon *alignableMuon,
                                                     AlignableExtras *alignableExtras)
    : theAlignables(),
      theAlignableTracker(alignableTracker),
      theAlignableMuon(alignableMuon),
      theAlignableExtras(alignableExtras) {}

//__________________________________________________________________________________________________
AlignmentParameterBuilder::AlignmentParameterBuilder(AlignableTracker *alignableTracker,
                                                     AlignableExtras *alignableExtras,
                                                     const edm::ParameterSet &pSet)
    : theAlignables(),
      theAlignableTracker(alignableTracker),
      theAlignableMuon(nullptr),
      theAlignableExtras(alignableExtras) {
  this->addAllSelections(pSet);
}

//__________________________________________________________________________________________________
AlignmentParameterBuilder::AlignmentParameterBuilder(AlignableTracker *alignableTracker,
                                                     AlignableMuon *alignableMuon,
                                                     AlignableExtras *alignableExtras,
                                                     const edm::ParameterSet &pSet)
    : theAlignables(),
      theAlignableTracker(alignableTracker),
      theAlignableMuon(alignableMuon),
      theAlignableExtras(alignableExtras) {
  this->addAllSelections(pSet);
}

const AlignableTracker *AlignmentParameterBuilder::alignableTracker() const { return theAlignableTracker; }

//__________________________________________________________________________________________________
void AlignmentParameterBuilder::addAllSelections(const edm::ParameterSet &pSet) {
  AlignmentParameterSelector selector(nullptr);
  std::vector<std::string> selsTypes(pSet.getParameter<std::vector<std::string> >("parameterTypes"));

  for (unsigned int i = 0; i < selsTypes.size(); ++i) {
    std::vector<std::string> selSetType(selector.decompose(selsTypes[i], ','));
    if (selSetType.size() != 2) {
      throw cms::Exception("BadConfig") << "AlignmentParameterBuilder"
                                        << "parameterTypes should contain 2 comma separated strings"
                                        << ", but found '" << selsTypes[i] << "'.";
    }
    this->addSelections(pSet.getParameter<edm::ParameterSet>(selSetType[0]),
                        AlignmentParametersFactory::parametersType(selSetType[1]));
  }
}

//__________________________________________________________________________________________________
unsigned int AlignmentParameterBuilder::addSelections(const edm::ParameterSet &pSet, ParametersType parType) {
  const unsigned int oldAliSize = theAlignables.size();

  AlignmentParameterSelector selector(theAlignableTracker, theAlignableMuon, theAlignableExtras);
  const unsigned int addedSets = selector.addSelections(pSet);

  const align::Alignables &alignables = selector.selectedAlignables();
  const std::vector<std::vector<char> > &paramSels = selector.selectedParameters();

  align::Alignables::const_iterator iAli = alignables.begin();
  std::vector<std::vector<char> >::const_iterator iParamSel = paramSels.begin();
  unsigned int nHigherLevel = 0;

  while (iAli != alignables.end() && iParamSel != paramSels.end()) {
    std::vector<bool> boolParSel;
    std::vector<char> parSel(*iParamSel);  // copy, since decodeParamSel may manipulate
    bool charSelIsGeneral = this->decodeParamSel(parSel, boolParSel);
    if (this->add(*iAli, parType, boolParSel))
      ++nHigherLevel;
    if (charSelIsGeneral)
      this->addFullParamSel((*iAli)->alignmentParameters(), parSel);

    ++iAli;
    ++iParamSel;
  }

  edm::LogInfo("Alignment") << "@SUB=AlignmentParameterBuilder::addSelections"
                            << " Added " << addedSets << " set(s) of alignables with "
                            << theAlignables.size() - oldAliSize << " alignables in total,"
                            << " of which " << nHigherLevel << " are higher level "
                            << "(using " << parametersTypeName(parType) << "AlignmentParameters).";

  return addedSets;
}

//__________________________________________________________________________________________________
bool AlignmentParameterBuilder::add(Alignable *alignable, ParametersType parType, const std::vector<bool> &sel) {
  AlignmentParameters *paras = AlignmentParametersFactory::createParameters(alignable, parType, sel);
  alignable->setAlignmentParameters(paras);
  theAlignables.push_back(alignable);

  const int aliTypeId = alignable->alignableObjectId();
  const bool isHigherLevel = (aliTypeId != align::AlignableDet && aliTypeId != align::AlignableDetUnit);
  return isHigherLevel;
}

//__________________________________________________________________________________________________
unsigned int AlignmentParameterBuilder::add(const align::Alignables &alignables,
                                            ParametersType parType,
                                            const std::vector<bool> &sel) {
  unsigned int nHigherLevel = 0;

  for (align::Alignables::const_iterator iAli = alignables.begin(); iAli != alignables.end(); ++iAli) {
    if (this->add(*iAli, parType, sel))
      ++nHigherLevel;
  }

  return nHigherLevel;
}

//__________________________________________________________________________________________________
void AlignmentParameterBuilder::fixAlignables(int n) {
  if (n < 1 || n > 3) {
    edm::LogError("BadArgument") << " n = " << n << " is not in [1,3]";
    return;
  }

  align::Alignables theNewAlignables;
  int i = 0;
  int imax = theAlignables.size();
  for (align::Alignables::const_iterator ia = theAlignables.begin(); ia != theAlignables.end(); ++ia) {
    i++;
    if (n == 1 && i > 1)
      theNewAlignables.push_back(*ia);
    else if (n == 2 && i > 1 && i < imax)
      theNewAlignables.push_back(*ia);
    else if (n == 3 && i > 2 && i < imax)
      theNewAlignables.push_back(*ia);
  }

  theAlignables = theNewAlignables;

  edm::LogInfo("Alignment") << "@SUB=AlignmentParameterBuilder::fixAlignables"
                            << "removing " << n << " alignables, so that " << theAlignables.size()
                            << " alignables left";
}

//__________________________________________________________________________________________________
bool AlignmentParameterBuilder::decodeParamSel(std::vector<char> &paramSelChar, std::vector<bool> &result) const {
  result.clear();
  // remove all spaces from paramSelChar - useful to group the parameters if they are many
  paramSelChar.erase(std::remove(paramSelChar.begin(), paramSelChar.end(), ' '), paramSelChar.end());

  bool anyNon01 = false;

  for (unsigned int pos = 0; pos < paramSelChar.size(); ++pos) {
    switch (paramSelChar[pos]) {
      default:
        anyNon01 = true;
        [[fallthrough]];
      case '1':
        result.push_back(true);
        break;
      case '0':
        result.push_back(false);
        break;
    }
  }

  return anyNon01;
}

//__________________________________________________________________________________________________
bool AlignmentParameterBuilder::addFullParamSel(AlignmentParameters *aliParams,
                                                const std::vector<char> &fullSel) const {
  if (!aliParams)
    return false;

  aliParams->setUserVariables(new SelectionUserVariables(fullSel));

  return true;
}