// #define VETO_THRESHOLD 3
// #define VETO_THRESHOLD 40
Float_t vetoThreshold64[6] = {8, 4, 4, 4, 3.5, 2.1};
-Float_t vetoThreshold80[6] = {100, 100, 100, 100, 100, 100};
+Float_t vetoThreshold80[6] = {7, 4, 6, 4, 3.5, 2.1};
#define LABORBUCH "RelevantRuns.csv"
// #define DATAPATH "/local/sstrohauer/data/data_15.10.2013_11:36/"
// #define DATAPATH "/local/sstrohauer/data/data_25.10.2013_14:50/"
// #define DATAPATH "/local/sstrohauer/data/data_4.11.2013_13:57/"
// #define DATAPATH "/local/sstrohauer/data/data_6.11.2013_17:50/"
+// #define DATAPATH "/local/sstrohauer/data/data_10.1.2014_21:20/"
+
#define DATAPATH "/local/sstrohauer/data/data_19.11.2013_14:11/"
#define USE_SEPARATE_NOISE_RUN true
#define DRAW_SEEDPIXEL_DISTRIBUTION false
#define DRAW_PIXEL_DISTRIBUTION_HISTOGRAMS false
// examples for fit functions: gaus, landau
-#define FIT_FUNC "gaus"
+#define FIT_FUNC "landau"
#define DRAW_FITS false
#define DRAW_NOISE_HISTOGRAM false
#define PRINT_EXPECTATION_VALUE_HITS_PER_FRAME false
-#define PRINT_HEADER true
+#define PRINT_ADDITIONAL_INFO false
+#define PRINT_HEADER false
#define xTitle "Gesammelte Ladung / ADC"
-#define yTitle "Haeufigkeit"
-#define graphTitle ""
+#define yTitle "Haeufigkeit / #frac{1}{2 ADC}"
+#define graphTitle "Seed Pixel Ladungsspektrum"
void PlotGraph(Float_t selTemp, Int_t selChip, std::string selSource, Int_t selMatrix6480, Int_t selSubmatrix, Float_t selRadDose)
{
histNtupleSum[nHistNtuple]->SetStats(kFALSE);
histNtupleSum[nHistNtuple]->SetLineStyle(2); // make summed graphs dashed
histNtupleVeto[nHistNtuple]->SetLineColor(nHistNtuple+1>4 ? nHistNtuple+2:nHistNtuple+1); // omit yellow line!
- histNtupleVeto[nHistNtuple]->SetStats(kFALSE);
+ histNtupleVeto[nHistNtuple]->SetStats(kTRUE);
histNtupleVeto[nHistNtuple]->SetLineStyle(1); // make summed graphs dotted
// histNtupleVeto[nHistNtuple]->SetLineStyle(3); // make summed graphs dotted
histPixelNumIsSeedSelection[nHistNtuple]->SetLineColor(nHistNtuple+1>4 ? nHistNtuple+2:nHistNtuple+1); // omit yellow line!
histPixelDistribution1[nHistNtuple]->SetLineColor(1); // black line
histPixelDistribution2[nHistNtuple]->SetLineColor(2); // red line
// TString legendEntry = Form("#splitline{Run %i: T_{set}=%.1f, T_{sens}=%.1f, Chip=%i, }{Source=%s, Matrix=%i, Submatrix=%i, RadDose=%.2E}", runNo, temperature, tempSens, chip, TString(source).Data(), matrix6480, submatrix, radDose);
- TString legendEntry = TString(Form("Run %i: ",runNo)) + (selTemp==-1 ? Form(" T_{set}=%.1f",temperature):"") + Form(" T_{sens}=%.1f",tempSens) + (selChip==-1 ? Form(" Chip=%i",chip):"") + (selSource=="-1" ? Form(" Source=%s",source.c_str()):"") + (selMatrix6480==-1 ? Form(" Matrix=%i",matrix6480):"") + (selSubmatrix==-1 ? Form(" Submatrix=%i",submatrix):"") + (selRadDose==-1 ? Form(" RadDose=%.2E",radDose):"");
+ TString legendEntry = TString(Form("Run %i: ",runNo)) + (selTemp==-1 ? Form(" T_{set}=%.1f",temperature):"") + Form(" T_{sens}=%.1f",tempSens) + (selChip==-1 ? Form(" Chip=%i",chip):"") + (selSource=="-1" ? Form(" Source=%s",source.c_str()):"") + (selMatrix6480==-1 ? Form(" Matrix=%i",matrix6480):"") + (selSubmatrix==-1 ? Form(" Submatrix=%i",submatrix):"") + (selRadDose==-1 ? TString(Form(" RadDose=%.2E ",radDose))+(chip==2 ? "rad":"n_{eq}"):"");
// TString legendEntry = TString(Form("test") + (selTemp==-1 ? Form(", T_{set}=%.1f",temperature):"");// + Form(", T_{sens}=%.1f",tempSens) + (selChip==-1 ? Form(", Chip=%i",chip):"") + (selSource=="-1" ? Form(", Source=%s",source):"") + (selMatrix6480==-1 ? Form(", Matrix=%i",matrix6480):"") + (selSubmatrix==-1 ? Form(", Submatrix=%i",submatrix):"") + (selRadDose==-1 ? Form(", RadDose=%.2E",radDose):"");
leg->AddEntry(histNtuple[nHistNtuple], legendEntry);
// leg->AddEntry((TObject*)0, "", "");
// histNtuple[nHistNtuple]->Draw();
// cout << "Hits per Frame: " << nentries/frame << endl;
+ if (PRINT_ADDITIONAL_INFO)
+ cout << "Total number of hits: " << nentries << endl;
//-------------------------------------------------
// calculate median of noise (with error estimation)
if (!runFinished_N)
continue;
- if ( (temperature_N == selTemp || selTemp == -1)
- && (chip_N == selChip || selChip == -1)
+ if ( (temperature_N == temperature)
+ && (chip_N == chip)
&& (source_N == "none")
- && (matrix6480_N == selMatrix6480 || selMatrix6480 == -1)
- && (radDose_N == selRadDose || selRadDose == -1))
+ && (matrix6480_N == matrix6480)
+ && (radDose_N == radDose))
{
noiseRunNo = runNo_N;
}
+// old version (obviously wrong):
+// if ( (temperature_N == selTemp || selTemp == -1)
+// && (chip_N == selChip || selChip == -1)
+// && (source_N == "none")
+// && (matrix6480_N == selMatrix6480 || selMatrix6480 == -1)
+// && (radDose_N == selRadDose || selRadDose == -1))
+// {
+// noiseRunNo = runNo_N;
+// }
}
catch(...){
cout << "Error while reading laboratory book (run number " << runNo << ")!\n";
if (noiseRunNo == 0)
cout << "No corresponding noise run found!\n";
+
+// cout << noiseRunNo << endl;
+
+
// open the file
TString path_N;
if (USE_SEPARATE_NOISE_RUN) {
if (TString(FIT_FUNC)=="gaus") {
(*it)->Fit(fitFunc, "N,Q,W", "", xValMax-10, xValMax+10);
posMax = fitFunc->GetParameter(1);
- fitMaxPosXError.push_back(fitFunc->GetParameter(2)); // sigma of the gaus fit
+ fitMaxPosXError.push_back(fitFunc->GetParError(1)); // error of the mean
fitFunc->DrawCopy("same");
}
else if (TString(FIT_FUNC)=="landau") {
Float_t maxFitMax = 1000;
(*it)->Fit(fitFunc, "N,Q,W", "", 20, RIGHT_BOUNDARY);
+ // Attention, attention: The parameter MPV from the ROOT-fit is NOT exactly the most probable value (only approximative)! It seems that one can calculate the real most probable value by the following formula:
+ // realMPV = MPV - 0.22278*Sigma
+ // see documentation: http://root.cern.ch/root/html/TMath.html#TMath:Landau
+ // The other way is simply to use the GetMaximum-function:
fitMax1 = fitFunc->GetMaximumX(0,RIGHT_BOUNDARY);
+ fitMaxPosXError.push_back(fitFunc->GetParError(1)); // error of the MPV. Does it make sense to choose the error calculated by ROOT? ROOT does not find the MPV at all...
+// cout << "error of the MPV = " << fitFunc->GetParError(1) << endl;
+ // the following fits are note required I think, because in my view a correct error estimation is the above fitFunc->GetParError(1).
fitFunc->DrawCopy("same");
(*it)->Fit(fitFunc, "N,Q,W", "", 20, fitMax1+10);
fitMax2 = fitFunc->GetMaximumX(0,RIGHT_BOUNDARY);
fitMax3 = fitFunc->GetMaximumX(0,RIGHT_BOUNDARY);
fitFunc->SetLineColor(kGreen);
fitFunc->DrawCopy("same");
- fitFunc->SetLineStyle(1); // normal for the following fits
+ fitFunc->SetLineStyle(1); // normal style for the following fits
- fitMaxPosXError.push_back(fitFunc->GetParameter(2)); // sigma of the landau fit
- cout << "yoho!\n";
- Int_t firstabove = (*it)->FindFirstBinAbove();
- Int_t lastabove = (*it)->FindLastBinAbove();
- cout << firstabove << endl;
- cout << lastabove << endl;
- Int_t integral = (*it)->Integral(firstabove,lastabove);
- cout << "integral = " << integral << endl;
- cout << "error = " << fitFunc->GetParameter(2) << endl;
- cout << "real error = " << fitFunc->GetParameter(2)/integral << endl << "end\n";
+// cout << "yoho!\n";
+// Int_t firstabove = (*it)->FindFirstBinAbove();
+// Int_t lastabove = (*it)->FindLastBinAbove();
+// cout << firstabove << endl;
+// cout << lastabove << endl;
+// Int_t integral = (*it)->Integral(firstabove,lastabove);
+// cout << "integral = " << integral << endl;
+// cout << "error of the MPV = " << fitFunc->GetParError(1) << endl;
+// cout << "sigma/sqrt(N) = " << fitFunc->GetParameter(2)/integral << endl << "end\n";
(*it)->GetXaxis()->SetRange((*it)->GetXaxis()->FindBin(fitMaxPosX[cnt]),(*it)->GetXaxis()->FindBin(RIGHT_BOUNDARY)); // look only for maxima with x greater than 20
Int_t xValMax = (*it)->GetBinCenter((*it)->GetMaximumBin());
(*it)->Fit(fitFunc, "N,Q,W", "", xValMax-15, xValMax+15);
- fitMaxPosXSumError.push_back(fitFunc->GetParameter(2)); // sigma of the gaus fit
+ fitMaxPosXSumError.push_back(fitFunc->GetParError(1)); // error of the mean
(*it)->GetXaxis()->UnZoom();
fitFunc->DrawCopy("same");
(*it)->GetXaxis()->SetRange((*it)->GetXaxis()->FindBin(fitMaxPosX[cnt]),(*it)->GetXaxis()->FindBin(RIGHT_BOUNDARY)); // look only for maxima with x greater than 20
Int_t xValMax = (*it)->GetBinCenter((*it)->GetMaximumBin());
(*it)->Fit(fitFunc, "N,Q,W", "", xValMax-15, xValMax+15);
- fitMaxPosXVetoError.push_back(fitFunc->GetParameter(2)); // sigma of the gaus fit
+ fitMaxPosXVetoError.push_back(fitFunc->GetParError(1)); // error of the mean
(*it)->GetXaxis()->UnZoom();
fitFunc->DrawCopy("same");
-
+ if (PRINT_ADDITIONAL_INFO) {
for (std::vector<TH1F*>::iterator it = histNtuple.begin(); it != histNtuple.end(); ++it){
Int_t cnt = std::distance(histNtupleSum.begin(), it) ;
- // get median and error estimation
- Double_t const probabilities[] = {0.30, 0.45, 0.50, 0.55, 0.70};
- Double_t quantiles[5];
- (*it)->GetQuantiles( 5, quantiles, probabilities);
- cout << "-2Quant -1Quant Median 1Quant 2Quant\n";
- cout << quantiles[0] << "\t" << quantiles[1] << "\t" << quantiles[2] << "\t" << quantiles[3] << "\t" << quantiles[4] << endl;
-// noiseLowest17percent.push_back((Float_t) quantiles[0]); // left error bar (lowest 17% noise values are below this quantile)
-// noiseMedian.push_back((Float_t) quantiles[1]); // median of the noise distribution
-// noiseHighest17percent.push_back((Float_t) quantiles[2]); // right error bar (highest 17% noise values are higher than this quantile)
-// cout << "q1: " << quantiles[1] - quantiles[0] << "\tq2: " << quantiles[1] << "\tq3: " << quantiles[2] - quantiles[1] << endl;
-
+ // get median and error estimation
+ Double_t const probabilities[] = {0.30, 0.45, 0.50, 0.55, 0.70};
+ Double_t quantiles[5];
+ (*it)->GetQuantiles( 5, quantiles, probabilities);
+ cout << "-2Quant -1Quant Median 1Quant 2Quant\n";
+ cout << quantiles[0] << "\t" << quantiles[1] << "\t" << quantiles[2] << "\t" << quantiles[3] << "\t" << quantiles[4] << endl;
+ // noiseLowest17percent.push_back((Float_t) quantiles[0]); // left error bar (lowest 17% noise values are below this quantile)
+ // noiseMedian.push_back((Float_t) quantiles[1]); // median of the noise distribution
+ // noiseHighest17percent.push_back((Float_t) quantiles[2]); // right error bar (highest 17% noise values are higher than this quantile)
+ // cout << "q1: " << quantiles[1] - quantiles[0] << "\tq2: " << quantiles[1] << "\tq3: " << quantiles[2] - quantiles[1] << endl;
+
}
-
+ }
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tColPeak":"")
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tColPkEr":"")
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tMPV":"")
+ + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tErrMPV":"")
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tErrL":"")
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tErrR":"")
+ (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? "\tSumColPeak":"")
// cout << runDetails[i] << Form("\t%.2f", fitMaxPosX[i]) << Form("\t%.2f", fitMaxPosXSum[i]) << endl;
Float_t noiseErrorMean = (noiseHighest17percent[i] - noiseLowest17percent[i])/2; // mean noise error
cout << runDetails[i]
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS ? Form("\t%.2f", fitMaxPosX[i]):"") // ColPeak/MPV
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.2f", fitMaxPosXError[i]):"") // gaus: sigma = error
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.2f", fitLandauErrorLeft[i]):"") // landau: ErrL
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.2f", fitLandauErrorRight[i]):"") // landau: ErrR
- << (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? Form("\t%.2f\t", fitMaxPosXSum[i]):"") // SumColPeak
- << (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? Form("\t%.2f\t", fitMaxPosXSumError[i]):"") // SumColPeakErr
- << (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? Form("\t%.2f", fitMaxPosXVeto[i]):"") // CalPeak
- << (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? Form("\t%.2f", fitMaxPosXVetoError[i]):"") // CalPkErr
- << Form("\t%.2f", noiseMedian[i]) // noise
- << Form("\t%.2f", noiseLowest17percent[i])
- << Form("\t%.2f", noiseHighest17percent[i])
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.2f", fitMaxPosX[i]/fitMaxPosXVeto[i]):"") // CCESeed = ColPeak/CalPeak
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.2f", sqrt( pow(1/fitMaxPosXVeto[i]*fitMaxPosXError[i], 2) + pow(fitMaxPosX[i]/pow(fitMaxPosXVeto[i], 2)*fitMaxPosXVetoError[i], 2) )):"") // CCESeedErr
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.2f", fitMaxPosXSum[i]/fitMaxPosXVeto[i]):"") // CCESum = SumColPeak/CalPeak
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.2f", sqrt( pow(1/fitMaxPosXVeto[i]*fitMaxPosXSumError[i], 2) + pow(fitMaxPosX[i]/pow((fitMaxPosXSum[i]), 2)*fitMaxPosXSumError[i], 2) )):"") // CCESeedErr
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.2f", fitMaxPosX[i]/noiseMedian[i]):"") // S/N-ratio
- << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.2f", sqrt( pow(1/noiseMedian[i]*fitMaxPosXError[i], 2) + pow(fitMaxPosX[i]/pow(noiseMedian[i], 2)*noiseErrorMean, 2) )) :"") // S/N error
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS ? Form("\t%.4f", fitMaxPosX[i]):"") // ColPeak/MPV
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.4f", fitMaxPosXError[i]):"") // gaus: sigma = error
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.4f", fitMaxPosXError[i]):"") // landau: MPVErr
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.4f", fitLandauErrorLeft[i]):"") // landau: ErrL
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.4f", fitLandauErrorRight[i]):"") // landau: ErrR
+ << (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? Form("\t%.4f\t", fitMaxPosXSum[i]):"") // SumColPeak
+ << (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? Form("\t%.4f\t", fitMaxPosXSumError[i]):"") // SumColPeakErr
+ << (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? Form("\t%.4f", fitMaxPosXVeto[i]):"") // CalPeak
+ << (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? Form("\t%.4f", fitMaxPosXVetoError[i]):"") // CalPkErr
+ << Form("\t%.4f", noiseMedian[i]) // noise
+ << Form("\t%.4f", noiseLowest17percent[i])
+ << Form("\t%.4f", noiseHighest17percent[i])
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.4f", fitMaxPosX[i]/fitMaxPosXVeto[i]):"") // CCESeed = ColPeak/CalPeak
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.4f", sqrt( pow(1/fitMaxPosXVeto[i]*fitMaxPosXError[i], 2) + pow(fitMaxPosX[i]/pow(fitMaxPosXVeto[i], 2)*fitMaxPosXVetoError[i], 2) )):"") // CCESeedErr
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_SUMMED_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.4f", fitMaxPosXSum[i]/fitMaxPosXVeto[i]):"") // CCESum = SumColPeak/CalPeak
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_SUMMED_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? Form("\t%.4f", sqrt( pow(1/fitMaxPosXVeto[i]*fitMaxPosXSumError[i], 2) + pow(fitMaxPosX[i]/pow((fitMaxPosXSum[i]), 2)*fitMaxPosXSumError[i], 2) )):"") // CCESeedErr
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.4f", fitMaxPosX[i]/noiseMedian[i]):"") // S/N-ratio
+ << (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? Form("\t%.4f", sqrt( pow(1/noiseMedian[i]*fitMaxPosXError[i], 2) + pow(fitMaxPosX[i]/pow(noiseMedian[i], 2)*noiseErrorMean, 2) )) :"") // S/N error
<< endl;
}
29080 6/13/2013 -70 -57.0 50000 5 none 64 3.00E+12 1.57 y
29081 6/13/2013 -70 -59.5 50000 5 none 80 3.00E+12 1.55 y
-29082 6/17/2013 20 21.5 50000 5 Sr90 64 3.00E+12 1.57 y
-29083 6/17/2013 20 21.0 50000 5 Sr90 80 3.00E+12 1.55 y
-29084 6/17/2013 -20 -14.0 50000 5 Sr90 64 3.00E+12 1.57 y
-29085 6/17/2013 -20 -15.0 50000 5 Sr90 80 3.00E+12 1.55 y
-29086 6/17/2013 -70 -54.0 50000 5 Sr90 64 3.00E+12 1.57 y
-29087 6/17/2013 -70 -54.5 50000 5 Sr90 80 3.00E+12 1.55 y
29088 6/20/2013 20 21.0 50000 6 Fe55,dimmed 64 1.00E+13 1.54 y 19 … 21 1.62 2.52
29111 10/18/2013 -70 -55.0 50000 1 Fe55,superdimmed 80 0.00E+00 3.10 y BULLSHIT (eingefroren)
29112 and 29113 were made by Benny and are in the DB (same as 29110 and 29111, and unfortunately also freezed)
-
+
+
+ the following measurements are for exploring the cause for the very high noise in the previous measurements
+29114 10/31/2013 20 21.0 50000 1 Cd109 64 0.00E+00 3.20 y measurement with Cd109 with the same parameters Vba, Vbo and Vdiode like all previous measurements 0.24 2.50
+29115 11/1/2013 20 21.0 50000 1 Cd109 64 0.00E+00 3.20 y 3.31 2.50
+29116 11/1/2013 20 21.0 50000 1 Cd109 64 0.00E+00 2.74 y Vbo, Vba, Vdiode 3.31 2.50
+29117 11/1/2013 20 21.0 50000 1 Cd109 64 0.00E+00 2.74 y 3.31 2.34
+29118 11/4/2013 20 21.0 50000 1 Cd109 64 0.00E+00 2.74 y 3.31 2.59
+
+29119 11/11/2013 20 21.0 50000 1 none 64 0.00E+00 3.31 y vdiode max 3.31 2.84
+29120 11/11/2013 20 21.0 50000 1 Cd109 64 0.00E+00 3.31 y vdiode max 3.31 2.84
+29121 11/11/2013 20 21.0 50000 1 none 64 0.00E+00 3.31 y vdiode min 3.31 2.57
+29122 11/11/2013 20 21.0 50000 1 Cd109 64 0.00E+00 3.31 y vdiode min 3.31 2.57
+
+
+
+
+29123 11/25/2013 20 23.0 1000 1 none 64 0.00E+00 3.20 y Varying Vba 0.00 2.50
+29124 11/25/2013 20 23.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y Varying Vba 0.00 2.50
+29125 11/25/2013 20 22.0 1000 1 none 64 0.00E+00 3.20 y Varying Vba 1.66 2.50
+29126 11/25/2013 20 21.5 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y Varying Vba 1.66 2.50
+29127 11/25/2013 20 23.0 1000 1 none 64 0.00E+00 3.20 y Varying Vba 3.31 2.50
+29128 11/25/2013 20 23.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y Varying Vba 3.31 2.50
+
+
+ The following measurements are to be made because I probably forgot to switch on the power supply of the sensor in run no 29082…29087.
+29129 1/21/2014 20 22.5 50000 5 Sr90 64 3.00E+12 1.57 y 1.86 2.52
+29130 1/21/2014 20 22.5 50000 5 Sr90 80 3.00E+12 1.55 y 1.69 2.52
+29131 1/22/2014 -20 -13.0 50000 5 Sr90 64 3.00E+12 1.57 y
+29132 1/22/2014 -20 -11.5 50000 5 Sr90 80 3.00E+12 1.55 y
+29133 1/22/2014 -30 -21.0 50000 5 Sr90 64 3.00E+12 1.57 y
+29134 1/22/2014 -30 -21.5 50000 5 Sr90 80 3.00E+12 1.55 y
+29135 -70 50000 5 Sr90 64 3.00E+12 1.57
+29136 -70 50000 5 Sr90 80 3.00E+12 1.55
29126 11/25/2013 20 21.5 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y Varying Vba 1.66 2.50
29127 11/25/2013 20 23.0 1000 1 none 64 0.00E+00 3.20 y Varying Vba 3.31 2.50
29128 11/25/2013 20 23.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y Varying Vba 3.31 2.50
+
+
+ The following measurements are to be made because I probably forgot to switch on the power supply of the sensor in run no 29082…29087.
+29129 1/21/2014 20 22.5 50000 5 Sr90 64 3.00E+12 1.57 y 1.86 2.52
+29130 1/21/2014 20 22.5 50000 5 Sr90 80 3.00E+12 1.55 y 1.69 2.52
+29131 1/22/2014 -20 -13.0 50000 5 Sr90 64 3.00E+12 1.57 y
+29132 1/22/2014 -20 -11.5 50000 5 Sr90 80 3.00E+12 1.55 y
+29133 1/22/2014 -30 -21.0 50000 5 Sr90 64 3.00E+12 1.57 y
+29134 1/22/2014 -30 -21.5 50000 5 Sr90 80 3.00E+12 1.55 y
+29135 -70 50000 5 Sr90 64 3.00E+12 1.57
+29136 -70 50000 5 Sr90 80 3.00E+12 1.55
}
if(clusterSum<3.*GetMeanNoise()){continue;} //Cut on whole cluster 3 times Mean(noise)
-
+// Or make a hard cut:
+// Double_t clusterSum=0;
+// for(Int_t i=0;i<25;i++)
+// {
+// clusterSum+=CLUSTER[i];// Charge of whole cluster
+//
+// }
+// if(clusterSum<3.*3){continue;} //Cut on whole cluster 3 times Mean(noise)
#define INPUTDATAPATH "/jspc12_F/Mi29/"
#define OUTPUTDATAPATH_MAINFOLDER "/local/sstrohauer/data/"
// use this #define for processing to a specific file
-#define OUTPUTDATAPATH_FIXED_PATH "/local/sstrohauer/data/data_19.11.2013_14:11/"
+// #define OUTPUTDATAPATH_FIXED_PATH "/local/sstrohauer/data/data_19.11.2013_14:11/"
using std::cout;
using std::endl;
+++ /dev/null
-set macro
-
-##### Color Palette by Color Scheme Designer
-##### Palette URL: http://colorschemedesigner.com/#3K40zsOsOK-K-
-
- blue_000 = "#A9BDE6" # = rgb(169,189,230)
- blue_025 = "#7297E6" # = rgb(114,151,230)
- blue_050 = "#1D4599" # = rgb(29,69,153)
- blue_075 = "#2F3F60" # = rgb(47,63,96)
- blue_100 = "#031A49" # = rgb(3,26,73)
-
- green_000 = "#A6EBB5" # = rgb(166,235,181)
- green_025 = "#67EB84" # = rgb(103,235,132)
- green_050 = "#11AD34" # = rgb(17,173,52)
- green_075 = "#2F6C3D" # = rgb(47,108,61)
- green_100 = "#025214" # = rgb(2,82,20)
-
- red_000 = "#F9B7B0" # = rgb(249,183,176)
- red_025 = "#F97A6D" # = rgb(249,122,109)
- red_050 = "#E62B17" # = rgb(230,43,23)
- red_075 = "#8F463F" # = rgb(143,70,63)
- red_100 = "#6D0D03" # = rgb(109,13,3)
-
- brown_000 = "#F9E0B0" # = rgb(249,224,176)
- brown_025 = "#F9C96D" # = rgb(249,201,109)
- brown_050 = "#E69F17" # = rgb(230,159,23)
- brown_075 = "#8F743F" # = rgb(143,116,63)
- brown_100 = "#6D4903" # = rgb(109,73,3)
-
- grid_color = "#d5e0c9"
- text_color = "#6a6a6a"
-
- my_font = "SVBasic Manual, 12"
- my_font_file = "~/local/share/fonts/defaults/LiberationMono-Regular.ttf"
- my_export_sz = "1024,768"
-
- my_line_width = "2"
- my_axis_width = "1.5"
- my_ps = "1.2"
- my_font_size = "14"
-
-# must convert font fo svg and ps
-#set term svg size @my_export_sz fname my_font fsize my_font_size enhanced dynamic rounded
-# set term png size @my_export_sz large font my_font
-# set term jpeg size @my_export_sz large font my_font
-#set term wxt enhanced font my_font
-
-set style data linespoints
-set style function lines
-set pointsize my_ps
-
-set style line 1 linecolor rgbcolor blue_025 linewidth @my_line_width pt 7
-set style line 2 linecolor rgbcolor green_025 linewidth @my_line_width pt 5
-set style line 3 linecolor rgbcolor red_025 linewidth @my_line_width pt 9
-set style line 4 linecolor rgbcolor brown_025 linewidth @my_line_width pt 13
-set style line 5 linecolor rgbcolor blue_050 linewidth @my_line_width pt 11
-set style line 6 linecolor rgbcolor green_050 linewidth @my_line_width pt 7
-set style line 7 linecolor rgbcolor red_050 linewidth @my_line_width pt 5
-set style line 8 linecolor rgbcolor brown_050 linewidth @my_line_width pt 9
-set style line 9 linecolor rgbcolor blue_075 linewidth @my_line_width pt 13
-set style line 10 linecolor rgbcolor green_075 linewidth @my_line_width pt 11
-set style line 11 linecolor rgbcolor red_075 linewidth @my_line_width pt 7
-set style line 12 linecolor rgbcolor brown_075 linewidth @my_line_width pt 5
-set style line 13 linecolor rgbcolor blue_100 linewidth @my_line_width pt 9
-set style line 14 linecolor rgbcolor green_100 linewidth @my_line_width pt 13
-set style line 15 linecolor rgbcolor red_100 linewidth @my_line_width pt 11
-set style line 16 linecolor rgbcolor brown_100 linewidth @my_line_width pt 7
-set style line 17 linecolor rgbcolor "#224499" linewidth @my_line_width pt 5
-
-#set style line 1 linecolor rgbcolor "#a0bae9" linewidth @my_line_width pt 7
-#set style line 2 linecolor rgbcolor "#ff7f7f" linewidth @my_line_width pt 5
-#set style line 3 linecolor rgbcolor "#80c65a" linewidth @my_line_width pt 9
-#set style line 4 linecolor rgbcolor "#ffcc7f" linewidth @my_line_width pt 13
-#set style line 5 linecolor rgbcolor "#dedc06" linewidth @my_line_width pt 11
-#set style line 6 linecolor rgbcolor "#7711ff" linewidth @my_line_width
-#set style line 7 linecolor rgbcolor "#ff0000" linewidth @my_line_width
-#set style line 8 linecolor rgbcolor "#008000" linewidth @my_line_width
-#set style line 9 linecolor rgbcolor "#ff9900" linewidth @my_line_width
-#set style line 10 linecolor rgbcolor "#aa9900" linewidth @my_line_width
-#set style line 11 linecolor rgbcolor "#990066" linewidth @my_line_width
-#set style line 12 linecolor rgbcolor "#990000" linewidth @my_line_width
-#set style line 13 linecolor rgbcolor "#003971" linewidth @my_line_width
-#set style line 14 linecolor rgbcolor "#76a4fb" linewidth @my_line_width
-#set style line 15 linecolor rgbcolor "#d5e0c9" linewidth @my_line_width
-#set style line 16 linecolor rgbcolor "#e5ecf9" linewidth @my_line_width
-#set style line 17 linecolor rgbcolor "#224499" linewidth @my_line_width
-
-## plot 1,2,3,4,5,6,7,8,9
-set style increment user
-set style arrow 1 filled
-
-## used for bar chart borders
-set style fill solid 0.5
-
-# Grey background
-#set object 1 rectangle from screen 0, screen 0 to screen 1, screen 1 behind fc rgbcolor "#cccccc"
-
-set size noratio
-set samples 300
-
-#set log x
-#set mxtics 10 # Makes logscale look good.
-
-set xtics textcolor rgb text_color font my_font
-set ytics textcolor rgb text_color font my_font
-set xlabel "X Label (unit)" textcolor rgb text_color font my_font
-set ylabel "Y Label (unit)" textcolor rgb text_color font my_font
-set label textcolor rgb text_color font my_font
-
-set title "Top Title" textcolor rgb text_color font my_font
-
-set xtics nomirror
-set ytics nomirror
-set border 3 back lw @my_axis_width lc rgb text_color
-# Bit plot splot
-# 1 bottom bottom left front
-# 2 left bottom left back
-# 4 top bottom right front
-# 8 right bottom right back
-# 16 no effect left vertical
-# 32 no effect back vertical
-# 64 no effect right vertical
-# 128 no effect front vertical
-# 256 no effect top left back
-# 512 no effect top right back
-# 1024 no effect top left front
-# 2048 no effect top right front
-
-set grid back lt 0 lc rgb grid_color # make dashed grid in the background
-#set key outside box width 2 height 2 enhanced spacing 2 # put legend outside the plot
-set key top right # legend in the upper right corner
\ No newline at end of file
+++ /dev/null
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-29024 20 21.5 1 Sr90 64 0 0 45.85 2.25 1.59 4.08 3.89 4.27 11.24
-29024 20 21.5 1 Sr90 64 0 1 35.28 0.57 0.02 2.84 2.74 2.98 12.41
-29024 20 21.5 1 Sr90 64 0 2 36.95 0.2 1.19 2.85 2.74 2.97 12.98
-29024 20 21.5 1 Sr90 64 0 3 36.65 0.14 0.67 2.83 2.71 2.95 12.96
-29024 20 21.5 1 Sr90 64 0 4 31.18 0.07 0.19 2.82 2.72 2.96 11.05
-29024 20 21.5 1 Sr90 64 0 5 25.68 0.09 0.43 2.44 2.35 2.57 10.54
-
-29025 20 21.5 1 Sr90 80 0 0 38.39 3.02 2.74 4.26 4.06 4.46 9.02
-29025 20 21.5 1 Sr90 80 0 1 29.5 0.78 0.17 2.81 2.69 2.97 10.5
-29025 20 21.5 1 Sr90 80 0 2 43.58 1.31 2.84 4.54 4.3 4.75 9.6
-29025 20 21.5 1 Sr90 80 0 3 32.7 0.01 0.65 2.54 2.39 2.76 12.87
-29025 20 21.5 1 Sr90 80 0 4 30.58 0.51 0.89 2.57 2.38 2.79 11.91
-29025 20 21.5 1 Sr90 80 0 5 23.24 0.14 0.11 2.26 2.14 2.4 10.29
-
-29026 -20 -14 1 Sr90 64 0 0 51.81 5.6 6.65 2.37 2.29 2.51 21.82
-29026 -20 -14 1 Sr90 64 0 1 38.11 1.14 0.13 1.99 1.93 2.1 19.19
-29026 -20 -14 1 Sr90 64 0 2 39.48 0.06 1.74 1.94 1.88 2.04 20.34
-29026 -20 -14 1 Sr90 64 0 3 38.94 0.32 1.07 1.92 1.86 2.04 20.32
-29026 -20 -14 1 Sr90 64 0 4 33.19 0.21 0.55 1.89 1.83 2 17.6
-29026 -20 -14 1 Sr90 64 0 5 27.01 0 1.1 1.78 1.73 1.89 15.15
-
-29027 -20 -14.5 1 Sr90 80 0 0 44.64 3.6 5.02 2.03 1.94 2.19 21.97
-29027 -20 -14.5 1 Sr90 80 0 1 31.96 1.92 0.45 1.6 1.54 1.74 19.96
-29027 -20 -14.5 1 Sr90 80 0 2 48.54 1.15 2.38 2.04 1.94 2.2 23.81
-29027 -20 -14.5 1 Sr90 80 0 3 34.28 0.11 1.15 1.58 1.51 1.71 21.74
-29027 -20 -14.5 1 Sr90 80 0 4 31.74 0.23 1.32 1.58 1.51 1.72 20.13
-29027 -20 -14.5 1 Sr90 80 0 5 24.2 0.17 0.34 1.48 1.42 1.62 16.35
-
-29028 -70 -53.5 1 Sr90 64 0 0 52.7 5.41 3.53 2.22 2.13 2.37 23.73
-29028 -70 -53.5 1 Sr90 64 0 1 37.34 0.54 0.82 1.67 1.62 1.81 22.29
-29028 -70 -53.5 1 Sr90 64 0 2 38.97 0.05 0.99 1.65 1.59 1.77 23.55
-29028 -70 -53.5 1 Sr90 64 0 3 38.63 0.41 0.26 1.64 1.58 1.76 23.61
-29028 -70 -53.5 1 Sr90 64 0 4 32.87 0.14 0.53 1.66 1.6 1.78 19.79
-29028 -70 -53.5 1 Sr90 64 0 5 26.79 0.05 1.19 1.56 1.5 1.69 17.18
-
-29029 -70 -54.5 1 Sr90 80 0 0 45.68 4.76 8.52 2.1 2.01 2.25 21.74
-29029 -70 -54.5 1 Sr90 80 0 1 31.69 0.69 0.61 1.5 1.44 1.64 21.09
-29029 -70 -54.5 1 Sr90 80 0 2 47.78 1.1 2.64 2.06 1.96 2.23 23.16
-29029 -70 -54.5 1 Sr90 80 0 3 33.99 0.3 0.79 1.49 1.44 1.63 22.74
-29029 -70 -54.5 1 Sr90 80 0 4 31.16 0.25 1.68 1.5 1.44 1.66 20.73
-29029 -70 -54.5 1 Sr90 80 0 5 24.31 0.14 0.28 1.43 1.37 1.57 17.05
-
-29046 20 21.6 2 Sr90 64 150000 0 42.2 1.25 3.42 4.42 4.22 4.63 9.55
-29046 20 21.6 2 Sr90 64 150000 1 33.48 0.19 0.53 3.13 3.03 3.25 10.7
-29046 20 21.6 2 Sr90 64 150000 2 36.19 0.22 0.66 3.1 3 3.22 11.67
-29046 20 21.6 2 Sr90 64 150000 3 35.4 0.12 1.12 3.06 2.97 3.17 11.57
-29046 20 21.6 2 Sr90 64 150000 4 31 0.21 0.01 2.99 2.89 3.11 10.36
-29046 20 21.6 2 Sr90 64 150000 5 25.48 0.15 0.65 2.62 2.53 2.73 9.74
-
-29047 20 21 2 Sr90 80 150000 0 37.48 0.26 2.53 4.22 4.05 4.39 8.89
-29047 20 21 2 Sr90 80 150000 1 29.19 0.8 0.01 2.76 2.66 2.89 10.57
-29047 20 21 2 Sr90 80 150000 2 43.54 0.2 2.79 4.38 4.16 4.56 9.95
-29047 20 21 2 Sr90 80 150000 3 32.38 0.13 0.29 2.67 2.57 2.79 12.14
-29047 20 21 2 Sr90 80 150000 4 30.51 0.25 1.11 2.64 2.48 2.81 11.54
-29047 20 21 2 Sr90 80 150000 5 23.19 0.24 0.3 2.24 2.14 2.36 10.37
-
-29048 -20 -15.5 2 Sr90 64 150000 0 49.19 2.82 3.53 2.44 2.36 2.59 20.13
-29048 -20 -15.5 2 Sr90 64 150000 1 35.88 0.07 0.93 2.09 2.03 2.2 17.16
-29048 -20 -15.5 2 Sr90 64 150000 2 38.53 0.54 1.84 2.04 1.98 2.15 18.86
-29048 -20 -15.5 2 Sr90 64 150000 3 38.07 0.06 1.61 2.01 1.95 2.12 18.93
-29048 -20 -15.5 2 Sr90 64 150000 4 32.51 0.11 0.23 1.97 1.91 2.08 16.49
-29048 -20 -15.5 2 Sr90 64 150000 5 26.71 0.06 1.06 1.87 1.81 1.98 14.31
-
-29049 -20 -16 2 Sr90 80 150000 0 42.62 2.54 5.63 1.97 1.89 2.13 21.66
-29049 -20 -16 2 Sr90 80 150000 1 30.89 1.04 0.27 1.62 1.56 1.75 19.07
-29049 -20 -16 2 Sr90 80 150000 2 47.72 1.47 2.44 1.97 1.88 2.15 24.22
-29049 -20 -16 2 Sr90 80 150000 3 33.8 0.22 1.42 1.61 1.55 1.73 21
-29049 -20 -16 2 Sr90 80 150000 4 31.33 0.25 1.19 1.61 1.55 1.75 19.45
-29049 -20 -16 2 Sr90 80 150000 5 24.03 0.2 0.43 1.51 1.45 1.64 15.96
-
-29050 -70 -54.5 2 Sr90 64 150000 0 51.01 4.9 5.5 2.16 2.08 2.32 23.6
-29050 -70 -54.5 2 Sr90 64 150000 1 35.65 1.03 0.47 1.76 1.7 1.87 20.25
-29050 -70 -54.5 2 Sr90 64 150000 2 37.42 0.06 1.25 1.74 1.68 1.84 21.54
-29050 -70 -54.5 2 Sr90 64 150000 3 37.41 0.16 1.81 1.72 1.66 1.83 21.79
-29050 -70 -54.5 2 Sr90 64 150000 4 32.61 0.18 0.39 1.68 1.62 1.8 19.4
-29050 -70 -54.5 2 Sr90 64 150000 5 26.86 0.05 0.94 1.6 1.54 1.7 16.84
-
-29051 -70 -54.5 2 Sr90 80 150000 0 43.9 5.07 4.02 1.96 1.88 2.12 22.34
-29051 -70 -54.5 2 Sr90 80 150000 1 29 0.62 0.18 1.58 1.51 1.7 18.37
-29051 -70 -54.5 2 Sr90 80 150000 2 47.21 1.29 2.54 1.94 1.86 2.11 24.29
-29051 -70 -54.5 2 Sr90 80 150000 3 32.68 0.06 0.93 1.58 1.51 1.71 20.66
-29051 -70 -54.5 2 Sr90 80 150000 4 30.5 0.23 0.92 1.61 1.53 1.74 18.97
-29051 -70 -54.5 2 Sr90 80 150000 5 24.1 0.22 0.38 1.42 1.36 1.54 16.95
-
-29064 20 22 4 Sr90 64 1000000000000 0 35.1 0.92 2.24 4.21 3.99 4.39 8.33
-29064 20 22 4 Sr90 64 1000000000000 1 31.95 0.2 0.32 2.66 2.54 2.8 12.02
-29064 20 22 4 Sr90 64 1000000000000 2 33 0.26 0.77 2.75 2.64 2.87 12.02
-29064 20 22 4 Sr90 64 1000000000000 3 32.46 0.13 0.51 2.78 2.63 2.91 11.7
-29064 20 22 4 Sr90 64 1000000000000 4 25.38 0.23 0.85 2.95 2.84 3.07 8.6
-29064 20 22 4 Sr90 64 1000000000000 5 23.11 0.5 0.71 2.36 2.27 2.46 9.81
-
-29065 20 21.5 4 Sr90 80 1000000000000 0 32.1 2.47 0.4 4.06 3.9 4.23 7.91
-29065 20 21.5 4 Sr90 80 1000000000000 1 27.35 0.22 0.1 2.66 2.54 2.79 10.29
-29065 20 21.5 4 Sr90 80 1000000000000 2 32.97 5.27 3.6 4.73 4.55 4.89 6.96
-29065 20 21.5 4 Sr90 80 1000000000000 3 27.39 0.1 0.63 2.75 2.65 2.88 9.96
-29065 20 21.5 4 Sr90 80 1000000000000 4 23.16 0.62 0.76 2.83 2.73 2.97 8.18
-29065 20 21.5 4 Sr90 80 1000000000000 5 20.34 0.3 0 2.33 2.25 2.43 8.74
-
-29066 -20 -13.5 4 Sr90 64 1000000000000 0 41.62 2.88 1.34 1.96 1.87 2.13 21.28
-29066 -20 -13.5 4 Sr90 64 1000000000000 1 35.52 0.57 0.28 1.55 1.49 1.7 22.86
-29066 -20 -13.5 4 Sr90 64 1000000000000 2 35.27 0.27 1.71 1.54 1.48 1.7 22.89
-29066 -20 -13.5 4 Sr90 64 1000000000000 3 34.85 0.17 0.75 1.54 1.48 1.7 22.6
-29066 -20 -13.5 4 Sr90 64 1000000000000 4 26.11 0.08 0.02 1.53 1.47 1.69 17.03
-29066 -20 -13.5 4 Sr90 64 1000000000000 5 24.71 0.36 0.83 1.45 1.39 1.59 17.03
-
-29067 -20 -14 4 Sr90 80 1000000000000 0 35.9 3.11 0.34 1.94 1.85 2.09 18.55
-29067 -20 -14 4 Sr90 80 1000000000000 1 29.46 0.86 0.19 1.55 1.48 1.68 19.06
-29067 -20 -14 4 Sr90 80 1000000000000 2 32.05 4.79 0.84 1.91 1.82 2.06 16.78
-29067 -20 -14 4 Sr90 80 1000000000000 3 27.22 0.16 1 1.53 1.47 1.66 17.74
-29067 -20 -14 4 Sr90 80 1000000000000 4 21.14 2.99 0.79 1.52 1.47 1.67 13.87
-29067 -20 -14 4 Sr90 80 1000000000000 5 20.69 0.22 0 1.45 1.4 1.56 14.3
-
-29068 -70 -53 4 Sr90 64 1000000000000 0 39.93 3.47 0.23 2.02 1.92 2.19 19.76
-29068 -70 -53 4 Sr90 64 1000000000000 1 35.13 0.43 0.36 1.5 1.43 1.65 23.48
-29068 -70 -53 4 Sr90 64 1000000000000 2 34.36 0.46 0.18 1.48 1.42 1.65 23.15
-29068 -70 -53 4 Sr90 64 1000000000000 3 33.18 0.07 0.69 1.48 1.42 1.65 22.42
-29068 -70 -53 4 Sr90 64 1000000000000 4 23.01 0.39 0.49 1.48 1.42 1.65 15.51
-29068 -70 -53 4 Sr90 64 1000000000000 5 23.79 0.24 0.46 1.4 1.34 1.53 17.01
-
-29069 -70 -54 4 Sr90 80 1000000000000 0 32.71 3.82 1.25 2.02 1.92 2.17 16.21
-29069 -70 -54 4 Sr90 80 1000000000000 1 28.96 0.36 0.38 1.48 1.42 1.63 19.6
-29069 -70 -54 4 Sr90 80 1000000000000 2 23.83 2.63 2.31 1.91 1.82 2.08 12.46
-29069 -70 -54 4 Sr90 80 1000000000000 3 23.99 0.06 0.29 1.46 1.4 1.62 16.45
-29069 -70 -54 4 Sr90 80 1000000000000 4 17.17 0 0.58 1.46 1.4 1.63 11.73
-29069 -70 -54 4 Sr90 80 1000000000000 5 19.09 0.87 0 1.39 1.34 1.51 13.7
+++ /dev/null
-\usepackage{sansmath}
-\usepackage[scaled=0.92]{helvet}
+++ /dev/null
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-
-
-# block index 0
-29024 20 21.5 1 Sr90 64 0 0 45.85 2.25 1.59 4.08 3.89 4.27 11.24
-29024 20 21.5 1 Sr90 64 0 1 35.28 0.57 0.02 2.84 2.74 2.98 12.41
-29024 20 21.5 1 Sr90 64 0 2 36.95 0.2 1.19 2.85 2.74 2.97 12.98
-29024 20 21.5 1 Sr90 64 0 3 36.65 0.14 0.67 2.83 2.71 2.95 12.96
-29024 20 21.5 1 Sr90 64 0 4 31.18 0.07 0.19 2.82 2.72 2.96 11.05
-29024 20 21.5 1 Sr90 64 0 5 25.68 0.09 0.43 2.44 2.35 2.57 10.54
-
-
-# block index 1
-29026 -20 -14 1 Sr90 64 0 0 51.81 5.6 6.65 2.37 2.29 2.51 21.82
-29026 -20 -14 1 Sr90 64 0 1 38.11 1.14 0.13 1.99 1.93 2.1 19.19
-29026 -20 -14 1 Sr90 64 0 2 39.48 0.06 1.74 1.94 1.88 2.04 20.34
-29026 -20 -14 1 Sr90 64 0 3 38.94 0.32 1.07 1.92 1.86 2.04 20.32
-29026 -20 -14 1 Sr90 64 0 4 33.19 0.21 0.55 1.89 1.83 2 17.6
-29026 -20 -14 1 Sr90 64 0 5 27.01 0 1.1 1.78 1.73 1.89 15.15
-
-
-# block index 2
-29028 -70 -53.5 1 Sr90 64 0 0 52.7 5.41 3.53 2.22 2.13 2.37 23.73
-29028 -70 -53.5 1 Sr90 64 0 1 37.34 0.54 0.82 1.67 1.62 1.81 22.29
-29028 -70 -53.5 1 Sr90 64 0 2 38.97 0.05 0.99 1.65 1.59 1.77 23.55
-29028 -70 -53.5 1 Sr90 64 0 3 38.63 0.41 0.26 1.64 1.58 1.76 23.61
-29028 -70 -53.5 1 Sr90 64 0 4 32.87 0.14 0.53 1.66 1.6 1.78 19.79
-29028 -70 -53.5 1 Sr90 64 0 5 26.79 0.05 1.19 1.56 1.5 1.69 17.18
+++ /dev/null
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-
-
-# block index 0
-29024 20 21.5 1 Sr90 64 0 0 45.85 2.25 1.59 4.08 3.89 4.27 11.24
-29024 20 21.5 1 Sr90 64 0 1 35.28 0.57 0.02 2.84 2.74 2.98 12.41
-29024 20 21.5 1 Sr90 64 0 2 36.95 0.2 1.19 2.85 2.74 2.97 12.98
-29024 20 21.5 1 Sr90 64 0 3 36.65 0.14 0.67 2.83 2.71 2.95 12.96
-29024 20 21.5 1 Sr90 64 0 4 31.18 0.07 0.19 2.82 2.72 2.96 11.05
-29024 20 21.5 1 Sr90 64 0 5 25.68 0.09 0.43 2.44 2.35 2.57 10.54
-
-
-# block index 1
-29026 -20 -14 1 Sr90 64 0 0 51.81 5.6 6.65 2.37 2.29 2.51 21.82
-29026 -20 -14 1 Sr90 64 0 1 38.11 1.14 0.13 1.99 1.93 2.1 19.19
-29026 -20 -14 1 Sr90 64 0 2 39.48 0.06 1.74 1.94 1.88 2.04 20.34
-29026 -20 -14 1 Sr90 64 0 3 38.94 0.32 1.07 1.92 1.86 2.04 20.32
-29026 -20 -14 1 Sr90 64 0 4 33.19 0.21 0.55 1.89 1.83 2 17.6
-29026 -20 -14 1 Sr90 64 0 5 27.01 0 1.1 1.78 1.73 1.89 15.15
-
-
-# block index 2
-29028 -70 -53.5 1 Sr90 64 0 0 52.7 5.41 3.53 2.22 2.13 2.37 23.73
-29028 -70 -53.5 1 Sr90 64 0 1 37.34 0.54 0.82 1.67 1.62 1.81 22.29
-29028 -70 -53.5 1 Sr90 64 0 2 38.97 0.05 0.99 1.65 1.59 1.77 23.55
-29028 -70 -53.5 1 Sr90 64 0 3 38.63 0.41 0.26 1.64 1.58 1.76 23.61
-29028 -70 -53.5 1 Sr90 64 0 4 32.87 0.14 0.53 1.66 1.6 1.78 19.79
-29028 -70 -53.5 1 Sr90 64 0 5 26.79 0.05 1.19 1.56 1.5 1.69 17.18
+++ /dev/null
-%!PS-Adobe-2.0 EPSF-2.0
-%%Title: ./output/Title.tex
-%%Creator: gnuplot 4.6 patchlevel 3
-%%CreationDate: Wed Jan 8 21:44:04 2014
-%%DocumentFonts:
-%%BoundingBox: 50 50 410 302
-%%EndComments
-%%BeginProlog
-/gnudict 256 dict def
-gnudict begin
-%
-% The following true/false flags may be edited by hand if desired.
-% The unit line width and grayscale image gamma correction may also be changed.
-%
-/Color true def
-/Blacktext false def
-/Solid true def
-/Dashlength 1 def
-/Landscape false def
-/Level1 false def
-/Rounded false def
-/ClipToBoundingBox false def
-/SuppressPDFMark false def
-/TransparentPatterns false def
-/gnulinewidth 5.000 def
-/userlinewidth gnulinewidth def
-/Gamma 1.0 def
-/BackgroundColor {-1.000 -1.000 -1.000} def
-%
-/vshift -73 def
-/dl1 {
- 10.0 Dashlength mul mul
- Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
-} def
-/dl2 {
- 10.0 Dashlength mul mul
- Rounded { currentlinewidth 0.75 mul add } if
-} def
-/hpt_ 31.5 def
-/vpt_ 31.5 def
-/hpt hpt_ def
-/vpt vpt_ def
-/doclip {
- ClipToBoundingBox {
- newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
- clip
- } if
-} def
-%
-% Gnuplot Prolog Version 4.6 (September 2012)
-%
-%/SuppressPDFMark true def
-%
-/M {moveto} bind def
-/L {lineto} bind def
-/R {rmoveto} bind def
-/V {rlineto} bind def
-/N {newpath moveto} bind def
-/Z {closepath} bind def
-/C {setrgbcolor} bind def
-/f {rlineto fill} bind def
-/g {setgray} bind def
-/Gshow {show} def % May be redefined later in the file to support UTF-8
-/vpt2 vpt 2 mul def
-/hpt2 hpt 2 mul def
-/Lshow {currentpoint stroke M 0 vshift R
- Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
- Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R
- Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
- /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
-/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
- {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
-/BL {stroke userlinewidth 2 mul setlinewidth
- Rounded {1 setlinejoin 1 setlinecap} if} def
-/AL {stroke userlinewidth 2 div setlinewidth
- Rounded {1 setlinejoin 1 setlinecap} if} def
-/UL {dup gnulinewidth mul /userlinewidth exch def
- dup 1 lt {pop 1} if 10 mul /udl exch def} def
-/PL {stroke userlinewidth setlinewidth
- Rounded {1 setlinejoin 1 setlinecap} if} def
-3.8 setmiterlimit
-% Default Line colors
-/LCw {1 1 1} def
-/LCb {0 0 0} def
-/LCa {0 0 0} def
-/LC0 {1 0 0} def
-/LC1 {0 1 0} def
-/LC2 {0 0 1} def
-/LC3 {1 0 1} def
-/LC4 {0 1 1} def
-/LC5 {1 1 0} def
-/LC6 {0 0 0} def
-/LC7 {1 0.3 0} def
-/LC8 {0.5 0.5 0.5} def
-% Default Line Types
-/LTw {PL [] 1 setgray} def
-/LTb {BL [] LCb DL} def
-/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
-/LT0 {PL [] LC0 DL} def
-/LT1 {PL [4 dl1 2 dl2] LC1 DL} def
-/LT2 {PL [2 dl1 3 dl2] LC2 DL} def
-/LT3 {PL [1 dl1 1.5 dl2] LC3 DL} def
-/LT4 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
-/LT5 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC5 DL} def
-/LT6 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC6 DL} def
-/LT7 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC7 DL} def
-/LT8 {PL [2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 4 dl2] LC8 DL} def
-/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
-/Dia {stroke [] 0 setdash 2 copy vpt add M
- hpt neg vpt neg V hpt vpt neg V
- hpt vpt V hpt neg vpt V closepath stroke
- Pnt} def
-/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
- currentpoint stroke M
- hpt neg vpt neg R hpt2 0 V stroke
- } def
-/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
- 0 vpt2 neg V hpt2 0 V 0 vpt2 V
- hpt2 neg 0 V closepath stroke
- Pnt} def
-/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
- hpt2 vpt2 neg V currentpoint stroke M
- hpt2 neg 0 R hpt2 vpt2 V stroke} def
-/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
- hpt neg vpt -1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt 1.62 mul V closepath stroke
- Pnt} def
-/Star {2 copy Pls Crs} def
-/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
- 0 vpt2 neg V hpt2 0 V 0 vpt2 V
- hpt2 neg 0 V closepath fill} def
-/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
- hpt neg vpt -1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt 1.62 mul V closepath fill} def
-/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
- hpt neg vpt 1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt -1.62 mul V closepath stroke
- Pnt} def
-/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
- hpt neg vpt 1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt -1.62 mul V closepath fill} def
-/DiaF {stroke [] 0 setdash vpt add M
- hpt neg vpt neg V hpt vpt neg V
- hpt vpt V hpt neg vpt V closepath fill} def
-/Pent {stroke [] 0 setdash 2 copy gsave
- translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
- closepath stroke grestore Pnt} def
-/PentF {stroke [] 0 setdash gsave
- translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
- closepath fill grestore} def
-/Circle {stroke [] 0 setdash 2 copy
- hpt 0 360 arc stroke Pnt} def
-/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
-/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
-/C1 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 90 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C2 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 90 180 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C3 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 180 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C4 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 180 270 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C5 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 90 arc
- 2 copy moveto
- 2 copy vpt 180 270 arc closepath fill
- vpt 0 360 arc} bind def
-/C6 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 90 270 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C7 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 270 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C8 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 270 360 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C9 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 270 450 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
- 2 copy moveto
- 2 copy vpt 90 180 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C11 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 180 arc closepath fill
- 2 copy moveto
- 2 copy vpt 270 360 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C12 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 180 360 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C13 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 0 90 arc closepath fill
- 2 copy moveto
- 2 copy vpt 180 360 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/C14 {BL [] 0 setdash 2 copy moveto
- 2 copy vpt 90 360 arc closepath fill
- vpt 0 360 arc} bind def
-/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
- vpt 0 360 arc closepath} bind def
-/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
- neg 0 rlineto closepath} bind def
-/Square {dup Rec} bind def
-/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
-/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
-/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
-/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
-/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
- exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
- 2 copy vpt Square fill Bsquare} bind def
-/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
-/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
- Bsquare} bind def
-/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
- Bsquare} bind def
-/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
-/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
- 2 copy vpt Square fill Bsquare} bind def
-/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
- 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
-/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
-/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
-/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
-/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
-/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
-/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
-/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
-/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
-/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
-/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
-/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
-/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
-/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
-/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
-/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
-/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
-/DiaE {stroke [] 0 setdash vpt add M
- hpt neg vpt neg V hpt vpt neg V
- hpt vpt V hpt neg vpt V closepath stroke} def
-/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
- 0 vpt2 neg V hpt2 0 V 0 vpt2 V
- hpt2 neg 0 V closepath stroke} def
-/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
- hpt neg vpt -1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt 1.62 mul V closepath stroke} def
-/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
- hpt neg vpt 1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt -1.62 mul V closepath stroke} def
-/PentE {stroke [] 0 setdash gsave
- translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
- closepath stroke grestore} def
-/CircE {stroke [] 0 setdash
- hpt 0 360 arc stroke} def
-/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
-/DiaW {stroke [] 0 setdash vpt add M
- hpt neg vpt neg V hpt vpt neg V
- hpt vpt V hpt neg vpt V Opaque stroke} def
-/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
- 0 vpt2 neg V hpt2 0 V 0 vpt2 V
- hpt2 neg 0 V Opaque stroke} def
-/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
- hpt neg vpt -1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt 1.62 mul V Opaque stroke} def
-/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
- hpt neg vpt 1.62 mul V
- hpt 2 mul 0 V
- hpt neg vpt -1.62 mul V Opaque stroke} def
-/PentW {stroke [] 0 setdash gsave
- translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
- Opaque stroke grestore} def
-/CircW {stroke [] 0 setdash
- hpt 0 360 arc Opaque stroke} def
-/BoxFill {gsave Rec 1 setgray fill grestore} def
-/Density {
- /Fillden exch def
- currentrgbcolor
- /ColB exch def /ColG exch def /ColR exch def
- /ColR ColR Fillden mul Fillden sub 1 add def
- /ColG ColG Fillden mul Fillden sub 1 add def
- /ColB ColB Fillden mul Fillden sub 1 add def
- ColR ColG ColB setrgbcolor} def
-/BoxColFill {gsave Rec PolyFill} def
-/PolyFill {gsave Density fill grestore grestore} def
-/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
-%
-% PostScript Level 1 Pattern Fill routine for rectangles
-% Usage: x y w h s a XX PatternFill
-% x,y = lower left corner of box to be filled
-% w,h = width and height of box
-% a = angle in degrees between lines and x-axis
-% XX = 0/1 for no/yes cross-hatch
-%
-/PatternFill {gsave /PFa [ 9 2 roll ] def
- PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
- PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
- TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
- clip
- currentlinewidth 0.5 mul setlinewidth
- /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
- 0 0 M PFa 5 get rotate PFs -2 div dup translate
- 0 1 PFs PFa 4 get div 1 add floor cvi
- {PFa 4 get mul 0 M 0 PFs V} for
- 0 PFa 6 get ne {
- 0 1 PFs PFa 4 get div 1 add floor cvi
- {PFa 4 get mul 0 2 1 roll M PFs 0 V} for
- } if
- stroke grestore} def
-%
-/languagelevel where
- {pop languagelevel} {1} ifelse
- 2 lt
- {/InterpretLevel1 true def}
- {/InterpretLevel1 Level1 def}
- ifelse
-%
-% PostScript level 2 pattern fill definitions
-%
-/Level2PatternFill {
-/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
- bind def
-/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke}
->> matrix makepattern
-/Pat1 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
- 0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
->> matrix makepattern
-/Pat2 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
- 8 8 L 8 0 L 0 0 L fill}
->> matrix makepattern
-/Pat3 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
- 0 12 M 12 0 L stroke}
->> matrix makepattern
-/Pat4 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
- 0 -4 M 12 8 L stroke}
->> matrix makepattern
-/Pat5 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
- 0 12 M 8 -4 L 4 12 M 10 0 L stroke}
->> matrix makepattern
-/Pat6 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
- 0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
->> matrix makepattern
-/Pat7 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
- 12 0 M -4 8 L 12 4 M 0 10 L stroke}
->> matrix makepattern
-/Pat8 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
- -4 0 M 12 8 L -4 4 M 8 10 L stroke}
->> matrix makepattern
-/Pat9 exch def
-/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
-/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
-/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
-/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
-/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
-/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
-/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
-} def
-%
-%
-%End of PostScript Level 2 code
-%
-/PatternBgnd {
- TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
-} def
-%
-% Substitute for Level 2 pattern fill codes with
-% grayscale if Level 2 support is not selected.
-%
-/Level1PatternFill {
-/Pattern1 {0.250 Density} bind def
-/Pattern2 {0.500 Density} bind def
-/Pattern3 {0.750 Density} bind def
-/Pattern4 {0.125 Density} bind def
-/Pattern5 {0.375 Density} bind def
-/Pattern6 {0.625 Density} bind def
-/Pattern7 {0.875 Density} bind def
-} def
-%
-% Now test for support of Level 2 code
-%
-Level1 {Level1PatternFill} {Level2PatternFill} ifelse
-%
-/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
-dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
-currentdict end definefont pop
-Level1 SuppressPDFMark or
-{} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
- userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
- /Title (./output/Title.tex)
- /Subject (gnuplot plot)
- /Creator (gnuplot 4.6 patchlevel 3)
- /Author (stefan)
-% /Producer (gnuplot)
-% /Keywords ()
- /CreationDate (Wed Jan 8 21:44:04 2014)
- /DOCINFO pdfmark
-end
-} ifelse
-end
-%%EndProlog
-%%Page: 1 1
-gnudict begin
-gsave
-doclip
-50 50 translate
-0.050 0.050 scale
-0 setgray
-newpath
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 704 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 704 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 1286 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 1286 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 1867 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 1867 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 2449 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 2449 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 3030 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 3030 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 3612 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 3612 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 4193 M
-4147 0 V
-1578 0 R
-132 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 4193 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 4775 M
-5857 0 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 4775 M
-63 0 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 946 704 M
-0 4071 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 2117 704 M
-0 4071 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 2117 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 3289 704 M
-0 4071 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 3289 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 4460 704 M
-0 4071 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 4460 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 5632 704 M
-0 3348 V
-0 660 R
-0 63 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 5632 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.000 UL
-LTa
-LCa setrgbcolor
-0.84 0.88 0.79 C 6803 704 M
-0 4071 V
-stroke
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 6803 704 M
-0 63 V
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C 946 4775 M
-946 704 L
-5857 0 V
-0 4071 R
--5857 0 R
-stroke
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C LTb
-LCb setrgbcolor
-LTb
-2.000 UP
-1.500 UL
-LTb
-0.00 0.00 0.00 C 2.000 UP
-3.000 UL
-LT0
-LC0 setrgbcolor
-0.45 0.59 0.90 C LCb setrgbcolor
-LT0
-LC0 setrgbcolor
-0.45 0.59 0.90 C 5885 4602 M
-654 0 V
-946 1065 M
-1171 340 V
-1172 166 V
-1171 -6 V
-5632 1009 L
-6803 861 L
-946 1065 CircleF
-2117 1405 CircleF
-3289 1571 CircleF
-4460 1565 CircleF
-5632 1009 CircleF
-6803 861 CircleF
-6212 4602 CircleF
-2.000 UP
-3.000 UL
-LT1
-LC1 setrgbcolor
-0.40 0.92 0.52 C LCb setrgbcolor
-LT1
-LC1 setrgbcolor
-0.40 0.92 0.52 C 5885 4382 M
-654 0 V
-946 4141 M
-2117 3376 L
-1172 335 V
-1171 -6 V
-5632 2914 L
-6803 2202 L
-946 4141 BoxF
-2117 3376 BoxF
-3289 3711 BoxF
-4460 3705 BoxF
-5632 2914 BoxF
-6803 2202 BoxF
-6212 4382 BoxF
-2.000 UP
-3.000 UL
-LT2
-LC2 setrgbcolor
-0.98 0.48 0.43 C LCb setrgbcolor
-LT2
-LC2 setrgbcolor
-0.98 0.48 0.43 C 5885 4162 M
-654 0 V
-946 4696 M
-2117 4278 L
-1172 366 V
-1171 18 V
-5632 3551 L
-6803 2792 L
-946 4696 TriUF
-2117 4278 TriUF
-3289 4644 TriUF
-4460 4662 TriUF
-5632 3551 TriUF
-6803 2792 TriUF
-6212 4162 TriUF
-2.000 UP
-1.500 UL
-LTb
-LCb setrgbcolor
-0.00 0.00 0.00 C stroke
-grestore
-end
-showpage
-%%Trailer
+++ /dev/null
-% GNUPLOT: LaTeX picture with Postscript
-\documentclass{minimal}
-% Set font size
-\makeatletter
-\def\@ptsize{1}
-\InputIfFileExists{size11.clo}{}{%
- \GenericError{(gnuplot) \space\space\space\@spaces}{%
- Gnuplot Error: File `size11.clo' not found! Could not set font size%
- }{See the gnuplot documentation for explanation.%
- }{For using a font size a file `size<fontsize>.clo' has to exist.
- Falling back ^^Jto default fontsize 10pt.}%
- \def\@ptsize{0}
- \input{size10.clo}%
-}%
-\makeatother
-\renewcommand*\rmdefault{phv}%
-% Load packages
-\usepackage{graphicx}
-\usepackage{color}
-\makeatletter
-% Select an appropriate default driver (from TeXLive graphics.cfg)
-\begingroup
- \chardef\x=0 %
- % check pdfTeX
- \@ifundefined{pdfoutput}{}{%
- \ifcase\pdfoutput
- \else
- \chardef\x=1 %
- \fi
- }%
- % check VTeX
- \@ifundefined{OpMode}{}{%
- \chardef\x=2 %
- }%
-\expandafter\endgroup
-\ifcase\x
- % default case
- \PassOptionsToPackage{dvips}{geometry}
-\or
- % pdfTeX is running in pdf mode
- \PassOptionsToPackage{pdftex}{geometry}
-\else
- % VTeX is running
- \PassOptionsToPackage{vtex}{geometry}
-\fi
-\makeatother
-% Set papersize
-\usepackage[papersize={360.00bp,252.00bp},text={360.00bp,252.00bp}]{geometry}
-% No page numbers and no paragraph indentation
-\pagestyle{empty}
-\setlength{\parindent}{0bp}%
-% Load configuration file
-\InputIfFileExists{gnuplot.cfg}{%
- \typeout{Using configuration file gnuplot.cfg}%
-}{%
- \typeout{No configuration file gnuplot.cfg found.}%
-}%
-\usepackage[pdftex, pdfauthor={Stefan Strohauer}, pdftitle={Title}, pdfsubject={MIMOSA-29 radiation hardness}, pdfkeywords={MIMOSA-29, radiation hardness}, pdfproducer={Latex with hyperref}, pdfcreator={pdflatex}]{hyperref}
-\usepackage[utf8]{inputenc}
-\usepackage[T1]{fontenc}
-\usepackage{siunitx}
-\sisetup{load-configurations = abbreviations,separate-uncertainty}
-\usepackage{amsmath}
-\usepackage[helvet]{sfmath}
-\begin{document}
-\begingroup
- \makeatletter
- \providecommand\color[2][]{%
- \GenericError{(gnuplot) \space\space\space\@spaces}{%
- Package color not loaded in conjunction with
- terminal option `colourtext'%
- }{See the gnuplot documentation for explanation.%
- }{Either use 'blacktext' in gnuplot or load the package
- color.sty in LaTeX.}%
- \renewcommand\color[2][]{}%
- }%
- \providecommand\includegraphics[2][]{%
- \GenericError{(gnuplot) \space\space\space\@spaces}{%
- Package graphicx or graphics not loaded%
- }{See the gnuplot documentation for explanation.%
- }{The gnuplot epslatex terminal needs graphicx.sty or graphics.sty.}%
- \renewcommand\includegraphics[2][]{}%
- }%
- \providecommand\rotatebox[2]{#2}%
- \@ifundefined{ifGPcolor}{%
- \newif\ifGPcolor
- \GPcolortrue
- }{}%
- \@ifundefined{ifGPblacktext}{%
- \newif\ifGPblacktext
- \GPblacktextfalse
- }{}%
- % define a \g@addto@macro without @ in the name:
- \let\gplgaddtomacro\g@addto@macro
- % define empty templates for all commands taking text:
- \gdef\gplbacktext{}%
- \gdef\gplfronttext{}%
- \makeatother
- \ifGPblacktext
- % no textcolor at all
- \def\colorrgb#1{}%
- \def\colorgray#1{}%
- \else
- % gray or color?
- \ifGPcolor
- \def\colorrgb#1{\color[rgb]{#1}}%
- \def\colorgray#1{\color[gray]{#1}}%
- \expandafter\def\csname LTw\endcsname{\color{white}}%
- \expandafter\def\csname LTb\endcsname{\color{black}}%
- \expandafter\def\csname LTa\endcsname{\color{black}}%
- \expandafter\def\csname LT0\endcsname{\color[rgb]{1,0,0}}%
- \expandafter\def\csname LT1\endcsname{\color[rgb]{0,1,0}}%
- \expandafter\def\csname LT2\endcsname{\color[rgb]{0,0,1}}%
- \expandafter\def\csname LT3\endcsname{\color[rgb]{1,0,1}}%
- \expandafter\def\csname LT4\endcsname{\color[rgb]{0,1,1}}%
- \expandafter\def\csname LT5\endcsname{\color[rgb]{1,1,0}}%
- \expandafter\def\csname LT6\endcsname{\color[rgb]{0,0,0}}%
- \expandafter\def\csname LT7\endcsname{\color[rgb]{1,0.3,0}}%
- \expandafter\def\csname LT8\endcsname{\color[rgb]{0.5,0.5,0.5}}%
- \else
- % gray
- \def\colorrgb#1{\color{black}}%
- \def\colorgray#1{\color[gray]{#1}}%
- \expandafter\def\csname LTw\endcsname{\color{white}}%
- \expandafter\def\csname LTb\endcsname{\color{black}}%
- \expandafter\def\csname LTa\endcsname{\color{black}}%
- \expandafter\def\csname LT0\endcsname{\color{black}}%
- \expandafter\def\csname LT1\endcsname{\color{black}}%
- \expandafter\def\csname LT2\endcsname{\color{black}}%
- \expandafter\def\csname LT3\endcsname{\color{black}}%
- \expandafter\def\csname LT4\endcsname{\color{black}}%
- \expandafter\def\csname LT5\endcsname{\color{black}}%
- \expandafter\def\csname LT6\endcsname{\color{black}}%
- \expandafter\def\csname LT7\endcsname{\color{black}}%
- \expandafter\def\csname LT8\endcsname{\color{black}}%
- \fi
- \fi
- \setlength{\unitlength}{0.0500bp}%
- \begin{picture}(7200.00,5040.00)%
- \gplgaddtomacro\gplbacktext{%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,704){\makebox(0,0)[r]{\strut{} 10}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,1286){\makebox(0,0)[r]{\strut{} 12}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,1867){\makebox(0,0)[r]{\strut{} 14}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,2449){\makebox(0,0)[r]{\strut{} 16}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,3030){\makebox(0,0)[r]{\strut{} 18}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,3612){\makebox(0,0)[r]{\strut{} 20}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,4193){\makebox(0,0)[r]{\strut{} 22}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(814,4775){\makebox(0,0)[r]{\strut{} 24}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(946,484){\makebox(0,0){\strut{} 0}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(2117,484){\makebox(0,0){\strut{} 1}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(3289,484){\makebox(0,0){\strut{} 2}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(4460,484){\makebox(0,0){\strut{} 3}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(5632,484){\makebox(0,0){\strut{} 4}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(6803,484){\makebox(0,0){\strut{} 5}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(176,2739){\rotatebox{-270}{\makebox(0,0){\strut{}S/N}}}%
- \colorrgb{0.00,0.00,0.00}%
- \put(3874,154){\makebox(0,0){\strut{}Rund ist der Fußball, aber \SI{50(3)e2}{\celsius \mV \ohm} $\frac{3}{5}$}}%
- \csname LTb\endcsname%
- \put(3874,4665){\makebox(0,0){\strut{}}}%
- }%
- \gplgaddtomacro\gplfronttext{%
- \csname LTb\endcsname%
- \put(5753,4602){\makebox(0,0)[r]{\strut{}\SI{20}{\celsius}}}%
- \csname LTb\endcsname%
- \put(5753,4382){\makebox(0,0)[r]{\strut{}-20 C}}%
- \csname LTb\endcsname%
- \put(5753,4162){\makebox(0,0)[r]{\strut{}-70 C}}%
- }%
- \gplbacktext
- \put(0,0){\includegraphics{./output/Title-inc}}%
- \gplfronttext
- \end{picture}%
-\endgroup
-\end{document}
+++ /dev/null
-#!/bin/bash
-
-file=$1
-
-tempfile="lastplot.tmp"
-
-
-colRunNo=1
-colT=2
-colTsens=3
-colChip=4
-colSource=5
-colMatrix=6
-colRadDose=7
-colSubMatrix=8
-colMPV=9
-colErrLeft=10
-colErrRight=11
-colNoise=12
-colNoiseLow17=13
-colNoiseHi17=14
-colSN=15
-
-colX=$colSubMatrix
-colY=$colSN
-groupBy=$colT
-pdfName="Title"
-plotTitle=''
-# xLabel='Matrix'
-# yLabel='S/N'
-xLabel='Rund ist der Fußball, aber \\SI{50(3)e2}{\\celsius \\mV \\ohm} $\\frac{3}{5}$'
-yLabel='S/N'
-# titles=('Submatrix 0' 'Submatrix 1' 'Submatrix 2' 'Submatrix 3' 'Submatrix 4' 'Submatrix 5')
-titles=("\\SI{20}{\celsius}" '-20 C' '-70 C')
-
-selRunNo=''
-# selT='20'
-selTsens=''
-selChip='1'
-selSource='S\S*'
-selMatrix='64'
-selRadDose=''
-selSubMatrix=''
-selMPV=''
-selErrLeft=''
-selErrRight=''
-selNoise=''
-selNoiseLow17=''
-selNoiseHi17=''
-selSN=''
-
-mkdir ./output
-mkdir ./plots
-
-##################
-# select only matching data points
-regexp='^' # new line character
-
-for i in selRunNo selT selTsens selChip selSource selMatrix selRadDose selSubMatrix selMPV selErrLeft selErrRight selNoise selNoiseLow17 selNoiseHi17 selSN
-do
- eval sel=\$$i # set variable sel to the string which is contained in the selected i
- if [ -z $sel ]; then
- pattern='\S*' # if length of string in selector is zero, accept every string in this column ('\S*' is a regular expression for any number of non-whitespace characters)
- else
- pattern=$sel # if length of string in selector is non-zero, accept only matching strings
- fi
-
- regexp=$regexp$pattern"\s+" # pattern plus at least one whitespace character
-done
-
-cat $file | egrep "^#" > $tempfile # overwrite tempfile
-cat $file | egrep "$regexp" >> $tempfile # append to tempfile
-
-
-
-##################
-# split matching data sets in groups (criterion: all elements of one group belong to the same submatrix)
-memory=""
-index="0"
-
-if [ -n "$groupBy" ]; then # is grouping non zero?
- tmpA=$(mktemp) # create temporary file :-)
- egrep "^#" $tempfile > $tmpA # write comment lines to temp file
- while IFS="" read line
- do
- newEntry=$(echo "$line" | cut -f $groupBy) # get entry of column $groupBy of line $line
- if [ "$memory" != "$newEntry" ]; then # number in column $groupBy does not match the previous one
- echo >> $tmpA # one blank line for splitting the data sets
- echo >> $tmpA # a second blank line for seperate linestyles
- echo "# block index $index" >> $tmpA
- index=$[ $index + 1 ]
- fi
- echo "$line" >> $tmpA
- memory=$newEntry
- done < <(egrep -v "^#" $tempfile | sort -k$groupBy,$groupBy -k$colT,$colT) # input for the while loop: sorted data file without comment lines. -k3,5 means "take cloumn 3 to 5 as one 'column to sort by'". normally one wants to take only one column as entity to sort by, so usally one chooses the option -k5,5 to sort by column 5. sort -k3,3 -k5,5 means sorting first by column 3 and then subsort by column 5.
- cat $tmpA > $tempfile # write final data file
- rm $tmpA
-fi
-
-
-##################
-# plotting
-echo
-echo
-echo making a plot with the following data:
-echo
-cat $tempfile # print tempfile for information
-echo
-echo
-echo starting gnuplot ...
-
-plotcmd="plot '"./$tempfile"' using $colX:$colY index 0 w lp title '${titles[0]}'"
-for i in $(seq 1 $[ $index - 1 ]); do
-plotcmd=$plotcmd", '' using $colX:$colY index $i w lp title '${titles[$i]}'"
-done
-
-echo
-echo "Plotting with plot-command:"
-echo $plotcmd
-echo
-
-cat <<EOF | gnuplot
-# 'phv,11' is helvetica, size 11. header sets the following string right before the \begin{document} in the .tex file
-# possibly easier by applying the "sed" command onto the resulting .tex file and replacing \begin{document} by \usepackage{abc}\nbegin{document}
-set terminal epslatex standalone solid color colortext 'phv,11' header \
- $(echo '"\\usepackage[pdftex,'\
- "pdfauthor={Stefan Strohauer},\
- pdftitle={$pdfName},\
- pdfsubject={MIMOSA-29 radiation hardness},\
- pdfkeywords={MIMOSA-29, radiation hardness},\
- pdfproducer={Latex with hyperref},\
- pdfcreator={pdflatex}]{hyperref}"\
- '\n\\usepackage[utf8]{inputenc}\
- \n\\usepackage[T1]{fontenc}\
- \n\\usepackage{siunitx}\
- \n\\sisetup{load-configurations = abbreviations,separate-uncertainty}\
- \n\\usepackage{amsmath}\
- \n\\usepackage[helvet]{sfmath}"')
-
-set output "./output/$pdfName.tex"
-
-set title "$plotTitle"
-set xlabel "$xLabel"
-set ylabel "$yLabel"
-
-# set logscale x
-# set mxtics 10 # Makes logscale look good. Really? It's a matter of taste...
-
-# set autoscale # scale axes automatically
-# set style line 1 lc rgb '#0060ad' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-# set style line 2 lc rgb 'blue' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-
-$plotcmd
-EOF
-
-sleep 0.5 # wait for gnuplot saving the file so that okular can open it
-
-pdflatex -output-directory=./output/ ./output/$pdfName.tex
-cp ./output/$pdfName.pdf ./plots/$pdfName.pdf
-
-echo open file with okular...
-okular --unique "./plots/$pdfName.pdf"& # open one instance of okular
-
-
-
-# old things:
-
-# latex test.tex
-# dvips -o test.ps test.dvi
-# ps2pdf test.ps test.pdf
-
-
-
-# set terminal epslatex solid color # with extra latex file
-# set output "graph1.tex" # with extra latex file
-# okular --unique "graph1.eps"& # open one instance of okular # with extra latex file
-
-
-
-# plotcmd="plot '"./$tempfile"' using $colX:$colY index 0 w lp ls 1"
-# for i in $(seq 1 $[ $index - 1 ]); do
-# plotcmd=$plotcmd", '' using $colX:$colY index $i w lp ls "$[ $i + 1 ]
-# done
-
-
-
-
-
-
-
-
-
-
-
-
-
-# Michas Version
-# cat <<EOF | gnuplot
-# set terminal postscript enhanced solid color
-# set output "|ps2pdf - $file.pdf"
-#
-# set title "$plotTitle"
-# set xlabel "$xLabel"
-# set ylabel "$yLabel"
-#
-# set autoscale # scale axes automatically
-# # set style line 1 lc rgb '#0060ad' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-# # set style line 2 lc rgb 'blue' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-#
-# $plotcmd
-# EOF
-#
-# sleep 0.5 # wait for gnuplot saving the file so that okular can open it
-#
-#
-# echo open file with okular...
-# okular --unique "$file.pdf"& # open one instance of okular
-
-
-
-
-
+++ /dev/null
-selRunNo
-selT
-selTsens
-selChip
-selSource
-selMatrix
-selRadDose
-selSubMatrix
-selMPV
-selErrLeft
-selErrRight
-selNoise
-selNoiseLow17
-selNoiseHi17
-selSN
+++ /dev/null
-1 1 0.1 0.2 0.5 0.4
-2 3 0.1 0.0 0.6 0.4
-3 4 0.1 -1 0.3 0.4
-4 4 0.4 0.2 0.9 0.5
-5 6 0.1 0.2 0.3 0.4
+++ /dev/null
-
-reset
-#set style line 1 lc rgb '#0060ad' lt 1 lw 0.3 pt 7 ps 5 # --- blue
-set style line 2 lc rgb 'blue' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-
-plot "data" using 1:3:4:5:6 with candlesticks
\ No newline at end of file
+++ /dev/null
-# Einfacher Plot
-#==========================================
-
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
-
-
-colRunNo = 1
-colT = 2
-colTsens = 3
-colChip = 4
-colSource = 5
-colMatrix = 6
-colRadDose = 7
-colSubMatrix = 8
-colMPV = 9
-colErrLeft = 10
-colErrRight = 11
-colNoise = 12
-colNoiseLow17 = 13
-colNoiseHi17 = 14
-colSN = 15
-
-selRunNo = -1
-selT = 20
-selTsens = -1
-selChip = 1
-selSource = '-1'
-selMatrix = -1
-selRadDose = -1
-selSubMatrix = 0
-
-selMPV = -1
-selErrLeft = -1
-selErrRight = -1
-selNoise = -1
-selNoiseLow17 = -1
-selNoiseHi17 = -1
-selSN = -1
-
-set title 'S/N ratio'
-set xlabel 'x'
-set ylabel 'S/N'
-
-#set xrange [-1:30]
-#set yrange [-10:250]
-set autoscale # scale axes automatically
-
-#print colSource
-#print selSource
-
-
-plot 'Sr90Auswertung.csv' using (column(colTsens)):( ((column(colRunNo) == selRunNo) || (selRunNo == -1)) \
- && ((column(colT) == selT) || (selT == -1)) \
- && ((column(colTsens) == selTsens) || (selTsens == -1)) \
- && ((column(colChip) == selChip) || (selChip == -1)) \
- && ((strcol(colSource) eq '-1') || (selSource eq '-1')) \
- && ((column(colMatrix) == selMatrix) || (selMatrix == -1)) \
- && ((column(colRadDose) == selRadDose) || (selRadDose == -1)) \
- && ((column(colSubMatrix) == selSubMatrix) || (selSubMatrix == -1)) \
- && ((column(colMPV) == selMPV) || (selMPV == -1)) \
- && ((column(colErrLeft) == selErrLeft) || (selErrLeft == -1)) \
- && ((column(colErrRight) == selErrRight) || (selErrRight == -1)) \
- && ((column(colNoise) == selNoise) || (selNoise == -1)) \
- && ((column(colNoiseLow17) == selNoiseLow17) || (selNoiseLow17 == -1)) \
- && ((column(colNoiseHi17) == selNoiseHi17) || (selNoiseHi17 == -1)) \
- && ((column(colSN) == selSN) || (selSN == -1)) \
- ? column(colSN) : 1/0) with points
-
-
-
-
-#plot 'Sr90Auswertung.csv' every ::1 using 2:15 with points
\ No newline at end of file
+++ /dev/null
-# for (int i=0;i<6;i++){PlotGraph(20,1,"Fe55,dimmed",64,i,-1);};
-# Submatrix, {0.30, 0.45, 0.50, 0.55, 0.70}, fitMaximum
-0 56.1977 68.3771 72.3286 76.6204 93.0051 66.54
-1 44.0672 50.0286 51.9974 54.1095 62.4898 49.24
-2 45.0728 52.7417 55.3718 58.1622 67.3485 50.52
-3 44.6763 51.7905 54.3085 56.9561 65.9129 48.45
-4 38.832 44.5131 46.6017 48.816 56.2371 41.22
-5 30.1511 35.536 37.495 39.4797 45.0684 44.07
+++ /dev/null
-
-reset
-#set style line 1 lc rgb '#0060ad' lt 1 lw 1 pt 7 ps 5 # --- blue
-set style line 2 lc rgb 'blue' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-
-set xrange [-1:6]
-set xlabel "Submatrix der 64µm-Matrix des unbestrahlten Mi-29 bei 20°C"
-set ylabel "CCE/ADC"
-
-set label "Box: 0.45-0.55 Quantil des Spektrums\nWhisker-Bars: 0.3-0.7 Quantil\nBlauer Punkt: Median\nSchwarzer Punkt: Maximum des Gaussfits ('traditioneller' CCE-Wert)" at 3, 90
-
-plot "data" using 1:3:2:6:5 with candlesticks whiskerbars,\
- '' using 1:7:7:7:7 with candlesticks lt -1 notitle,\
- '' using 1:4:4:4:4 with candlesticks lt 2 notitle
+++ /dev/null
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-29024 20 21.5 1 Sr90 64 0 0 45.85 2.25 1.59 4.08 3.89 4.27 11.24
-29024 20 21.5 1 Sr90 64 0 1 35.28 0.57 0.02 2.84 2.74 2.98 12.41
-29024 20 21.5 1 Sr90 64 0 2 36.95 0.2 1.19 2.85 2.74 2.97 12.98
-29024 20 21.5 1 Sr90 64 0 3 36.65 0.14 0.67 2.83 2.71 2.95 12.96
-29024 20 21.5 1 Sr90 64 0 4 31.18 0.07 0.19 2.82 2.72 2.96 11.05
-29024 20 21.5 1 Sr90 64 0 5 25.68 0.09 0.43 2.44 2.35 2.57 10.54
-
-29025 20 21.5 1 Sr90 80 0 0 38.39 3.02 2.74 4.26 4.06 4.46 9.02
-29025 20 21.5 1 Sr90 80 0 1 29.5 0.78 0.17 2.81 2.69 2.97 10.5
-29025 20 21.5 1 Sr90 80 0 2 43.58 1.31 2.84 4.54 4.3 4.75 9.6
-29025 20 21.5 1 Sr90 80 0 3 32.7 0.01 0.65 2.54 2.39 2.76 12.87
-29025 20 21.5 1 Sr90 80 0 4 30.58 0.51 0.89 2.57 2.38 2.79 11.91
-29025 20 21.5 1 Sr90 80 0 5 23.24 0.14 0.11 2.26 2.14 2.4 10.29
-
-29026 -20 -14 1 Sr90 64 0 0 51.81 5.6 6.65 2.37 2.29 2.51 21.82
-29026 -20 -14 1 Sr90 64 0 1 38.11 1.14 0.13 1.99 1.93 2.1 19.19
-29026 -20 -14 1 Sr90 64 0 2 39.48 0.06 1.74 1.94 1.88 2.04 20.34
-29026 -20 -14 1 Sr90 64 0 3 38.94 0.32 1.07 1.92 1.86 2.04 20.32
-29026 -20 -14 1 Sr90 64 0 4 33.19 0.21 0.55 1.89 1.83 2 17.6
-29026 -20 -14 1 Sr90 64 0 5 27.01 0 1.1 1.78 1.73 1.89 15.15
-
-29027 -20 -14.5 1 Sr90 80 0 0 44.64 3.6 5.02 2.03 1.94 2.19 21.97
-29027 -20 -14.5 1 Sr90 80 0 1 31.96 1.92 0.45 1.6 1.54 1.74 19.96
-29027 -20 -14.5 1 Sr90 80 0 2 48.54 1.15 2.38 2.04 1.94 2.2 23.81
-29027 -20 -14.5 1 Sr90 80 0 3 34.28 0.11 1.15 1.58 1.51 1.71 21.74
-29027 -20 -14.5 1 Sr90 80 0 4 31.74 0.23 1.32 1.58 1.51 1.72 20.13
-29027 -20 -14.5 1 Sr90 80 0 5 24.2 0.17 0.34 1.48 1.42 1.62 16.35
-
-29028 -70 -53.5 1 Sr90 64 0 0 52.7 5.41 3.53 2.22 2.13 2.37 23.73
-29028 -70 -53.5 1 Sr90 64 0 1 37.34 0.54 0.82 1.67 1.62 1.81 22.29
-29028 -70 -53.5 1 Sr90 64 0 2 38.97 0.05 0.99 1.65 1.59 1.77 23.55
-29028 -70 -53.5 1 Sr90 64 0 3 38.63 0.41 0.26 1.64 1.58 1.76 23.61
-29028 -70 -53.5 1 Sr90 64 0 4 32.87 0.14 0.53 1.66 1.6 1.78 19.79
-29028 -70 -53.5 1 Sr90 64 0 5 26.79 0.05 1.19 1.56 1.5 1.69 17.18
-
-29029 -70 -54.5 1 Sr90 80 0 0 45.68 4.76 8.52 2.1 2.01 2.25 21.74
-29029 -70 -54.5 1 Sr90 80 0 1 31.69 0.69 0.61 1.5 1.44 1.64 21.09
-29029 -70 -54.5 1 Sr90 80 0 2 47.78 1.1 2.64 2.06 1.96 2.23 23.16
-29029 -70 -54.5 1 Sr90 80 0 3 33.99 0.3 0.79 1.49 1.44 1.63 22.74
-29029 -70 -54.5 1 Sr90 80 0 4 31.16 0.25 1.68 1.5 1.44 1.66 20.73
-29029 -70 -54.5 1 Sr90 80 0 5 24.31 0.14 0.28 1.43 1.37 1.57 17.05
-
-29046 20 21.6 2 Sr90 64 150000 0 42.2 1.25 3.42 4.42 4.22 4.63 9.55
-29046 20 21.6 2 Sr90 64 150000 1 33.48 0.19 0.53 3.13 3.03 3.25 10.7
-29046 20 21.6 2 Sr90 64 150000 2 36.19 0.22 0.66 3.1 3 3.22 11.67
-29046 20 21.6 2 Sr90 64 150000 3 35.4 0.12 1.12 3.06 2.97 3.17 11.57
-29046 20 21.6 2 Sr90 64 150000 4 31 0.21 0.01 2.99 2.89 3.11 10.36
-29046 20 21.6 2 Sr90 64 150000 5 25.48 0.15 0.65 2.62 2.53 2.73 9.74
-
-29047 20 21 2 Sr90 80 150000 0 37.48 0.26 2.53 4.22 4.05 4.39 8.89
-29047 20 21 2 Sr90 80 150000 1 29.19 0.8 0.01 2.76 2.66 2.89 10.57
-29047 20 21 2 Sr90 80 150000 2 43.54 0.2 2.79 4.38 4.16 4.56 9.95
-29047 20 21 2 Sr90 80 150000 3 32.38 0.13 0.29 2.67 2.57 2.79 12.14
-29047 20 21 2 Sr90 80 150000 4 30.51 0.25 1.11 2.64 2.48 2.81 11.54
-29047 20 21 2 Sr90 80 150000 5 23.19 0.24 0.3 2.24 2.14 2.36 10.37
-
-29048 -20 -15.5 2 Sr90 64 150000 0 49.19 2.82 3.53 2.44 2.36 2.59 20.13
-29048 -20 -15.5 2 Sr90 64 150000 1 35.88 0.07 0.93 2.09 2.03 2.2 17.16
-29048 -20 -15.5 2 Sr90 64 150000 2 38.53 0.54 1.84 2.04 1.98 2.15 18.86
-29048 -20 -15.5 2 Sr90 64 150000 3 38.07 0.06 1.61 2.01 1.95 2.12 18.93
-29048 -20 -15.5 2 Sr90 64 150000 4 32.51 0.11 0.23 1.97 1.91 2.08 16.49
-29048 -20 -15.5 2 Sr90 64 150000 5 26.71 0.06 1.06 1.87 1.81 1.98 14.31
-
-29049 -20 -16 2 Sr90 80 150000 0 42.62 2.54 5.63 1.97 1.89 2.13 21.66
-29049 -20 -16 2 Sr90 80 150000 1 30.89 1.04 0.27 1.62 1.56 1.75 19.07
-29049 -20 -16 2 Sr90 80 150000 2 47.72 1.47 2.44 1.97 1.88 2.15 24.22
-29049 -20 -16 2 Sr90 80 150000 3 33.8 0.22 1.42 1.61 1.55 1.73 21
-29049 -20 -16 2 Sr90 80 150000 4 31.33 0.25 1.19 1.61 1.55 1.75 19.45
-29049 -20 -16 2 Sr90 80 150000 5 24.03 0.2 0.43 1.51 1.45 1.64 15.96
-
-29050 -70 -54.5 2 Sr90 64 150000 0 51.01 4.9 5.5 2.16 2.08 2.32 23.6
-29050 -70 -54.5 2 Sr90 64 150000 1 35.65 1.03 0.47 1.76 1.7 1.87 20.25
-29050 -70 -54.5 2 Sr90 64 150000 2 37.42 0.06 1.25 1.74 1.68 1.84 21.54
-29050 -70 -54.5 2 Sr90 64 150000 3 37.41 0.16 1.81 1.72 1.66 1.83 21.79
-29050 -70 -54.5 2 Sr90 64 150000 4 32.61 0.18 0.39 1.68 1.62 1.8 19.4
-29050 -70 -54.5 2 Sr90 64 150000 5 26.86 0.05 0.94 1.6 1.54 1.7 16.84
-
-29051 -70 -54.5 2 Sr90 80 150000 0 43.9 5.07 4.02 1.96 1.88 2.12 22.34
-29051 -70 -54.5 2 Sr90 80 150000 1 29 0.62 0.18 1.58 1.51 1.7 18.37
-29051 -70 -54.5 2 Sr90 80 150000 2 47.21 1.29 2.54 1.94 1.86 2.11 24.29
-29051 -70 -54.5 2 Sr90 80 150000 3 32.68 0.06 0.93 1.58 1.51 1.71 20.66
-29051 -70 -54.5 2 Sr90 80 150000 4 30.5 0.23 0.92 1.61 1.53 1.74 18.97
-29051 -70 -54.5 2 Sr90 80 150000 5 24.1 0.22 0.38 1.42 1.36 1.54 16.95
-
-29064 20 22 4 Sr90 64 1000000000000 0 35.1 0.92 2.24 4.21 3.99 4.39 8.33
-29064 20 22 4 Sr90 64 1000000000000 1 31.95 0.2 0.32 2.66 2.54 2.8 12.02
-29064 20 22 4 Sr90 64 1000000000000 2 33 0.26 0.77 2.75 2.64 2.87 12.02
-29064 20 22 4 Sr90 64 1000000000000 3 32.46 0.13 0.51 2.78 2.63 2.91 11.7
-29064 20 22 4 Sr90 64 1000000000000 4 25.38 0.23 0.85 2.95 2.84 3.07 8.6
-29064 20 22 4 Sr90 64 1000000000000 5 23.11 0.5 0.71 2.36 2.27 2.46 9.81
-
-29065 20 21.5 4 Sr90 80 1000000000000 0 32.1 2.47 0.4 4.06 3.9 4.23 7.91
-29065 20 21.5 4 Sr90 80 1000000000000 1 27.35 0.22 0.1 2.66 2.54 2.79 10.29
-29065 20 21.5 4 Sr90 80 1000000000000 2 32.97 5.27 3.6 4.73 4.55 4.89 6.96
-29065 20 21.5 4 Sr90 80 1000000000000 3 27.39 0.1 0.63 2.75 2.65 2.88 9.96
-29065 20 21.5 4 Sr90 80 1000000000000 4 23.16 0.62 0.76 2.83 2.73 2.97 8.18
-29065 20 21.5 4 Sr90 80 1000000000000 5 20.34 0.3 0 2.33 2.25 2.43 8.74
-
-29066 -20 -13.5 4 Sr90 64 1000000000000 0 41.62 2.88 1.34 1.96 1.87 2.13 21.28
-29066 -20 -13.5 4 Sr90 64 1000000000000 1 35.52 0.57 0.28 1.55 1.49 1.7 22.86
-29066 -20 -13.5 4 Sr90 64 1000000000000 2 35.27 0.27 1.71 1.54 1.48 1.7 22.89
-29066 -20 -13.5 4 Sr90 64 1000000000000 3 34.85 0.17 0.75 1.54 1.48 1.7 22.6
-29066 -20 -13.5 4 Sr90 64 1000000000000 4 26.11 0.08 0.02 1.53 1.47 1.69 17.03
-29066 -20 -13.5 4 Sr90 64 1000000000000 5 24.71 0.36 0.83 1.45 1.39 1.59 17.03
-
-29067 -20 -14 4 Sr90 80 1000000000000 0 35.9 3.11 0.34 1.94 1.85 2.09 18.55
-29067 -20 -14 4 Sr90 80 1000000000000 1 29.46 0.86 0.19 1.55 1.48 1.68 19.06
-29067 -20 -14 4 Sr90 80 1000000000000 2 32.05 4.79 0.84 1.91 1.82 2.06 16.78
-29067 -20 -14 4 Sr90 80 1000000000000 3 27.22 0.16 1 1.53 1.47 1.66 17.74
-29067 -20 -14 4 Sr90 80 1000000000000 4 21.14 2.99 0.79 1.52 1.47 1.67 13.87
-29067 -20 -14 4 Sr90 80 1000000000000 5 20.69 0.22 0 1.45 1.4 1.56 14.3
-
-29068 -70 -53 4 Sr90 64 1000000000000 0 39.93 3.47 0.23 2.02 1.92 2.19 19.76
-29068 -70 -53 4 Sr90 64 1000000000000 1 35.13 0.43 0.36 1.5 1.43 1.65 23.48
-29068 -70 -53 4 Sr90 64 1000000000000 2 34.36 0.46 0.18 1.48 1.42 1.65 23.15
-29068 -70 -53 4 Sr90 64 1000000000000 3 33.18 0.07 0.69 1.48 1.42 1.65 22.42
-29068 -70 -53 4 Sr90 64 1000000000000 4 23.01 0.39 0.49 1.48 1.42 1.65 15.51
-29068 -70 -53 4 Sr90 64 1000000000000 5 23.79 0.24 0.46 1.4 1.34 1.53 17.01
-
-29069 -70 -54 4 Sr90 80 1000000000000 0 32.71 3.82 1.25 2.02 1.92 2.17 16.21
-29069 -70 -54 4 Sr90 80 1000000000000 1 28.96 0.36 0.38 1.48 1.42 1.63 19.6
-29069 -70 -54 4 Sr90 80 1000000000000 2 23.83 2.63 2.31 1.91 1.82 2.08 12.46
-29069 -70 -54 4 Sr90 80 1000000000000 3 23.99 0.06 0.29 1.46 1.4 1.62 16.45
-29069 -70 -54 4 Sr90 80 1000000000000 4 17.17 0 0.58 1.46 1.4 1.63 11.73
-29069 -70 -54 4 Sr90 80 1000000000000 5 19.09 0.87 0 1.39 1.34 1.51 13.7
+++ /dev/null
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-
-
-# block index 0
-29024 20 21.5 1 Sr90 64 0 0 45.85 2.25 1.59 4.08 3.89 4.27 11.24
-29024 20 21.5 1 Sr90 64 0 1 35.28 0.57 0.02 2.84 2.74 2.98 12.41
-29024 20 21.5 1 Sr90 64 0 2 36.95 0.2 1.19 2.85 2.74 2.97 12.98
-29024 20 21.5 1 Sr90 64 0 3 36.65 0.14 0.67 2.83 2.71 2.95 12.96
-29024 20 21.5 1 Sr90 64 0 4 31.18 0.07 0.19 2.82 2.72 2.96 11.05
-29024 20 21.5 1 Sr90 64 0 5 25.68 0.09 0.43 2.44 2.35 2.57 10.54
-
-
-# block index 1
-29026 -20 -14 1 Sr90 64 0 0 51.81 5.6 6.65 2.37 2.29 2.51 21.82
-29026 -20 -14 1 Sr90 64 0 1 38.11 1.14 0.13 1.99 1.93 2.1 19.19
-29026 -20 -14 1 Sr90 64 0 2 39.48 0.06 1.74 1.94 1.88 2.04 20.34
-29026 -20 -14 1 Sr90 64 0 3 38.94 0.32 1.07 1.92 1.86 2.04 20.32
-29026 -20 -14 1 Sr90 64 0 4 33.19 0.21 0.55 1.89 1.83 2 17.6
-29026 -20 -14 1 Sr90 64 0 5 27.01 0 1.1 1.78 1.73 1.89 15.15
-
-
-# block index 2
-29028 -70 -53.5 1 Sr90 64 0 0 52.7 5.41 3.53 2.22 2.13 2.37 23.73
-29028 -70 -53.5 1 Sr90 64 0 1 37.34 0.54 0.82 1.67 1.62 1.81 22.29
-29028 -70 -53.5 1 Sr90 64 0 2 38.97 0.05 0.99 1.65 1.59 1.77 23.55
-29028 -70 -53.5 1 Sr90 64 0 3 38.63 0.41 0.26 1.64 1.58 1.76 23.61
-29028 -70 -53.5 1 Sr90 64 0 4 32.87 0.14 0.53 1.66 1.6 1.78 19.79
-29028 -70 -53.5 1 Sr90 64 0 5 26.79 0.05 1.19 1.56 1.5 1.69 17.18
+++ /dev/null
-#!/bin/bash
-
-file=$1
-
-tempfile="lastplot.tmp"
-
-
-colRunNo=1
-colT=2
-colTsens=3
-colChip=4
-colSource=5
-colMatrix=6
-colRadDose=7
-colSubMatrix=8
-colMPV=9
-colErrLeft=10
-colErrRight=11
-colNoise=12
-colNoiseLow17=13
-colNoiseHi17=14
-colSN=15
-
-colX=$colSubMatrix
-colY=$colSN
-groupBy=$colT
-plotTitle='S/N-ratio'
-xLabel='Matrix'
-yLabel='S/N'
-# titles=('Submatrix 0' 'Submatrix 1' 'Submatrix 2' 'Submatrix 3' 'Submatrix 4' 'Submatrix 5')
-titles=('20 C' '-20 C' '-70 C')
-
-selRunNo=''
-# selT='20'
-selTsens=''
-selChip='1'
-selSource='S\S*'
-selMatrix='64'
-selRadDose=''
-selSubMatrix=''
-selMPV=''
-selErrLeft=''
-selErrRight=''
-selNoise=''
-selNoiseLow17=''
-selNoiseHi17=''
-selSN=''
-
-##################
-# select only matching data points
-regexp='^' # new line character
-
-for i in selRunNo selT selTsens selChip selSource selMatrix selRadDose selSubMatrix selMPV selErrLeft selErrRight selNoise selNoiseLow17 selNoiseHi17 selSN
-do
- eval sel=\$$i # set variable sel to the string which is contained in the selected i
- if [ -z $sel ]; then
- pattern='\S*' # if length of string in selector is zero, accept every string in this column ('\S*' is a regular expression for any number of non-whitespace characters)
- else
- pattern=$sel # if length of string in selector is non-zero, accept only matching strings
- fi
-
- regexp=$regexp$pattern"\s+" # pattern plus at least one whitespace character
-done
-
-cat $file | egrep "^#" > $tempfile # overwrite tempfile
-cat $file | egrep "$regexp" >> $tempfile # append to tempfile
-
-
-
-##################
-# split matching data sets in groups (criterion: all elements of one group belong to the same submatrix)
-memory=""
-index="0"
-
-if [ -n "$groupBy" ]; then # is grouping non zero?
- tmpA=$(mktemp) # create temporary file :-)
- egrep "^#" $tempfile > $tmpA # write comment lines to temp file
- while IFS="" read line
- do
- newEntry=$(echo "$line" | cut -f $groupBy) # get entry of column $groupBy of line $line
- if [ "$memory" != "$newEntry" ]; then # number in column $groupBy does not match the previous one
- echo >> $tmpA # one blank line for splitting the data sets
- echo >> $tmpA # a second blank line for seperate linestyles
- echo "# block index $index" >> $tmpA
- index=$[ $index + 1 ]
- fi
- echo "$line" >> $tmpA
- memory=$newEntry
- done < <(egrep -v "^#" $tempfile | sort -k$groupBy,$groupBy -k$colT,$colT) # input for the while loop: sorted data file without comment lines. -k3,5 means "take cloumn 3 to 5 as one 'column to sort by'". normally one wants to take only one column as entity to sort by, so usally one chooses the option -k5,5 to sort by column 5. sort -k3,3 -k5,5 means sorting first by column 3 and then subsort by column 5.
- cat $tmpA > $tempfile # write final data file
- rm $tmpA
-fi
-
-
-##################
-# plotting
-echo
-echo
-echo making a plot with the following data:
-echo
-cat $tempfile # print tempfile for information
-echo
-echo
-echo starting gnuplot ...
-
-plotcmd="plot '"./$tempfile"' using $colX:$colY index 0 w lp title '${titles[0]}'"
-for i in $(seq 1 $[ $index - 1 ]); do
-plotcmd=$plotcmd", '' using $colX:$colY index $i w lp title '${titles[$i]}'"
-done
-
-
-cat <<EOF | gnuplot
-set terminal postscript enhanced solid color
-set output "|ps2pdf - $file.pdf"
-set title "$plotTitle"
-set xlabel "$xLabel"
-set ylabel "$yLabel"
-
-set autoscale # scale axes automatically
-# set style line 1 lc rgb '#0060ad' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-# set style line 2 lc rgb 'blue' lt 1 lw 2 pt 7 ps 1.5 # --- blue
-
-$plotcmd
-
-EOF
-
-sleep 0.1
-
-echo open file with okular...
-okular "$file.pdf"&
-
-
-
-
-# plotcmd="plot '"./$tempfile"' using $colX:$colY index 0 w lp ls 1"
-# for i in $(seq 1 $[ $index - 1 ]); do
-# plotcmd=$plotcmd", '' using $colX:$colY index $i w lp ls "$[ $i + 1 ]
-# done
-
-
-
-
+++ /dev/null
-selRunNo
-selT
-selTsens
-selChip
-selSource
-selMatrix
-selRadDose
-selSubMatrix
-selMPV
-selErrLeft
-selErrRight
-selNoise
-selNoiseLow17
-selNoiseHi17
-selSN
+++ /dev/null
-# Einfacher Plot
-#==========================================
-
-# runNo T Tsens Chip Source Matrix RadDose Submtrx MPV ErrLeft ErrRight Noise NoiseLow17 NoiseHi17 S/N Edited by hand Comment
-# 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
-
-
-colRunNo = 1
-colT = 2
-colTsens = 3
-colChip = 4
-colSource = 5
-colMatrix = 6
-colRadDose = 7
-colSubMatrix = 8
-colMPV = 9
-colErrLeft = 10
-colErrRight = 11
-colNoise = 12
-colNoiseLow17 = 13
-colNoiseHi17 = 14
-colSN = 15
-
-selRunNo = -1
-selT = 20
-selTsens = -1
-selChip = 1
-selSource = '-1'
-selMatrix = -1
-selRadDose = -1
-selSubMatrix = 0
-
-selMPV = -1
-selErrLeft = -1
-selErrRight = -1
-selNoise = -1
-selNoiseLow17 = -1
-selNoiseHi17 = -1
-selSN = -1
-
-set title 'S/N ratio'
-set xlabel 'x'
-set ylabel 'S/N'
-
-#set xrange [-1:30]
-#set yrange [-10:250]
-set autoscale # scale axes automatically
-
-#print colSource
-#print selSource
-
-
-plot 'Sr90Auswertung.csv' using (column(colTsens)):( ((column(colRunNo) == selRunNo) || (selRunNo == -1)) \
- && ((column(colT) == selT) || (selT == -1)) \
- && ((column(colTsens) == selTsens) || (selTsens == -1)) \
- && ((column(colChip) == selChip) || (selChip == -1)) \
- && ((strcol(colSource) eq '-1') || (selSource eq '-1')) \
- && ((column(colMatrix) == selMatrix) || (selMatrix == -1)) \
- && ((column(colRadDose) == selRadDose) || (selRadDose == -1)) \
- && ((column(colSubMatrix) == selSubMatrix) || (selSubMatrix == -1)) \
- && ((column(colMPV) == selMPV) || (selMPV == -1)) \
- && ((column(colErrLeft) == selErrLeft) || (selErrLeft == -1)) \
- && ((column(colErrRight) == selErrRight) || (selErrRight == -1)) \
- && ((column(colNoise) == selNoise) || (selNoise == -1)) \
- && ((column(colNoiseLow17) == selNoiseLow17) || (selNoiseLow17 == -1)) \
- && ((column(colNoiseHi17) == selNoiseHi17) || (selNoiseHi17 == -1)) \
- && ((column(colSN) == selSN) || (selSN == -1)) \
- ? column(colSN) : 1/0) with points
-
-
-
-
-#plot 'Sr90Auswertung.csv' every ::1 using 2:15 with points
\ No newline at end of file
jspc29 / radhard.git / commitdiff
