MillePedeFileReader Class Reference

#include <MillePedeFileReader.h>

Public Member Functions
std::array< double, 6 > const &	getTXobs () const
std::array< double, 6 > const &	getTXobsErr () const
std::array< double, 6 > const &	getTYobs () const
std::array< double, 6 > const &	getTYobsErr () const
std::array< double, 6 > const &	getTZobs () const
std::array< double, 6 > const &	getTZobsErr () const
std::array< double, 6 > const &	getXobs () const
std::array< double, 6 > const &	getXobsErr () const
std::array< double, 6 > const &	getYobs () const
std::array< double, 6 > const &	getYobsErr () const
std::array< double, 6 > const &	getZobs () const
std::array< double, 6 > const &	getZobsErr () const
	MillePedeFileReader (const edm::ParameterSet &)
void	read ()
bool	storeAlignments ()
	~MillePedeFileReader ()

Private Member Functions
void	readMillePedeLogFile ()
void	readMillePedeResultFile ()

Private Attributes
double	Cutoffs [6]
bool	Error = false
bool	HitErrorMax = false
bool	HitMax = false
double	maxErrorCut_
double	maxMoveCut_
std::string	millePedeLogFile_
std::string	millePedeResFile_
bool	Movements = false
int	Nrec = 0
bool	PedeSuccess = false
double	sigCut_
bool	Significant = false
double	tXcut_
std::array< double, 6 >	tXobs = {{0.,0.,0.,0.,0.,0.}}
std::array< double, 6 >	tXobsErr = {{0.,0.,0.,0.,0.,0.}}
double	tYcut_
std::array< double, 6 >	tYobs = {{0.,0.,0.,0.,0.,0.}}
std::array< double, 6 >	tYobsErr = {{0.,0.,0.,0.,0.,0.}}
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
double	Xcut_
std::array< double, 6 >	Xobs = {{0.,0.,0.,0.,0.,0.}}
std::array< double, 6 >	XobsErr = {{0.,0.,0.,0.,0.,0.}}
double	Ycut_
std::array< double, 6 >	Yobs = {{0.,0.,0.,0.,0.,0.}}
std::array< double, 6 >	YobsErr = {{0.,0.,0.,0.,0.,0.}}
double	Zcut_
std::array< double, 6 >	Zobs = {{0.,0.,0.,0.,0.,0.}}
std::array< double, 6 >	ZobsErr = {{0.,0.,0.,0.,0.,0.}}

Detailed Description

Definition at line 16 of file MillePedeFileReader.h.

Constructor & Destructor Documentation

MillePedeFileReader::MillePedeFileReader ( const edm::ParameterSet &  config )

explicit

Definition at line 103 of file MillePedeFileReader.h.

                                                   :
  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"))
{
}

MillePedeFileReader::~MillePedeFileReader ( )

inline

Definition at line 22 of file MillePedeFileReader.h.

{}

Member Function Documentation

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

inline

Definition at line 29 of file MillePedeFileReader.h.

References tXobs.

{ return tXobs;    }

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

inline

Definition at line 30 of file MillePedeFileReader.h.

References tXobsErr.

{ return tXobsErr; }

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

inline

Definition at line 34 of file MillePedeFileReader.h.

References tYobs.

{ return tYobs;    }

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

inline

Definition at line 35 of file MillePedeFileReader.h.

References tYobsErr.

{ return tYobsErr; }

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

inline

Definition at line 39 of file MillePedeFileReader.h.

References tZobs.

{ return tZobs;    }

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

inline

Definition at line 40 of file MillePedeFileReader.h.

References tZobsErr.

{ return tZobsErr; }

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

inline

Definition at line 27 of file MillePedeFileReader.h.

References Xobs.

{ return Xobs;     }

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

inline

Definition at line 28 of file MillePedeFileReader.h.

References XobsErr.

{ return XobsErr;  }

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

inline

Definition at line 32 of file MillePedeFileReader.h.

References Yobs.

{ return Yobs;     }

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

inline

Definition at line 33 of file MillePedeFileReader.h.

References YobsErr.

