/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Renaud/analyze_prod.cpp

Go to the documentation of this file.

// compile with:
// gcc `root-config --cflags` `root-config --glibs` analyze_prod.cpp -o analyze_prod


#include <map>
#include "TKey.h"
#include <TROOT.h>
#include <TStyle.h>
#include <TRint.h>
#include <TKey.h>
#include <TTree.h>
#include <TCanvas.h>
#include "TH1.h"
#include "TF1.h"
#include <TH1I.h>
#include <TH2I.h>
#include <TFile.h>
#include <stdlib.h>
#include <TSystem.h>
#include <iostream>
#include "TApplication.h"
#include <sstream>
//#include "Log.hh"
#include <Rtypes.h>
#include <TCanvas.h>
#include "TColor.h"
#include "TFunction.h"
#include "TAttLine.h"
#include "TPaveStats.h"
//#include "MicroException.hh"
#include <string>
#include "Riostream.h"

using namespace std;


int main(int argc, char**argv)
//analyze_prod(const Int_t prod)
{       
        gROOT->Reset();

//      gROOT->SetStyle("Plain");
        gStyle->SetFrameFillStyle(0);
        gStyle->SetPalette(1);
        gStyle->SetOptStat(111111);
        gStyle->SetFrameFillStyle(0);
        gStyle->SetOptStat("");
        gStyle->SetOptFit(0111);
        gStyle->SetStatX(0.9);
        gStyle->SetStatY(0.9);
        gStyle->SetStatH(0.1);
        gStyle->SetStatW(0.1);

//      gStyle->SetOptStat("") ;
//      gStyle->SetOptFit(0000);        
        
        gStyle->SetOptStat("");
        gStyle->SetOptFit(10001);
        
        if (argc !=2)
        {
                cout << "usage: " << argv[0] << " prod_run_number" << endl ;
                exit(-1) ;
        }
        Int_t prod = atoi(argv[1]);
        
        const Int_t NbOfChannels = 64 ;
        
        const TString directory = Form("/gpfs/LAPP-DATA/LC/Detecteurs/MicroMegas/data/MICROROC/PROD_%d/RESULTS/", prod);
        cout << directory << endl ;
        
        Int_t start, stop ;
        switch (prod)
        {       case 1:
                        start = 1 ;
                        stop = 341 ;
                        //stop = 30 ;
                        break ;
                case 2:
                        start = 342 ;
                        stop = 1136 ;
                        break ;
                default:
                        cout << "This production run does not exist yet !" << endl ;
                        exit(-1) ;
        }
        
        TApplication theApp("App", &argc, argv);                                        // see $ROOTSYS/test/hworld.cxx
        //TApplication *theApp  = NULL;                         // from Jean
        //theApp  = new TRint("App",NULL,NULL );        // from Jean
        
//      for (Int_t chip=start;chip<=stop;chip++)
//      {       
//              const TString name_file_gain = directory +  Form("chip_%d_gain_OK.txt", chip) ;
//              const TString name_file_gain = Form("%schip_%d_gain_OK.txt", directory, chip) ;
//              cout << name_file_gain << endl ;
//      }
        
        // Variable definition
        Double_t n, gain ;      // to read text file
        Double_t pedestal_inflex, pedestal_inflex_err, pedestal_noise, pedestal_noise_err ; // to read text file
        Double_t data_inflex, data_inflex_err, data_noise, data_noise_err ;     // to read text file
        TH1F * h_gain = new TH1F ("h_gain", "",100, 0.12, 0.16) ; // histo to store result
        h_gain->SetXTitle("");
        TH1F * h_ped = new TH1F ("h_ped", "", 1024, 0, 1024) ;
        h_ped->SetXTitle("Threshold [DACu]");
        h_ped->SetYTitle("count");
        TH1F * h_data = new TH1F ("h_data", "", 1024, 0, 1024) ;
        h_data->SetXTitle("Threshold [DACu]");
        h_data->SetYTitle("count");
        ifstream file_gain, file_data, file_pedestal ;
        for (Int_t chip=start;chip<=stop;chip++)
        {       const TString name_file_gain = directory + Form("chip_%d_gain_OK.txt", chip) ;
                file_gain.open(name_file_gain.Data()) ;
                if (!file_gain.good())
                {       cout << "Chip " << chip << " has been rejected during test" << endl ;
                        //TString name_file = Form("chip_%d_gain_BAD.txt", chip) ;
                        //if (!fp.good()) cout << "Error opening " << name_file << endl ;
                        //break ;
                }
                else
                {       const TString name_file_pedestal = directory + Form("chip_%d_pedestal.txt", chip) ;
                        file_pedestal.open(name_file_pedestal.Data()) ;
                        const TString name_file_data = directory + Form("chip_%d_data.txt", chip) ;
                        file_data.open(name_file_data.Data()) ;
                        for (Int_t ch=0;ch<NbOfChannels;ch++)
                        {       file_gain >> gain ;
                                if (!file_gain.good()) break;
                                file_pedestal >> n >> pedestal_inflex >> pedestal_inflex_err >> pedestal_noise >> pedestal_noise_err ;
                                file_data >> n >> data_inflex >> data_inflex_err >> data_noise >> data_noise_err ;
                                h_gain->Fill(gain) ;
                                h_ped->Fill(pedestal_inflex) ;
                                h_data->Fill(data_inflex) ;
                        }
                }
                file_gain.close();      
                file_pedestal.close();  
                file_data.close();      
                
                
        }
        
        TCanvas * c1 = new TCanvas("c1","MICROROC gain");
        h_gain->SetLineColor(kBlue+1);
        h_gain->Draw("");
        h_gain->Fit("gaus","RQ", "SAME", 0.13 , 0.15);
        h_gain->GetFunction("gaus")->SetLineColor(kRed);
        h_gain->GetXaxis()->SetTitle("Conversion ratio [fCV/DACu]");
        h_gain->GetYaxis()->SetTitle("Count");
        
        
        ////////////////////////////////////////////////////////////////////////////////////////
        TCanvas * c2 = new TCanvas("c2","MICROROC data and pedestal");
        c2->cd() ;

        h_ped->Draw("");
        //h_ped->SetLineColor(kBlue+1);
        h_ped->Fit("gaus","RQ", "SAME", 160, 210);
        h_ped->GetFunction("gaus")->SetLineColor(kBlue+1);

        h_data->Draw("SAMES");
        //h_data->SetLineColor(kGreen+1);
        h_data->Fit("gaus","RQ", "SAME", 500, 580);
        h_data->GetFunction("gaus")->SetLineColor(kOrange+5);

        c2->Update(); // to create TGraph !!! COMPULSORY !!

        TPaveStats * p1 = (TPaveStats*) h_ped->GetListOfFunctions()->FindObject("stats");
        h_ped->GetListOfFunctions()->Print();
//      p1->Print();
        h_ped->GetListOfFunctions()->Remove(p1);
        h_ped->SetStats(0);
        p1->SetTextColor(kBlue+1);
        p1->SetX1NDC(0.75);
        p1->SetY1NDC(0.8);
        p1->SetY2NDC(0.9);
        p1->Draw();


        // NEED SAMES (!=SAME) in Drawin histo !!!!! http://root.cern.ch/root/html/THistPainter.html#THistPainter:Paint find sames in text
        //h_data->GetListOfFunctions()->Print();
        TPaveStats * p2 = (TPaveStats*) h_data->GetListOfFunctions()->FindObject("stats");
        p2->Print();
        h_data->GetListOfFunctions()->Remove(p2);
        h_data->SetStats(0);
        p2->SetTextColor(kOrange+5);
        p2->SetX1NDC(0.75);
        p2->SetX2NDC(0.9);
        p2->SetY1NDC(0.7);
        p2->SetY2NDC(0.8);
        p2->Draw();



        theApp.Run(kTRUE) ;
        //theApp->Run();
        //delete theApp;
        return 0 ;

}

---