/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Alignment/MillePedeAlignmentAlgorithm/src/PedeSteerer.h

Go to the documentation of this file.

#ifndef PEDESTEERER_H
#define PEDESTEERER_H

#include <vector>
#include <map> 
#include <set> 
#include <string>
// forward ofstream:
#include <iosfwd> 

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

class Alignable;
class AlignableTracker;
class AlignableMuon;
class AlignableExtras;
class AlignmentParameterStore;
class PedeLabelerBase;

/***************************************
****************************************/
class PedeSteerer
{
 public:
  PedeSteerer(AlignableTracker *aliTracker, AlignableMuon *aliMuon, AlignableExtras *aliExtras,
              AlignmentParameterStore *store,
              const PedeLabelerBase *labels, const edm::ParameterSet &config,
              const std::string &defaultDir, bool noSteerFiles);
  ~PedeSteerer();
    
  bool isNoHiera(const Alignable* ali) const;

  void buildSubSteer(AlignableTracker *aliTracker, AlignableMuon *aliMuon, AlignableExtras *aliExtras);
  std::string buildMasterSteer(const std::vector<std::string> &binaryFiles);
  int runPede(const std::string &masterSteer) const;
  void correctToReferenceSystem();
  bool isCorrectToRefSystem(const std::vector<Alignable*> &coordDefiners) const;


  double cmsToPedeFactor(unsigned int parNum) const;
  int parameterSign() const { return myParameterSign; }
  const std::string& directory() const { return myDirectory;}

 private:
  typedef std::map<const Alignable*,std::vector<float> > AlignablePresigmasMap;

  bool checkParameterChoices(const std::vector<Alignable*> &alis) const;
  unsigned int buildNoHierarchyCollection(const std::vector<Alignable*> &alis);
  std::pair<unsigned int, unsigned int> fixParameters(const std::vector<Alignable*> &alignables,
                                                      const std::string &file);
  int fixParameter(Alignable *ali, unsigned int iRunRange,
                   unsigned int iParam, char selector, std::ofstream* &filePtr,
                   const std::string &fileName);
  std::vector<Alignable*> selectCoordinateAlis(const std::vector<Alignable*> &alignables) const;
  void defineCoordinates(const std::vector<Alignable*> &alis, Alignable *aliMaster,
                         const std::string &fileName);

  unsigned int hierarchyConstraints(const std::vector<Alignable*> &alis, const std::string &file);
  void hierarchyConstraint(const Alignable *ali, const std::vector<Alignable*> &components,
                           std::ofstream &file) const;

  unsigned int presigmas(const std::vector<edm::ParameterSet> &cffPresi,
                         const std::string &fileName, const std::vector<Alignable*> &alis,
                         AlignableTracker *aliTracker, AlignableMuon *aliMuon, AlignableExtras *aliExtras);
  unsigned int presigmasFile(const std::string &fileName, const std::vector<Alignable*> &alis,
                             const AlignablePresigmasMap &aliPresisMap); 
  std::string fileName(const std::string &addendum) const;
  std::ofstream* createSteerFile(const std::string &name, bool addToList);

  // data members
  const AlignmentParameterStore *myParameterStore; 
  const PedeLabelerBase         *myLabels; 

  edm::ParameterSet myConfig;
  std::string myDirectory; 
  bool myNoSteerFiles; 
  bool myIsSteerFileDebug; 
  int myParameterSign; 
  double theMinHieraConstrCoeff; 
  unsigned int theMinHieraParPerConstr; 

  std::vector<std::string> mySteeringFiles; 

  std::set<const Alignable*> myNoHieraCollection; 
  Alignable *theCoordMaster;                      
  std::vector<Alignable*> theCoordDefiners;      
};

#endif