Geant4 Cross Reference |
1 { 1 2 3 TFile *input_file_1 = new TFile("test35a.roo 4 TFile *input_file_2 = new TFile("test35b.roo 5 6 TCanvas *c1 = new TCanvas("c1", "test35", 20 7 c1->SetGridx(); 8 c1->SetGridy(); 9 c1->SetLogx(); 10 c1->SetLogy(); 11 12 // histogram for energy spectra 13 int n = 41; 14 float bin[41]; 15 16 for (int i = 0; i < n; i++) { 17 bin[i] =pow(10,(-2+0.1*i)); 18 } 19 // 20 TH1 *h_1 = new TH1D("unbiased","Source spect 21 TH1 *h_2 = new TH1D("biased","Source spectru 22 input_file_1->cd(); 23 input_file_1->ls(); 24 // get the tuple t1 25 double energy, weight; 26 TTree *t1 = (TTree *) input_file_1->Get("MyT 27 t1->SetBranchAddress("Energy", &energy); 28 t1->SetBranchAddress("Weight", &weight); 29 cout <<t1->GetEntries() << endl; 30 for (int i = 0; i < t1->GetEntries(); i++) { 31 t1.GetEntry(i); 32 // cout << energy << " " << weight << e 33 h_1->Fill(energy,weight); 34 } 35 input_file_2->cd(); 36 TTree *t2 = (TTree *) input_file_2->Get("MyT 37 t2->SetBranchAddress("Energy", &energy); 38 t2->SetBranchAddress("Weight", &weight); 39 cout <<t2->GetEntries() << endl; 40 for (int i = 0; i < t1->GetEntries(); i++) { 41 t2.GetEntry(i); 42 h_2->Fill(energy,weight); 43 } 44 // h_2->SetFillColor(kRed); 45 h_2->SetLineStyle(kDashed); 46 h_2->SetLineColor(kBlue); 47 h_2->Draw(); 48 h_1->Draw("same") ; 49 c1->Update(); 50 c1->Print("./test35.png"); 51 52 input_file_1->Close(); 53 input_file_2->Close(); 54 } 55