{ return YobsErr;  }

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

inline

Definition at line 37 of file MillePedeFileReader.h.

References Zobs.

{ return Zobs;     }

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

inline

Definition at line 38 of file MillePedeFileReader.h.

References ZobsErr.

{ return ZobsErr;  }

void MillePedeFileReader::read

(

)

Definition at line 120 of file MillePedeFileReader.h.

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

       {
  readMillePedeLogFile();
  readMillePedeResultFile();
}

void MillePedeFileReader::readMillePedeLogFile ( )

private

Definition at line 137 of file MillePedeFileReader.h.

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) {
          PedeSuccess = false;
          Movements   = false;
          Error       = false;
          Significant = false;
          updateDB   = false;
        }
      }
    }

  } else {
    edm::LogError("MillePedeFileReader") << "Could not read millepede log-file.";

    PedeSuccess = false;
    Movements   = false;
    Error       = false;
    Significant = false;
    updateDB   = false;
    Nrec = 0;
  }
}

void MillePedeFileReader::readMillePedeResultFile ( )

private

Definition at line 177 of file MillePedeFileReader.h.

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

{
  std::ifstream resFile;
  resFile.open(millePedeResFile_.c_str());

  if (resFile.is_open()) {
    edm::LogInfo("MillePedeFileReader") << "Reading millepede result-file";
    double Multiplier[6] = {10000.,10000.,10000.,1000000.,1000000.,1000000.};

    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*/) {
        PedeSuccess = true;

        int alignable      = std::stoi(tokens[0]);
        int alignableIndex = alignable % 10 - 1;

        double ObsMove = std::stof(tokens[3]) * Multiplier[alignableIndex];
        double ObsErr  = std::stof(tokens[4]) * Multiplier[alignableIndex];

        int det = -1;

        if (alignable >= 60 && alignable <= 69) {
          det = 2; // TPBHalfBarrel (x+)
        } else if (alignable >= 8780 && alignable <= 8789) {
          det = 3; // TPBHalfBarrel (x-)
        } else if (alignable >= 17520 && alignable <= 17529) {
          det = 4; // TPEHalfCylinder (x+,z+)
        } else if (alignable >= 22380 && alignable <= 22389) {
          det = 5; // TPEHalfCylinder (x-,z+)
        } else if (alignable >= 27260 && alignable <= 27269) {
          det = 0; // TPEHalfCylinder (x+,z-)
        } else if (alignable >= 32120 && alignable <= 32129) {
          det = 1; //TPEHalfCylinder (x-,z-)
        } else {
          continue;
        }

        if (alignableIndex == 0 && det >= 0 && det <= 5) {
          Xobs[det] = ObsMove;
          XobsErr[det] = ObsErr;
        } else if (alignableIndex == 1 && det >= 0 && det <= 5) {
          Yobs[det] = ObsMove;
          YobsErr[det] = ObsErr;
        } else if (alignableIndex == 2 && det >= 0 && det <= 5) {
          Zobs[det] = ObsMove;
          ZobsErr[det] = ObsErr;
        } else if (alignableIndex == 3 && det >= 0 && det <= 5) {
          tXobs[det] = ObsMove;
          tXobsErr[det] = ObsErr;
        } else if (alignableIndex == 4 && det >= 0 && det <= 5) {
          tYobs[det] = ObsMove;
          tYobsErr[det] = ObsErr;
        } else if (alignableIndex == 5 && det >= 0 && det <= 5) {
          tZobs[det] = ObsMove;
          tZobsErr[det] = ObsErr;
        }

        if (abs(ObsMove) > maxMoveCut_) {
          Movements   = false;
          Error       = false;
          Significant = false;
          updateDB   = false;
          HitMax      = false;
          continue;

        } else if (abs(ObsMove) > Cutoffs[alignableIndex]) {
          Movements = true;

          if (abs(ObsErr) > maxErrorCut_) {
            Error       = false;
            Significant = false;
            updateDB   = false;
            HitErrorMax = true;
            continue;
          } else {
            Error = true;
            if (abs(ObsMove/ObsErr) > sigCut_) {
              Significant = true;
            }
          }
        }
        updateDB = true;
      }
    }
  } else {
    edm::LogError("MillePedeFileReader") << "Could not read millepede result-file.";

    PedeSuccess = false;
    Movements   = false;
    Error       = false;
    Significant = false;
    updateDB   = false;
    Nrec = 0;
  }
}

