/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Renaud/ROOT/linearityoffset.C

Go to the documentation of this file.

//#include "Riostream.h"

linearityoffset()
{
        //.X show.C(1,2)
        
        TString directory = "/gpfs/LAPP-DATA/LC/Detecteurs/MicroMegas/data/MICROROC/MR1/TEST_DAC/RESULTS/";
        
        const Int_t NbOfRun = 16 ; // one run for each pedestal of channel 7
//      const Int_t NbOfPoints = 16 ; // pied: 0 to 15
        const Int_t NbOfChannels = 64 ;

        Double_t n, a, b, c, d ;

        ///////////////////////////
        // Plot noise vs t
        ///////////////////////////
//      Float_t temperature[]={25.5, 47.0, 68.5, 87.0, 105.0};
        Double_t x[NbOfRun];
        for (Int_t i=0;i<NbOfRun;i++)
        {
                x[i]=i ;
        }
        
        Double_t data_inflex[NbOfChannels][NbOfRun] ; // first dim is for x / y
        
        Double_t mean_inflex[NbOfRun] ;
        // Gather y data from files
        for (Int_t run=0;run<NbOfRun;run++) // x axis
        {       //for (Int_t ped=0;ped<NbOfPoints;ped++)
                //{
                        mean_inflex [run] = 0 ;
                        TString name_file = directory + "run" + Form("%d", 2*run+1) + "_pied" + Form("%d", ped) + ".txt" ;
                        cout << name_file << endl ;
                        ifstream fp ;
                        fp.open(name_file.Data()) ; 
                
                        if (!fp.good()) cout << "pb !";
                        for (Int_t ch=0;ch<NbOfChannels;ch++)
                        {               fp >> n >> a >> b >> c >> d ;
                                        //cout << n << " " << a << " " << b << " " << c << " " << d <<  endl ;
                                        if (!fp.good()) break;
                                        // fill data
                                        data_inflex [ch][run] = a ; // DAC units 
                                        // fill mean for all channels
                                        mean_inflex [run] += data_inflex [ch][run] ;
                        }
                        fp.close();     
                        mean_inflex [run] /= NbOfChannels ;
                //}
        }
        
//      cout << "Temperature" ;
//      for (Int_t run=0;run<NbOfPoints;run++) // x axis
//      {       cout << "\t" << temperature[run] ;
//      }
//      cout << endl << "Offset" << "\t" ;
//      for (Int_t run=0;run<NbOfPoints;run++) // x axis
//      {       cout << "\t" << data_noise[1][run] ;
//      }
//      cout << endl ;
        
//      exit(-1);
        
        //----------------------------------------------------------------
        //                       BUILD GRAPH 
        //----------------------------------------------------------------      
        TCanvas * c1 = new TCanvas("c1","MICROROC offset dac");
        TMultiGraph * mg = new TMultiGraph("mg","");
        TGraph * gr[NbOfChannels];
        for(Int_t igraph=0;igraph<NbOfRun;igraph++)     // one graph per temperature
        {       gr [igraph] = new TGraph(NbOfPoints, x, mean_inflex [igraph] ); // TGraph (n, x, y)
                gr [igraph]->SetLineColor(kBlue+igraph) ;
                gr [igraph]->SetLineWidth(1);
                gr [igraph]->SetMarkerColor(kBlue+igraph);
                gr [igraph]->SetMarkerStyle(8);
                gr [igraph]->SetMarkerSize(1);
                mg->Add(gr [igraph]);
        }
        
        mg->Draw("ALP");
        mg->GetXaxis()->SetTitle("Pedestal [DACu]");
        mg->GetXaxis()->SetRangeUser(0.,110.);
        mg->GetYaxis()->SetTitle("Inflexion Point [DACu]");
//      mg->GetYaxis()->SetRangeUser(2.4,2.5);
        mg->GetYaxis()->SetTitleOffset(1.1);
        mg->Draw("ALP"); // to refresh
//      Legend
//      TString temp_string = Form ("%f", x[1]) ;
//      temp_string 
        TLegend * leg = new TLegend(0.1,0.75,0.25,0.9) ;
        for(Int_t igraph=0;igraph<NbOfRun;igraph++)     // one graph per temperature
        {       
                TString temp_string = Form ("%5.1f #circC", temperature[igraph]) ;
                leg->AddEntry(gr [igraph], temp_string, "lp");
        }
        leg->SetFillStyle(1);
        leg->SetFillColor(0);
        leg->Draw();
        
        return 0 ;

}

---