]> jspc29.x-matter.uni-frankfurt.de Git - radhard.git/commitdiff
Dennis Version ChargeSpektrum
authorDennis Doering <doering@physik.uni-frankfurt.de>
Wed, 17 Jun 2015 12:01:29 +0000 (14:01 +0200)
committerDennis Doering <doering@physik.uni-frankfurt.de>
Wed, 17 Jun 2015 12:01:29 +0000 (14:01 +0200)
MABS_run_analyzer/ChargeSpektrum.c

index 20687607fc18ba46f91d634bf8f0daf696406a49..90d46b86504d9ea984b46629729269ce7b4e9925 100644 (file)
@@ -176,7 +176,7 @@ void ChargeSpektrum(TString runnumber = "")
                 {
             //         gROOT->SetBatch(kTRUE);
                     //           runs[runi]->plotSeed();
-                    runs[runi]->plotSeedThresholdCalibrated();
+          //          runs[runi]->plotSeedThresholdCalibrated();
 //                     runs[runi]->plotSeedThreshold();
 //                     runs[runi]->plotSeed();
             //         runs[runi]->plotSum();
@@ -198,6 +198,7 @@ void ChargeSpektrum(TString runnumber = "")
 //     writeObservableToFile("seed threshold calibrated");
     writeObservableToFile("seed threshold");
        writeObservableToFile("seed threshold calibrated");
+       writeObservableToFile("sum threshold");
 }
 
 Int_t ReadRunList(std::vector<int>* runlist)
@@ -284,7 +285,8 @@ Bool_t writeObservableToFile()
     for(Int_t runi=0;runi<numberRuns;runi++) // loop over runs read from file
     {
                TString labbookTempSens=Form("%d",runs[runi]->labbook.tempSens);
-        headerInfo+=runs[runi]->runcode+"\t"+labbookTempSens+"\t";
+               TString labbookFrames=Form("%d",runs[runi]->labbook.frames_foundDB);
+        headerInfo+=runs[runi]->runcode+"\t"+labbookTempSens+labbookFrames+"\t\t";
     }
     TH1F* plothistogrampointer = *runs[0]->plothistogrampointer;
     TString runnumberListe=Form("%s_",plothistogrampointer->GetName());
@@ -297,7 +299,7 @@ Bool_t writeObservableToFile()
             {
                 plothistogrampointer = *runs[runi]->plothistogrampointer;
                 runnumberListe+=Form("%d_",runs[runi]->labbook.runnumber);
-                header+=Form("%d\t%d_%s\t", runs[runi]->labbook.runnumber, runs[runi]->labbook.runnumber, plothistogrampointer->GetName());
+                header+=Form("%d\t%d_%s\t\t", runs[runi]->labbook.runnumber, runs[runi]->labbook.runnumber, plothistogrampointer->GetName());
             }
         }
        }
@@ -305,17 +307,19 @@ Bool_t writeObservableToFile()
     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 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;
-                    outline=Form("%.1f\t%.1f\t",plothistogrampointer->GetBinCenter(bini),plothistogrampointer->GetBinContent(bini));
+                                       Float_t histoNorm=plothistogrampointer->GetBinContent(bini)/runs[runi]->labbook.frames_foundDB*10000000;
+                    outline=Form("%.1f\t%.1f\t%.1f\t",plothistogrampointer->GetBinCenter(bini),plothistogrampointer->GetBinContent(bini),histoNorm);
                     *fout<<outline;
                 }
             }