Changeset 197 in JEM-EUSO

Timestamp:

Sep 27, 2013, 9:36:51 AM (11 years ago)

Author:

dagoret

Message:

bug due to a comment left corrected

Location:

equalization_gain/trunk

Files:

• MySetup.sh (modified) (1 diff)
• README (modified) (1 diff)
• src/PlotScurvesAll-diffFit.cc (modified) (4 diffs)
• src/PlotScurvesAll-diffFit.d (deleted)
• src/PlotScurvesAll-diffFit.h (modified) (3 diffs)
• src/PlotScurvesAllMain (modified) (previous)
• src/PlotScurvesAllMain.d (deleted)

equalization_gain/trunk/MySetup.sh

@@ -1 +1 @@
 #!/bin/sh
 
-export srcDIR="/Users/dagoret-campagnesylvie/MacOsX/EUSO-BALLOON/GainEqualisation/analysis/src"
-export dataDIR="/Users/dagoret-campagnesylvie/MacOsX/EUSO-BALLOON/GainEqualisation/analysis/data/"
+export srcDIR="/Users/dagoret-campagnesylvie/MacOsX/EUSO-BALLOON/GainEqualisation/SVN/analysis/trunk/src"
+export dataDIR="/Users/dagoret-campagnesylvie/MacOsX/EUSO-BALLOON/GainEqualisation/SVN/analysis/trunk/data/"

equalization_gain/trunk/README

@@ -5 +5 @@
 *
 **********************************************************
+
+Repository
+===========
+
+to see the projects for JEM-EUSO just dou :
+
+==> svn ls https://svn.lal.in2p3.fr/projects/JEM-EUSO
+
+==> svn ls https://svn.lal.in2p3.fr/projects/JEM-EUSO/equalization_gain
+branches/
+tags/
+trunk/
+
+
 The Gain Equalisation files are inside the analysis directory.
 

equalization_gain/trunk/src/PlotScurvesAll-diffFit.cc

@@ -216 +216 @@
   ofstream out(kFnameOut.Data()); 
 
