diff --git a/macros/CCEScan/voltageFit_spline.C b/macros/CCEScan/voltageFit_spline.C index 0ee80d6..e90c165 100755 --- a/macros/CCEScan/voltageFit_spline.C +++ b/macros/CCEScan/voltageFit_spline.C @@ -177,7 +177,7 @@ TString dn_data(Form("%s/data/ST/Aging" ,diskVar)); - TString fn_data(Form("%s/CCEScan.root",//Michele for TESTING + TString fn_data(Form("%s/CCEScan_TEST.root",//Michele for TESTING dn_data.Data())); TString fn_dummy(Form("%s/check",dn_data.Data())); @@ -321,7 +321,7 @@ TMultiGraph* mg_volt = new TMultiGraph("mg_volt","Voltage multigraph"); - // Extract the ratio values for the sampling point + // Extract the ratio values for the sampling point TVectorD* vp_ratio = (TVectorD*)f_ratio->Get(Form("v_ratio_%s",lay.Data())); TVectorD* vp_ratio_err = (TVectorD*)f_ratio->Get(Form("v_ratio_err_%s",lay.Data())); TVectorD* vp_secRatio = (TVectorD*)f_ratio->Get(Form("v_sector_%s",lay.Data())); @@ -346,22 +346,22 @@ ratio_err = STTool::TTratio_err; } - /* Elena - // if (false && det.EqualTo("TT")) { - */ - // Using ratio per sector - if (true) { - int ratioInd = getClosestIndex(v_secRatio,(double)sector); - if (!(ratioInd<0) && (int)v_secRatio(ratioInd)==sector && true) { - Info("voltageFit_spline","Using specific ratio for sector %d",sector); - ratio_val = v_ratio(ratioInd); - ratio_err = v_ratio_err(ratioInd); + // Elena - Cancelled by Michele because of what said in l.159 in the Radiation Damage paper + if (false && det.EqualTo("TT")) { + + // Using ratio per sector + if (true) { + int ratioInd = getClosestIndex(v_secRatio,(double)sector);//Gets the index in v_secRatio that has a double that is closest to sector + if (!(ratioInd<0) && (int)v_secRatio(ratioInd)==sector && true) {//if this is not senseless and the sector is correct + Info("voltageFit_spline","Using specific ratio for sector %d",sector); + ratio_val = v_ratio(ratioInd); + ratio_err = v_ratio_err(ratioInd); + } } } - // Define fit range double fitMax = 510.0; double fitMin = 280.0;