{
    const Int_t NThrows = 100000;
    //Histogram for real data
    TH1D *eventbin1 = new TH1D("eventbin1","",100,0,50);
    TH1D *eventbin2 = new TH1D("eventbin2","",100,0,50);
    TH1D *eventbin3 = new TH1D("eventbin3","",100,0,50);
    TH1D *his_event_average = new TH1D("event_average", "", 100, 0, 50);
    Double_t event1;
    Double_t event2;
    Double_t event3;
    Double_t event_aver;

    Double_t event1mean = 12.9522;
    Double_t event1rms = TMath::Sqrt(event1mean);
    Double_t event2mean = 14.0023;
    Double_t event2rms = TMath::Sqrt(event2mean);
    Double_t event3mean = 15.9277;
    Double_t event3rms = TMath::Sqrt(event3mean);

    for (Int_t ithrow = 0; ithrow < NThrows; ++ ithrow){
         event1 = gRandom->Gaus(event1mean, event1rms);
         eventbin1 -> Fill(event1);
         event2 = gRandom->Gaus(event2mean, event2rms);
         eventbin2 -> Fill(event2);
         event3 = gRandom->Gaus(event3mean, event3rms);
         eventbin3 -> Fill(event3);

         event_aver = sum(event1, event2, event3);
         event_aver -> Scale(1.0/3.0);
         his_event_average -> Fill(event_aver);
}
    
     new TCanvas;
    gStyle->SetOptStat(1);
    eventbin1->Draw("hist");
    eventbin1 -> SetTitle("Number of event at E = 0.55 GeV");
    eventbin1 ->GetXaxis()->SetTitle("Number of event");
  //  gPad->Print("bin1.pdf");
    cout<<"for event1mean "<<eventbin1->GetMean()<<" rms "<<eventbin1->GetRMS()<<endl;
 
    new TCanvas;
    gStyle->SetOptStat(1);
    eventbin2->Draw("hist");
    eventbin2 -> SetTitle("Number of event at E = 0.65 GeV");
    eventbin2 ->GetXaxis()->SetTitle("Number of event");
  //  gPad->Print("bin1.pdf");
    cout<<"for event2mean "<<eventbin2->GetMean()<<" rms "<<eventbin2->GetRMS()<<endl;

    
    new TCanvas;
    gStyle->SetOptStat(1);
    eventbin3->Draw("hist");
    eventbin3 -> SetTitle("Number of event at E = 0.75 GeV");
    eventbin3 ->GetXaxis()->SetTitle("Number of event");
  //  gPad->Print("bin1.pdf");
    cout<<"for event3mean "<<eventbin3->GetMean()<<" rms "<<eventbin3->GetRMS()<<endl;

    new TCanvas;
    gStyle->SetOptStat(1);
    his_event_average->Draw("hist");
    his_event_average -> SetTitle("Number of event at E = 0.6 GeV");
    his_event_average ->GetXaxis()->SetTitle("Number of event");
  //  gPad->Print("bin1.pdf");
    cout<<"for event_aver_mean "<<his_event_average->GetMean()<<" rms "<<his_event_average->GetRMS()<<endl;

}