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