#include <MillePedeFileReader.h>

Public Member Functions
const std::array< double, 6 > &	getTXobs () const

const std::array< double, 6 > &	getTXobsErr () const

const std::array< double, 6 > &	getTYobs () const

const std::array< double, 6 > &	getTYobsErr () const

const std::array< double, 6 > &	getTZobs () const

const std::array< double, 6 > &	getTZobsErr () const

const std::array< double, 6 > &	getXobs () const

const std::array< double, 6 > &	getXobsErr () const

const std::array< double, 6 > &	getYobs () const

const std::array< double, 6 > &	getYobsErr () const

const std::array< double, 6 > &	getZobs () const

const std::array< double, 6 > &	getZobsErr () const

	MillePedeFileReader (const edm::ParameterSet &, const std::shared_ptr< const PedeLabelerBase > &)

void	read ()

bool	storeAlignments ()

virtual	~MillePedeFileReader ()=default

Private Types
enum	PclHLS : int { PclHLS::NotInPCL = -1, PclHLS::TPBHalfBarrelXplus = 2, PclHLS::TPBHalfBarrelXminus = 3, PclHLS::TPEHalfCylinderXplusZplus = 4, PclHLS::TPEHalfCylinderXminusZplus = 5, PclHLS::TPEHalfCylinderXplusZminus = 0, PclHLS::TPEHalfCylinderXminusZminus = 1 }

Private Member Functions
PclHLS	getHLS (const Alignable *)

void	readMillePedeLogFile ()

void	readMillePedeResultFile ()

Private Attributes
const std::array< double, 6 >	cutoffs_

const double	maxErrorCut_

const double	maxMoveCut_

const std::string	millePedeLogFile_

const std::string	millePedeResFile_

int	Nrec_ {0}

const std::shared_ptr< const PedeLabelerBase >	pedeLabeler_

const double	sigCut_

const double	tXcut_

std::array< double, 6 >	tXobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	tXobsErr_ = {{0.,0.,0.,0.,0.,0.}}

const double	tYcut_

std::array< double, 6 >	tYobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	tYobsErr_ = {{0.,0.,0.,0.,0.,0.}}

const double	tZcut_

std::array< double, 6 >	tZobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	tZobsErr_ = {{0.,0.,0.,0.,0.,0.}}

bool	updateDB_ {false}

const double	Xcut_

std::array< double, 6 >	Xobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	XobsErr_ = {{0.,0.,0.,0.,0.,0.}}

const double	Ycut_

std::array< double, 6 >	Yobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	YobsErr_ = {{0.,0.,0.,0.,0.,0.}}

const double	Zcut_

std::array< double, 6 >	Zobs_ = {{0.,0.,0.,0.,0.,0.}}

std::array< double, 6 >	ZobsErr_ = {{0.,0.,0.,0.,0.,0.}}

Static Private Attributes
static std::array< double, 6 >	multiplier_

Detailed Description

Definition at line 15 of file MillePedeFileReader.h.

Member Enumeration Documentation

enum MillePedeFileReader::PclHLS : int

strongprivate

Enumerator
NotInPCL
TPBHalfBarrelXplus
TPBHalfBarrelXminus
TPEHalfCylinderXplusZplus
TPEHalfCylinderXminusZplus
TPEHalfCylinderXplusZminus
TPEHalfCylinderXminusZminus

Definition at line 46 of file MillePedeFileReader.h.

                     : int { NotInPCL = -1,
                             TPBHalfBarrelXplus = 2,
                             TPBHalfBarrelXminus = 3,
                             TPEHalfCylinderXplusZplus = 4,
                             TPEHalfCylinderXminusZplus = 5,
                             TPEHalfCylinderXplusZminus = 0,
                             TPEHalfCylinderXminusZminus = 1};

Constructor & Destructor Documentation

MillePedeFileReader::MillePedeFileReader	(	const edm::ParameterSet &	config,
		const std::shared_ptr< const PedeLabelerBase > &	pedeLabeler
	)

explicit

Definition at line 20 of file MillePedeFileReader.cc.

                                                                                :
   pedeLabeler_(pedeLabeler),
   millePedeLogFile_(config.getParameter<std::string>("millePedeLogFile")),
   millePedeResFile_(config.getParameter<std::string>("millePedeResFile")),
 
   sigCut_     (config.getParameter<double>("sigCut")),
   Xcut_       (config.getParameter<double>("Xcut")),
   tXcut_      (config.getParameter<double>("tXcut")),
   Ycut_       (config.getParameter<double>("Ycut")),
   tYcut_      (config.getParameter<double>("tYcut")),
   Zcut_       (config.getParameter<double>("Zcut")),
   tZcut_      (config.getParameter<double>("tZcut")),
   maxMoveCut_ (config.getParameter<double>("maxMoveCut")),
   maxErrorCut_(config.getParameter<double>("maxErrorCut"))
 {
 }

