InvMatrixUtils.h

Go to the documentation of this file.

#ifndef InvMatrixUtils_h
#define InvMatrixUtils_h

#include <string>
#include <map>
#include "TObject.h"
#include "TF1.h"
#include "TH2.h"
#include "TProfile.h"
#include "TCanvas.h"
#include "CLHEP/Geometry/Point3D.h"
#include "CLHEP/Matrix/GenMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "TFile.h"

//#include "ConfigParser.h"
#include "Calibration/Tools/interface/InvMatrixCommonDefs.h"

void setStyle();

TCanvas* getGlobalCanvas(std::string name = "Inv MatrixCanvas");

TFile* getGlobalTFile(std::string name = "Inv MatrixTFile.root");
//TFile * getGlobalTFile (const char* name) ;

int saveGlobalTFile(std::string name = "Inv MatrixFile.root");

CLHEP::HepMatrix* getSavedMatrix(const std::string& name);

HepGeom::Point3D<Float_t> TBposition(const Float_t amplit[7][7],
                                     const Float_t beamEne,
                                     const Float_t w0 = 4.0,
                                     const Float_t x0 = 8.9,  //mm
                                     const Float_t a0 = 6.2,
                                     const Float_t sideX = 24.06,   //mm
                                     const Float_t sideY = 22.02);  //mm

double get5x5(const Float_t energy[7][7]);

double get3x3(const Float_t energy[7][7]);

int parseConfigFile(const TString& config);

int xtalFromEtaPhi(const int& myEta, const int& myPhi);

int xtalFromiEtaiPhi(const int& iEta, const int& iPhi);

int etaFromXtal(const int& xtal);

int phiFromXtal(const int& xtal);

int ietaFromXtal(const int& xtal);

int iphiFromXtal(const int& xtal);

int extract(std::vector<int>* output, const std::string& dati);

int writeCalibTxt(const CLHEP::HepMatrix& AmplitudeMatrix,
                  const CLHEP::HepMatrix& SigmaMatrix,
                  const CLHEP::HepMatrix& StatisticMatrix,
                  std::string fileName = "calibOutput.txt");

int writeCMSSWCoeff(const CLHEP::HepMatrix& amplMatrix,
                    double calibThres,
                    float ERef,
                    const CLHEP::HepMatrix& sigmaMatrix,
                    const CLHEP::HepMatrix& statisticMatrix,
                    std::string fileName = "calibOutput.txt",
                    std::string genTag = "CAL_GENTAG",
                    std::string method = "CAL_METHOD",
                    std::string version = "CAL_VERSION",
                    std::string type = "CAL_TYPE");

int writeCMSSWCoeff(const CLHEP::HepMatrix& amplMatrix,
                    double calibThres,
                    int etaRef,
                    int phiRef,
                    const CLHEP::HepMatrix& sigmaMatrix,
                    const CLHEP::HepMatrix& statisticMatrix,
                    std::string fileName = "calibOutput.txt",
                    std::string genTag = "CAL_GENTAG",
                    std::string method = "CAL_METHOD",
                    std::string version = "CAL_VERSION",
                    std::string type = "CAL_TYPE");

int translateCoeff(const CLHEP::HepMatrix& calibcoeff,
                   const CLHEP::HepMatrix& sigmaMatrix,
                   const CLHEP::HepMatrix& statisticMatrix,
                   std::string SMnumber = "1",
                   double calibThres = 0.01,
                   std::string fileName = "calibOutput.txt",
                   std::string genTag = "CAL_GENTAG",
                   std::string method = "CAL_METHOD",
                   std::string version = "CAL_VERSION",
                   std::string type = "CAL_TYPE");

int readCMSSWcoeff(CLHEP::HepMatrix& calibcoeff, const std::string& inputFileName, double defaultVal = 1.);

int readCMSSWcoeffForComparison(CLHEP::HepMatrix& calibcoeff, const std::string& inputFileName);

TH1D* smartProfile(TH2F* strip, double width);

TH1D* smartGausProfile(TH2F* strip, double width);

TH1D* smartError(TH1D* strip);

double effectiveSigma(TH1F& histogram, int vSteps = 100);

std::pair<int, int> findSupport(TH1F& histogram, double thres = 0.);

void mtrTransfer(double output[SCMaxEta][SCMaxPhi], CLHEP::HepMatrix* input, double Default);

template <class Type>
void mtrReset(Type superModules[SCMaxEta][SCMaxPhi], const Type val) {
  for (int e = 0; e < SCMaxEta; ++e)
    for (int p = 0; p < SCMaxPhi; ++p) {
      superModules[e][p] = val;
    }
}

double etaCorrE1E9(int eta);
double etaCorrE1E49(int eta);
double etaCorrE1E25(int eta);

#endif