#include <MillePedeFileReader.h>

Public Member Functions
const AlignPCLThresholds::threshold_map	getThresholdMap () const

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 > &, const std::shared_ptr< const AlignPCLThresholds > &)

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 *)

std::string	getStringFromHLS (PclHLS HLS)

void	readMillePedeLogFile ()

void	readMillePedeResultFile ()

Private Attributes
const std::string	millePedeLogFile_

const std::string	millePedeResFile_

int	Nrec_ {0}

const std::shared_ptr< const PedeLabelerBase >	pedeLabeler_

const std::shared_ptr< const AlignPCLThresholds >	theThresholds_

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

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

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

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

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}

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

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

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

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

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 16 of file MillePedeFileReader.h.

Member Enumeration Documentation

enum MillePedeFileReader::PclHLS : int

strongprivate

Enumerator
NotInPCL
TPBHalfBarrelXplus
TPBHalfBarrelXminus
TPEHalfCylinderXplusZplus
TPEHalfCylinderXminusZplus
TPEHalfCylinderXplusZminus
TPEHalfCylinderXminusZminus

Definition at line 51 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,
		const std::shared_ptr< const AlignPCLThresholds > &	theThresholds
	)

explicit

Definition at line 20 of file MillePedeFileReader.cc.

References read().

                                                                             :
   pedeLabeler_(pedeLabeler),
   theThresholds_(theThresholds),
   millePedeLogFile_(config.getParameter<std::string>("millePedeLogFile")),
   millePedeResFile_(config.getParameter<std::string>("millePedeResFile"))
 {
 }

virtual MillePedeFileReader::~MillePedeFileReader ( )

virtualdefault

Member Function Documentation

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

private

Definition at line 208 of file MillePedeFileReader.cc.

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

Referenced by readMillePedeResultFile().

                                    {
   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;
   }
 }

std::string MillePedeFileReader::getStringFromHLS ( MillePedeFileReader::PclHLS HLS )

private

Definition at line 257 of file MillePedeFileReader.cc.

References constexpr, Exception, multiplier_, TPBHalfBarrelXminus, TPBHalfBarrelXplus, TPEHalfCylinderXminusZminus, TPEHalfCylinderXminusZplus, TPEHalfCylinderXplusZminus, and TPEHalfCylinderXplusZplus.

                                                                   {
   switch (HLS)
     {
     case PclHLS::TPBHalfBarrelXminus         : return "TPBHalfBarrelXminus";
     case PclHLS::TPBHalfBarrelXplus          : return "TPBHalfBarrelXplus";
     case PclHLS::TPEHalfCylinderXminusZminus : return "TPEHalfCylinderXminusZminus";
     case PclHLS::TPEHalfCylinderXplusZminus  : return "TPEHalfCylinderXplusZminus";
     case PclHLS::TPEHalfCylinderXminusZplus  : return "TPEHalfCylinderXminusZplus";
     case PclHLS::TPEHalfCylinderXplusZplus   : return "TPEHalfCylinderXplusZplus";
     default: 
       throw cms::Exception("LogicError")
         << "@SUB=MillePedeFileReader::getStringFromHLS\n"
         << "Found an alignable structure not possible to map in the default AlignPCLThresholds partitions";
     }
 }

const AlignPCLThresholds::threshold_map MillePedeFileReader::getThresholdMap ( ) const

inline

Definition at line 45 of file MillePedeFileReader.h.

References theThresholds_.

45 {return theThresholds_.get()->getThreshold_Map (); }

MillePedeFileReader::theThresholds_

const std::shared_ptr< const AlignPCLThresholds > theThresholds_

Definition: MillePedeFileReader.h:74

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

inline

Definition at line 32 of file MillePedeFileReader.h.

References tXobs_.

32 { return tXobs_; }

MillePedeFileReader::tXobs_

std::array< double, 6 > tXobs_

Definition: MillePedeFileReader.h:93

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

inline

Definition at line 33 of file MillePedeFileReader.h.

References tXobsErr_.

33 { return tXobsErr_; }

MillePedeFileReader::tXobsErr_

std::array< double, 6 > tXobsErr_

Definition: MillePedeFileReader.h:94

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

inline

Definition at line 37 of file MillePedeFileReader.h.

References tYobs_.

37 { return tYobs_; }

MillePedeFileReader::tYobs_

std::array< double, 6 > tYobs_

Definition: MillePedeFileReader.h:98

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

inline

Definition at line 38 of file MillePedeFileReader.h.

References tYobsErr_.

38 { return tYobsErr_; }

MillePedeFileReader::tYobsErr_

std::array< double, 6 > tYobsErr_

Definition: MillePedeFileReader.h:99

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

inline

Definition at line 42 of file MillePedeFileReader.h.

References tZobs_.

42 { return tZobs_; }

MillePedeFileReader::tZobs_

std::array< double, 6 > tZobs_

Definition: MillePedeFileReader.h:103

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

inline

Definition at line 43 of file MillePedeFileReader.h.

