LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 3 Activity

Merged in master (pull request #11)

Browse Source 
Master
This commit is contained in:
Zaher Salman 2019-08-31 08:31:26 +00:00 committed by Andreas Suter
commit a1924c2756
2 changed files with 70 additions and 56 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/classes/PMsrHandler.cpp
View File

@ -2337,6 +2337,11 @@ Int_t PMsrHandler::WriteMsrFile(const Char_t *filename, std::map<UInt_t, TString
fout << "logy" << std::endl;
}
// lifetimecorrection
if (fPlots[i].fLifeTimeCorrection) {
fout << "lifetimecorrection" << std::endl;
}
// rrf_packing
if (fPlots[i].fRRFPacking) {
fout << "rrf_packing " << fPlots[i].fRRFPacking << std::endl;
@ -4797,6 +4802,8 @@ Bool_t PMsrHandler::HandlePlotEntry(PMsrLines &lines)
param.fLogX = true;
} else if (iter1->fLine.Contains("logy", TString::kIgnoreCase)) {
param.fLogY = true;
} else if (iter1->fLine.Contains("lifetimecorrection", TString::kIgnoreCase)) {
param.fLifeTimeCorrection = true;
} else if (iter1->fLine.Contains("view_packing", TString::kIgnoreCase)) {
tokens = iter1->fLine.Tokenize(" \t");
if (!tokens) {

109
src/external/MuSRFitGUI/MSR.pm vendored
View File

@ -299,32 +299,36 @@ sub CreateMSRUni {
}
$FRANGE_Line = "fit TINI TFIN";
$PAC_Line = "packing BINNING";
$PAC_Line = "packing BINNING";
$Single_RUN =
"$RUN_Line\n$Type_Line\n$Alpha_Line$Hist_Lines\n$T0DataBg$MAP_Line\n$FRANGE_Line\n$PAC_Line\n\n";
$Single_RUN =
"$RUN_Line\n$Alpha_Line$Hist_Lines\n$T0DataBg$MAP_Line\n";
# This is from CreateRUNBlk
#$DEBUG =~s /T0LINE/$T0_Line/;
#$DEBUG =~s /BGLINE/$Bg_Line/;
#$DEBUG =~s /DATALINE/$Data_Line/;
# Now add the appropriate values of fit range and packing
my $Range_Min = 8;
my $Range_Max = 0;
my $k = 0;
# Now add the appropriate values of fit range and packing
# If there no multiple ranges/packing then simply use GLOBAL block
if ($#TiVals == 0) {
$GLOBAL_Block = "GLOBAL\n$Type_Line\n$FRANGE_Line\n$PAC_Line\n\n";
} else {
$GLOBAL_Block = "GLOBAL\n$Type_Line\n\n";
$Single_RUN = $Single_RUN."$FRANGE_Line\n$PAC_Line\n";
}
my $Range_Min = 8;
my $Range_Max = 0;
my $k = 0;
foreach my $Ti (@TiVals) {
my $Tf = $TfVals[$k];
my $BIN = $BINVals[$k];
$RUN_Block = $RUN_Block . $Single_RUN;
$RUN_Block =~ s/TINI/$Ti/g;
$RUN_Block =~ s/TFIN/$Tf/g;
$RUN_Block =~ s/BINNING/$BIN/g;
$RUN_Block =~ s/TINI/$Ti/g;
$RUN_Block =~ s/TFIN/$Tf/g;
$RUN_Block =~ s/BINNING/$BIN/g;
$GLOBAL_Block =~ s/TINI/$Ti/g;
$GLOBAL_Block =~ s/TFIN/$Tf/g;
$GLOBAL_Block =~ s/BINNING/$BIN/g;
# For multiple ranges use this
if ( $Ti < $Range_Min ) { $Range_Min = $Ti; }
# For multiple ranges use this
if ( $Ti < $Range_Min ) { $Range_Min = $Ti; }
if ( $Tf > $Range_Max ) { $Range_Max = $Tf; }
$RUNS_Line = "$RUNS_Line " . $Range_Order;
++$k;
++$Range_Order;
@ -404,17 +408,12 @@ sub CreateMSRUni {
if ($Shared) {
$NoBg_Line = "norm 1\n";
}
else {
$NoBg_Line = "norm $PCount\n";
}
# Optional - add lifetime correction for SingleHist fits
if ( $All{"ltc"} eq "y" ) {
$NoBg_Line = $NoBg_Line . "lifetimecorrection\n";
}
}
# If you encounter NBg in the parameters list make sure
# to fill this line for the RUN block.
else {
$NoBg_Line = "norm $PCount\n";
}
}
# If you encounter NBg in the parameters list make sure
# to fill this line for the RUN block.
elsif ( $Param_ORG eq "NBg" ) {
if ($Shared) {
$NoBg_Line = $NoBg_Line . "backgr.fit 2\n";
@ -484,34 +483,38 @@ sub CreateMSRUni {
$Data_Line=$Data_Line."\nt0 ".$All{"t0$Hist"};
}
$FRANGE_Line = "fit TINI TFIN";
$PAC_Line = "packing BINNING";
$FRANGE_Line = "fit TINI TFIN";
$PAC_Line = "packing BINNING";
$Single_RUN = $EMPTY;
$Tmp_Hist_Line = $Hist_Lines;
$Tmp_Hist_Line =~ s/HIST/$Hist/g;
$Single_RUN = $Single_RUN
. "$RUN_Line\n$Type_Line\n$NoBg_Line$Tmp_Hist_Line\n$Data_Line\n$MAP_Line\n$FRANGE_Line\n$PAC_Line\n\n";
$Tmp_Hist_Line = $Hist_Lines;
$Tmp_Hist_Line =~ s/HIST/$Hist/g;
$Single_RUN = "$RUN_Line\n$NoBg_Line$Tmp_Hist_Line\n$Data_Line\n$MAP_Line\n";
# This is from CreateRUNBlk
#$DEBUG =~s /T0LINE/$T0_Line/;
#$DEBUG =~s /BGLINE/$Bg_Line/;
#$DEBUG =~s /DATALINE/$Data_Line/;
# Now add the appropriate values of fit range and packing
# If there no multiple ranges/packing then simply use GLOBAL block
if ($#TiVals == 0) {
$GLOBAL_Block = "GLOBAL\n$Type_Line\n$FRANGE_Line\n$PAC_Line\n\n";
} else {
$GLOBAL_Block = "GLOBAL\n$Type_Line\n\n";
$Single_RUN = $Single_RUN."$FRANGE_Line\n$PAC_Line\n";
}
# Now add the appropriate values of fit range and packing
my $Range_Min = 8;
my $Range_Max = 0;
my $k = 0;
my $Range_Min = 8;
my $Range_Max = 0;
my $k = 0;
foreach my $Ti (@TiVals) {
my $Tf = $TfVals[$k];
my $BIN = $BINVals[$k];
$RUN_Block = $RUN_Block . $Single_RUN;
$RUN_Block =~ s/TINI/$Ti/g;
$RUN_Block =~ s/TFIN/$Tf/g;
$RUN_Block =~ s/BINNING/$BIN/g;
$RUN_Block =~ s/TINI/$Ti/g;
$RUN_Block =~ s/TFIN/$Tf/g;
$RUN_Block =~ s/BINNING/$BIN/g;
$GLOBAL_Block =~ s/TINI/$Ti/g;
$GLOBAL_Block =~ s/TFIN/$Tf/g;
$GLOBAL_Block =~ s/BINNING/$BIN/g;
# For multiple ranges use this
if ( $Ti < $Range_Min ) { $Range_Min = $Ti; }
# For multiple ranges use this
if ( $Ti < $Range_Min ) { $Range_Min = $Ti; }
if ( $Tf > $Range_Max ) { $Range_Max = $Tf; }
$RUNS_Line = "$RUNS_Line " . $Range_Order;
++$k;
@ -546,7 +549,7 @@ sub CreateMSRUni {
$PCount $Param $value $error $error $minvalue $maxvalue";
}
$Full_T_Block = "THEORY\n###############################################################\n$Full_T_Block\n";
$Full_T_Block = "THEORY\n$Full_T_Block\n";
# COMMAND Block
$COMMANDS_Block = "COMMANDS\nFITMINTYPE\nSAVE\n";
@ -558,6 +561,8 @@ sub CreateMSRUni {
my $logxy = $EMPTY;
if ( $All{"logx"} eq "y" ) { $logxy = $logxy . "logx\n"; }
if ( $All{"logy"} eq "y" ) { $logxy = $logxy . "logy\n"; }
# Optional - add lifetime correction for SingleHist fits
if ( $All{"ltc"} eq "y" ) { $logxy = $logxy . "lifetimecorrection"; }
# Check if a plot range is defined (i.e. different from fit)
$PRANGE_Line = "use_fit_ranges";
@ -573,7 +578,7 @@ sub CreateMSRUni {
if ( $All{"ViewBin"}!=0 ) { $VIEWBIN_Line = "view_packing ".$All{"ViewBin"};}
my $RRFBlock=MSR::CreateRRFBlk(\%All);
$PLOT_Block = "PLOT $PLT\nruns $RUNS_Line\n$PRANGE_Line\n$VIEWBIN_Line\n$RRFBlock\n$logxy";
$PLOT_Block = "PLOT $PLT\n$logxy\nruns $RUNS_Line\n$PRANGE_Line\n$VIEWBIN_Line\n$RRFBlock";
#END - PLOT Block
# FFT Block
@ -594,6 +599,8 @@ $Full_T_Block
###############################################################
$FUNCTIONS_Block
###############################################################
$GLOBAL_Block
###############################################################
$RUN_Block
###############################################################
$COMMANDS_Block