RandomScripts/nucHW1.C at master · rr521111/RandomScripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
void nucHW1(){

    Double_t A = 0;
    Double_t Z = 0;
    Double_t BE = 0;
    vector<TString> types_line;
    vector<vector<TString>> type_table;
    //TH3D *h3 = new TH3D("A vs Z vs BE 3d", "A vs Z vs BE 3d", 270, 0, 271, 110, 0, 111, 1000, 1, -1);
    TGraph2DErrors *t3 = new TGraph2DErrors(2438);
    TH2D *h2 = new TH2D("A vs Z vs BE", "A vs Z vs BE", 270, 0, 271, 110, 0, 111);
    TH2D *h2fit = new TH2D("Fit Result A vs Z vs BE", "Fit Result A vs Z vs BE", 270, 0, 271, 110, 0, 111);
    TH2D *h2res = new TH2D("Fit Residual A vs Z vs BE", "Fit Residual A vs Z vs BE", 270, 0, 271, 110, 0, 111);
    TH2D *h2rescent = new TH2D("Fit Residual  A vs Z vs BE Percent Diff", "Fit Residual A vs Z vs BE Percent Diff", 270, 0, 271, 110, 0, 111);

    ifstream infile;
    infile.open("./AME2012.txt");

    string ins;

    bool isemptyline = false;
    Int_t point = 0;

    while (!infile.eof()) {
        //splits the file at line breaks
        getline(infile, ins, '\n');
        stringstream line(ins);
        types_line.clear();
        //splits the line at commas
        while (line.good()) {
            string substr;
            getline(line, substr, ',');
            if(substr == ""){
                isemptyline = true;
                break;
            }
            types_line.push_back(substr);
        }

        if(isemptyline == true){
            break;
        }

        A = stod(types_line[2].Data());
        Z = stod(types_line[1].Data());
        BE = stod(types_line[6].Data())*A;

        //h3->Fill(A, Z, BE);
        t3->SetPoint(point, A, Z, BE);
        t3->SetPointError(point, 0, 0, stod(types_line[7].Data())*A);
        h2->Fill(A, Z, BE);
        point++;
        //cout << types_line[0] << " " << types_line[1] << " " << types_line[6] << "+-" << types_line[7] << " " << endl;

    }

    infile.close();

    TCanvas *c1 = new TCanvas();
    c1->Divide(2,2);
    c1->cd(1);
    TF2* f2 = new TF2("f2", "[0]*x + [1]*(x^(2/3)) + [2]*y*(y-1)/(x^(1/3)) + [3]*((y-x/2)^2)/x + [4]*(sqrt(x)^(-1))");
    t3->Fit("f2");
    gStyle->SetOptStat(0);
    //gPad->SetLogz(1);
    h2->Draw("colz");

    infile.open("./AME2012.txt");

    isemptyline = false;

    while (!infile.eof()) {
        //splits the file at line breaks
        getline(infile, ins, '\n');
        stringstream line(ins);
        types_line.clear();
        //splits the line at commas
        while (line.good()) {
            string substr;
            getline(line, substr, ',');
            if(substr == ""){
                isemptyline = true;
                break;
            }
            types_line.push_back(substr);

        }

        if(isemptyline == true){
            break;
        }

        A = stod(types_line[2].Data());
        Z = stod(types_line[1].Data());
        BE = stod(types_line[6].Data())*A;

        h2fit->Fill(A, Z, f2->Eval(A, Z));
        h2res->Fill(A, Z, BE - f2->Eval(A, Z));
        h2rescent->Fill(A, Z, abs((BE - f2->Eval(A, Z))/(BE + f2->Eval(A, Z))));
    }

    infile.close();

    c1->cd(2);
    //gPad->SetLogz(1);
    h2fit->Draw("colz");
    c1->cd(3);
    //gPad->SetLogz(1);
    h2res->Draw("colz");
    c1->cd(4);
    //gPad->SetLogz(1);
    h2rescent->Draw("colz");

    return;
}