References tZobsErr_.

43 { return tZobsErr_; }

MillePedeFileReader::tZobsErr_

std::array< double, 6 > tZobsErr_

Definition: MillePedeFileReader.h:104

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

inline

Definition at line 30 of file MillePedeFileReader.h.

References Xobs_.

30 { return Xobs_; }

MillePedeFileReader::Xobs_

std::array< double, 6 > Xobs_

Definition: MillePedeFileReader.h:91

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

inline

Definition at line 31 of file MillePedeFileReader.h.

References XobsErr_.

31 { return XobsErr_; }

MillePedeFileReader::XobsErr_

std::array< double, 6 > XobsErr_

Definition: MillePedeFileReader.h:92

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

inline

Definition at line 35 of file MillePedeFileReader.h.

References Yobs_.

35 { return Yobs_; }

MillePedeFileReader::Yobs_

std::array< double, 6 > Yobs_

Definition: MillePedeFileReader.h:96

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

inline

Definition at line 36 of file MillePedeFileReader.h.

References YobsErr_.

36 { return YobsErr_; }

MillePedeFileReader::YobsErr_

std::array< double, 6 > YobsErr_

Definition: MillePedeFileReader.h:97

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

inline

Definition at line 40 of file MillePedeFileReader.h.

References Zobs_.

40 { return Zobs_; }

MillePedeFileReader::Zobs_

std::array< double, 6 > Zobs_

Definition: MillePedeFileReader.h:101

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

inline

Definition at line 41 of file MillePedeFileReader.h.

References ZobsErr_.

41 { return ZobsErr_; }

MillePedeFileReader::ZobsErr_

std::array< double, 6 > ZobsErr_

Definition: MillePedeFileReader.h:102

void MillePedeFileReader::read ( )

Definition at line 31 of file MillePedeFileReader.cc.

References storeAlignments().

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

        {
   readMillePedeLogFile();
   readMillePedeResultFile();
 }

void MillePedeFileReader::readMillePedeLogFile ( )

private

Definition at line 48 of file MillePedeFileReader.cc.

References geometryCSVtoXML::line, heppy_check::logFile, readMillePedeResultFile(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by storeAlignments().

 {
   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_ < theThresholds_->getNrecords() ) {
       edm::LogInfo("MillePedeFileReader")<<"Number of records used "<<theThresholds_->getNrecords()<<std::endl;
           updateDB_   = false;
         }
       }
     }
 
   } else {
     edm::LogError("MillePedeFileReader") << "Could not read millepede log-file.";
 
     updateDB_   = false;
     Nrec_ = 0;
   }
 }

void MillePedeFileReader::readMillePedeResultFile ( )

private

Definition at line 81 of file MillePedeFileReader.cc.

References funct::abs(), getHLS(), geometryCSVtoXML::line, AlCaHLTBitMon_QueryRunRegistry::string, AlignPCLThresholds::theta_X, AlignPCLThresholds::theta_Y, AlignPCLThresholds::theta_Z, AlignPCLThresholds::X, AlignPCLThresholds::Y, and AlignPCLThresholds::Z.

