00001
00011 #include <iostream>
00012 #include "TH1F.h"
00013 #include "TF1.h"
00014 #include "TCanvas.h"
00015 #include "TPaveText.h"
00016 #include "NueAna/Shwfit.h"
00017 #include "AnalysisNtuples/ANtpDefaultValue.h"
00018 #include "MessageService/MsgService.h"
00019
00020 CVSID("$Id: Shwfit.cxx,v 1.21 2008/07/28 17:56:42 boehm Exp $");
00021
00022 ClassImp(Shwfit)
00023
00024 Shwfit::Shwfit():
00025 par_a(ANtpDefaultValue::kFloat),
00026 par_b(ANtpDefaultValue::kFloat),
00027 par_e0(ANtpDefaultValue::kFloat),
00028 chisq(ANtpDefaultValue::kFloat),
00029 shwmax(ANtpDefaultValue::kFloat),
00030 shwmaxplane(ANtpDefaultValue::kInt),
00031 conv(ANtpDefaultValue::kInt),
00032 chisq_ndf(ANtpDefaultValue::kFloat),
00033 e0_pe_ratio(ANtpDefaultValue::kFloat),
00034 caldet_comp(ANtpDefaultValue::kFloat),
00035 max_pe_plane(ANtpDefaultValue::kFloat),
00036 shwmaxplane_diff(ANtpDefaultValue::kFloat),
00037 hiPhStripCountM4(ANtpDefaultValue::kInt),
00038 hiPhPlaneCountM4(ANtpDefaultValue::kInt),
00039 hiPhStripCountM2(ANtpDefaultValue::kInt),
00040 hiPhPlaneCountM2(ANtpDefaultValue::kInt),
00041 hiPhStripCount(ANtpDefaultValue::kInt),
00042 hiPhPlaneCount(ANtpDefaultValue::kInt),
00043 hiPhStripCountP2(ANtpDefaultValue::kInt),
00044 hiPhPlaneCountP2(ANtpDefaultValue::kInt),
00045 hiPhStripCountP4(ANtpDefaultValue::kInt),
00046 hiPhPlaneCountP4(ANtpDefaultValue::kInt),
00047
00048 Beta_Maxwell(ANtpDefaultValue::kFloat),
00049 Energy_Maxwell(ANtpDefaultValue::kFloat),
00050 chisq_Maxwell(ANtpDefaultValue::kFloat),
00051 ndf_Maxwell(ANtpDefaultValue::kFloat),
00052 Beta_Maxwell3(ANtpDefaultValue::kFloat),
00053 Energy_Maxwell3(ANtpDefaultValue::kFloat),
00054 chisq_Maxwell3(ANtpDefaultValue::kFloat),
00055 ndf_Maxwell3(ANtpDefaultValue::kFloat),
00056 trans_u_mean(ANtpDefaultValue::kFloat),
00057 trans_u_sigma(ANtpDefaultValue::kFloat),
00058 trans_u_chisq(ANtpDefaultValue::kFloat),
00059 trans_u_ndf(ANtpDefaultValue::kFloat),
00060 trans_v_mean(ANtpDefaultValue::kFloat),
00061 trans_v_sigma(ANtpDefaultValue::kFloat),
00062 trans_v_chisq(ANtpDefaultValue::kFloat),
00063 trans_v_ndf(ANtpDefaultValue::kFloat),
00064 E_ratio_half(ANtpDefaultValue::kFloat),
00065 n_ratio_half(ANtpDefaultValue::kFloat),
00066 E_ratio_2(ANtpDefaultValue::kFloat),
00067 n_ratio_2(ANtpDefaultValue::kFloat),
00068 pos_E_split(ANtpDefaultValue::kFloat),
00069 pos_n_split(ANtpDefaultValue::kFloat),
00070
00071
00072 contPlaneCount(ANtpDefaultValue::kInt),
00073 contPlaneCount015(ANtpDefaultValue::kInt),
00074 contPlaneCount030(ANtpDefaultValue::kInt),
00075 contPlaneCount050(ANtpDefaultValue::kInt),
00076 contPlaneCount075(ANtpDefaultValue::kInt),
00077 contPlaneCount100(ANtpDefaultValue::kInt),
00078 contPlaneCount200(ANtpDefaultValue::kInt),
00079
00080 u_asym_peak(ANtpDefaultValue::kFloat),
00081 u_asym_vert(ANtpDefaultValue::kFloat),
00082 u_molrad_peak(ANtpDefaultValue::kFloat),
00083 u_molrad_vert(ANtpDefaultValue::kFloat),
00084 u_mean(ANtpDefaultValue::kFloat),
00085 u_rms(ANtpDefaultValue::kFloat),
00086 u_skew(ANtpDefaultValue::kFloat),
00087 u_kurt(ANtpDefaultValue::kFloat),
00088 v_asym_peak(ANtpDefaultValue::kFloat),
00089 v_asym_vert(ANtpDefaultValue::kFloat),
00090 v_molrad_peak(ANtpDefaultValue::kFloat),
00091 v_molrad_vert(ANtpDefaultValue::kFloat),
00092 v_mean(ANtpDefaultValue::kFloat),
00093 v_rms(ANtpDefaultValue::kFloat),
00094 v_skew(ANtpDefaultValue::kFloat),
00095 v_kurt(ANtpDefaultValue::kFloat),
00096 uv_asym_peak(ANtpDefaultValue::kFloat),
00097 uv_asym_vert(ANtpDefaultValue::kFloat),
00098 uv_molrad_peak(ANtpDefaultValue::kFloat),
00099 uv_molrad_vert(ANtpDefaultValue::kFloat),
00100 uv_mean(ANtpDefaultValue::kFloat),
00101 uv_rms(ANtpDefaultValue::kFloat),
00102 uv_skew(ANtpDefaultValue::kFloat),
00103 uv_kurt(ANtpDefaultValue::kFloat),
00104 uv_ratio(ANtpDefaultValue::kFloat),
00105 vtxEnergy(ANtpDefaultValue::kFloat),
00106 energyPlane0(ANtpDefaultValue::kFloat),
00107 energyPlane1(ANtpDefaultValue::kFloat),
00108 energyPlane2(ANtpDefaultValue::kFloat),
00109 USlope(ANtpDefaultValue::kFloat),
00110 VSlope(ANtpDefaultValue::kFloat),
00111 UOffset(ANtpDefaultValue::kFloat),
00112 VOffset(ANtpDefaultValue::kFloat),
00113 UFitquality(ANtpDefaultValue::kFloat),
00114 VFitquality(ANtpDefaultValue::kFloat),
00115 UWDiff(ANtpDefaultValue::kFloat),
00116 VWDiff(ANtpDefaultValue::kFloat),
00117 USlopeMom(ANtpDefaultValue::kFloat),
00118 VSlopeMom(ANtpDefaultValue::kFloat),
00119 slopefix(ANtpDefaultValue::kFloat),
00120 UBeamLike(ANtpDefaultValue::kFloat),
00121 VBeamLike(ANtpDefaultValue::kFloat),
00122 UVSlope(ANtpDefaultValue::kFloat),
00123 ULongE(ANtpDefaultValue::kFloat),
00124 VLongE(ANtpDefaultValue::kFloat),
00125 LongE(ANtpDefaultValue::kFloat),
00126 complexity(ANtpDefaultValue::kFloat),
00127 wcomplexity(ANtpDefaultValue::kFloat),
00128 u_asym_peak_9s_2pe(ANtpDefaultValue::kFloat),
00129 u_asym_vert_9s_2pe(ANtpDefaultValue::kFloat),
00130 u_molrad_peak_9s_2pe(ANtpDefaultValue::kFloat),
00131 u_molrad_vert_9s_2pe(ANtpDefaultValue::kFloat),
00132 u_mean_9s_2pe(ANtpDefaultValue::kFloat),
00133 u_rms_9s_2pe(ANtpDefaultValue::kFloat),
00134 u_skew_9s_2pe(ANtpDefaultValue::kFloat),
00135 u_kurt_9s_2pe(ANtpDefaultValue::kFloat),
00136 v_asym_peak_9s_2pe(ANtpDefaultValue::kFloat),
00137 v_asym_vert_9s_2pe(ANtpDefaultValue::kFloat),
00138 v_molrad_peak_9s_2pe(ANtpDefaultValue::kFloat),
00139 v_molrad_vert_9s_2pe(ANtpDefaultValue::kFloat),
00140 v_mean_9s_2pe(ANtpDefaultValue::kFloat),
00141 v_rms_9s_2pe(ANtpDefaultValue::kFloat),
00142 v_skew_9s_2pe(ANtpDefaultValue::kFloat),
00143 v_kurt_9s_2pe(ANtpDefaultValue::kFloat),
00144 uv_asym_peak_9s_2pe(ANtpDefaultValue::kFloat),
00145 uv_asym_vert_9s_2pe(ANtpDefaultValue::kFloat),
00146 uv_molrad_peak_9s_2pe(ANtpDefaultValue::kFloat),
00147 uv_molrad_vert_9s_2pe(ANtpDefaultValue::kFloat),
00148 uv_mean_9s_2pe(ANtpDefaultValue::kFloat),
00149 uv_rms_9s_2pe(ANtpDefaultValue::kFloat),
00150 uv_skew_9s_2pe(ANtpDefaultValue::kFloat),
00151 uv_kurt_9s_2pe(ANtpDefaultValue::kFloat),
00152 uv_ratio_9s_2pe(ANtpDefaultValue::kFloat),
00153
00154 u_asym_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00155 u_asym_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00156 u_molrad_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00157 u_molrad_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00158 u_mean_9s_2pe_dw(ANtpDefaultValue::kFloat),
00159 u_rms_9s_2pe_dw(ANtpDefaultValue::kFloat),
00160 u_skew_9s_2pe_dw(ANtpDefaultValue::kFloat),
00161 u_kurt_9s_2pe_dw(ANtpDefaultValue::kFloat),
00162 v_asym_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00163 v_asym_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00164 v_molrad_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00165 v_molrad_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00166 v_mean_9s_2pe_dw(ANtpDefaultValue::kFloat),
00167 v_rms_9s_2pe_dw(ANtpDefaultValue::kFloat),
00168 v_skew_9s_2pe_dw(ANtpDefaultValue::kFloat),
00169 v_kurt_9s_2pe_dw(ANtpDefaultValue::kFloat),
00170 uv_asym_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00171 uv_asym_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00172 uv_molrad_peak_9s_2pe_dw(ANtpDefaultValue::kFloat),
00173 uv_molrad_vert_9s_2pe_dw(ANtpDefaultValue::kFloat),
00174 uv_mean_9s_2pe_dw(ANtpDefaultValue::kFloat),
00175 uv_rms_9s_2pe_dw(ANtpDefaultValue::kFloat),
00176 uv_skew_9s_2pe_dw(ANtpDefaultValue::kFloat),
00177 uv_kurt_9s_2pe_dw(ANtpDefaultValue::kFloat),
00178 uv_ratio_9s_2pe_dw(ANtpDefaultValue::kFloat),
00179
00180 lenepl(0),
00181 ph_hist(0),
00182 tenestu(0),
00183 tenestv(0),
00184 tenestu_9s_2pe(0),
00185 tenestv_9s_2pe(0),
00186 tenestu_9s_2pe_dw(0),
00187 tenestv_9s_2pe_dw(0),
00188 efit(0),
00189 hfit(0),
00190 efit_maxwell(0),
00191 efit_maxwell3(0),
00192 ufit(0),
00193 vfit(0),
00194 info1(0)
00195 {}
00196
00197 Shwfit::Shwfit(const Shwfit &s):
00198 TObject(),
00199 par_a(s.par_a),
00200 par_b(s.par_b),
00201 par_e0(s.par_e0),
00202 chisq(s.chisq),
00203 shwmax(s.shwmax),
00204 shwmaxplane(s.shwmaxplane),
00205 conv(s.conv),
00206 chisq_ndf(s.chisq_ndf),
00207 e0_pe_ratio(s.e0_pe_ratio),
00208 caldet_comp(s.caldet_comp),
00209 max_pe_plane(s.max_pe_plane),
00210 shwmaxplane_diff(s.shwmaxplane_diff),
00211 hiPhStripCountM4(s.hiPhStripCountM4),
00212 hiPhPlaneCountM4(s.hiPhPlaneCountM4),
00213 hiPhStripCountM2(s.hiPhStripCountM2),
00214 hiPhPlaneCountM2(s.hiPhPlaneCountM2),
00215 hiPhStripCount(s.hiPhStripCount),
00216 hiPhPlaneCount(s.hiPhPlaneCount),
00217 hiPhStripCountP2(s.hiPhStripCountP2),
00218 hiPhPlaneCountP2(s.hiPhPlaneCountP2),
00219 hiPhStripCountP4(s.hiPhStripCountP4),
00220 hiPhPlaneCountP4(s.hiPhPlaneCountP4),
00221
00222 Beta_Maxwell(s.Beta_Maxwell),
00223 Energy_Maxwell(s.Energy_Maxwell),
00224 chisq_Maxwell(s.chisq_Maxwell),
00225 ndf_Maxwell(s.ndf_Maxwell),
00226 Beta_Maxwell3(s.Beta_Maxwell3),
00227 Energy_Maxwell3(s.Energy_Maxwell3),
00228 chisq_Maxwell3(s.chisq_Maxwell3),
00229 ndf_Maxwell3(s.ndf_Maxwell3),
00230 trans_u_mean(s.trans_u_mean),
00231 trans_u_sigma(s.trans_u_sigma),
00232 trans_u_chisq(s.trans_u_chisq),
00233 trans_u_ndf(s.trans_u_ndf),
00234 trans_v_mean(s.trans_v_mean),
00235 trans_v_sigma(s.trans_v_sigma),
00236 trans_v_chisq(s.trans_v_chisq),
00237 trans_v_ndf(s.trans_v_ndf),
00238 E_ratio_half(s.E_ratio_half),
00239 n_ratio_half(s.n_ratio_half),
00240 E_ratio_2(s.E_ratio_2),
00241 n_ratio_2(s.n_ratio_2),
00242 pos_E_split(s.pos_E_split),
00243 pos_n_split(s.pos_n_split),
00244
00245
00246 contPlaneCount(s.contPlaneCount),
00247 contPlaneCount015(s.contPlaneCount015),
00248 contPlaneCount030(s.contPlaneCount030),
00249 contPlaneCount050(s.contPlaneCount050),
00250 contPlaneCount075(s.contPlaneCount075),
00251 contPlaneCount100(s.contPlaneCount100),
00252 contPlaneCount200(s.contPlaneCount200),
00253
00254 u_asym_peak(s.u_asym_peak),
00255 u_asym_vert(s.u_asym_vert),
00256 u_molrad_peak(s.u_molrad_peak),
00257 u_molrad_vert(s.u_molrad_vert),
00258 u_mean(s.u_mean),
00259 u_rms(s.u_rms),
00260 u_skew(s.u_skew),
00261 u_kurt(s.u_kurt),
00262 v_asym_peak(s.v_asym_peak),
00263 v_asym_vert(s.v_asym_vert),
00264 v_molrad_peak(s.v_molrad_peak),
00265 v_molrad_vert(s.v_molrad_vert),
00266 v_mean(s.v_mean),
00267 v_rms(s.v_rms),
00268 v_skew(s.v_skew),
00269 v_kurt(s.v_kurt),
00270 uv_asym_peak(s.uv_asym_peak),
00271 uv_asym_vert(s.uv_asym_vert),
00272 uv_molrad_peak(s.uv_molrad_peak),
00273 uv_molrad_vert(s.uv_molrad_vert),
00274 uv_mean(s.uv_mean),
00275 uv_rms(s.uv_rms),
00276 uv_skew(s.uv_skew),
00277 uv_kurt(s.uv_kurt),
00278 uv_ratio(s.uv_ratio),
00279 vtxEnergy(s.vtxEnergy),
00280 energyPlane0(s.energyPlane0),
00281 energyPlane1(s.energyPlane1),
00282 energyPlane2(s.energyPlane2),
00283 USlope(s.USlope),
00284 VSlope(s.VSlope),
00285 UOffset(s.UOffset),
00286 VOffset(s.VOffset),
00287 UFitquality(s.UFitquality),
00288 VFitquality(s.VFitquality),
00289 UWDiff(s.UWDiff),
00290 VWDiff(s.VWDiff),
00291 USlopeMom(s.USlopeMom),
00292 VSlopeMom(s.VSlopeMom),
00293 slopefix(s.slopefix),
00294 UBeamLike(s.UBeamLike),
00295 VBeamLike(s.VBeamLike),
00296 UVSlope(s.UVSlope),
00297 ULongE(s.ULongE),
00298 VLongE(s.VLongE),
00299 LongE(s.LongE),
00300 complexity(s.complexity),
00301 wcomplexity(s.wcomplexity),
00302 u_asym_peak_9s_2pe(s.u_asym_peak_9s_2pe),
00303 u_asym_vert_9s_2pe(s.u_asym_vert_9s_2pe),
00304 u_molrad_peak_9s_2pe(s.u_molrad_peak_9s_2pe),
00305 u_molrad_vert_9s_2pe(s.u_molrad_vert_9s_2pe),
00306 u_mean_9s_2pe(s.u_mean_9s_2pe),
00307 u_rms_9s_2pe(s.u_rms_9s_2pe),
00308 u_skew_9s_2pe(s.u_skew_9s_2pe),
00309 u_kurt_9s_2pe(s.u_kurt_9s_2pe),
00310 v_asym_peak_9s_2pe(s.v_asym_peak_9s_2pe),
00311 v_asym_vert_9s_2pe(s.v_asym_vert_9s_2pe),
00312 v_molrad_peak_9s_2pe(s.v_molrad_peak_9s_2pe),
00313 v_molrad_vert_9s_2pe(s.v_molrad_vert_9s_2pe),
00314 v_mean_9s_2pe(s.v_mean_9s_2pe),
00315 v_rms_9s_2pe(s.v_rms_9s_2pe),
00316 v_skew_9s_2pe(s.v_skew_9s_2pe),
00317 v_kurt_9s_2pe(s.v_kurt_9s_2pe),
00318 uv_asym_peak_9s_2pe(s.uv_asym_peak_9s_2pe),
00319 uv_asym_vert_9s_2pe(s.uv_asym_vert_9s_2pe),
00320 uv_molrad_peak_9s_2pe(s.uv_molrad_peak_9s_2pe),
00321 uv_molrad_vert_9s_2pe(s.uv_molrad_vert_9s_2pe),
00322 uv_mean_9s_2pe(s.uv_mean_9s_2pe),
00323 uv_rms_9s_2pe(s.uv_rms_9s_2pe),
00324 uv_skew_9s_2pe(s.uv_skew_9s_2pe),
00325 uv_kurt_9s_2pe(s.uv_kurt_9s_2pe),
00326 uv_ratio_9s_2pe(s.uv_ratio_9s_2pe),
00327
00328 u_asym_peak_9s_2pe_dw(s.u_asym_peak_9s_2pe_dw),
00329 u_asym_vert_9s_2pe_dw(s.u_asym_vert_9s_2pe_dw),
00330 u_molrad_peak_9s_2pe_dw(s.u_molrad_peak_9s_2pe_dw),
00331 u_molrad_vert_9s_2pe_dw(s.u_molrad_vert_9s_2pe_dw),
00332 u_mean_9s_2pe_dw(s.u_mean_9s_2pe_dw),
00333 u_rms_9s_2pe_dw(s.u_rms_9s_2pe_dw),
00334 u_skew_9s_2pe_dw(s.u_skew_9s_2pe_dw),
00335 u_kurt_9s_2pe_dw(s.u_kurt_9s_2pe_dw),
00336 v_asym_peak_9s_2pe_dw(s.v_asym_peak_9s_2pe_dw),
00337 v_asym_vert_9s_2pe_dw(s.v_asym_vert_9s_2pe_dw),
00338 v_molrad_peak_9s_2pe_dw(s.v_molrad_peak_9s_2pe_dw),
00339 v_molrad_vert_9s_2pe_dw(s.v_molrad_vert_9s_2pe_dw),
00340 v_mean_9s_2pe_dw(s.v_mean_9s_2pe_dw),
00341 v_rms_9s_2pe_dw(s.v_rms_9s_2pe_dw),
00342 v_skew_9s_2pe_dw(s.v_skew_9s_2pe_dw),
00343 v_kurt_9s_2pe_dw(s.v_kurt_9s_2pe_dw),
00344 uv_asym_peak_9s_2pe_dw(s.uv_asym_peak_9s_2pe_dw),
00345 uv_asym_vert_9s_2pe_dw(s.uv_asym_vert_9s_2pe_dw),
00346 uv_molrad_peak_9s_2pe_dw(s.uv_molrad_peak_9s_2pe_dw),
00347 uv_molrad_vert_9s_2pe_dw(s.uv_molrad_vert_9s_2pe_dw),
00348 uv_mean_9s_2pe_dw(s.uv_mean_9s_2pe_dw),
00349 uv_rms_9s_2pe_dw(s.uv_rms_9s_2pe_dw),
00350 uv_skew_9s_2pe_dw(s.uv_skew_9s_2pe_dw),
00351 uv_kurt_9s_2pe_dw(s.uv_kurt_9s_2pe_dw),
00352 uv_ratio_9s_2pe_dw(s.uv_ratio_9s_2pe_dw)
00353
00354 {
00355 if(s.lenepl!=0){
00356 lenepl=new TH1F(*(s.lenepl));
00357 }
00358 else{
00359 lenepl=0;
00360 }
00361 if(s.tenestu!=0){
00362 tenestu=new TH1F(*(s.tenestu));
00363 }
00364 else{
00365 tenestu=0;
00366 }
00367 if(s.tenestv!=0){
00368 tenestv=new TH1F(*(s.tenestv));
00369 }
00370 else{
00371 tenestv=0;
00372 }
00373 if(s.efit!=0){
00374 efit=new TF1(*(s.efit));
00375 }
00376 else{
00377 efit=0;
00378 }
00379 if(s.hfit!=0){
00380 hfit=new TF1(*(s.hfit));
00381 }
00382 else{
00383 hfit=0;
00384 }
00385 if(s.info1!=0){
00386 info1=new TPaveText(*(s.info1));
00387 }
00388 else{
00389 info1=0;
00390 }
00391
00392 if(s.efit_maxwell!=0) { efit_maxwell=new TF1(*(s.efit_maxwell)); }
00393 else{efit_maxwell=0; }
00394 if(s.efit_maxwell3!=0){ efit_maxwell3=new TF1(*(s.efit_maxwell3));}
00395 else{ efit_maxwell3=0; }
00396 if(s.ufit!=0){ ufit=new TF1(*(s.ufit)); }
00397 else{ ufit=0; }
00398 if(s.vfit!=0){ vfit=new TF1(*(s.vfit)); }
00399 else{ vfit=0; }
00400
00401 if(s.tenestv_9s_2pe!=0){ tenestv_9s_2pe=new TH1F(*(s.tenestv_9s_2pe));}
00402 else{ tenestv_9s_2pe=0; }
00403 if(s.tenestv_9s_2pe_dw!=0){
00404 tenestv_9s_2pe_dw=new TH1F(*(s.tenestv_9s_2pe_dw));}
00405 else{ tenestv_9s_2pe_dw=0; }
00406
00407 if(s.tenestu_9s_2pe!=0){ tenestu_9s_2pe=new TH1F(*(s.tenestu_9s_2pe));}
00408 else{ tenestu_9s_2pe=0; }
00409 if(s.tenestu_9s_2pe_dw!=0){
00410 tenestu_9s_2pe_dw=new TH1F(*(s.tenestu_9s_2pe_dw));}
00411 else{ tenestu_9s_2pe_dw=0; }
00412
00413 if(s.ph_hist!=0){
00414 ph_hist = new TH1F(*(s.ph_hist));
00415 }
00416 else{ ph_hist=0; }
00417 }
00418
00419
00420 Shwfit::~Shwfit()
00421 {
00422 if(lenepl!=0){
00423 delete lenepl;
00424 lenepl=0;
00425 }
00426 if(tenestu!=0){
00427 delete tenestu;
00428 tenestu=0;
00429 }
00430 if(tenestv!=0){
00431 delete tenestv;
00432 tenestv=0;
00433 }
00434
00435 if(efit!=0){
00436 delete efit;
00437 efit=0;
00438 }
00439 if(hfit!=0){
00440 delete hfit;
00441 hfit=0;
00442 }
00443 if(info1!=0){
00444 delete info1;
00445 info1=0;
00446 }
00447
00448 if(efit_maxwell!=0){ delete efit_maxwell; efit_maxwell=0; }
00449 if(efit_maxwell3!=0){ delete efit_maxwell3; efit_maxwell3=0; }
00450 if(ufit!=0){ delete ufit; ufit=0; }
00451 if(vfit!=0){ delete vfit; vfit=0; }
00452
00453 if(tenestv_9s_2pe!=0) { delete tenestv_9s_2pe; tenestv_9s_2pe=0;}
00454 if(tenestv_9s_2pe_dw!=0){ delete tenestv_9s_2pe_dw; tenestv_9s_2pe_dw=0;}
00455 if(tenestu_9s_2pe!=0) { delete tenestu_9s_2pe; tenestu_9s_2pe=0;}
00456 if(tenestu_9s_2pe_dw!=0){ delete tenestu_9s_2pe_dw; tenestu_9s_2pe_dw=0;}
00457
00458 if(ph_hist !=0) { delete ph_hist; ph_hist = 0; }
00459
00460 }
00461
00462
00463 void Shwfit::Draw(TPad *pad)
00464 {
00465
00466 pad->cd(1);
00467 Print("test");
00468
00469
00470 if(info1) info1->Delete();
00471 info1= new TPaveText(.05,.05,0.95,0.95);
00472
00473 info1->Clear();
00474 info1->SetFillColor(kWhite);
00475 info1->SetBorderSize(1);
00476 info1->SetTextAlign(13);
00477 info1->SetTextFont(42);
00478 char temp[100];
00479
00480 sprintf(temp," ");
00481 info1->AddText(temp);
00482
00483 sprintf(temp," Long Fit:");
00484 info1->AddText(temp);
00485 sprintf(temp," a = %2.2f, b = %2.2f, E_0 = %2.2f ",par_a,par_b,par_e0);
00486 info1->AddText(temp);
00487 sprintf(temp," chisq/ndf = %2.2f, E0/MEU ratio = %2.2f, Conv = %d, ContPlane = %d ",chisq_ndf,e0_pe_ratio,conv,contPlaneCount);
00488 info1->AddText(temp);
00489
00490 sprintf(temp,"U Trans: U RMS = %2.2f, U kurt = %2.2f, U MolRad = %2.2f", u_rms, u_kurt, u_molrad_vert);
00491 info1->AddText(temp);
00492
00493 sprintf(temp,"V Trans: V RMS = %2.2f, V kurt = %2.2f, V MolRad = %2.2f", v_rms, v_kurt, v_molrad_vert);
00494 info1->AddText(temp);
00495
00496 sprintf(temp,"UV Trans: UV RMS = %2.2f, UV kurt = %2.2f, UV MolRad = %2.2f", uv_rms, uv_kurt, uv_molrad_vert);
00497 info1->AddText(temp);
00498
00499 info1->Draw();
00500
00501 pad->cd(3);
00502 lenepl->Draw();
00503 if(conv==1) efit->Draw("sames");
00504 pad->cd(2);
00505
00506 tenestu->Draw();
00507 pad->cd(4);
00508 tenestv->Draw();
00509 pad->cd();
00510 pad->Modified();
00511
00512
00513 }
00514 void Shwfit::Draw(Option_t *)
00515 {
00517 }
00518
00519
00520 void Shwfit::Print(Option_t *) const
00521 {
00522 MSG("Shwfit::Print",Msg::kDebug) << "par_a = " << par_a
00523 << " par_b = " << par_b
00524 << " par_e0 = " << par_e0
00525 << " chisq_ndf = " << chisq_ndf
00526 << " e0_pe_ratio = " << e0_pe_ratio
00527 << " conv = " << conv
00528 << endl;
00529 }
00530
00531 void Shwfit::Clear(Option_t *)
00532 {
00533 Reset();
00534 if(lenepl!=0) { delete lenepl; lenepl=0; }
00535 if(tenestu!=0){ delete tenestu; tenestu=0; }
00536 if(tenestv!=0){ delete tenestv; tenestv=0; }
00537 if(efit!=0){ delete efit; efit=0; }
00538 if(hfit!=0){ delete hfit; hfit=0; }
00539 if(info1!=0){ delete info1; info1=0; }
00540 if(tenestv_9s_2pe!=0) { delete tenestv_9s_2pe; tenestv_9s_2pe=0;}
00541 if(tenestv_9s_2pe_dw!=0){ delete tenestv_9s_2pe_dw; tenestv_9s_2pe_dw=0;}
00542 if(tenestu_9s_2pe!=0) { delete tenestu_9s_2pe; tenestu_9s_2pe=0;}
00543 if(tenestu_9s_2pe_dw!=0){ delete tenestu_9s_2pe_dw; tenestu_9s_2pe_dw=0;}
00544
00545 if(efit_maxwell!=0){ delete efit_maxwell; efit_maxwell=0; }
00546 if(efit_maxwell3!=0){ delete efit_maxwell3; efit_maxwell3=0; }
00547 if(ufit!=0){ delete ufit; ufit=0; }
00548 if(vfit!=0){ delete vfit; vfit=0; }
00549 if(ph_hist !=0) { delete ph_hist; ph_hist = 0; }
00550
00551 }
00552
00553 void Shwfit::Reset()
00554 {
00555 par_a= ANtpDefaultValue::kFloat;
00556 par_b= ANtpDefaultValue::kFloat;
00557 par_e0=ANtpDefaultValue::kFloat;
00558 chisq= ANtpDefaultValue::kFloat;
00559 shwmax=ANtpDefaultValue::kFloat;
00560 shwmaxplane= ANtpDefaultValue::kInt;
00561 conv= ANtpDefaultValue::kInt;
00562 chisq_ndf=ANtpDefaultValue::kFloat;
00563 e0_pe_ratio=ANtpDefaultValue::kFloat;
00564 caldet_comp=ANtpDefaultValue::kFloat;
00565 max_pe_plane=ANtpDefaultValue::kFloat;
00566 shwmaxplane_diff=ANtpDefaultValue::kFloat;
00567
00568 Beta_Maxwell= ANtpDefaultValue::kFloat;
00569 Energy_Maxwell= ANtpDefaultValue::kFloat;
00570 chisq_Maxwell= ANtpDefaultValue::kFloat;
00571 ndf_Maxwell= ANtpDefaultValue::kFloat;
00572 Beta_Maxwell3= ANtpDefaultValue::kFloat;
00573 Energy_Maxwell3= ANtpDefaultValue::kFloat;
00574 chisq_Maxwell3= ANtpDefaultValue::kFloat;
00575 ndf_Maxwell3= ANtpDefaultValue::kFloat;
00576 trans_u_mean= ANtpDefaultValue::kFloat;
00577 trans_u_sigma= ANtpDefaultValue::kFloat;
00578 trans_u_chisq= ANtpDefaultValue::kFloat;
00579 trans_u_ndf= ANtpDefaultValue::kFloat;
00580 trans_v_mean= ANtpDefaultValue::kFloat;
00581 trans_v_sigma= ANtpDefaultValue::kFloat;
00582 trans_v_chisq= ANtpDefaultValue::kFloat;
00583 trans_v_ndf= ANtpDefaultValue::kFloat;
00584 E_ratio_half= ANtpDefaultValue::kFloat;
00585 n_ratio_half= ANtpDefaultValue::kFloat;
00586 E_ratio_2= ANtpDefaultValue::kFloat;
00587 n_ratio_2= ANtpDefaultValue::kFloat;
00588 pos_E_split= ANtpDefaultValue::kFloat;
00589 pos_n_split= ANtpDefaultValue::kFloat;
00590
00591 hiPhStripCountM4=ANtpDefaultValue::kInt;
00592 hiPhPlaneCountM4=ANtpDefaultValue::kInt;
00593 hiPhStripCountM2=ANtpDefaultValue::kInt;
00594 hiPhPlaneCountM2=ANtpDefaultValue::kInt;
00595 hiPhStripCount=ANtpDefaultValue::kInt;
00596 hiPhPlaneCount=ANtpDefaultValue::kInt;
00597 hiPhStripCountP2=ANtpDefaultValue::kInt;
00598 hiPhPlaneCountP2=ANtpDefaultValue::kInt;
00599 hiPhStripCountP4=ANtpDefaultValue::kInt;
00600 hiPhPlaneCountP4=ANtpDefaultValue::kInt;
00601 contPlaneCount=ANtpDefaultValue::kInt;
00602 contPlaneCount015=ANtpDefaultValue::kInt;
00603 contPlaneCount030=ANtpDefaultValue::kInt;
00604 contPlaneCount050=ANtpDefaultValue::kInt;
00605 contPlaneCount075=ANtpDefaultValue::kInt;
00606 contPlaneCount100=ANtpDefaultValue::kInt;
00607 contPlaneCount200=ANtpDefaultValue::kInt;
00608 u_asym_peak=ANtpDefaultValue::kFloat;
00609 u_asym_vert=ANtpDefaultValue::kFloat;
00610 u_molrad_peak=ANtpDefaultValue::kFloat;
00611 u_molrad_vert=ANtpDefaultValue::kFloat;
00612 u_mean=ANtpDefaultValue::kFloat;
00613 u_rms= ANtpDefaultValue::kFloat;
00614 u_skew=ANtpDefaultValue::kFloat;
00615 u_kurt=ANtpDefaultValue::kFloat;
00616 v_asym_peak=ANtpDefaultValue::kFloat;
00617 v_asym_vert=ANtpDefaultValue::kFloat;
00618 v_molrad_peak=ANtpDefaultValue::kFloat;
00619 v_molrad_vert=ANtpDefaultValue::kFloat;
00620 v_mean=ANtpDefaultValue::kFloat;
00621 v_rms= ANtpDefaultValue::kFloat;
00622 v_skew=ANtpDefaultValue::kFloat;
00623 v_kurt=ANtpDefaultValue::kFloat;
00624 uv_asym_peak=ANtpDefaultValue::kFloat;
00625 uv_asym_vert=ANtpDefaultValue::kFloat;
00626 uv_molrad_peak=ANtpDefaultValue::kFloat;
00627 uv_molrad_vert=ANtpDefaultValue::kFloat;
00628 uv_mean= ANtpDefaultValue::kFloat;
00629 uv_rms = ANtpDefaultValue::kFloat;
00630 uv_skew= ANtpDefaultValue::kFloat;
00631 uv_kurt= ANtpDefaultValue::kFloat;
00632 uv_ratio=ANtpDefaultValue::kFloat;
00633 vtxEnergy=ANtpDefaultValue::kFloat;
00634 energyPlane0 = ANtpDefaultValue::kFloat;
00635 energyPlane1 = ANtpDefaultValue::kFloat;
00636 energyPlane2 = ANtpDefaultValue::kFloat;
00637 USlope = ANtpDefaultValue::kFloat;
00638 VSlope = ANtpDefaultValue::kFloat;
00639 UOffset = ANtpDefaultValue::kFloat;
00640 VOffset = ANtpDefaultValue::kFloat;
00641 UFitquality = ANtpDefaultValue::kFloat;
00642 VFitquality = ANtpDefaultValue::kFloat;
00643 UWDiff = ANtpDefaultValue::kFloat;
00644 VWDiff = ANtpDefaultValue::kFloat;
00645 USlopeMom = ANtpDefaultValue::kFloat;
00646 VSlopeMom = ANtpDefaultValue::kFloat;
00647 slopefix = ANtpDefaultValue::kFloat;
00648 UBeamLike = ANtpDefaultValue::kFloat;
00649 VBeamLike = ANtpDefaultValue::kFloat;
00650 UVSlope = ANtpDefaultValue::kFloat;
00651 ULongE = ANtpDefaultValue::kFloat;
00652 VLongE = ANtpDefaultValue::kFloat;
00653 LongE = ANtpDefaultValue::kFloat;
00654 complexity = ANtpDefaultValue::kFloat;
00655 wcomplexity = ANtpDefaultValue::kFloat;
00656
00657 u_asym_peak_9s_2pe = ANtpDefaultValue::kFloat;
00658 u_asym_vert_9s_2pe = ANtpDefaultValue::kFloat;
00659 u_molrad_peak_9s_2pe = ANtpDefaultValue::kFloat;
00660 u_molrad_vert_9s_2pe = ANtpDefaultValue::kFloat;
00661 u_mean_9s_2pe = ANtpDefaultValue::kFloat;
00662 u_rms_9s_2pe = ANtpDefaultValue::kFloat;
00663 u_skew_9s_2pe = ANtpDefaultValue::kFloat;
00664 u_kurt_9s_2pe = ANtpDefaultValue::kFloat;
00665 v_asym_peak_9s_2pe = ANtpDefaultValue::kFloat;
00666 v_asym_vert_9s_2pe = ANtpDefaultValue::kFloat;
00667 v_molrad_peak_9s_2pe = ANtpDefaultValue::kFloat;
00668 v_molrad_vert_9s_2pe = ANtpDefaultValue::kFloat;
00669 v_mean_9s_2pe = ANtpDefaultValue::kFloat;
00670 v_rms_9s_2pe = ANtpDefaultValue::kFloat;
00671 v_skew_9s_2pe = ANtpDefaultValue::kFloat;
00672 v_kurt_9s_2pe = ANtpDefaultValue::kFloat;
00673 uv_asym_peak_9s_2pe = ANtpDefaultValue::kFloat;
00674 uv_asym_vert_9s_2pe = ANtpDefaultValue::kFloat;
00675 uv_molrad_peak_9s_2pe = ANtpDefaultValue::kFloat;
00676 uv_molrad_vert_9s_2pe = ANtpDefaultValue::kFloat;
00677 uv_mean_9s_2pe = ANtpDefaultValue::kFloat;
00678 uv_rms_9s_2pe = ANtpDefaultValue::kFloat;
00679 uv_skew_9s_2pe = ANtpDefaultValue::kFloat;
00680 uv_kurt_9s_2pe = ANtpDefaultValue::kFloat;
00681 uv_ratio_9s_2pe = ANtpDefaultValue::kFloat;
00682
00683 u_asym_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00684 u_asym_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00685 u_molrad_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00686 u_molrad_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00687 u_mean_9s_2pe_dw = ANtpDefaultValue::kFloat;
00688 u_rms_9s_2pe_dw = ANtpDefaultValue::kFloat;
00689 u_skew_9s_2pe_dw = ANtpDefaultValue::kFloat;
00690 u_kurt_9s_2pe_dw = ANtpDefaultValue::kFloat;
00691 v_asym_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00692 v_asym_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00693 v_molrad_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00694 v_molrad_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00695 v_mean_9s_2pe_dw = ANtpDefaultValue::kFloat;
00696 v_rms_9s_2pe_dw = ANtpDefaultValue::kFloat;
00697 v_skew_9s_2pe_dw = ANtpDefaultValue::kFloat;
00698 v_kurt_9s_2pe_dw = ANtpDefaultValue::kFloat;
00699 uv_asym_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00700 uv_asym_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00701 uv_molrad_peak_9s_2pe_dw = ANtpDefaultValue::kFloat;
00702 uv_molrad_vert_9s_2pe_dw = ANtpDefaultValue::kFloat;
00703 uv_mean_9s_2pe_dw = ANtpDefaultValue::kFloat;
00704 uv_rms_9s_2pe_dw = ANtpDefaultValue::kFloat;
00705 uv_skew_9s_2pe_dw = ANtpDefaultValue::kFloat;
00706 uv_kurt_9s_2pe_dw = ANtpDefaultValue::kFloat;
00707 uv_ratio_9s_2pe_dw = ANtpDefaultValue::kFloat;
00708
00709
00710
00711
00712
00713
00714
00715
00716
00717 }
00718
