24 |
#include <TCanvas.h> |
#include <TCanvas.h> |
25 |
#include <TStyle.h> |
#include <TStyle.h> |
26 |
#include <TH1F.h> |
#include <TH1F.h> |
27 |
|
#include <TF1.h> |
28 |
|
#include <TGraph.h> |
29 |
#include <TDecompSVD.h> |
#include <TDecompSVD.h> |
30 |
// |
// |
31 |
|
|
43 |
// |
// |
44 |
PamLevel2 *L2; |
PamLevel2 *L2; |
45 |
Bool_t debug; |
Bool_t debug; |
46 |
TFile *file; |
TFile *lfile; |
47 |
|
TFile *ffile; |
48 |
Int_t nbin; |
Int_t nbin; |
49 |
TArrayF *brig; |
TArrayF *brig; |
50 |
TArrayF *brigm; |
TArrayF *brigm; |
51 |
TArrayF *qplmean[17]; |
TArrayF *qplmean[17]; |
52 |
TMatrixD *hmat[17]; |
TMatrixD *hmat[17]; |
53 |
TMatrixF *hfmat[17]; |
TMatrixD *hfmat[17]; |
54 |
TMatrixD *fqplmean[17]; |
TMatrixD *fqplmean[17]; |
55 |
Int_t N; |
Int_t N; |
56 |
Int_t NC; |
Int_t NC; |
71 |
Bool_t dofull; |
Bool_t dofull; |
72 |
Bool_t dolong; |
Bool_t dolong; |
73 |
Int_t degfre; |
Int_t degfre; |
74 |
|
Int_t fdegfre; |
75 |
// |
// |
76 |
Float_t longtzeta; ///< longitudinal covariance parameter used for event selection |
Float_t longtzeta; ///< longitudinal covariance parameter used for event selection |
77 |
Float_t fulltzeta; ///< full calorimeter covariance parameter used for event selection |
Float_t fulltzeta; ///< full calorimeter covariance parameter used for event selection |
78 |
// |
// |
79 |
|
Float_t negfulltzeta; |
80 |
|
Float_t posfulltzeta; |
81 |
|
Float_t minsvalue; |
82 |
|
Float_t maxsvalue; |
83 |
|
Float_t aveposvar; |
84 |
|
Float_t avenegvar; |
85 |
|
Int_t numneg; |
86 |
|
Int_t numpos; |
87 |
|
// |
88 |
|
|
89 |
public: |
public: |
90 |
// |
// |
98 |
Float_t GetNormLongTZeta(); ///< Returns longitudinal covariance parameter used for event selection |
Float_t GetNormLongTZeta(); ///< Returns longitudinal covariance parameter used for event selection |
99 |
Float_t GetNormFullTZeta(); ///< Returns full calorimeter covariance parameter used for event selection |
Float_t GetNormFullTZeta(); ///< Returns full calorimeter covariance parameter used for event selection |
100 |
// |
// |
101 |
|
|
102 |
|
Float_t GetFullPositive(){Process(); return posfulltzeta;}; |
103 |
|
Float_t GetFullNegative(){Process(); return negfulltzeta;}; |
104 |
|
Int_t GetFullNumPositive(){Process(); return numpos;}; |
105 |
|
Int_t GetFullNumNegative(){Process(); return numneg;}; |
106 |
|
Float_t GetFullMaxVar(){Process(); return maxsvalue;}; |
107 |
|
Float_t GetFullMinVar(){Process(); return minsvalue;}; |
108 |
|
Float_t GetFullAveragePosVar(){Process(); return aveposvar;}; |
109 |
|
Float_t GetFullAverageNegVar(){Process(); return avenegvar;}; |
110 |
|
|
111 |
|
// |
112 |
void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT); |
void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT); |
113 |
void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used. |
void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used. |
114 |
Int_t GetDegreeOfFreedom(){return degfre;}; |
Int_t GetDegreeOfFreedom(){return degfre;}; |
115 |
|
Int_t GetLongDegreeOfFreedom(){return degfre;}; |
116 |
|
Int_t GetFullDegreeOfFreedom(){return fdegfre;}; |
117 |
// |
// |
118 |
TArrayF *LoadLongAverage(Float_t rig); |
TArrayF *LoadLongAverage(Float_t rig); |
119 |
TMatrixD *LoadCovarianceMatrix(Float_t rig); |
TMatrixD *LoadCovarianceMatrix(Float_t rig); |
125 |
Float_t GetFullHmatrixAt(Int_t i, Int_t j, Float_t rig); |
Float_t GetFullHmatrixAt(Int_t i, Int_t j, Float_t rig); |
126 |
Float_t GetFullAverageAt(Int_t plane, Int_t strip, Float_t rig, Int_t rigbin); |
Float_t GetFullAverageAt(Int_t plane, Int_t strip, Float_t rig, Int_t rigbin); |
127 |
Float_t GetFullHmatrixAt(Int_t i, Int_t j, Float_t rig, Int_t rigbin); |
Float_t GetFullHmatrixAt(Int_t i, Int_t j, Float_t rig, Int_t rigbin); |
128 |
|
Float_t GetFullHmatrixAt(Int_t i, Int_t j, Float_t rig, Int_t rigbin, Int_t mtherig); |
129 |
// |
// |
130 |
Bool_t Open(TString matrixfile); |
Bool_t Open(TString matrixfile); |
131 |
|
Bool_t Open(TString longmatrixfile,TString fullmatrixfile); |
132 |
Bool_t LoadBin(); |
Bool_t LoadBin(); |
133 |
|
Bool_t LoadBin(Bool_t full); |
134 |
Bool_t LoadLong(); |
Bool_t LoadLong(); |
135 |
Bool_t LoadFull(); |
Bool_t LoadFull(); |
136 |
Bool_t LoadMatrices(); |
Bool_t LoadMatrices(); |
137 |
Bool_t LoadFullMatrices(); |
Bool_t LoadFullMatrices(); |
138 |
TFile *GetFile(){return file;}; |
TFile *GetFile(){return lfile;}; |
139 |
|
TFile *GetLongFile(){return lfile;}; |
140 |
|
TFile *GetFullFile(){return ffile;}; |
141 |
// |
// |
142 |
TMatrixD *LoadFullAverage(Int_t rigbin); |
TMatrixD *LoadFullAverage(Int_t rigbin); |
143 |
TMatrixD *LoadFullNAverage(Int_t rigbin); |
TMatrixD *LoadFullNAverage(Int_t rigbin); |
165 |
void WriteInvertedFullMatrix(TMatrixD matrix, Int_t bin); |
void WriteInvertedFullMatrix(TMatrixD matrix, Int_t bin); |
166 |
void CloseMatrixFile(); |
void CloseMatrixFile(); |
167 |
// |
// |
168 |
void CalculateLongTZeta(){ dofull = false; dolong=true; }; |
void CalculateLongTZeta(){ dolong=true; }; |
169 |
void CalculateFullTZeta(){ dolong = false; dofull=true; }; |
void CalculateFullTZeta(){ dofull=true; }; |
170 |
|
void CalculateLongTZeta(Bool_t bo){ dolong=bo; }; |
171 |
|
void CalculateFullTZeta(Bool_t bo){ dofull=bo; }; |
172 |
void SetDebug(Bool_t d){ debug=d; }; |
void SetDebug(Bool_t d){ debug=d; }; |
173 |
void SetNoWpreSampler(Int_t n); |
void SetNoWpreSampler(Int_t n); |
174 |
void SetNoWcalo(Int_t n); |
void SetNoWcalo(Int_t n); |