/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/Alignment/MillePedeAlignmentAlgorithm/src/MomentumDependentPedeLabeler.h

Go to the documentation of this file.

#ifndef MomentumDEPENDENTPEDELABELER_H
#define MomentumDEPENDENTPEDELABELER_H

#include <vector>
#include <map> 

#include <Alignment/MillePedeAlignmentAlgorithm/interface/PedeLabelerBase.h>

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

/***************************************
****************************************/
class MomentumDependentPedeLabeler : public PedeLabelerBase
{
 public:

  MomentumDependentPedeLabeler(const PedeLabelerBase::TopLevelAlignables& alignables,
                               const edm::ParameterSet& config);
  ~MomentumDependentPedeLabeler();
  
  unsigned int alignableLabel(Alignable *alignable) const;
  unsigned int alignableLabelFromParamAndInstance(Alignable *alignable,
                                                  unsigned int param,
                                                  unsigned int instance) const;
  unsigned int lasBeamLabel(unsigned int lasBeamId) const;
  unsigned int parameterLabel(unsigned int aliLabel, unsigned int parNum) const;
  unsigned int parameterLabel(Alignable *alignable, unsigned int parNum,
                              const AlignmentAlgorithmBase::EventInfo &eventInfo,
                              const TrajectoryStateOnSurface &tsos) const;
  bool hasSplitParameters(Alignable *alignable) const;
  unsigned int numberOfParameterInstances(Alignable *alignable, int param=-1) const;
  unsigned int maxNumberOfParameterInstances() const { return theMaxNumberOfParameterInstances; }

  unsigned int paramNumFromLabel(unsigned int paramLabel) const;
  unsigned int alignableLabelFromLabel(unsigned int label) const;
  Alignable* alignableFromLabel(unsigned int label) const;
  unsigned int lasBeamIdFromLabel(unsigned int label) const;

 private:
  typedef std::map <Alignable*, unsigned int> AlignableToIdMap;
  typedef AlignableToIdMap::value_type AlignableToIdPair;
  typedef std::pair <float,float> MomentumRange;
  typedef std::vector <MomentumRange> MomentumRangeVector;
  typedef std::map <unsigned int, MomentumRangeVector> MomentumRangeParamMap;
  typedef std::map <Alignable*, MomentumRangeParamMap> AlignableToMomentumRangeMap;
  typedef AlignableToMomentumRangeMap::value_type AlignableToMomentumRangePair;
  typedef std::map <unsigned int, Alignable*> IdToAlignableMap;
  typedef std::map <unsigned int, unsigned int> UintUintMap;
  
  std::vector<std::string> decompose(const std::string &s, std::string::value_type delimiter) const;
  std::vector<unsigned int> convertParamSel(const std::string &selString) const;
  unsigned int buildMomentumDependencyMap(AlignableTracker *aliTracker, AlignableMuon* aliMuon,
                                          AlignableExtras *extras,
                                          const edm::ParameterSet &config);

  unsigned int buildMap(const std::vector<Alignable*> &alis);
  unsigned int buildReverseMap();

  // data members
  AlignableToIdMap            theAlignableToIdMap; 
  AlignableToMomentumRangeMap theAlignableToMomentumRangeMap; 
  IdToAlignableMap            theIdToAlignableMap; 
  UintUintMap                 theLasBeamToLabelMap;  
  UintUintMap                 theLabelToLasBeamMap; 
  const MomentumRange         theOpenMomentumRange;
  unsigned int                theMaxNumberOfParameterInstances;
};

#endif