bool MillePedeFileReader::storeAlignments

(

)

Definition at line 126 of file MillePedeFileReader.h.

Referenced by PCLTrackerAlProducer::finish().

                  {
  return updateDB;
}

Member Data Documentation

double MillePedeFileReader::Cutoffs[6]

private

Initial value:

= {  Xcut_,  Ycut_,  Zcut_,
                          tXcut_, tYcut_, tZcut_}

Definition at line 64 of file MillePedeFileReader.h.

bool MillePedeFileReader::Error = false

private

Definition at line 69 of file MillePedeFileReader.h.

bool MillePedeFileReader::HitErrorMax = false

private

Definition at line 73 of file MillePedeFileReader.h.

bool MillePedeFileReader::HitMax = false

private

Definition at line 72 of file MillePedeFileReader.h.

double MillePedeFileReader::maxErrorCut_

private

Definition at line 62 of file MillePedeFileReader.h.

double MillePedeFileReader::maxMoveCut_

private

Definition at line 62 of file MillePedeFileReader.h.

std::string MillePedeFileReader::millePedeLogFile_

private

Definition at line 52 of file MillePedeFileReader.h.

std::string MillePedeFileReader::millePedeResFile_

private

Definition at line 53 of file MillePedeFileReader.h.

bool MillePedeFileReader::Movements = false

private

Definition at line 68 of file MillePedeFileReader.h.

int MillePedeFileReader::Nrec = 0

private

Definition at line 75 of file MillePedeFileReader.h.

bool MillePedeFileReader::PedeSuccess = false

private

Definition at line 67 of file MillePedeFileReader.h.

double MillePedeFileReader::sigCut_

private

Definition at line 56 of file MillePedeFileReader.h.

bool MillePedeFileReader::Significant = false

private

Definition at line 70 of file MillePedeFileReader.h.

double MillePedeFileReader::tXcut_

private

Definition at line 58 of file MillePedeFileReader.h.

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

private

Definition at line 81 of file MillePedeFileReader.h.

Referenced by getTXobs().

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

private

Definition at line 82 of file MillePedeFileReader.h.

Referenced by getTXobsErr().

double MillePedeFileReader::tYcut_

private

Definition at line 59 of file MillePedeFileReader.h.

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

private

Definition at line 86 of file MillePedeFileReader.h.

Referenced by getTYobs().

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

private

Definition at line 87 of file MillePedeFileReader.h.

Referenced by getTYobsErr().

double MillePedeFileReader::tZcut_

private

Definition at line 60 of file MillePedeFileReader.h.

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

private

Definition at line 91 of file MillePedeFileReader.h.

Referenced by getTZobs().

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

private

Definition at line 92 of file MillePedeFileReader.h.

Referenced by getTZobsErr().

bool MillePedeFileReader::updateDB = false

private

Definition at line 71 of file MillePedeFileReader.h.

double MillePedeFileReader::Xcut_

private

Definition at line 58 of file MillePedeFileReader.h.

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

private

Definition at line 79 of file MillePedeFileReader.h.

Referenced by getXobs().

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

private

Definition at line 80 of file MillePedeFileReader.h.

Referenced by getXobsErr().

double MillePedeFileReader::Ycut_

private

Definition at line 59 of file MillePedeFileReader.h.

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

private

Definition at line 84 of file MillePedeFileReader.h.

Referenced by getYobs().

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

private

Definition at line 85 of file MillePedeFileReader.h.

Referenced by getYobsErr().

double MillePedeFileReader::Zcut_

private

Definition at line 60 of file MillePedeFileReader.h.

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

private

Definition at line 89 of file MillePedeFileReader.h.

Referenced by getZobs().

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

private

Definition at line 90 of file MillePedeFileReader.h.

Referenced by getZobsErr().