virtual MillePedeFileReader::~MillePedeFileReader ( )

virtualdefault

Member Function Documentation

MillePedeFileReader::PclHLS MillePedeFileReader::getHLS ( const Alignable * alignable )

private

Definition at line 179 of file MillePedeFileReader.cc.

References Alignable::alignableObjectId(), Exception, Alignable::id(), align::TPBHalfBarrel, and align::TPEHalfCylinder.

                                    {
   if (!alignable) return PclHLS::NotInPCL;
 
   const auto& tns = pedeLabeler_->alignableTracker()->trackerNameSpace();
 
   switch (alignable->alignableObjectId()) {
   case align::TPBHalfBarrel:
     switch (tns.tpb().halfBarrelNumber(alignable->id())) {
     case 1: return PclHLS::TPBHalfBarrelXminus;
     case 2: return PclHLS::TPBHalfBarrelXplus;
     default:
       throw cms::Exception("LogicError")
         << "@SUB=MillePedeFileReader::getHLS\n"
         << "Found a pixel half-barrel number that should not exist: "
         << tns.tpb().halfBarrelNumber(alignable->id());
     }
   case align::TPEHalfCylinder:
     switch (tns.tpe().endcapNumber(alignable->id())) {
     case 1:
       switch (tns.tpe().halfCylinderNumber(alignable->id())) {
       case 1: return PclHLS::TPEHalfCylinderXminusZminus;
       case 2: return PclHLS::TPEHalfCylinderXplusZminus;
       default:
         throw cms::Exception("LogicError")
           << "@SUB=MillePedeFileReader::getHLS\n"
           << "Found a pixel half-cylinder number that should not exist: "
           << tns.tpe().halfCylinderNumber(alignable->id());
       }
     case 2:
       switch (tns.tpe().halfCylinderNumber(alignable->id())) {
       case 1: return PclHLS::TPEHalfCylinderXminusZplus;
       case 2: return PclHLS::TPEHalfCylinderXplusZplus;
       default:
         throw cms::Exception("LogicError")
           << "@SUB=MillePedeFileReader::getHLS\n"
           << "Found a pixel half-cylinder number that should not exist: "
           << tns.tpe().halfCylinderNumber(alignable->id());
       }
     default:
       throw cms::Exception("LogicError")
         << "@SUB=MillePedeFileReader::getHLS\n"
         << "Found a pixel endcap number that should not exist: "
         << tns.tpe().endcapNumber(alignable->id());
     }
   default: return PclHLS::NotInPCL;
   }
 }

const std::array<double, 6>& MillePedeFileReader::getTXobs ( ) const

inline

Definition at line 29 of file MillePedeFileReader.h.

References tXobs_.

29 { return tXobs_; }

MillePedeFileReader::tXobs_

std::array< double, 6 > tXobs_

Definition: MillePedeFileReader.h:96

const std::array<double, 6>& MillePedeFileReader::getTXobsErr ( ) const

inline

Definition at line 30 of file MillePedeFileReader.h.

References tXobsErr_.

30 { return tXobsErr_; }

MillePedeFileReader::tXobsErr_

std::array< double, 6 > tXobsErr_

Definition: MillePedeFileReader.h:97

const std::array<double, 6>& MillePedeFileReader::getTYobs ( ) const

inline

Definition at line 34 of file MillePedeFileReader.h.

References tYobs_.

34 { return tYobs_; }

MillePedeFileReader::tYobs_

std::array< double, 6 > tYobs_

Definition: MillePedeFileReader.h:101

const std::array<double, 6>& MillePedeFileReader::getTYobsErr ( ) const

inline

Definition at line 35 of file MillePedeFileReader.h.

References tYobsErr_.

35 { return tYobsErr_; }

MillePedeFileReader::tYobsErr_

std::array< double, 6 > tYobsErr_

Definition: MillePedeFileReader.h:102

const std::array<double, 6>& MillePedeFileReader::getTZobs ( ) const

inline

Definition at line 39 of file MillePedeFileReader.h.

References tZobs_.

39 { return tZobs_; }

MillePedeFileReader::tZobs_

std::array< double, 6 > tZobs_

Definition: MillePedeFileReader.h:106

const std::array<double, 6>& MillePedeFileReader::getTZobsErr ( ) const

inline

Definition at line 40 of file MillePedeFileReader.h.

References tZobsErr_.

40 { return tZobsErr_; }

MillePedeFileReader::tZobsErr_

std::array< double, 6 > tZobsErr_

Definition: MillePedeFileReader.h:107

const std::array<double, 6>& MillePedeFileReader::getXobs ( ) const

inline

Definition at line 27 of file MillePedeFileReader.h.

