MAJ des fichiers analyse des datas

.gitignore

@@ -62,6 +62,7 @@ _deps
 
 *.pcm
 *.root
+*.txt
 
 
 

.vscode/settings.json

 {
-    "C_Cpp.errorSquiggles": "disabled"
+    "C_Cpp.errorSquiggles": "disabled",
+    "files.associations": {
+        "new": "cpp",
+        "string_view": "cpp"
+    }
 }
\ No newline at end of file

README.md

@@ -38,7 +38,7 @@ Personnaly, I used the vrml.mac but you can create another one. Just to remember
 ## Commit #1 le 09/01/2024 [AIFIRA2023Sim.0.0.0]
 - Simulation issue de la version TP_Simulation.0.9.0 -> Commit initial avec changement des noms des fichiers et des variables
 
-## Commit #1 le 10/01/2024 [AIFIRA2023Sim.0.1.0]
+## Commit #2 le 10/01/2024 [AIFIRA2023Sim.0.1.0]
 1 Allègement du code issus de TP_Simulation afin de ne garder que l'essentiel 
 - Suppression des fichiers relatifs aux fibres optiques
 - Suppression des fichiers liées à d'autres géométries

Resultats/Analyse/Data/Analyse.cc

+#include <iostream>
+#include <fstream>
+#include <sstream>
+#include <vector>
+#include "Analyse.hh"
+
+using namespace std;
+
+void Analyse()
+{
+    gStyle->SetOptStat(0);
+    path = "/mnt/hgfs/shared/Simulations/AIFIRA2023_Simulation/Resultats/Analyse/Data/CAMERA_AIFIRA/ZnS-1mm/Linearite/10ms/Linearite_";
+    //path = "/mnt/hgfs/shared/Simulations/AIFIRA2023_Simulation/Resultats/Analyse/Data/CAMERA_AIFIRA/ZnS-1mm/Degradation/Film_degradation_";
+    //path = "/mnt/hgfs/shared/Simulations/AIFIRA2023_Simulation/Resultats/Analyse/Data/ORCA/EJ262-0.1mm/irradiation_5ms/";
+    //path = "/mnt/hgfs/shared/Simulations/AIFIRA2023_Simulation/Resultats/Analyse/Data/ORCA/ZnS-1mm/ZnS_1mm_integration_camera_1s";
+    //path = "/mnt/hgfs/shared/Simulations/AIFIRA2023_Simulation/Resultats/Analyse/Data/ORCA/ZnS-1mm/ZnS_1mm_temps_integration_5s_temps_irradiation_10ms";
+
+    N_files = 1980;
+    BitCodage = 8;
+    Rebin =2; // ATTENTION : Pour l'analyse ORCA, il faut penser à changer le rebinning effectué pour la CAMERA AIFIRA
+
+    //BdF_Analyse(path, N_files);
+    // TCanvas *c1 = new TCanvas;
+    // c1->SetGrayscale();
+    // BdF->Draw("colz");
+
+
+    // for (int i=0; i<=10; i++)
+    // {
+    //     filename = path + to_string(i) + ".txt";
+         filename = path + "1797.txt";
+    //     //filename = path + "EJ262_0.1mm_temps_integration_500ms_temps_irradiation_10ms.txt";
+         Signal_Analyse(filename);
+         Draw_Results(Rebin);
+    //     // string png = to_string(i) + ".png";
+    //     // c1->SaveAs(png.c_str());
+    //      //c1->SaveAs("EJ262_0.1mm_temps_integration_500ms_temps_irradiation_1ms_1.png");
+    //     // c1->Modified();
+    //     // c1->Update();
+    // }
+
+   //Analyse_Degradation(path, 14, 5417);
+
+
+
+    // TFile *file_BdF = new TFile("BdF_ORCA_500ms_CLEAN.root"); // BdF ORCA
+    // BdF_Mean = (TH2D *)file_BdF->Get("BdF");
+
+    // string name = "Signal";
+
+    // Signal = Plot_from_CSV(filename.c_str(), name.c_str());
+    // Signal_true = (TH2D *)Signal->Clone("Signal_true");
+    // // Signal_true->SetName("Signal_true");
+    // ComputeBdFSubstraction(Signal_true, BdF_Mean);
+    // Signal_true->Draw("colz");
+    // // new TCanvas;
+    // // BdF_Mean->Draw("colz");
+
+    vector<float> Value;
+    float Mean = 0;
+
+    for (int i = 1; i <= 500; i++)
+    {
+        for (int j = 1; j <= 500; j++)
+        {
+            //if (i<800 && j >800 && j<1000)
+            //
+                Value.push_back(Signal_true->GetBinContent(i, j));
+                //Value.push_back(BdF_Mean->GetBinContent(i, j));
+            //}
+        }
+    }
+
+    // // for (int i=0; i<Value.size(); i++)
+
+    Mean = CalculMoyenne(Value);
+    float EcartType = CalculEcartType(Value);
+    cout << "n = " << Value.size() << endl;
+    cout << "Mean = " << Mean << endl;
+    cout << "sigma = " << EcartType << endl;
+
+}
\ No newline at end of file

Resultats/Analyse/Data/Analyse_bis.cc

+#include <iostream>
+#include <fstream>
+#include <sstream>
+#include <vector>
+#include "Analyse.hh"
+
+using namespace std;
+
+void Analyse_bis()
+{
+
+    TFile* file1 = new TFile("BdF_ORCA_4ms_CLEAN.root");
+    TH2D* a = (TH2D*)file1->Get("BdF");
+    //a->Draw("colz");
+
+    //new TCanvas;
+    TFile* file = new TFile("BdF_ORCA_500ms_CLEAN.root");
+    TH2D* h = (TH2D*)file->Get("BdF");
+    //h->Draw("colz");
+
+    vector<float> Value;
+    float Mean = 0;
+
+    TH2D* final = (TH2D *)a->Clone("BdF");
+    
+    
+
+    for (int i = 1; i <= h->GetNbinsX(); i++)
+    {
+        for (int j = 1; j <= h->GetNbinsY(); j++)
+        {
+            final->SetBinContent(i, j, 1.452*a->GetBinContent(i,j) - h->GetBinContent(i,j));
+            //final->SetBinContent(i, j, 1.452*a->GetBinContent(i,j));
+
+
+            if (j <100 || j>1900)
+            {
+                Value.push_back(final->GetBinContent(i, j));
+            }
+               //Value.push_back(final->GetBinContent(i, j));
+            
+        }
+
+    }
+
+    //final->Draw("colz");
+    //final->SaveAs("BdF_ORCA_500ms_CLEAN_bis.root", "recreate");
+
+    Mean = CalculMoyenne(Value);
+    float EcartType = CalculEcartType(Value);
+    cout << "n = " << Value.size() << endl;
+    cout << "Mean = " << Mean << endl;
+    cout << "sigma = " << EcartType << endl;
+
+}
\ No newline at end of file