Changeset 169 in JEM-EUSO

Timestamp:

May 31, 2013, 1:44:36 PM (11 years ago)

Author:

moretto

Message:

analysis ne permet pas d'avoir les courbes d'efficacité

File:

• esaf_lal/branches/camille/camille_work/KIT_phase1_study/macros/Analysis.C (modified) (6 diffs)

esaf_lal/branches/camille/camille_work/KIT_phase1_study/macros/Analysis.C

@@ -26 +26 @@
 void Analysis()
 {
+    //Déclaration d'un tableau de TProfile permettant de stocker les TProfiles créé par l'analyse des fichiers
     TProfile *hprof_multi[2];
     
-    //Initialisation & process
+    //Initialisation & process du premier fichier
     init(inputfile);
     process_profile();
+    TCanvas *c3 = new TCanvas("c3", "TProfile of the GTU distribution for 100 events for 30 ns dead time");
     hprof_multi[1] = hprof;
+    hprof_multi[1]->Draw();
     
+    //Initialisation & process du second fichier
     init(inputfile3);
     process_profile();
+    TCanvas *c4 = new TCanvas("c4", "TProfile of the GTU distribution for 100 events for 5 ns dead time");
     hprof_multi[2] = hprof;
+    hprof_multi[2]->Draw();
     
-    TCanvas *c4 = new TCanvas("c4", "Efficiency vs GTU");
+    
+    TCanvas *c4 = new TCanvas("c5", "Efficiency vs GTU");
     TProfile *hprof_efficiency = new TProfile("hprof_efficiency", "Efficiency vs GTU", 128, 0, 128);
-    
-    for (Int_t i=1; i<128; i++)
+  
+    for (Int_t i=0; i<128; i++)
     {
         Float_t bin1 = hprof_multi[1]->GetBinContent(i);
@@ -60 +67 @@
 
 //______________________________________________________________________________
-void process_profile ()
+void process_profile()
+//ici on récupÚre les données de temps pour les 100 événements afin de réaliser un TProfile
 {
     //Déclaration des histogrammes et autres
     TH1F *histo_total = new TH1F("histo_total", "Histogram of photons read by the electronics for 100 showers", 128, 0, 128);
     TH1F *histo_multi[100];
-    TProfile *hprof = new TProfile("hprof", "Average number of photons per GTU", 128, 0, 128);
+    hprof = new TProfile("hprof", "Average number of photons per GTU", 128, 0, 128);
+    
+    histo_total->Reset();
+    hprof->Reset();
     
     //Boucle sur l'ensemble des événements pour créer une distribution totale
@@ -88 +99 @@
         
     }
-
-    TCanvas *c1 = new TCanvas("c1", "Distribution totale");
-    histo_total->GetXaxis()->SetTitle("GTU");
-    histo_total->Draw();
+    
+    /* TCanvas *c1 = new TCanvas("c1", "Distribution totale");
+     histo_total->GetXaxis()->SetTitle("GTU");
+     histo_total->Draw(); */
     
     //Boucle sur l'ensemble des événements pour tracer séparément les histogrammes
-    TCanvas *c2 = new TCanvas("c2", "Collection of GTU histograms for 100 events");
-    c2->Divide(10,10);
+    //TCanvas *c2 = new TCanvas("c2", "Collection of GTU histograms for 100 events");
+    //c2->Divide(10,10);
     
     for (Int_t i=0; i<nb_entries; i++)
+    {
         histo_multi[i] = new TH1F("histo_multi", "Single shower histograms", 128, 0, 128);
+        histo_multi[i]->Reset();
+    }
     
     for (Int_t  i=0; i<100 ; i++)
@@ -118 +132 @@
         }
         
-        c2->cd(i+1);
-        histo_multi[i]->Draw();
+        //c2->cd(i+1);
+        //histo_multi[i]->Draw();
         
     }
@@ -126 +140 @@
     for (Int_t j=0; j<nb_entries; j++)
     {
-        for (Int_t i=1; i<=128; i++)
+        for (Int_t i=0; i<128; i++)
         {
             hprof->Fill(i, histo_multi[j]->GetBinContent(i));
@@ -133 +147 @@
     }
     
-    TCanvas *c3 = new TCanvas("c3", "TProfile of the GTU distribution for 100 events");
-    hprof->Draw();
-
 }