References Xobs_.

27 { return Xobs_; }

MillePedeFileReader::Xobs_

std::array< double, 6 > Xobs_

Definition: MillePedeFileReader.h:94

const std::array<double, 6>& MillePedeFileReader::getXobsErr ( ) const

inline

Definition at line 28 of file MillePedeFileReader.h.

References XobsErr_.

28 { return XobsErr_; }

MillePedeFileReader::XobsErr_

std::array< double, 6 > XobsErr_

Definition: MillePedeFileReader.h:95

const std::array<double, 6>& MillePedeFileReader::getYobs ( ) const

inline

Definition at line 32 of file MillePedeFileReader.h.

References Yobs_.

32 { return Yobs_; }

MillePedeFileReader::Yobs_

std::array< double, 6 > Yobs_

Definition: MillePedeFileReader.h:99

const std::array<double, 6>& MillePedeFileReader::getYobsErr ( ) const

inline

Definition at line 33 of file MillePedeFileReader.h.

References YobsErr_.

33 { return YobsErr_; }

MillePedeFileReader::YobsErr_

std::array< double, 6 > YobsErr_

Definition: MillePedeFileReader.h:100

const std::array<double, 6>& MillePedeFileReader::getZobs ( ) const

inline

Definition at line 37 of file MillePedeFileReader.h.

References Zobs_.

37 { return Zobs_; }

MillePedeFileReader::Zobs_

std::array< double, 6 > Zobs_

Definition: MillePedeFileReader.h:104

const std::array<double, 6>& MillePedeFileReader::getZobsErr ( ) const

inline

Definition at line 38 of file MillePedeFileReader.h.

References ZobsErr_.

38 { return ZobsErr_; }

MillePedeFileReader::ZobsErr_

std::array< double, 6 > ZobsErr_

Definition: MillePedeFileReader.h:105

void MillePedeFileReader::read ( )

Definition at line 39 of file MillePedeFileReader.cc.

