Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
/**
 * \file IntegratedCalibrationBase.cc
 *
 *  \author Gero Flucke
 *  \date August 2012
 *  $Revision: 1.77 $
 *  $Date: 2011/09/06 13:46:08 $
 *  (last update by $Author: mussgill $)
 */

#include "Alignment/CommonAlignmentAlgorithm/interface/IntegratedCalibrationBase.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"

// Already included in header:
//#include <vector>
//#include <utility>

//============================================================================
IntegratedCalibrationBase::IntegratedCalibrationBase(const edm::ParameterSet &cfg)
    : name_(cfg.getParameter<std::string>("calibrationName")) {}

//============================================================================
std::vector<IntegratedCalibrationBase::Values> IntegratedCalibrationBase::derivatives(
    const TrackingRecHit &hit,
    const TrajectoryStateOnSurface &tsos,
    const edm::EventSetup &setup,
    const EventInfo &eventInfo) const {
  // Prepare result vector, initialised all with 0.:
  std::vector<Values> result(this->numParameters(), Values(0., 0.));

  // Get non-zero derivatives and their index:
  std::vector<ValuesIndexPair> derivsIndexPairs;
  const unsigned int numNonZero = this->derivatives(derivsIndexPairs, hit, tsos, setup, eventInfo);

  // Put non-zero values into result:
  for (unsigned int i = 0; i < numNonZero; ++i) {
    const ValuesIndexPair &valuesIndex = derivsIndexPairs[i];
    result[valuesIndex.second] = valuesIndex.first;
  }

  return result;
}