diff --git a/macros/CCEScan/ratioFit_spline.C b/macros/CCEScan/ratioFit_spline.C index 3d921e6..6fb6059 100755 --- a/macros/CCEScan/ratioFit_spline.C +++ b/macros/CCEScan/ratioFit_spline.C @@ -93,10 +93,10 @@ const char* diskVar = std::getenv ("DISK"); const char* homeVar = std::getenv ("CCEHOME"); - int fill = 1944; + int fill = 1616; TString det("TT"); TString lay("TTaU"); - int val = 3; + int val = 5; bool notSave = false; for (int i=1; iGet(Form("v_sector_%s", lay.Data())); TVectorD v_vdepl = *(TVectorD*)f_vdepl->Get(Form("v_vdepl_%s", lay.Data())); - + f_vdepl->Close(); - + Double_t a_tt_excl[18] = {0}; /*{ 2625,2626,2627,2628,2629,2630, 2631,2632,2633,2634,2635,2636, @@ -259,7 +249,7 @@ double height_val = v_input(0); double height_err = v_input(1); - + if (true || det.EqualTo("IT")) { height_val /= 0.1306; // transform back to the amplitude height_err /= 0.1306; // transform back to the amplitude @@ -268,8 +258,6 @@ TMath::Power(v_input(3)/v_input(2),2)+ 2*v_input(8)/(height_val*v_input(2))); } - - Printf(" %8.2f V: %8.4f+/-%8.4f",volt,height_val,height_err); @@ -436,7 +424,7 @@ lay.Data(),h_ratio->GetMean(),h_ratio->GetRMS()); TString dn_output(Form("%s/data/db",homeVar)); - TString fn_output(Form("%s/ratio.root", + TString fn_output(Form("%s/ratio_int.root", dn_output.Data())); TVectorD vr_ratio = getROOTVector(v_ratio); @@ -445,9 +433,9 @@ TFile* f_output = new TFile(fn_output.Data(),"UPDATE"); - TString vn_ratio = TString::Format("v_ratio_%s",lay.Data()); - TString vn_ratio_err = TString::Format("v_ratio_err_%s",lay.Data()); - TString vn_secRatio = TString::Format("v_sector_%s",lay.Data()); + TString vn_ratio = TString::Format("v_ratio_%s_val%d",lay.Data(),val); + TString vn_ratio_err = TString::Format("v_ratio_err_%s_val%d",lay.Data(),val); + TString vn_secRatio = TString::Format("v_sector_%s_val%d",lay.Data(),val); f_output->Delete(Form("%s;*",vn_ratio.Data())); f_output->Delete(Form("%s;*",vn_ratio_err.Data()));