Geant4_10
plottest35.C
Go to the documentation of this file.
1 {
2 
3  TFile *input_file_1 = new TFile("test35a.root");
4  TFile *input_file_2 = new TFile("test35b.root");
5 
6  TCanvas *c1 = new TCanvas("c1", "test35", 200, 10, 700, 500);
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 spectrum",40,bin);
21  TH1 *h_2 = new TH1D("biased","Source spectrum",40,bin);
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("MyTuple");
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 << endl;
33  h_1->Fill(energy,weight);
34  }
35  input_file_2->cd();
36  TTree *t2 = (TTree *) input_file_2->Get("MyTuple");
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 }
TTree * t1
Definition: plottest35.C:26
TH1 * h_1
Definition: plottest35.C:20
float bin[41]
Definition: plottest35.C:14
TCanvas * c1
Definition: plotHisto.C:7
double weight
Definition: plottest35.C:25
TH1 * h_2
Definition: plottest35.C:21
Char_t n[5]
double energy
Definition: plottest35.C:25
TTree * t2
Definition: plottest35.C:36
tuple input_file_1
Definition: plottest35.py:8
TFile * input_file_2
Definition: plottest35.C:4