| 1 | sensitivity(){
|
|---|
| 2 | TGraph dwarf;
|
|---|
| 3 | dwarf.SetName("DWARF");
|
|---|
| 4 | dwarf.SetTitle("Differential Flux Sensitivity");
|
|---|
| 5 | TLegend leg(0.75, 0.75, 0.95, 0.95); //x1, y1, x2, y2 in %
|
|---|
| 6 | //dwarf.SetPoint(dwarf.GetN()+1, x1, y1);
|
|---|
| 7 | //dwarf.SetPoint(dwarf.GetN()+1, x2, y2);
|
|---|
| 8 | dwarf.SetPoint(dwarf.GetN(), 0.001, 6E-12);
|
|---|
| 9 | dwarf.SetPoint(dwarf.GetN(), 0.25, 6E-12);
|
|---|
| 10 | dwarf.SetPoint(dwarf.GetN(), 0.31, 4.38E-12);
|
|---|
| 11 | dwarf.SetPoint(dwarf.GetN(), 0.39, 2.1E-12);
|
|---|
| 12 | dwarf.SetPoint(dwarf.GetN(), 0.45, 1.54E-12);
|
|---|
| 13 | dwarf.SetPoint(dwarf.GetN(), 0.5, 1.4E-12);
|
|---|
| 14 | dwarf.SetPoint(dwarf.GetN(), 1, 7.5E-13);
|
|---|
| 15 | dwarf.SetPoint(dwarf.GetN(), 10, 2.25E-13);
|
|---|
| 16 | dwarf.SetPoint(dwarf.GetN(), 50, 3.13E-13);
|
|---|
| 17 | dwarf.SetLineColor(kBlue);
|
|---|
| 18 | dwarf.SetLineWidth(5);
|
|---|
| 19 |
|
|---|
| 20 | TGraph *gclone = dwarf.DrawClone("AC");
|
|---|
| 21 |
|
|---|
| 22 | TH1 *h = gclone->GetHistogram();
|
|---|
| 23 |
|
|---|
| 24 | TAxis *axex = h->GetXaxis();
|
|---|
| 25 | TAxis *axey = h->GetYaxis();
|
|---|
| 26 | //axex->SetLabelOffset(x);
|
|---|
| 27 | //axey->SetLabelOffset(y);
|
|---|
| 28 | //axex->SetLabelSize(x);
|
|---|
| 29 | //axey->SetLabelSize(y);
|
|---|
| 30 | axex->SetTitleOffset(1.30);
|
|---|
| 31 | axey->SetTitleOffset(1.2);
|
|---|
| 32 | //axex->SetTitleSize(x);
|
|---|
| 33 | //axey->SetTitleSize(y);
|
|---|
| 34 |
|
|---|
| 35 | axex->SetMoreLogLabels();
|
|---|
| 36 | //axey->SetMoreLogLabels();
|
|---|
| 37 |
|
|---|
| 38 | axex->SetTitle("E_0 [TeV]");
|
|---|
| 39 | axey->SetTitle("F(E<E_0) [cm^{-2}s^{-1}]");
|
|---|
| 40 |
|
|---|
| 41 | axex->SetRangeUser(0.005, 100);
|
|---|
| 42 | axey->SetRangeUser(1E-13, 2E-10);
|
|---|
| 43 | gclone->RemovePoint(0);
|
|---|
| 44 | leg.AddEntry("DWARF", "DWARF", "l");
|
|---|
| 45 |
|
|---|
| 46 |
|
|---|
| 47 | TGraph hegras;
|
|---|
| 48 | hegras.SetName("HEGRAS");
|
|---|
| 49 | hegras.SetPoint(hegras.GetN(), 0.5, 1.7E-12);
|
|---|
| 50 | hegras.SetPoint(hegras.GetN(), 1, 1E-12);
|
|---|
| 51 | hegras.SetPoint(hegras.GetN(), 10, 2.5E-13);
|
|---|
| 52 | hegras.SetPoint(hegras.GetN(), 50, 2.5E-13);
|
|---|
| 53 | hegras.SetLineColor(kGreen);
|
|---|
| 54 | hegras.SetLineWidth(5);
|
|---|
| 55 | hegras.DrawClone("C");
|
|---|
| 56 | leg.AddEntry("HEGRAS", "HEGRA CT System", "l");
|
|---|
| 57 |
|
|---|
| 58 | //TGraph hegra2;
|
|---|
| 59 | //hegra.SetName("HEGRA2");
|
|---|
| 60 | //hegra2.SetPoint(hegra2.GetN(), 0.4, 1.7E-12); //1.2E-12
|
|---|
| 61 | //hegra2.SetLineColor(kBlue);
|
|---|
| 62 | //hegra2.SetMarkerStyle(kFullDotMedium);
|
|---|
| 63 | //hegra2.DrawClone("P");
|
|---|
| 64 | //leg.AddEntry("HEGRA2", "HEGRA CT System 2", "p");
|
|---|
| 65 |
|
|---|
| 66 | TF1* pow2=new TF1("pow2","1.7e-12/pow(0.4,-0.73)*pow(x,-0.73)",0.4,7);
|
|---|
| 67 | pow2->SetLineColor(kGreen);
|
|---|
| 68 | pow2->SetLineStyle(kDashed);
|
|---|
| 69 | pow2->Draw("same");
|
|---|
| 70 | leg.AddEntry("pow2", "HEGRA CT System 2", "l");
|
|---|
| 71 |
|
|---|
| 72 | //TGraph hegra;
|
|---|
| 73 | //hegra.SetName("HEGRA");
|
|---|
| 74 | //hegra.SetPoint(hegra.GetN(), 0.75, 7.1E-12); //1E-11
|
|---|
| 75 | //hegra.SetLineColor(kBlue);
|
|---|
| 76 | //hegra.SetMarkerStyle(kFullDotMedium);
|
|---|
| 77 | //hegra.SetMarkerSize(5);
|
|---|
| 78 | //hegra.DrawClone("P");
|
|---|
| 79 | //leg.AddEntry("HEGRA", "1 HEGRA CT", "p");
|
|---|
| 80 |
|
|---|
| 81 | TF1* pow3=new TF1("pow3","7.1e-12/pow(0.75,-0.73)*pow(x,-0.73)",0.75,7);
|
|---|
| 82 | pow3->SetLineColor(kBlue);
|
|---|
| 83 | pow3->SetLineStyle(kDashed);
|
|---|
| 84 | pow3->Draw("same");
|
|---|
| 85 | leg.AddEntry("pow3", "1 single HEGRA CT", "l");
|
|---|
| 86 |
|
|---|
| 87 | TGraph whipple;
|
|---|
| 88 | whipple.SetName("Whipple");
|
|---|
| 89 | whipple.SetPoint(whipple.GetN(), 0.25, 1E-11);
|
|---|
| 90 | whipple.SetPoint(whipple.GetN(), 1, 2.7E-12);
|
|---|
| 91 | whipple.SetPoint(whipple.GetN(), 3, 1E-12);
|
|---|
| 92 | whipple.SetPoint(whipple.GetN(), 10, 1E-12);
|
|---|
| 93 | whipple.SetLineColor(kRed);
|
|---|
| 94 | whipple.SetLineWidth(5);
|
|---|
| 95 | whipple.DrawClone("C");
|
|---|
| 96 | leg.AddEntry("Whipple", "Whipple", "l");
|
|---|
| 97 |
|
|---|
| 98 | //TGraph whipple2;
|
|---|
| 99 | //whipple2.SetName("Whipple2");
|
|---|
| 100 | //whipple2.SetPoint(whipple2.GetN(), 0.3, 9.2E-12); //1.3E-11
|
|---|
| 101 | //whipple2.SetLineColor(kBlue);
|
|---|
| 102 | //whipple2.SetMarkerStyle(kFullDotMedium);
|
|---|
| 103 | //whipple2.DrawClone("P");
|
|---|
| 104 | //leg.AddEntry("Whipple2", "Whipple2", "p");
|
|---|
| 105 |
|
|---|
| 106 | TF1* pow1=new TF1("pow1","1.3e-11/pow(0.3,-0.95)*pow(x,-0.95)",0.3,3.5);
|
|---|
| 107 | pow1->SetLineColor(kRed);
|
|---|
| 108 | pow1->SetLineStyle(kDashed);
|
|---|
| 109 | pow1->Draw("same");
|
|---|
| 110 | leg.AddEntry("pow1", "Whipple 2", "l");
|
|---|
| 111 |
|
|---|
| 112 | TGraph magic_loi;
|
|---|
| 113 | magic_loi.SetName("magic_loi");
|
|---|
| 114 | magic_loi.SetPoint(magic_loi.GetN(), 0.01, 1E-10);
|
|---|
| 115 | magic_loi.SetPoint(magic_loi.GetN(), 0.045, 1E-11);
|
|---|
| 116 | magic_loi.SetPoint(magic_loi.GetN(), 0.095, 4.5E-12);
|
|---|
| 117 | magic_loi.SetPoint(magic_loi.GetN(), 0.55, 1E-12);
|
|---|
| 118 | magic_loi.SetPoint(magic_loi.GetN(), 1, 6E-13);
|
|---|
| 119 | magic_loi.SetLineColor(kCyan);
|
|---|
| 120 | magic_loi.SetLineWidth(5);
|
|---|
| 121 | magic_loi.SetLineStyle(kDashed);
|
|---|
| 122 | magic_loi.DrawClone("C");
|
|---|
| 123 | leg.AddEntry("magic_loi", "MAGIC LoI", "l");
|
|---|
| 124 |
|
|---|
| 125 | TGraph magic;
|
|---|
| 126 | magic.SetName("magic");
|
|---|
| 127 | magic.SetPoint(magic.GetN(), 0.036, 1.2E-10);
|
|---|
| 128 | magic.SetPoint(magic.GetN(), 0.045, 8E-11);
|
|---|
| 129 | magic.SetPoint(magic.GetN(), 0.063, 5E-11);
|
|---|
| 130 | magic.SetPoint(magic.GetN(), 0.095, 2.3E-11);
|
|---|
| 131 | magic.SetPoint(magic.GetN(), 0.12, 1.3E-11);
|
|---|
| 132 | magic.SetPoint(magic.GetN(), 0.15, 7E-12);
|
|---|
| 133 | magic.SetPoint(magic.GetN(), 0.18, 5E-12);
|
|---|
| 134 | magic.SetPoint(magic.GetN(), 0.25, 3E-12);
|
|---|
| 135 | magic.SetPoint(magic.GetN(), 0.31, 2.5E-12);
|
|---|
| 136 | magic.SetPoint(magic.GetN(), 0.45, 1.1E-12);
|
|---|
| 137 | magic.SetPoint(magic.GetN(), 0.5, 1.05E-12);
|
|---|
| 138 | magic.SetPoint(magic.GetN(), 0.6, 9E-13);
|
|---|
| 139 | magic.SetPoint(magic.GetN(), 0.8, 6E-13);
|
|---|
| 140 | magic.SetPoint(magic.GetN(), 0.9, 6E-13);
|
|---|
| 141 | magic.SetPoint(magic.GetN(), 1, 4E-13);
|
|---|
| 142 | magic.SetLineColor(kCyan);
|
|---|
| 143 | magic.SetLineWidth(5);
|
|---|
| 144 | magic.DrawClone("C");
|
|---|
| 145 | leg.AddEntry("magic", "MAGIC", "l");
|
|---|
| 146 |
|
|---|
| 147 | leg.DrawClone();
|
|---|
| 148 | gPad->SetGridx();
|
|---|
| 149 | gPad->SetGridx();
|
|---|
| 150 | gPad->SetLogx();
|
|---|
| 151 | gPad->SetLogy();
|
|---|
| 152 | gPad->SetFillColor(kWhite);
|
|---|
| 153 | gPad->SetBottomMargin(0.13);
|
|---|
| 154 | gPad->SetLeftMargin(0.11);
|
|---|
| 155 |
|
|---|
| 156 | } |
|---|