Referenced by edmIntegrityCheck.PublishToFileSystem::get(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::goto(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches(), and MillePedeAlignmentAlgorithm::storeAlignments().

        {
   readMillePedeLogFile();
   readMillePedeResultFile();
 }

void MillePedeFileReader::readMillePedeLogFile ( )

private

Definition at line 56 of file MillePedeFileReader.cc.

References geometryCSVtoXML::line, heppy_check::logFile, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   std::ifstream logFile;
   logFile.open(millePedeLogFile_.c_str());
 
   if (logFile.is_open()) {
     edm::LogInfo("MillePedeFileReader") << "Reading millepede log-file";
     std::string line;
 
     while (getline(logFile, line)) {
       std::string Nrec_string = "NREC =";
 
       if (line.find(Nrec_string) != std::string::npos) {
         std::istringstream iss(line);
         std::string trash;
         iss >> trash >> trash >> Nrec_;
 
         if (Nrec_ < 25000) {
           updateDB_   = false;
         }
       }
     }
 
   } else {
     edm::LogError("MillePedeFileReader") << "Could not read millepede log-file.";
 
     updateDB_   = false;
     Nrec_ = 0;
   }
 }

void MillePedeFileReader::readMillePedeResultFile ( )

private

Definition at line 88 of file MillePedeFileReader.cc.

References funct::abs(), geometryCSVtoXML::line, AlCaHLTBitMon_QueryRunRegistry::string, and unpackBuffers-CaloStage2::token.

 {
   updateDB_ = false;
   std::ifstream resFile;
   resFile.open(millePedeResFile_.c_str());
 
   if (resFile.is_open()) {
     edm::LogInfo("MillePedeFileReader") << "Reading millepede result-file";
 
     std::string line;
     getline(resFile, line); // drop first line
 
     while (getline(resFile, line)) {
       std::istringstream iss(line);
 
       std::vector<std::string> tokens;
       std::string token;
       while (iss >> token) {
         tokens.push_back(token);
       }
 
       if (tokens.size() > 4 /*3*/) {
 
         auto alignableLabel = std::stoul(tokens[0]);
         auto alignableIndex = alignableLabel % 10 - 1;
         const auto alignable = pedeLabeler_->alignableFromLabel(alignableLabel);
 
         double ObsMove = std::stof(tokens[3]) * multiplier_[alignableIndex];
         double ObsErr  = std::stof(tokens[4]) * multiplier_[alignableIndex];
 
         auto det = getHLS(alignable);
 
         if (det != PclHLS::NotInPCL) {
           switch (alignableIndex) {
           case 0:
             Xobs_[static_cast<int>(det)] = ObsMove;
             XobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           case 1:
             Yobs_[static_cast<int>(det)] = ObsMove;
             YobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           case 2:
             Zobs_[static_cast<int>(det)] = ObsMove;
             ZobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           case 3:
             tXobs_[static_cast<int>(det)] = ObsMove;
             tXobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           case 4:
             tYobs_[static_cast<int>(det)] = ObsMove;
             tYobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           case 5:
             tZobs_[static_cast<int>(det)] = ObsMove;
             tZobsErr_[static_cast<int>(det)] = ObsErr;
             break;
           }
         } else {
           continue;
         }
 
         if (std::abs(ObsMove) > maxMoveCut_) {
           updateDB_    = false;
           break;
 
         } else if (std::abs(ObsMove) > cutoffs_[alignableIndex]) {
 
           if (std::abs(ObsErr) > maxErrorCut_) {
             updateDB_    = false;
             break;
           } else {
             if (std::abs(ObsMove/ObsErr) < sigCut_) {
               continue;
             }
           }
           updateDB_ = true;
         }
       }
     }
   } else {
     edm::LogError("MillePedeFileReader") << "Could not read millepede result-file.";
 
     updateDB_   = false;
     Nrec_ = 0;
   }
 }

bool MillePedeFileReader::storeAlignments ( )

Definition at line 45 of file MillePedeFileReader.cc.

Referenced by MillePedeAlignmentAlgorithm::storeAlignments().

                   {
   return updateDB_;
 }

Member Data Documentation

const std::array<double, 6> MillePedeFileReader::cutoffs_

private

Initial value:

= {{ Xcut_, Ycut_, Zcut_,

tXcut_, tYcut_, tZcut_}}

Definition at line 88 of file MillePedeFileReader.h.

const double MillePedeFileReader::maxErrorCut_

private

Definition at line 78 of file MillePedeFileReader.h.

const double MillePedeFileReader::maxMoveCut_

private

Definition at line 78 of file MillePedeFileReader.h.

const std::string MillePedeFileReader::millePedeLogFile_

private

Definition at line 68 of file MillePedeFileReader.h.

const std::string MillePedeFileReader::millePedeResFile_

private

Definition at line 69 of file MillePedeFileReader.h.

std::array< double, 6 > MillePedeFileReader::multiplier_

staticprivate

Initial value:

= {{ 10000.,      
                                                            10000.,      
                                                            10000.,      
                                                            1000000.,    
                                                            1000000.,    
                                                            1000000. }}

Definition at line 81 of file MillePedeFileReader.h.

int MillePedeFileReader::Nrec_ {0}

private

Definition at line 92 of file MillePedeFileReader.h.

const std::shared_ptr<const PedeLabelerBase> MillePedeFileReader::pedeLabeler_

private

Definition at line 65 of file MillePedeFileReader.h.

const double MillePedeFileReader::sigCut_

private

Definition at line 72 of file MillePedeFileReader.h.

const double MillePedeFileReader::tXcut_

private

Definition at line 74 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::tXobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 96 of file MillePedeFileReader.h.

Referenced by getTXobs().

std::array<double, 6> MillePedeFileReader::tXobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 97 of file MillePedeFileReader.h.

Referenced by getTXobsErr().

const double MillePedeFileReader::tYcut_

private

Definition at line 75 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::tYobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 101 of file MillePedeFileReader.h.

Referenced by getTYobs().

std::array<double, 6> MillePedeFileReader::tYobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 102 of file MillePedeFileReader.h.

Referenced by getTYobsErr().

const double MillePedeFileReader::tZcut_

private

Definition at line 76 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::tZobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 106 of file MillePedeFileReader.h.

Referenced by getTZobs().

std::array<double, 6> MillePedeFileReader::tZobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 107 of file MillePedeFileReader.h.

Referenced by getTZobsErr().

bool MillePedeFileReader::updateDB_ {false}

private

Definition at line 91 of file MillePedeFileReader.h.

const double MillePedeFileReader::Xcut_

private

Definition at line 74 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::Xobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 94 of file MillePedeFileReader.h.

Referenced by getXobs().

std::array<double, 6> MillePedeFileReader::XobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 95 of file MillePedeFileReader.h.

Referenced by getXobsErr().

const double MillePedeFileReader::Ycut_

private

Definition at line 75 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::Yobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 99 of file MillePedeFileReader.h.

Referenced by getYobs().

std::array<double, 6> MillePedeFileReader::YobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 100 of file MillePedeFileReader.h.

Referenced by getYobsErr().

const double MillePedeFileReader::Zcut_

private

Definition at line 76 of file MillePedeFileReader.h.

std::array<double, 6> MillePedeFileReader::Zobs_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 104 of file MillePedeFileReader.h.

Referenced by getZobs().

std::array<double, 6> MillePedeFileReader::ZobsErr_ = {{0.,0.,0.,0.,0.,0.}}

private

Definition at line 105 of file MillePedeFileReader.h.

Referenced by getZobsErr().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation