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

added new option to musrt0: the key 'T' will set the t0 cursor to the maximum of the histogram

Browse Source
This commit is contained in:
nemu 2009-11-12 13:29:42 +00:00
parent 4071ab40bd
commit 7ccd021cfc
2 changed files with 49 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
src/classes/PMusrT0.cpp
View File

@ -107,8 +107,14 @@ cout << endl << "run Name = " << rawRunData->GetRunName()->Data() << ", runNo =
fHisto->SetMarkerSize(0.5); fHisto->SetMarkerSize(0.5);
fHisto->SetMarkerColor(TColor::GetColor(0,0,0)); // black fHisto->SetMarkerColor(TColor::GetColor(0,0,0)); // black
Double_t maxVal = 0.0;
fT0Estimated = 0;
for (UInt_t i=0; i<rawRunData->GetDataBin(histoNo-1)->size(); i++) { for (UInt_t i=0; i<rawRunData->GetDataBin(histoNo-1)->size(); i++) {
fHisto->SetBinContent(i+1, rawRunData->GetDataBin(histoNo-1)->at(i)); fHisto->SetBinContent(i+1, rawRunData->GetDataBin(histoNo-1)->at(i));
if (rawRunData->GetDataBin(histoNo-1)->at(i) > maxVal) {
maxVal = rawRunData->GetDataBin(histoNo-1)->at(i);
fT0Estimated = i;
}
} }
// generate canvas etc // generate canvas etc
@ -206,6 +212,8 @@ void PMusrT0::HandleCmdKey(Int_t event, Int_t x, Int_t y, TObject *selected)
Done(0); Done(0);
} else if (x == 'u') { // unzoom to the original range } else if (x == 'u') { // unzoom to the original range
UnZoom(); UnZoom();
} else if (x == 'T') { // set estimated t0 channel
SetEstimatedT0Channel();
} else if (x == 't') { // set t0 channel } else if (x == 't') { // set t0 channel
SetT0Channel(); SetT0Channel();
} else if (x == 'b') { // set first background channel } else if (x == 'b') { // set first background channel
@ -421,6 +429,45 @@ cout << endl << ">> PMusrT0::SetT0Channel(): binx = " << binx << endl;
fMainCanvas->Update(); fMainCanvas->Update();
} }
//--------------------------------------------------------------------------
// SetEstimatedT0Channel
//--------------------------------------------------------------------------
/**
* <p>
*
*/
void PMusrT0::SetEstimatedT0Channel()
{
// set t0 bin in msr-Handler
UInt_t idx = 0;
switch(fDetectorTag) {
case DETECTOR_TAG_FORWARD:
idx = fAddRunNo * fAddRunOffset / 2;
break;
case DETECTOR_TAG_BACKWARD:
idx = 1 + fAddRunNo * fAddRunOffset / 2;
break;
case DETECTOR_TAG_LEFT:
idx = 2 + fAddRunNo * fAddRunOffset / 2;
break;
case DETECTOR_TAG_RIGHT:
idx = 3 + fAddRunNo * fAddRunOffset / 2;
break;
default:
break;
}
fMsrHandler->SetMsrT0Entry(fRunNo, idx, fT0Estimated);
Double_t x = fHisto->GetXaxis()->GetBinCenter(fT0Estimated);
// shift line to the proper position
fT0Line->SetX1(x);
fT0Line->SetX2(x);
fMainCanvas->Modified(); // needed that Update is actually working
fMainCanvas->Update();
}
//-------------------------------------------------------------------------- //--------------------------------------------------------------------------
// SetDataFirstChannel // SetDataFirstChannel
//-------------------------------------------------------------------------- //--------------------------------------------------------------------------

2
src/include/PMusrT0.h
View File

@ -79,6 +79,7 @@ class PMusrT0 : public TObject, public TQObject
Int_t fDetectorTag; Int_t fDetectorTag;
Int_t fAddRunNo; Int_t fAddRunNo;
Int_t fAddRunOffset; Int_t fAddRunOffset;
Int_t fT0Estimated;
// canvas related variables // canvas related variables
TCanvas *fMainCanvas; TCanvas *fMainCanvas;
@ -101,6 +102,7 @@ class PMusrT0 : public TObject, public TQObject
void InitDataAndBkg(); void InitDataAndBkg();
void SetT0Channel(); void SetT0Channel();
void SetEstimatedT0Channel();
void SetDataFirstChannel(); void SetDataFirstChannel();
void SetDataLastChannel(); void SetDataLastChannel();
void SetBkgFirstChannel(); void SetBkgFirstChannel();