#define COLDONE 10
#define NUMPIXELS 25
-#define RIGHT_BOUNDARY 1000
-#define VETO_BINSIZE 1000
+// #define RIGHT_BOUNDARY 350
+#define RIGHT_BOUNDARY 384 // 2**8*3 = 384
+// #define SEED_AND_SUM_BINS 200
+#define SEED_AND_SUM_BINS RIGHT_BOUNDARY/2
+#define VETO_BINS RIGHT_BOUNDARY/2
-#define VETO_THRESHOLD 5
+// #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};
#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_15.11.2013_16:3/"
+#define DATAPATH "/local/sstrohauer/data/data_19.11.2013_14:11/"
#define USE_SEPARATE_NOISE_RUN true
#define DRAW_SINGLE_HISTOGRAMS true
-#define DRAW_SUMMED_HISTOGRAMS true
-#define DRAW_VETO_HISTOGRAMS true
+#define DRAW_SUMMED_HISTOGRAMS false
+#define DRAW_VETO_HISTOGRAMS false
#define DRAW_SEEDPIXEL_DISTRIBUTION false
#define DRAW_PIXEL_DISTRIBUTION_HISTOGRAMS false
// examples for fit functions: gaus, landau
#define FIT_FUNC "gaus"
-#define DRAW_FITS true
+#define DRAW_FITS false
#define DRAW_NOISE_HISTOGRAM false
#define PRINT_EXPECTATION_VALUE_HITS_PER_FRAME false
-#define PRINT_HEADER false
+#define PRINT_HEADER true
+
+#define xTitle "Gesammelte Ladung / ADC"
+#define yTitle "Haeufigkeit"
+#define graphTitle ""
void PlotGraph(Float_t selTemp, Int_t selChip, std::string selSource, Int_t selMatrix6480, Int_t selSubmatrix, Float_t selRadDose)
{
hitsNtuple->SetBranchAddress("counts", &frameCounts, &frameCountsBranch);
// create histogram
- histNtuple.push_back(new TH1F(Form("hist%i",nHistNtuple), "Histogram title", 200, 0, RIGHT_BOUNDARY));
- histNtupleSum.push_back(new TH1F(Form("histSum%i",nHistNtuple), "Histogram title", 200, 0, RIGHT_BOUNDARY));
- histNtupleVeto.push_back(new TH1F(Form("histVeto%i",nHistNtuple), "Histogram title", VETO_BINSIZE, 0, RIGHT_BOUNDARY));
+ histNtuple.push_back(new TH1F(Form("hist%i",nHistNtuple), "Histogram title", SEED_AND_SUM_BINS, 0, RIGHT_BOUNDARY));
+ histNtupleSum.push_back(new TH1F(Form("histSum%i",nHistNtuple), "Histogram title", SEED_AND_SUM_BINS, 0, RIGHT_BOUNDARY));
+ histNtupleVeto.push_back(new TH1F(Form("histVeto%i",nHistNtuple), "Histogram title", VETO_BINS, 0, RIGHT_BOUNDARY));
histPixelNumIsSeed.push_back(new TH1F(Form("histPixelNumIsSeed_%i",nHistNtuple), "Histogram title", totalPixelNum, 0, totalPixelNum));
histPixelNumIsSeedSelection.push_back(new TH1F(Form("histPixelNumIsSeedSelection_%i",nHistNtuple), "Histogram title", totalPixelNum, 0, totalPixelNum));
histPixelDistribution1.push_back(new TH1F(Form("histPixelDistribution1_%i",nHistNtuple), "Histogram title", totalPixelNum, 0, totalPixelNum));
// histogram with the single pixel and with "veto trigger"
// for the "histNtupleVeto" take only hits where only the seed pixel contains charge
- Double_t notSeedSum = 0;
- for (Int_t i=0; i<12; i++)
- notSeedSum += pixel[i];
- for (Int_t i=13; i<NUMPIXELS; i++)
- notSeedSum += pixel[i];
- if (TMath::Abs(notSeedSum) > VETO_THRESHOLD)
+ Bool_t skipHit = false;
+ Float_t* vetoThreshold;
+ if (matrix6480 == 64)
+ vetoThreshold = vetoThreshold64;
+ else if (matrix6480 == 80)
+ vetoThreshold = vetoThreshold80;
+
+ for (Int_t i=0; i<12; i++) {
+ if (pixel[i] > vetoThreshold[submatrix])
+ skipHit = true;
+ }
+ for (Int_t i=13; i<NUMPIXELS; i++) {
+ if (pixel[i] > vetoThreshold[submatrix])
+ skipHit = true;
+ }
+ if (skipHit)
continue;
-// meanNotSeedSum+=notSeedSum;
-// nhits++;
+
+
+ // alternative shitty vetotrigger 1
+// Int_t numVetoPixel = 0;
+// Float_t* vetoThreshold;
+// if (matrix6480 == 64)
+// vetoThreshold = vetoThreshold64;
+// else if (matrix6480 == 80)
+// vetoThreshold = vetoThreshold80;
+//
+// for (Int_t i=0; i<12; i++) {
+// if (pixel[i] > vetoThreshold[submatrix])
+// numVetoPixel++;
+// }
+// for (Int_t i=13; i<NUMPIXELS; i++) {
+// if (pixel[i] > vetoThreshold[submatrix])
+// numVetoPixel++;
+// }
+// if (numVetoPixel > 0)
+// continue;
+
+ // alternatative shitty veto trigger 2 (watch out! Summation over all pixels except the seed is
+ // nonsense: remember the pedestal of ~ -1 for every pixel)
+ // Double_t notSeedSum = 0;
+// for (Int_t i=0; i<12; i++)
+// notSeedSum += pixel[i];
+// for (Int_t i=13; i<NUMPIXELS; i++)
+// notSeedSum += pixel[i];
+// if (TMath::Abs(notSeedSum) > VETO_THRESHOLD)
+// continue;
+
histNtupleVeto[nHistNtuple]->Fill(pixel[12]); // histogram with the single pixel
}
histNtuple[nHistNtuple]->SetLineColor(nHistNtuple+1);
+ histNtuple[nHistNtuple]->SetLineColor(nHistNtuple+1>4 ? nHistNtuple+2:nHistNtuple+1); // omit yellow line!
histNtuple[nHistNtuple]->SetStats(kFALSE);
- histNtupleSum[nHistNtuple]->SetLineColor(nHistNtuple+1);
+ histNtupleSum[nHistNtuple]->SetLineColor(nHistNtuple+1>4 ? nHistNtuple+2:nHistNtuple+1); // omit yellow line!
histNtupleSum[nHistNtuple]->SetStats(kFALSE);
histNtupleSum[nHistNtuple]->SetLineStyle(2); // make summed graphs dashed
- histNtupleVeto[nHistNtuple]->SetLineColor(nHistNtuple+1);
+ histNtupleVeto[nHistNtuple]->SetLineColor(nHistNtuple+1>4 ? nHistNtuple+2:nHistNtuple+1); // omit yellow line!
histNtupleVeto[nHistNtuple]->SetStats(kFALSE);
- histNtupleVeto[nHistNtuple]->SetLineStyle(3); // make summed graphs dotted
- histPixelNumIsSeedSelection[nHistNtuple]->SetLineColor(nHistNtuple+1);
+ 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);
// cout << "Hits per Frame: " << nentries/frame << endl;
-
//-------------------------------------------------
// calculate median of noise (with error estimation)
Int_t runNo_N; // run number
TCanvas* c1 = new TCanvas();
- c1->SetTitle("Testtitel");
TF1* fitFunc = new TF1("fitFunc",FIT_FUNC,0,RIGHT_BOUNDARY);
// plot every single pixel histogram in one canvas
if (DRAW_SINGLE_HISTOGRAMS){
(*it)->SetMaximum(1.05*maxY);
(*it)->Draw((cnt==0)?"":"same");
+ (*it)->SetTitle(graphTitle);
+ (*it)->GetXaxis()->SetTitle(xTitle);
+ (*it)->GetYaxis()->SetTitle(yTitle);
+
if (DRAW_FITS) {
Float_t posMax = 0;
(*it)->GetXaxis()->SetRange((*it)->GetXaxis()->FindBin(20),(*it)->GetXaxis()->FindBin(RIGHT_BOUNDARY)); // look only for maxima with x greater than 20
fitFunc->DrawCopy("same");
fitFunc->SetLineStyle(1); // normal 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";
+
+
+
+
+
+
+
+
+
+
+
if (DRAW_SUMMED_HISTOGRAMS){
for (std::vector<TH1F*>::iterator it = histNtupleSum.begin(); it != histNtupleSum.end(); ++it){
Int_t cnt = std::distance(histNtupleSum.begin(), it) ;
- (*it)->SetMaximum(1.05*maxY);
if (!DRAW_SINGLE_HISTOGRAMS)
(*it)->Draw((cnt==0)?"":"same");
else
}
}
}
+ // add legend
+ TString legendEntry = TString("Same Parameters: ") + (selTemp!=-1 ? Form(" T_{set}=%.1f",selTemp):"") + (selChip!=-1 ? Form(" Chip=%i",selChip):"") + (selSource!="-1" ? Form(" Source=%s",selSource.c_str()):"") + (selMatrix6480!=-1 ? Form(" Matrix=%i",selMatrix6480):"") + (selSubmatrix!=-1 ? Form(" Submatrix=%i",selSubmatrix):"") + (selRadDose!=-1 ? Form(" RadDose=%.2E",selRadDose):"");
+ leg->AddEntry((TObject*) 0, legendEntry, "");
+ leg->SetTextSize(0.02);
+ leg->Draw();
+
// plot single pixel histograms with veto trigger in the same histogram
if (DRAW_VETO_HISTOGRAMS){
+ TCanvas* c4 = new TCanvas();
for (std::vector<TH1F*>::iterator it = histNtupleVeto.begin(); it != histNtupleVeto.end(); ++it){
Int_t cnt = std::distance(histNtupleVeto.begin(), it) ;
- (*it)->SetMaximum(1.05*maxY);
- if (!DRAW_SINGLE_HISTOGRAMS && !DRAW_SUMMED_HISTOGRAMS)
- (*it)->Draw((cnt==0)?"":"same");
- else
- (*it)->Draw("same");
+ (*it)->Draw((cnt==0)?"":"same");
if (DRAW_FITS) {
// (*it)->GetXaxis()->SetRange((*it)->GetXaxis()->FindBin(fitMaxPosX[cnt] + 25),(*it)->GetXaxis()->FindBin(RIGHT_BOUNDARY)); // look only for maxima with x greater than 25 more right than the peak of the seed pixel
}
}
}
- // add legend
- TString legendEntry = TString("Same Parameters: ") + (selTemp!=-1 ? Form(" T_{set}=%.1f",selTemp):"") + (selChip!=-1 ? Form(" Chip=%i",selChip):"") + (selSource!="-1" ? Form(" Source=%s",selSource.c_str()):"") + (selMatrix6480!=-1 ? Form(" Matrix=%i",selMatrix6480):"") + (selSubmatrix!=-1 ? Form(" Submatrix=%i",selSubmatrix):"") + (selRadDose!=-1 ? Form(" RadDose=%.2E",selRadDose):"");
- leg->AddEntry((TObject*) 0, legendEntry, "");
- leg->SetTextSize(0.02);
- leg->Draw();
// plot information about the seed pixel distribution in another histogram
}
+
+
+
+
+
+
+
+
+ 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;
+
+ }
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
// print all information
TString runDetailsHeader = TString("runNo\tT\tTSens\tChip\tSource\tMatrix\tRadDose\t\tSubmtrx")
+ (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tColPeak":"")
+ (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? "\tSumColPeak":"")
+ (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? "\tSumColErr":"")
+ (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? "\tCalPeak":"")
- + (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? "\tCalPkErr":"")
+ + (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? "\tCalPkEr":"")
+ "\tNoise"
+ "\tNseLo17"
+ "\tNseHi17"
<< endl;
}
-
gettimeofday(&end, 0);
diff = (int) (end.tv_sec-start.tv_sec);
if (PRINT_HEADER)
-// Int_t runNo = 29012;
-// Int_t submatrix = 0;
-//
-//
-//
-//
-// // Variables used to store the data
-// Int_t totalSize = 0; // Sum of data size (in bytes) of all events
-// Int_t eventSize = 0; // Size of the current event
-// TBranch *eventSizeBranch = 0; // Pointer to the event.fEventsize branch
-//
-// // open the file
-// TString path = TString(OUTPUTDATAPATH) + Form("%i/%i_0_%i.root", runNo, runNo, submatrix);
-// TFile *f = TFile::Open(path);
-//
-// if (f == 0) {
-// // if we cannot open the file, print an error message and return immediatly
-// printf("Error: cannot open %s\n", path);
-// return;
-// }
-//
-// // get a pointer to the tree
-// TNtuple* hitNtuple = (TNtuple*) f->Get("hit");
-//
-// hitNtuple->Draw("p13");
-// // hitNtuple->Draw("p10","","same");
-//
-// delete f;
-//
-// Int_t runNo = 29016;
-// Int_t submatrix = 0;
-// TString path = TString(OUTPUTDATAPATH) + Form("%i/%i_0_%i.root", runNo, runNo, submatrix);
-// TFile *f = TFile::Open(path);
-// if (f == 0) {
-// // if we cannot open the file, print an error message and return immediatly
-// printf("Error: cannot open %s\n", path);
-// return;
+// // plot single pixel histograms with veto trigger in the same histogram
+// if (DRAW_VETO_HISTOGRAMS){
+// for (std::vector<TH1F*>::iterator it = histNtupleVeto.begin(); it != histNtupleVeto.end(); ++it){
+// Int_t cnt = std::distance(histNtupleVeto.begin(), it) ;
+// (*it)->SetMaximum(1.05*maxY);
+// if (!DRAW_SINGLE_HISTOGRAMS && !DRAW_SUMMED_HISTOGRAMS)
+// (*it)->Draw((cnt==0)?"":"same");
+// else
+// (*it)->Draw("same");
+//
+// if (DRAW_FITS) {
+// // (*it)->GetXaxis()->SetRange((*it)->GetXaxis()->FindBin(fitMaxPosX[cnt] + 25),(*it)->GetXaxis()->FindBin(RIGHT_BOUNDARY)); // look only for maxima with x greater than 25 more right than the peak of the seed pixel
+// (*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
+//
+// (*it)->GetXaxis()->UnZoom();
+// fitFunc->DrawCopy("same");
+// // cout << "Maximum bin: " << (*it)->GetBinCenter((*it)->GetMaximumBin()) << endl;
+// // cout << "Maximum of gaussian fit: " << fitFunc->GetParameter(1) << "\n\n";
+// fitMaxPosXVeto.push_back(fitFunc->GetParameter(1));
+// }
+// }
// }
-// hitNtuple = (TNtuple*) f->Get("hit");
-// hitNtuple->Draw("p13","","same");
-
-// TH1F *hist = new TH1F("hHitNtuple", "p13!! ist das :)"); // dynamic allocation (delete at the end)
-// Long64_t nentries = hitNtuple->GetEntries();
-// for (Long64_t i; i=0; i<nentries) {
-// hist->Fill()
-// }
-// delete hist;
-
-
-
-
-
-
+// // add legend
+// TString legendEntry = TString("Same Parameters: ") + (selTemp!=-1 ? Form(" T_{set}=%.1f",selTemp):"") + (selChip!=-1 ? Form(" Chip=%i",selChip):"") + (selSource!="-1" ? Form(" Source=%s",selSource.c_str()):"") + (selMatrix6480!=-1 ? Form(" Matrix=%i",selMatrix6480):"") + (selSubmatrix!=-1 ? Form(" Submatrix=%i",selSubmatrix):"") + (selRadDose!=-1 ? Form(" RadDose=%.2E",selRadDose):"");
+// leg->AddEntry((TObject*) 0, legendEntry, "");
+// leg->SetTextSize(0.02);
+// leg->Draw();
-// {
-// TFile f1("essai1.root");
-// TH1F *h1=f1.Get("2");
-// TFile f2("essai2.root");
-// TH1F *h2=f2.Get("2");
-// TFile f3("essai3.root");
-// TH1F *h3=f3.Get("2");
-//
-// h2->SetLineColor(kRed);
-// h3->SetLineColor(kGreen);
-//
-// h1->Draw("l");
-// h2->Draw("same");
-// h3->Draw("same");
-// }
-
-
-
-// Int_t runNo = 29012;
-// Int_t submatrix = 0;
-//
-//
-//
-//
-// // Variables used to store the data
-// Int_t totalSize = 0; // Sum of data size (in bytes) of all events
-// Int_t eventSize = 0; // Size of the current event
-// TBranch *eventSizeBranch = 0; // Pointer to the event.fEventsize branch
+// // print all information
+// TString runDetailsHeader = TString("runNo\tT\tTSens\tChip\tSource\tMatrix\tRadDose\t\tSubmtrx")
+// + (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" ? "\tErrL":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tErrR":"")
+// + (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? "\tSumColPeak":"")
+// + (DRAW_SUMMED_HISTOGRAMS && DRAW_FITS ? "\tSumColErr":"")
+// + (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? "\tCalPeak":"")
+// + (DRAW_VETO_HISTOGRAMS && DRAW_FITS ? "\tCalPkEr":"")
+// + "\tNoise"
+// + "\tNseLo17"
+// + "\tNseHi17"
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tCCESeed":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tCCESdEr":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tCCESum":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="gaus" ? "\tCCESmEr":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tS/N":"")
+// + (DRAW_SINGLE_HISTOGRAMS && DRAW_FITS && TString(FIT_FUNC)=="landau" ? "\tS/NErr":"")
+// + "\n";
+//
+// if (PRINT_HEADER)
+// cout << endl << runDetailsHeader;
//
-// // open the file
-// TString path = TString(OUTPUTDATAPATH) + Form("%i/%i_0_%i.root", runNo, runNo, submatrix);
-// TFile *f = TFile::Open(path);
-//
-// if (f == 0) {
-// // if we cannot open the file, print an error message and return immediatly
-// printf("Error: cannot open %s\n", path);
-// return;
+// for (UInt_t i=0; i<runDetails.size(); ++i){
+// // 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
+// << endl;
// }
-//
-// // get a pointer to the tree
-// TNtuple* hitNtuple = (TNtuple*) f->Get("hit");
-//
-// hitNtuple->Draw("p13");
-// // hitNtuple->Draw("p10","","same");
-//
-// delete f;
-//
-// Int_t runNo = 29016;
-// Int_t submatrix = 0;
-// TString path = TString(OUTPUTDATAPATH) + Form("%i/%i_0_%i.root", runNo, runNo, submatrix);
-// TFile *f = TFile::Open(path);
-// if (f == 0) {
-// // if we cannot open the file, print an error message and return immediatly
-// printf("Error: cannot open %s\n", path);
-// return;
-// }
-// hitNtuple = (TNtuple*) f->Get("hit");
-// hitNtuple->Draw("p13","","same");
-
Runnr. Date Temperature T-sensor Events Chip Source 64/80 Radiation dose Vbo Done? Allowed T-range Comment Vba Vdiode
29012 5/17/2013 20 21.5 50000 1 Fe55 64 0.00E+00 3.20 y 19 … 21 All the yellow marked measurements were made with a baseline in the range between 500 and 3500 for -20°C and were not modified for the other highlighted measurements. 0.24 2.50
-29013 5/17/2013 20 21.5 50000 1 Fe55 80 0.00E+00 3.12 y 0.24 2.50
+29013 5/17/2013 20 21.5 50000 1 Fe55 80 0.00E+00 3.12 y 0.44 2.50
29014 5/17/2013 -20 -15.5 50000 1 Fe55 64 0.00E+00 3.19 y -12 … -19 0.24 2.50
-29015 5/17/2013 -20 -16.0 50000 1 Fe55 80 0.00E+00 3.12 y 0.24 2.50
+29015 5/17/2013 -20 -16.0 50000 1 Fe55 80 0.00E+00 3.12 y 0.44 2.50
29016 5/18/2013 -70 -58.0 50000 1 Fe55 64 0.00E+00 3.18 y -30 … -60 0.24 2.50
-29017 5/18/2013 -70 -58.0 50000 1 Fe55 80 0.00E+00 3.10 y 0.24 2.50
+29017 5/18/2013 -70 -58.0 50000 1 Fe55 80 0.00E+00 3.10 y 0.44 2.50
29018 5/16/2013 20 20.0 50000 1 none 64 0.00E+00 3.20 y
29019 5/16/2013 20 20.5 50000 1 none 80 0.00E+00 3.12 y
29033 5/28/2013 -20 -16.5 50000 1 Fe55,dimmed 80 0.00E+00 3.12 y same as run 29015, but dimmed
-29034 5/29/2013 20 21.5 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y 19 … 21
-29035 5/29/2013 20 21.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y
+29034 5/29/2013 20 21.5 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y 19 … 21 0.00 2.53
+29035 5/29/2013 20 21.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y 0.00 2.53
29036 5/29/2013 -20 -16.0 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y -12 … -19
29037 5/30/2013 -20 -16.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y
29038 5/30/2013 -70 -57.0 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y -30 … -60
29051 6/2/2013 -70 -54.5 50000 2 Sr90 80 1.50E+05 3.15 y
-29052 6/7/2013 20 21.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y 19 … 21
-29053 6/7/2013 20 21.0 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y
+29052 6/7/2013 20 21.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y 19 … 21 1.62 2.52
+29053 6/7/2013 20 21.0 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y 1.72 2.52
29054 6/7/2013 -20 -15.5 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y -12 … -19
29055 6/7/2013 -20 -16.5 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y
29056 6/7/2013 -70 -58.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y -30 … -60
29069 6/10/2013 -70 -54.0 50000 4 Sr90 80 1.00E+12 1.52 y
-29070 6/14/2013 20 21.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y 19 … 21
-29071 6/14/2013 20 21.0 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y
+29070 6/14/2013 20 21.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y 19 … 21 1.86 2.52
+29071 6/14/2013 20 21.0 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y 1.69 2.52
29072 6/14/2013 -20 -16.0 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y -12 … -19
29073 6/14/2013 -20 -16.5 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y
29074 6/14/2013 -70 -58.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y -30 … -60
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
-29089 6/20/2013 20 21.0 50000 6 Fe55,dimmed 80 1.00E+13 1.53 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
+29089 6/20/2013 20 21.0 50000 6 Fe55,dimmed 80 1.00E+13 1.53 y 1.71 2.52
29090 6/20/2013 -20 -16.5 50000 6 Fe55,dimmed 64 1.00E+13 1.54 y -12 … -19
29091 6/20/2013 -20 -16.5 50000 6 Fe55,dimmed 80 1.00E+13 1.53 y
29092 6/20/2013 -70 -58.5 50000 6 Fe55,dimmed 64 1.00E+13 1.54 y -30 … -60
End of 1st measurement phase.
-29106 10/17/2013 20 41415.0 50000 1 Fe55,superdimmed 64 0.00E+00 43891.00 y
-29107 10/17/2013 20 41415.0 50000 1 Fe55,superdimmed 80 0.00E+00 41611.00 y
-29108 10/18/2013 -20 -16.0 50000 1 Fe55,superdimmed 64 0.00E+00 43525.00 y BULLSHIT (eingefroren)
-29109 10/18/2013 -20 -16.5 50000 1 Fe55,superdimmed 80 0.00E+00 41611.00 y BULLSHIT (eingefroren)
-29110 10/18/2013 -70 -48.5 50000 1 Fe55,superdimmed 64 0.00E+00 43160.00 y BULLSHIT (eingefroren)
-29111 10/18/2013 -70 41415.0 50000 1 Fe55,superdimmed 80 0.00E+00 41550.00 y BULLSHIT (eingefroren)
+29106 10/17/2013 20 21.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y
+29107 10/17/2013 20 21.0 50000 1 Fe55,superdimmed 80 0.00E+00 3.12 y
+29108 10/18/2013 -20 -15.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.19 y BULLSHIT (eingefroren)
+29109 10/18/2013 -20 -16.0 50000 1 Fe55,superdimmed 80 0.00E+00 3.12 y BULLSHIT (eingefroren)
+29110 10/18/2013 -70 -54.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.18 y BULLSHIT (eingefroren)
+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 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
+
29012 5/17/2013 20 21.5 50000 1 Fe55 64 0.00E+00 3.20 y 19 … 21 All the yellow marked measurements were made with a baseline in the range between 500 and 3500 for -20°C and were not modified for the other highlighted measurements. 0.24 2.50
-29013 5/17/2013 20 21.5 50000 1 Fe55 80 0.00E+00 3.12 y 0.24 2.50
+29013 5/17/2013 20 21.5 50000 1 Fe55 80 0.00E+00 3.12 y 0.44 2.50
29014 5/17/2013 -20 -15.5 50000 1 Fe55 64 0.00E+00 3.19 y -12 … -19 0.24 2.50
-29015 5/17/2013 -20 -16.0 50000 1 Fe55 80 0.00E+00 3.12 y 0.24 2.50
+29015 5/17/2013 -20 -16.0 50000 1 Fe55 80 0.00E+00 3.12 y 0.44 2.50
29016 5/18/2013 -70 -58.0 50000 1 Fe55 64 0.00E+00 3.18 y -30 … -60 0.24 2.50
-29017 5/18/2013 -70 -58.0 50000 1 Fe55 80 0.00E+00 3.10 y 0.24 2.50
+29017 5/18/2013 -70 -58.0 50000 1 Fe55 80 0.00E+00 3.10 y 0.44 2.50
29018 5/16/2013 20 20.0 50000 1 none 64 0.00E+00 3.20 y
29019 5/16/2013 20 20.5 50000 1 none 80 0.00E+00 3.12 y
29033 5/28/2013 -20 -16.5 50000 1 Fe55,dimmed 80 0.00E+00 3.12 y same as run 29015, but dimmed
-29034 5/29/2013 20 21.5 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y 19 … 21
-29035 5/29/2013 20 21.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y
+29034 5/29/2013 20 21.5 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y 19 … 21 0.00 2.53
+29035 5/29/2013 20 21.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y 0.00 2.53
29036 5/29/2013 -20 -16.0 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y -12 … -19
29037 5/30/2013 -20 -16.0 50000 2 Fe55,dimmed 80 1.50E+05 3.15 y
29038 5/30/2013 -70 -57.0 50000 2 Fe55,dimmed 64 1.50E+05 3.15 y -30 … -60
29051 6/2/2013 -70 -54.5 50000 2 Sr90 80 1.50E+05 3.15 y
-29052 6/7/2013 20 21.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y 19 … 21
-29053 6/7/2013 20 21.0 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y
+29052 6/7/2013 20 21.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y 19 … 21 1.62 2.52
+29053 6/7/2013 20 21.0 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y 1.72 2.52
29054 6/7/2013 -20 -15.5 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y -12 … -19
29055 6/7/2013 -20 -16.5 50000 4 Fe55,dimmed 80 1.00E+12 1.52 y
29056 6/7/2013 -70 -58.0 50000 4 Fe55,dimmed 64 1.00E+12 1.52 y -30 … -60
29069 6/10/2013 -70 -54.0 50000 4 Sr90 80 1.00E+12 1.52 y
-29070 6/14/2013 20 21.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y 19 … 21
-29071 6/14/2013 20 21.0 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y
+29070 6/14/2013 20 21.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y 19 … 21 1.86 2.52
+29071 6/14/2013 20 21.0 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y 1.69 2.52
29072 6/14/2013 -20 -16.0 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y -12 … -19
29073 6/14/2013 -20 -16.5 50000 5 Fe55,dimmed 80 3.00E+12 1.55 y
29074 6/14/2013 -70 -58.5 50000 5 Fe55,dimmed 64 3.00E+12 1.57 y -30 … -60
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
-29089 6/20/2013 20 21.0 50000 6 Fe55,dimmed 80 1.00E+13 1.53 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
+29089 6/20/2013 20 21.0 50000 6 Fe55,dimmed 80 1.00E+13 1.53 y 1.71 2.52
29090 6/20/2013 -20 -16.5 50000 6 Fe55,dimmed 64 1.00E+13 1.54 y -12 … -19
29091 6/20/2013 -20 -16.5 50000 6 Fe55,dimmed 80 1.00E+13 1.53 y
29092 6/20/2013 -70 -58.5 50000 6 Fe55,dimmed 64 1.00E+13 1.54 y -30 … -60
End of 1st measurement phase.
-29106 10/17/2013 20 41415.0 50000 1 Fe55,superdimmed 64 0.00E+00 43891.00 y
-29107 10/17/2013 20 41415.0 50000 1 Fe55,superdimmed 80 0.00E+00 41611.00 y
-29108 10/18/2013 -20 -16.0 50000 1 Fe55,superdimmed 64 0.00E+00 43525.00 y BULLSHIT (eingefroren)
-29109 10/18/2013 -20 -16.5 50000 1 Fe55,superdimmed 80 0.00E+00 41611.00 y BULLSHIT (eingefroren)
-29110 10/18/2013 -70 -48.5 50000 1 Fe55,superdimmed 64 0.00E+00 43160.00 y BULLSHIT (eingefroren)
-29111 10/18/2013 -70 41415.0 50000 1 Fe55,superdimmed 80 0.00E+00 41550.00 y BULLSHIT (eingefroren)
+29106 10/17/2013 20 21.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.20 y
+29107 10/17/2013 20 21.0 50000 1 Fe55,superdimmed 80 0.00E+00 3.12 y
+29108 10/18/2013 -20 -15.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.19 y BULLSHIT (eingefroren)
+29109 10/18/2013 -20 -16.0 50000 1 Fe55,superdimmed 80 0.00E+00 3.12 y BULLSHIT (eingefroren)
+29110 10/18/2013 -70 -54.0 50000 1 Fe55,superdimmed 64 0.00E+00 3.18 y BULLSHIT (eingefroren)
+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 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
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
#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_4.11.2013_13:57/"
+#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