-  // first loop on channels
+  // first loop on channels to fill histonames
   //----------------------
   for(ch=0;ch<kNch;ch++){
@@ -799 +799 @@
 //--------------------------------------------------------------------------------------------------------
 void Analysis::PlotScurve(
-        /*      Int_t    kDrawCh=31,
+                Int_t    kDrawCh=31,
                 TString  fname=kFname,
                 TString  fname_ped=kFnamePed, //ASIC_D ped_file
@@ -814 +814 @@
 {
 
+  // reset all variables
+  Init();
+
+  //OpenGainFile(fname_gain);
+
+  FillHistoName();
+
+  ReadGainFile(fname_gain,kBoard);
+}
+//--------------------------------------------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------------------------------------------
+// Initialize the program by reseting counters or variables and create objects
+//--------------------------------------------------------------------------------------------------------
+void Analysis::Init()
+{
+
   if(input!=0) delete input;
   if(input_ped!=0) delete input_ped;
@@ -831 +849 @@
 
 
+  // create canvas
+
+  c1                = new TCanvas(kCanvasName1,kCanvasTitle1,200,20,1100,600); //smoothed S-curves
+  c1->Divide(11,6);
+  c2                = new TCanvas(kCanvasName2,kCanvasTitle2,300,50,1100,600); //differentiated histograms of smoothed s-curves
+  c2->Divide(11,6);
+  c5                = new TCanvas(kCanvasName5,kCanvasTitle5,400,70,1100,600);
+  c5->Divide(11,6);
+  c3                = new TCanvas(kCanvasName3,kCanvasTitle3,500,100,1200,300);
+  c3->Divide(4,1);
+  c4                = new TCanvas(kCanvasName4,kCanvasTitle4,900,400,300,300);
+
+
+  // create histos and function
+  hist_scurve       = new TH1D[kNch]; // smoothed S-curves
+  hist_scurve_ped   = new TH1D[kNch]; // S-curves for pedestal runs
+  hist_diff         = new TH1D[kNch]; // differentiated histograms of smoothed s-curves
+  hist_diff_fit     = new TH1D[kNch]; // differentiated histograms of fitted s-curves
+  fit_sdiff         = new TF1[kNch];  // function fitted on differentiated histograms of fitted s-curves
+
+
+ for(ch=0;ch<kNch;ch++){
+    dac_half[ch]=0;
+    dac_ped[ch]=0;
+    dac_spe[ch]=0;
+    dac_sfit[ch]=0;
+    gain_org[ch]=0;
+    mean_diff[ch]=0;
+    sfit_chiS[ch]=0;
+    flag_fit[ch]=0;
+ }
+
+
   
 }
 //--------------------------------------------------------------------------------------------------------
 
+//--------------------------------------------------------------------------------------------------------
+// Initialize the program by reseting counters or variables and create objects
+//--------------------------------------------------------------------------------------------------------
+void Analysis::OpenGainFile(const TString fname_gain)
+{
+  TString input_gain(fname_gain);
+  ifstream in_gain(input_gain.Data());
+  if(!in_gain.good()){
+    cerr << "File " << input_gain << " open failed!" << endl;
+    exit(-1);
+  }
+  in_gain.close();
+  in_gain.clear();
+
+}
+//--------------------------------------------------------------------------------------------------------
+
+//--------------------------------------------------------------------------------------------------------
+// Fill Histogram name
+//--------------------------------------------------------------------------------------------------------
+void Analysis::FillHistoName()
+{
+  
+// first loop on channels to fill histonames
+  //----------------------
+  for(ch=0;ch<kNch;ch++){
+    
+    name_hist[ch]           ="hist";
+    name_hist[ch]          +=ch;
+    name_hist_ped[ch]       ="hist_ped";
+    name_hist_ped[ch]      +=ch;
+    name_hist_diff[ch]      ="hist_diff";
+    name_hist_diff[ch]     +=ch;
+    name_hist_diff_fit[ch]  ="hist_diff_fit";
+    name_hist_diff_fit[ch] +=ch;
+    name_fit_diff[ch]       ="fit_diff";
+    name_fit_diff[ch]      +=ch;
+  }
+
+}
+//--------------------------------------------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------------------------------------------
+// Read the gain
+//--------------------------------------------------------------------------------------------------------
+
+void Analysis::ReadGainFile(const TString fname_gain,Bool_t kBoard)
+{
+  TString input_gain(fname_gain);
+  ifstream in_gain(input_gain.Data());
+  if(!in_gain.good()){
+    cerr << "File " << input_gain << " open failed!" << endl;
+    exit(-1);
+  }
+ // read gain
+  Double_t da;
+  Int_t count=0;
+  ch=0;
+  while(in_gain >> da){
+    if(!kBoard && (count>=128 && count<192)){
+      gain_org[ch]=da;
+      cerr << "!kBoard, count, gain_org[" << ch << "] "
+           << !kBoard << " " << count << " " << gain_org[ch] << endl;
+      ch++;
+    }
+    if(kBoard && (count>=192 && count<256)){
+      gain_org[ch]=da;
+      cerr << "kBoard, count, gain_org[" << ch << "] "
+           << kBoard << " " << count << " " << gain_org[ch] << endl;
+      ch++;
+    }
+    count++;
+  }// end if first loop on all channels      
+  in_gain.close();
+  in_gain.clear();
+ 
+
+}
+//--------------------------------------------------------------------------------------------------------

equalization_gain/trunk/src/PlotScurvesAll-diffFit.h

@@ -6 +6 @@
 #define _PLOTSCURVESALLDIFFFIT_H_
 
+#include <iostream>
+#include "TString.h"
 #include "PlotScurvesAll-diffFitConst.h"
 
+
+class TCanvas;
+class TH1D;
+class TF1;
 
 void PlotScurve(
@@ -61 +67 @@
 private:
 
+ void Init();
+ void OpenGainFile(const TString);
+ void FillHistoName();
+ void ReadGainFile(const TString,Bool_t);
+
 // common variables to the analysis
      TString *input;
@@ -83 +94 @@
      Double_t flag_fit[kNch];
 
+     //canvas pointers
+
+     TCanvas *c1; //smoothed S-curves
+     TCanvas *c2; //differentiated histograms of smoothed s-curves
+     TCanvas *c5;
+     TCanvas *c3;            
+     TCanvas *c4;
+
+     //histograms
+
+     TString  name_hist[kNch];
+     TString  name_hist_ped[kNch];
+     TString  name_hist_diff[kNch];
+     TString  name_hist_diff_fit[kNch];
+     TString  name_fit_diff[kNch];
+     TH1D    *hist_scurve; // smoothed S-curves
+     TH1D    *hist_scurve_ped; // S-curves for pedestal runs
+     TH1D    *hist_diff; // differentiated histograms of smoothed s-curves
+     TH1D    *hist_diff_fit; // differentiated histograms of fitted s-curves
+     TF1     *fit_sdiff;  // function fitted on differentiated histograms of fitted s-curves
+
+
+     
+
 };