ROOT学习——显示晶须定义的蜡烛图示例(candleplotwhiskers.C)

本文主要是介绍ROOT学习——显示晶须定义的蜡烛图示例(candleplotwhiskers.C)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

蜡烛图又称日本线、K线、阴阳线、棒线等，目前常用的说法百是“K线”以下统称K线。K线源于日本德川幕府时代（1603～1867年）的米市交易，用来计算米价每天的涨跌，后来把它引入股票市场价格走势的分析中，目前已成为股票、外汇技术分析中的重要方法。是技术分析的一种，最早为日本人于十九世纪所度创，被当时日本米市的商人用来记录米市的行知情与价格波动，包括开市价、收市价、最高价及最低价，阳烛代表当日升市，阴烛代表跌市。这种图表分析法在当时的中国以至整个东南亚地区均尤为流行。由于用这种方法绘制出来的图表形状颇似一根根蜡烛道，加上这些蜡烛有黑白之分，因而也叫阴阳线图表。通过K线图，人们能够把每日或某一周期的市况表现完全记录下来，股价经过一段时间的盘档后，在图上即形成一种特殊区域或形版态，不同的形态显示出不同意义。可以从这些形态的变化中摸索出一些有规律的东西出来。权K线图形态可分为反转形态、整理形态及缺口和趋向线等。后K线图因其细腻独到的标画方式而被引入到股市及期货市场。

Q1 (-25%): -0.675525 Median: 0.00168511 Q3 (+25%): 0.676189FCN=11.7941 FROM MIGRAD    STATUS=CONVERGED     138 CALLS         139 TOTALEDM=1.08103e-12    STRATEGY= 1      ERROR MATRIX ACCURATE EXT PARAMETER                                   STEP         FIRST   NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE 1  Constant     3.72831e+03   1.82107e+02   4.21637e-02   1.36245e-082  Mean        -1.10774e-01   7.20921e-02   8.42717e-06  -1.49401e-043  Sigma        9.59469e-01   2.51499e-02   2.82317e-06   5.37129e-04FCN=8.04689 FROM MIGRAD    STATUS=CONVERGED     139 CALLS         140 TOTALEDM=6.08152e-08    STRATEGY= 1      ERROR MATRIX ACCURATE EXT PARAMETER                                   STEP         FIRST   NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE 1  Constant     4.07186e+03   2.44580e+02   3.87026e-02   3.94837e-072  Mean        -2.92179e-02   8.71835e-02   7.07837e-06  -2.13598e-023  Sigma        1.00971e+00   2.88317e-02   2.27510e-06  -1.25158e-01FCN=6.93238 FROM MIGRAD    STATUS=CONVERGED     106 CALLS         107 TOTALEDM=6.56099e-08    STRATEGY= 1      ERROR MATRIX ACCURATE EXT PARAMETER                                   STEP         FIRST   NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE 1  Constant     3.97002e+03   2.57596e+01   2.40433e-02  -9.00146e-062  Mean        -1.89022e-04   1.18284e-02   1.62701e-05   7.29142e-033  Sigma        1.02465e+00   3.35469e-02   1.82218e-05   7.40149e-03

下面给出具体代码步骤：

创建一个画板，将其纵向分为两块：

auto c1 = new TCanvas("c1","Candle Presets",700,800);
c1->Divide(1,2);

定义一个随机数，新建二维直方图和一维直方图：

auto rng = new TRandom();
auto h1 = new TH2I("h1","Gaus",100,-5,5,1,0,1);
auto h2 = new TH1I("h2","Gaus",100,-5,5);

设置直方图x、y轴标题：

h1->GetXaxis()->SetTitle("Standard deviation #sigma");
h2->GetXaxis()->SetTitle("Standard deviation #sigma");
h2->GetYaxis()->SetTitle("dN/d#sigma");

生成随机数填充直方图：

float myRand;
for (int i = 0; i < 100000; i++) {myRand = rng->Gaus(0,1);h1->Fill(myRand,0);h2->Fill(myRand);
}

定义分位数数组以及概率密度参数数组：

Double_t *q = new Double_t[3];
Double_t *p = new Double_t[3];
q[0] = 0.; q[1] = 0.; q[2] = 0.;
p[0] = 0.25; p[1] = 0.5; p[2] = 0.75;

获取h2分位数值：

