ReadPGInfo Class Reference

#include <ReadPGInfo.h>

Public Member Functions
char *	getId (int, int, int)
TMatrixD	giveQC (int, int, int)
TMatrixD	giveQCCal (int, int, int)
TMatrixD	giveR (int, int, int)
TMatrixD	giveSurvey (int, int, int)
	ReadPGInfo (const char *name)
	~ReadPGInfo ()

Private Attributes
TFile *	rootFile

Detailed Description

Date:

2010/02/25 11:33:32

Revision:

1.3

Author

Luca Scodellaro Luca..nosp@m.Scod.nosp@m.ellar.nosp@m.o@ce.nosp@m.rn.ch

Definition at line 16 of file ReadPGInfo.h.

Constructor & Destructor Documentation

ReadPGInfo::ReadPGInfo

(

const char *

name

)

Definition at line 321 of file ReadPGInfo.cc.

References rootFile.

                                       {
  TDirectory *dirSave = gDirectory;
  rootFile = new TFile(name);
  dirSave->cd();
}

ReadPGInfo::~ReadPGInfo

(

)

Definition at line 327 of file ReadPGInfo.cc.

References rootFile.

{ delete rootFile; }

Member Function Documentation

char * ReadPGInfo::getId	(	int	wheel,
		int	station,
		int	sector
	)

Definition at line 329 of file ReadPGInfo.cc.

References chambers, counter, position, and TOTALCHAMBERS.

Referenced by giveQC(), giveR(), and giveSurvey().

                                                          {
  for (int counter = 0; counter < TOTALCHAMBERS; ++counter) {
    if (wheel == position[counter][0] && sector == position[counter][1] && station == position[counter][2])
      return chambers[counter];
  }
  return nullptr;
}

TMatrixD ReadPGInfo::giveQC	(	int	wheel,
		int	station,
		int	sector
	)

Definition at line 371 of file ReadPGInfo.cc.

References relativeConstraints::empty, getId(), and rootFile.

Referenced by giveQCCal().

                                                              {
  TMatrixD *empty = new TMatrixD(0, 0);
  char *id = getId(wheel, station, sector);
  if (id == nullptr)
    return *empty;
  TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
  TDirectoryFile *myQC = (TDirectoryFile *)myDir->Get("QCW");
  TMatrixD *QC;
  if (myQC == nullptr) {
    QC = new TMatrixD(0, 0);
  } else {
    QC = (TMatrixD *)myQC->Get("matrix");
  }
  return *QC;
}

TMatrixD ReadPGInfo::giveQCCal	(	int	wheel,
		int	station,
		int	sector
	)

Definition at line 348 of file ReadPGInfo.cc.

References c, relativeConstraints::error, giveQC(), and edmScanValgrind::maxCount.

Referenced by DTMuonMillepede::getbcsMatrix(), DTMuonMillepede::getbqcMatrix(), DTMuonMillepede::getCcsMatrix(), and DTMuonMillepede::getCqcMatrix().

                                                                 {
  TMatrixD *mat = new TMatrixD(0, 0);
  TMatrixD qc = giveQC(wheel, station, sector);
  if (qc.GetNrows() == 0)
    return *mat;
  mat->ResizeTo(12, 2);
  int maxCount = 12;
  if (station == 4)
    maxCount = 8;
  for (int c = 0; c < maxCount; ++c) {
    float error;
    if (qc(c, 1) == 0 || qc(c, 3) == 0) {
      (*mat)(c, 0) = (qc(c, 0) + qc(c, 2)) / 2.0;
      (*mat)(c, 1) = 500;
    } else {
      error = 1.0 / (1.0 / (qc(c, 1) * qc(c, 1)) + 1.0 / (qc(c, 3) * qc(c, 3)));
      (*mat)(c, 0) = (qc(c, 0) / (qc(c, 1) * qc(c, 1)) + qc(c, 2) / (qc(c, 3) * qc(c, 3))) * error;
      (*mat)(c, 1) = TMath::Sqrt(error);
    }
  }
  return *mat;
}

TMatrixD ReadPGInfo::giveR	(	int	wheel,
		int	station,
		int	sector
	)

Definition at line 337 of file ReadPGInfo.cc.

References relativeConstraints::empty, getId(), dttmaxenums::R, and rootFile.

                                                             {
  TMatrixD *empty = new TMatrixD(0, 0);
  char *id = getId(wheel, station, sector);
  if (id == nullptr)
    return *empty;
  TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
  TDirectoryFile *myR = (TDirectoryFile *)myDir->Get("R");
  TMatrixD *R = (TMatrixD *)myR->Get("matrix");
  return *R;
}

TMatrixD ReadPGInfo::giveSurvey	(	int	wheel,
		int	station,
		int	sector
	)

Definition at line 387 of file ReadPGInfo.cc.

References relativeConstraints::empty, getId(), and rootFile.

Referenced by DTMuonMillepede::getbcsMatrix(), and DTMuonMillepede::getbsurveyMatrix().

                                                                  {
  TMatrixD *empty = new TMatrixD(0, 0);
  char *id = getId(wheel, station, sector);
  if (id == nullptr)
    return *empty;
  TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
  TDirectoryFile *mySur = (TDirectoryFile *)myDir->Get("Survey");
  TMatrixD *Survey = (TMatrixD *)mySur->Get("matrix_layer");
  return *Survey;
}

Member Data Documentation

TFile* ReadPGInfo::rootFile

private

Definition at line 27 of file ReadPGInfo.h.

Referenced by giveQC(), giveR(), giveSurvey(), ReadPGInfo(), and ~ReadPGInfo().