#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]

double	maxErrorCut_

double	maxMoveCut_

std::string	millePedeLogFile_

std::string	millePedeResFile_

int	Nrec = 0

double	sigCut_

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

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file MillePedeFileReader.cc.

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

18 {}

Member Function Documentation

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

inline

Definition at line 25 of file MillePedeFileReader.h.

References tXobs.

25 { return tXobs; }

MillePedeFileReader::tXobs

std::array< double, 6 > tXobs

Definition: MillePedeFileReader.h:68

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

inline

Definition at line 26 of file MillePedeFileReader.h.

References tXobsErr.

26 { return tXobsErr; }

MillePedeFileReader::tXobsErr

std::array< double, 6 > tXobsErr

Definition: MillePedeFileReader.h:69

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

inline

Definition at line 30 of file MillePedeFileReader.h.

References tYobs.

30 { return tYobs; }

MillePedeFileReader::tYobs

std::array< double, 6 > tYobs

Definition: MillePedeFileReader.h:73

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

inline

Definition at line 31 of file MillePedeFileReader.h.

References tYobsErr.

31 { return tYobsErr; }

MillePedeFileReader::tYobsErr

std::array< double, 6 > tYobsErr

Definition: MillePedeFileReader.h:74

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

inline

Definition at line 35 of file MillePedeFileReader.h.

References tZobs.

35 { return tZobs; }

MillePedeFileReader::tZobs

std::array< double, 6 > tZobs

Definition: MillePedeFileReader.h:78

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

inline

Definition at line 36 of file MillePedeFileReader.h.

References tZobsErr.

36 { return tZobsErr; }

MillePedeFileReader::tZobsErr

std::array< double, 6 > tZobsErr

Definition: MillePedeFileReader.h:79

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

inline

Definition at line 23 of file MillePedeFileReader.h.

References Xobs.

23 { return Xobs; }

MillePedeFileReader::Xobs

std::array< double, 6 > Xobs

Definition: MillePedeFileReader.h:66

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

inline

Definition at line 24 of file MillePedeFileReader.h.

References XobsErr.

24 { return XobsErr; }

MillePedeFileReader::XobsErr

std::array< double, 6 > XobsErr

Definition: MillePedeFileReader.h:67

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

inline

Definition at line 28 of file MillePedeFileReader.h.

References Yobs.

28 { return Yobs; }

MillePedeFileReader::Yobs

std::array< double, 6 > Yobs

Definition: MillePedeFileReader.h:71

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

inline

Definition at line 29 of file MillePedeFileReader.h.

References YobsErr.

29 { return YobsErr; }

MillePedeFileReader::YobsErr

std::array< double, 6 > YobsErr

Definition: MillePedeFileReader.h:72

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

inline

Definition at line 33 of file MillePedeFileReader.h.

References Zobs.

33 { return Zobs; }

MillePedeFileReader::Zobs

std::array< double, 6 > Zobs

Definition: MillePedeFileReader.h:76

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

inline

Definition at line 34 of file MillePedeFileReader.h.

References ZobsErr.

34 { return ZobsErr; }

MillePedeFileReader::ZobsErr

std::array< double, 6 > ZobsErr

Definition: MillePedeFileReader.h:77

void MillePedeFileReader::read ( )

Definition at line 35 of file MillePedeFileReader.cc.

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 52 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 84 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";
     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*/) {
 
         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 (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 41 of file MillePedeFileReader.cc.

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

double MillePedeFileReader::maxErrorCut_

private

Definition at line 58 of file MillePedeFileReader.h.

double MillePedeFileReader::maxMoveCut_

private

Definition at line 58 of file MillePedeFileReader.h.

std::string MillePedeFileReader::millePedeLogFile_

private

Definition at line 48 of file MillePedeFileReader.h.

std::string MillePedeFileReader::millePedeResFile_

private

Definition at line 49 of file MillePedeFileReader.h.

int MillePedeFileReader::Nrec = 0

private

Definition at line 64 of file MillePedeFileReader.h.

double MillePedeFileReader::sigCut_

private

Definition at line 52 of file MillePedeFileReader.h.

double MillePedeFileReader::tXcut_

private

Definition at line 54 of file MillePedeFileReader.h.

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

private

Definition at line 68 of file MillePedeFileReader.h.

Referenced by getTXobs().

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

private

Definition at line 69 of file MillePedeFileReader.h.

Referenced by getTXobsErr().

double MillePedeFileReader::tYcut_

private

Definition at line 55 of file MillePedeFileReader.h.

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

private

Definition at line 73 of file MillePedeFileReader.h.

Referenced by getTYobs().

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

private

Definition at line 74 of file MillePedeFileReader.h.

Referenced by getTYobsErr().

double MillePedeFileReader::tZcut_

private

Definition at line 56 of file MillePedeFileReader.h.

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

private

Definition at line 78 of file MillePedeFileReader.h.

Referenced by getTZobs().

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

private

Definition at line 79 of file MillePedeFileReader.h.

Referenced by getTZobsErr().

bool MillePedeFileReader::updateDB = false

private

Definition at line 63 of file MillePedeFileReader.h.

double MillePedeFileReader::Xcut_

private

Definition at line 54 of file MillePedeFileReader.h.

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

private

Definition at line 66 of file MillePedeFileReader.h.

Referenced by getXobs().

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

private

Definition at line 67 of file MillePedeFileReader.h.

Referenced by getXobsErr().

double MillePedeFileReader::Ycut_

private

Definition at line 55 of file MillePedeFileReader.h.

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

private

Definition at line 71 of file MillePedeFileReader.h.

Referenced by getYobs().

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

private

Definition at line 72 of file MillePedeFileReader.h.

Referenced by getYobsErr().

double MillePedeFileReader::Zcut_

private

Definition at line 56 of file MillePedeFileReader.h.

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

private

Definition at line 76 of file MillePedeFileReader.h.

Referenced by getZobs().

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

private

Definition at line 77 of file MillePedeFileReader.h.

Referenced by getZobsErr().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation