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

Go to the documentation of this file.

//#include "Riostream.h"

noise_vs_t()
{
        //.X show.C(1,2)
        //cout << chip_num << " " << ch_num << endl ;
        TString directory = "/gpfs/LAPP-DATA/LC/Detecteurs/MicroMegas/data/MICROROC/MR1/TEST_TEMPERATURE/RESULTS/";
        
        const Int_t NbOfPoints = 9 ; // x
        const Int_t NbOfChannels = 64 ;
        Int_t pied = 7 ;

        Double_t n, a, b, c, d ;

        ///////////////////////////
        // Plot noise vs t
        ///////////////////////////
        Double_t temperature[]={25.5, 36.7, 47.0, 55.8, 68.5, 78.0, 87.0, 94.5, 105.0};
        Double_t data_noise[NbOfChannels][NbOfPoints] ; // first dim is for x / y
        Double_t data_inflex[NbOfChannels][NbOfPoints] ; // first dim is for x / y
        
        Double_t mean_noise[NbOfPoints] ;
        Double_t mean_inflex[NbOfPoints] ;
        // Gather y data from files
        for (Int_t run=0;run<NbOfPoints;run++) // x axis
        {       
                mean_noise [run] = 0 ;
                mean_inflex [run] = 0 ;
                TString name_file = directory + "run" + Form("%d", run+1) + "_pied" + Form("%d", pied) + ".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 
                                data_noise [ch][run] = c * 0.13 ; // DAC unit to fC http://lappwiki01.in2p3.fr/Linear_Collider_LAPP/index.php/MicrorocProd
                                // fill mean
                                mean_inflex [run] += data_inflex [ch][run] ;
                                mean_noise [run] += data_noise [ch][run] ;
                }
                fp.close();     
                mean_inflex [run] /= NbOfChannels ;
                mean_noise [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 ;
        
        //----------------------------------------------------------------
        //                       BUILD GRAPH FOR NOISE
        //----------------------------------------------------------------      
        TCanvas * c1 = new TCanvas("c1","MICROROC offset dac");
        TMultiGraph * mg = new TMultiGraph("mg","");
        TGraph * gr[NbOfChannels];
        for(Int_t igraph=0;igraph<NbOfChannels;igraph++)
        {       gr [igraph] = new TGraph(NbOfPoints, temperature, data_noise[igraph] ); // TGraph (n, x, y)
                gr [igraph]->SetLineColor(kBlue-10+igraph%8) ;
                cout << igraph%32 << " " ;
                gr [igraph]->SetLineWidth(1);
                gr [igraph]->SetMarkerColor(kBlue-10+igraph%8);
                gr [igraph]->SetMarkerStyle(8);
                gr [igraph]->SetMarkerSize(1);
                mg->Add(gr [igraph]);
        }
        // Add mean
        TGraph * mean = new TGraph(NbOfPoints, temperature, mean_noise) ;
        mean->SetLineColor(kRed) ;
        mean->SetLineWidth(5);
        mean->SetMarkerColor(kRed);
        mg->Add(mean);

        // Add simulation
        Double_t sim_temp[] = {20., 40., 60., 80., 100.};
        Double_t sim_noise[] ={0.196, 0.202, 0.208, 0.214, 0.221};
        TGraph * sim = new TGraph(5, sim_temp, sim_noise) ;
        sim->SetLineColor(kGreen) ;
        sim->SetLineWidth(3);
        sim->SetMarkerColor(kGreen);
        mg->Add(sim);
        
        mg->Draw("ALP");
        mg->GetXaxis()->SetTitle("Temperature [#circC]");
        mg->GetXaxis()->SetRangeUser(0.,110.);
        mg->GetYaxis()->SetTitle("Noise #sigma [fC]");
//      mg->GetYaxis()->SetRangeUser(2.4,2.5);
        mg->GetYaxis()->SetTitleOffset(1.1);
        mg->Draw("ALP"); // to refresh
//      
        TLegend * leg = new TLegend(0.7,0.75,0.9,0.9) ;
        leg->AddEntry(gr [7],"All channels","lp");      //p: marker only --  l: line
        leg->AddEntry(mean,"Mean","lp");
        leg->AddEntry(sim,"Simulation","lp");
        leg->SetFillStyle(0);
        leg->Draw();
        
        
        //----------------------------------------------------------------
        //                       BUILD GRAPH FOR TEMPERATURE
        //----------------------------------------------------------------      
        TCanvas * c2 = new TCanvas("c2","MICROROC offset dac");
        TMultiGraph * mg2 = new TMultiGraph("mg2","");
        TGraph * gr[NbOfChannels];
        for(Int_t igraph=0;igraph<NbOfChannels;igraph++)
        {       gr [igraph] = new TGraph(NbOfPoints, temperature, data_inflex[igraph] ); // TGraph (n, x, y)
                gr [igraph]->SetLineColor(kBlue-10+igraph%8) ;
                cout << igraph%32 << " " ;
                gr [igraph]->SetLineWidth(1);
                gr [igraph]->SetMarkerColor(kBlue-10+igraph%8);
                gr [igraph]->SetMarkerStyle(8);
                gr [igraph]->SetMarkerSize(1);
                mg2->Add(gr [igraph]);
        }
        // Add mean
        TGraph * mean = new TGraph(NbOfPoints, temperature, mean_inflex) ;
        mean->SetLineColor(kRed) ;
        mean->SetLineWidth(5);
        mean->SetMarkerColor(kRed);
        mg2->Add(mean);

        
        mg2->Draw("ALP");
        mg2->GetXaxis()->SetTitle("Temperature [#circC]");
        mg2->GetXaxis()->SetRangeUser(0.,110.);
        mg2->GetYaxis()->SetTitle("Pedestal [DACu]");
        mg2->GetYaxis()->SetRangeUser(200,217.5);
        mg2->GetYaxis()->SetTitleOffset(1.1);
        mg2->Draw("ALP"); // to refresh
//      
        TLegend * leg2 = new TLegend(0.7,0.8,0.9,0.9) ;
        leg2->AddEntry(gr [7],"All channels","lp");     //p: marker only --  l: line
        leg2->AddEntry(mean,"Mean","lp");
        leg2->SetFillStyle(0);
        leg2->Draw();
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
//      
//      
        
//      Double_t data[2][NbOfGraph][NbOfPoints] ; // first dim is for x / y


//      for(Int_t i=0;i<NbOfPoints;i++)
//      {       data[0][0][i] = temperature_measured[i] ; // x
//              data[1][0][i] = voltage_measured[i] ;     // y
//              data[0][1][i] = temperature_simulated[i] ;
//              data[1][1][i] = voltage_simulated[i] ;
//              cout << data[0][0][i] << " " ;
//      }
//      cout << endl ;
//      
//      
        
 ////   Pedestal: Scurve and derive
//      TFile * f_ped = new TFile(name_file_ped) ;
//      TString * h_ped_name_h = new TString() ;
//      h_ped_name_h->Form("hscurve_%d", ch_num) ;
//      TString * h_ped_name_d = new TString() ;
//      h_ped_name_d->Form("deriv_%d", ch_num) ;
//      TH1I * h_ped_h = (TH1I*) f_ped->Get(h_ped_name_h->Data()) ;
//      TH1I * h_ped_d = (TH1I*) f_ped->Get(h_ped_name_d->Data()) ;
//      h_ped_h->GetXaxis()->SetRangeUser(100,500);
//      h_ped_d->GetXaxis()->SetRangeUser(100,500);
//      
//      TCanvas * c1 = new TCanvas("c1","Pedestals",0, 0, 800, 300);
//      c1->Divide(2,1);
//      c1->cd(1);
//      h_ped_h->Draw();
//      c1->cd(2);
//      h_ped_d->Draw();
//      
//      //      Data: Scurve and derive
//      TFile * f_dat = new TFile(name_file_data) ;
//      TString * h_dat_name_h = new TString() ;
//      h_dat_name_h->Form("hscurve_%d", ch_num) ;
//      TString * h_dat_name_d = new TString() ;
//      h_dat_name_d->Form("deriv_%d", ch_num) ;
//      TH1I * h_dat_h = (TH1I*) f_dat->Get(h_dat_name_h->Data()) ;
//      TH1I * h_dat_d = (TH1I*) f_dat->Get(h_dat_name_d->Data()) ;
//      h_dat_h->GetXaxis()->SetRangeUser(350,750);
//      h_dat_d->GetXaxis()->SetRangeUser(350,750);
//      
//      TCanvas * c2 = new TCanvas("c2","Data",0, 350, 800, 300);
//      c2->Divide(2,1);
//      c2->cd(1);
//      h_dat_h->Draw();
//      c2->cd(2);
//      h_dat_d->Draw();
        
        return 0 ;

}

---