Excelfunktion angefangen

diff --git a/MABS_run_analyzer/ChargeSpektrum.c b/MABS_run_analyzer/ChargeSpektrum.c
index 316b814fa1e7061dca71d8be307d0c788aa4d89e..452f09a884177b7058d234fea574871841755823 100644 (file)
--- a/MABS_run_analyzer/ChargeSpektrum.c
+++ b/MABS_run_analyzer/ChargeSpektrum.c
@@ -143,16 +143,16 @@ void ChargeSpektrum(TString runnumber = "")
                     //         runs[runi]->plotNoise();
                     if (!isBatch)
                         gROOT->SetBatch(kFALSE);
-                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->histogram->calibrated);
-                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->histogramthreshold->calibrated);
-                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->dennismegacut->calibrated);
-                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->histogramfixedthreshold->calibrated);
+                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->histogram);
+                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->histogramthreshold);
+                    runs[runi]->compareHistogramClassVector.push_back(runs[runi]->dennismegacut);
                     runs[runi]->plotCompareHistograms();
-                    runs[runi]->plotAllHistograms(runs[runi]->histogram->calibrated);
-                    runs[runi]->plotAllHistograms(runs[runi]->dennismegacut->calibrated);
-                    runs[runi]->plot1DHistogram(runs[runi]->histogram->calibrated, runs[runi]->histogram->Seed, "landau");
-                    runs[runi]->plot1DHistogram(runs[runi]->histogram->calibrated, runs[runi]->histogram->calibrated->Seed, "landau");
-                    compareHistogramClassVector.push_back(runs[runi]->histogram->calibrated);
+                    // runs[runi]->plotAllHistograms(runs[runi]->histogram->calibrated);
+                    // runs[runi]->plotAllHistograms(runs[runi]->dennismegacut->calibrated);
+                    // runs[runi]->plot1DHistogram(runs[runi]->histogram->calibrated, runs[runi]->histogram->Seed, "landau");
+                    // runs[runi]->plot1DHistogram(runs[runi]->histogram->calibrated, runs[runi]->histogram->calibrated->Seed, "landau");
+                     compareHistogramClassVector.push_back(runs[runi]->histogram);
+                                        compareHistogramClassVector.push_back(runs[runi]->dennismegacut);
 //                     runs[runi]->integrateSr90Spectra(&runs[runi]->histogramfixedthresholdCalibrated, runs[runi]->histogramfixedthresholdCalibrated.Seed, 0); 
 //                     runs[runi]->plot1DHistogram(&runs[runi]->histogramfixedthresholdCalibrated, runs[runi]->histogramfixedthresholdCalibrated.Seed, "landau");
 //                     runs[runi]->plot1DHistogram(&runs[runi]->histogramthresholdCalibrated, runs[runi]->histogramthresholdCalibrated.Seed, "landau");
@@ -171,8 +171,8 @@ void ChargeSpektrum(TString runnumber = "")
     }
     plotAllRuns();
 //     plotAllRuns("seed threshold calibrated");
-//     setCustomPath("Excel/");
-//     writeObservableToFile("seed threshold calibrated");
+       setCustomPath("Excel/");
+ //    writeDataToExcel();
 //    writeObservableToFile("seed threshold");
 //     writeObservableToFile("seed threshold");
 //     writeObservableToFile("sum threshold");

diff --git a/MABS_run_analyzer/ChargeSpektrumFunctions.c b/MABS_run_analyzer/ChargeSpektrumFunctions.c
index 9362888fdf24baa970739861501cf044b7c783bf..40a43debe47ffb9e9bb605348ca676d035018b2b 100644 (file)
--- a/MABS_run_analyzer/ChargeSpektrumFunctions.c
+++ b/MABS_run_analyzer/ChargeSpektrumFunctions.c
@@ -177,7 +177,62 @@ Bool_t writeObservableToFile()
     
     return 0;
 }
+Bool_t writeDataToExcel()
+{      
+    system("mkdir "+ runs[0]->savepathresults + " -p");
+    TString headerInfo = "";
+    for(Int_t runi=0;runi<numberRuns;runi++) // loop over runs read from file
+    {
+        headerInfo+=runs[runi]->runcode+"\t\t\t";
+    }
+    TH1F* plothistogrampointer = runs[0]->plothistogrampointer;
+    TString runnumberListe="";
+    TString header="";
+       cout<<"komme ich hierher3"<<endl;
+       for (UInt_t histogrami=0; histogrami < compareHistogramClassVector.size(); histogrami++)
+        {
+            HistogramType* curhistogramclassp = compareHistogramClassVector.at(histogrami);
+            TH1F* curhistogramclone;
+            curhistogramclone = (TH1F*) curhistogramclassp->Seed->Clone();
+                       runnumberListe+=Form("%d_",curhistogramclassp->labbook.runnumber);
+         //   header+=Form("%d\t%s\t\t\t", runs[runi]->labbook.runnumber, curhistogramclone->GetTitle());
+                       
+               }
 
+    // TString filename= runs[0]->savepathresults + "/" + ownpath + runnumberListe + Form("%s",plothistogrampointer->GetName())+"histogram.dat";
+       // fstream* fout = new fstream(filename,ios::out);
+    // *fout << headerInfo << endl;
+    // *fout << header << endl;
+       
+       
+    // for(Int_t bini=0;bini<plothistogrampointer->GetNbinsX();bini++)
+    // {
+        // for(Int_t runi=0;runi<numberRuns;runi++) // loop over runs read from file
+        // {
+            // if (runs[runi] != nullptr)
+            // {
+                // if (!runs[runi]->error)
+                // {
+                    // plothistogrampointer = runs[runi]->plothistogrampointer;
+                    // TString outline;
+                                       // Double_t binContentNorm=plothistogrampointer->GetBinContent(bini)/runs[runi]->labbook.frames_foundDB*10000000;
+                                       // Double_t binSumme=0.0;
+                                       // for(Int_t binSummei=bini;binSummei<plothistogrampointer->GetNbinsX();binSummei++)//Sum up from bini beginning (cut condition for detection eff.)
+                                       // {
+                                               // Double_t binContentNormInSumme=plothistogrampointer->GetBinContent(binSummei)/runs[runi]->labbook.frames_foundDB*10000000;
+                                               // binSumme=binSumme+binContentNormInSumme;
+                                       // }
+                                       // outline=Form("%.1f\t%.1f\t%.1f\t%.1f\t",plothistogrampointer->GetBinCenter(bini),plothistogrampointer->GetBinContent(bini),binContentNorm,binSumme);
+                    // *fout<<outline;
+                // }
+            // }
+        // }
+        // *fout<<endl;
+    // }
+  //  fout->close();
+    
+    return 0;
+}
 
 
 Bool_t plotAllRuns()

diff --git a/MABS_run_analyzer/Run.c b/MABS_run_analyzer/Run.c
index f5a3304a25e8cbe5053f6eefb07db50a238f5c85..62ed3221ebfaf867c5baf01e10303332507ea737 100644 (file)
--- a/MABS_run_analyzer/Run.c
+++ b/MABS_run_analyzer/Run.c
@@ -816,7 +816,7 @@ Bool_t Run::binSeedSumVeto()
                             histogramfixedthreshold->Veto->Fill(processed->fFrameInfo.p[12][hiti]);    // histogram with the single pixel
                     }
                     
-                    if (processed->fFrameInfo.pixel[hiti]>200)
+                    if (processed->fFrameInfo.pixel[hiti]>800)
                     {
                         dennismegacut->numberofhits++;