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