h2->GetQuantiles(3,q,p);
//其中GetQuantitiles函数：
//输入参数
//    -此一维直方图（TH1F，D等）。 也可以是TProfile
//    -nprobSum数组q的最大大小和数组probSum的大小（如果给定）
//    -将计算分位数的位置的probSum数组。
//      如果probSum为null，则probSum将在内部进行计算，其大小=箱数+ h中的1。 它将对应于在直方图的最低边缘（分位数= 0）和bin的所有较高边缘计算的分位数。
//      如果probSum不为null，则假定至少包含nprobSum值。
//输出
//    -返回值nq（<= nprobSum）以及计算的分位数
//    -用nq分位数填充的数组q

输出获得的分位数：

cout << "Q1 (-25%): " << q[0] << " Median: " << q[1] << " Q3 (+25%): " << q[2] << endl;

定义3个一维直方图：

double iqr = q[2]-q[0];
auto mygaus_1_middle = new TF1("mygaus_1_middle","gaus",q[0],q[2]);
auto mygaus_1_left   = new TF1("mygaus_1_left","gaus",q[0]-1.5*iqr,q[0]);
mygaus_1_left->SetLineColor(kGreen);
auto mygaus_1_right  = new TF1("mygaus_1_right","gaus",q[2],q[2]+1.5*iqr);
mygaus_1_right->SetLineColor(kGreen);

绘制h1：

c1->cd(1); //定位到画板上半部分
h1->SetLineWidth(3); //设置线宽为3 
h1->SetFillStyle(0); //设置填充类型
h1->Draw("candley2 scat"); //用candley2类型绘制和散点图类型绘制

绘制h2：

c1->cd(2);
h2->Draw("");

对h2区间范围内进行拟合：

h2->Fit("mygaus_1_left","R"); //会自动输出拟合参数

不改变坐标轴对拟合函数进行绘制：

mygaus_1_left->Draw("same");

绘制图形线(q[0]-1.5*iqr的竖线)：

auto l3 = new TLine(q[0]-1.5*iqr,0,q[0]-1.5*iqr,mygaus_1_left->Eval(q[0]-1.5*iqr)); //从(x1,y1)到(x2,y2)
l3->SetLineColor(kGreen); //设置线的颜色 
l3->SetLineWidth(2); //设置线的宽度
l3->Draw(""); //绘制

绘制图形线(q[0]的竖线)：

auto l1 = new TLine(q[0]        ,0,q[0]        ,mygaus_1_left->Eval(q[0]));
l1->SetLineWidth(2);      
l1->SetLineColor(kGreen); 
l1->Draw("");

同理，拟合绘制右区间：

h2->Fit("mygaus_1_right","R","");
mygaus_1_right->Draw("same");
auto l4 = new TLine(q[2]+1.5*iqr,0,q[2]+1.5*iqr,mygaus_1_left->Eval(q[2]+1.5*iqr));
l4->SetLineColor(kGreen); 
l4->SetLineWidth(2);      
l4->Draw("");
auto l5 = new TLine(q[2]        ,0,q[2]        ,mygaus_1_right->Eval(q[2]));
l5->SetLineWidth(2);      
l5->SetLineColor(kGreen); 
l5->Draw("");

拟合绘制中间区间：

h2->Fit("mygaus_1_middle","R");
mygaus_1_middle->Draw("same");

向图像中添加文本(原则上，也可以通过h2-> Integral()计算这些值)：

auto t = new TText(); t->SetTextFont(42); //设置本文格式
t->DrawText(0,mygaus_1_middle->Eval(0)/2,"50%"); //x坐标，y坐标，本文
t->DrawText(-1.5,mygaus_1_middle->Eval(-1.5)/2,"24.65%");
t->DrawText(+1,mygaus_1_middle->Eval(+1.5)/2,"24.65%");
t->DrawText(q[0]-1.5*iqr,1000,Form("%.3f",q[0]-1.5*iqr))->SetTextAngle(90); //设置文本角度
t->DrawText(q[2]+1.5*iqr,1000,Form("%.3f",q[2]+1.5*iqr))->SetTextAngle(90);
t->DrawText(q[0],1000,Form("%.3f",q[0]))->SetTextAngle(90);
t->DrawText(q[2],1000,Form("%.3f",q[2]))->SetTextAngle(90);

代码地址：https://github.com/root-project/root/blob/master/tutorials/hist/candleplotwhiskers.C

这篇关于ROOT学习——显示晶须定义的蜡烛图示例(candleplotwhiskers.C)的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！