Referenced by readMillePedeLogFile().

 {
   
   // cutoffs by coordinate and by alignable 
   std::map<std::string,std::array<float, 6 > > cutoffs_; 
   std::map<std::string,std::array<float, 6 > > significances_;
   std::map<std::string,std::array<float, 6 > > thresholds_;
   std::map<std::string,std::array<float, 6 > > errors_;
 
   std::vector<std::string> alignables_ = theThresholds_->getAlignableList();
   for (auto &ali : alignables_){
     cutoffs_[ali]       = theThresholds_->getCut(ali);
     significances_[ali] = theThresholds_->getSigCut(ali);
     thresholds_[ali]    = theThresholds_->getMaxMoveCut(ali);
     errors_[ali]        = theThresholds_->getMaxErrorCut(ali);
   }
 
   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);
     int detIndex = static_cast<int>(det);
     auto coord = static_cast<AlignPCLThresholds::coordType>(alignableIndex); 
     std::string detLabel = getStringFromHLS(det);
 
         if (det != PclHLS::NotInPCL) {
           switch (coord) {
           case AlignPCLThresholds::X:
             Xobs_[detIndex] = ObsMove;
             XobsErr_[detIndex] = ObsErr;
             break;
           case AlignPCLThresholds::Y:
             Yobs_[detIndex] = ObsMove;
             YobsErr_[detIndex] = ObsErr;
             break;
           case AlignPCLThresholds::Z:
             Zobs_[detIndex] = ObsMove;
             ZobsErr_[detIndex] = ObsErr;
             break;
           case AlignPCLThresholds::theta_X:
             tXobs_[detIndex] = ObsMove;
             tXobsErr_[detIndex] = ObsErr;
             break;
           case AlignPCLThresholds::theta_Y:
             tYobs_[detIndex] = ObsMove;
             tYobsErr_[detIndex] = ObsErr;
             break;
           case AlignPCLThresholds::theta_Z:
             tZobs_[detIndex] = ObsMove;
             tZobsErr_[detIndex] = ObsErr;
             break;
       default:
         edm::LogError("MillePedeFileReader") << "Currently not able to handle DOF " << coord 
                          << std::endl;
         break;
           }
         } else {
           continue;
         }
 
     edm::LogVerbatim("MillePedeFileReader")<<" alignableLabel: "<< alignableLabel <<" with alignableIndex "<<alignableIndex <<" detIndex"<< detIndex <<"\n"
                            <<" i.e. detLabel: "<< detLabel <<" ("<< coord <<")\n"
                            <<" has movement: "<< ObsMove <<" +/- "<< ObsErr <<"\n"
                            <<" cutoff (cutoffs_["<< detLabel <<"]["<< coord <<"]): "<<  cutoffs_[detLabel][alignableIndex] <<"\n"  
                            <<" significance (significances_["<< detLabel <<"]["<< coord <<"]): "<<  significances_[detLabel][alignableIndex] <<"\n"
                            <<" error thresolds (errors_["<< detLabel <<"]["<< coord <<"]): "<< errors_[detLabel][alignableIndex] <<"\n"
                            <<" max movement (thresholds_["<< detLabel <<"]["<< coord <<"]): "<< thresholds_[detLabel][alignableIndex] <<"\n"
                            <<"============="<< std::endl;
 
         if (std::abs(ObsMove) > thresholds_[detLabel][alignableIndex]) {
       edm::LogWarning("MillePedeFileReader")<<"Aborting payload creation."
                         <<" Exceeding maximum thresholds for movement: "<<std::abs(ObsMove)<<" for"<< detLabel <<"("<<coord<<")" ;    
           updateDB_    = false;
           break;
 
         } else if (std::abs(ObsMove) > cutoffs_[detLabel][alignableIndex]) {
 
           if (std::abs(ObsErr) > errors_[detLabel][alignableIndex]) {
         edm::LogWarning("MillePedeFileReader")<<"Aborting payload creation." 
                           <<" Exceeding maximum thresholds for error: "<<std::abs(ObsErr)<<" for"<< detLabel <<"("<<coord<<")" ;          
             updateDB_    = false;
             break;
           } else {
             if (std::abs(ObsMove/ObsErr) < significances_[detLabel][alignableIndex]) {
               continue;
             }
           }
           updateDB_ = true;
       edm::LogInfo("MillePedeFileReader")<<"This correction: "<<ObsMove<<"+/-" <<ObsErr<<" for "<< detLabel <<"("<<coord<<") will trigger a new Tracker Alignment payload!";
         }
       }
     }
   } else {
     edm::LogError("MillePedeFileReader") << "Could not read millepede result-file.";
 
     updateDB_   = false;
     Nrec_ = 0;
   }
 }

bool MillePedeFileReader::storeAlignments ( )

Definition at line 37 of file MillePedeFileReader.cc.

References readMillePedeLogFile().

Referenced by read(), and MillePedeAlignmentAlgorithm::storeAlignments().

                   {
   return updateDB_;
 }

Member Data Documentation

const std::string MillePedeFileReader::millePedeLogFile_

private

Definition at line 77 of file MillePedeFileReader.h.

const std::string MillePedeFileReader::millePedeResFile_

private

Definition at line 78 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.

Referenced by getStringFromHLS().

int MillePedeFileReader::Nrec_ {0}

private

Definition at line 89 of file MillePedeFileReader.h.

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

private

Definition at line 71 of file MillePedeFileReader.h.

const std::shared_ptr<const AlignPCLThresholds> MillePedeFileReader::theThresholds_

private

Definition at line 74 of file MillePedeFileReader.h.

Referenced by getThresholdMap().

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

private

Definition at line 93 of file MillePedeFileReader.h.

Referenced by getTXobs().

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

private

Definition at line 94 of file MillePedeFileReader.h.

Referenced by getTXobsErr().

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

private

Definition at line 98 of file MillePedeFileReader.h.

Referenced by getTYobs().

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

private

Definition at line 99 of file MillePedeFileReader.h.

Referenced by getTYobsErr().

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

private

Definition at line 103 of file MillePedeFileReader.h.

Referenced by getTZobs().

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

private

Definition at line 104 of file MillePedeFileReader.h.

Referenced by getTZobsErr().

bool MillePedeFileReader::updateDB_ {false}

private

Definition at line 88 of file MillePedeFileReader.h.

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

private

Definition at line 91 of file MillePedeFileReader.h.

Referenced by getXobs().

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

private

Definition at line 92 of file MillePedeFileReader.h.

Referenced by getXobsErr().

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

private

Definition at line 96 of file MillePedeFileReader.h.

Referenced by getYobs().

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

private

Definition at line 97 of file MillePedeFileReader.h.

Referenced by getYobsErr().

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

private

Definition at line 101 of file MillePedeFileReader.h.

Referenced by getZobs().

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

private

Definition at line 102 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