3 // read ALICE INEL data from file
7 cout << "================================" << endl;
9 cout << "read ALICE INEL data from file" <<endl;
10 cout << "Number of bins: " << binsInelAlice <<endl;
11 cout << "Filename: " << filenameInelAlice <<endl;
13 cout << "================================" << endl;
16 ifstream fileInelAlice;
17 fileInelAlice.open(filenameInelAlice);
20 while(!fileInelAlice.eof()) {
21 if(i == binsInelAlice) break;
22 fileInelAlice >> centerPtInelAlice[i] >> inelAlice[i] >> statInelAlice[i] >> systInelAlice[i];
23 //systInelAlice[i] = 0.15*inelAlice[i]+1e-7; // sys error has to be provided, currently 10%
25 // the width of the pt bins, currently hardwired....
26 widthPtInelAlice[i] = 0.05;
27 if (centerPtInelAlice[i] > 1) widthPtInelAlice[i] = 0.1;
28 if (centerPtInelAlice[i] > 2) widthPtInelAlice[i] = 0.2;
29 if (centerPtInelAlice[i] > 4) widthPtInelAlice[i] = 0.5;
30 if (centerPtInelAlice[i] > 7) widthPtInelAlice[i] = 1.0;
32 ptInelAlice[i] = centerPtInelAlice[i];
33 errPtInelAlice[i] = widthPtInelAlice[i] / 2.0;
34 lowPtInelAlice[i] = centerPtInelAlice[i] - errPtInelAlice[i];
35 highPtInelAlice[i] = centerPtInelAlice[i] + errPtInelAlice[i];
37 lowStatInelAlice[i] = inelAlice[i] - statInelAlice[i];
38 highStatInelAlice[i] = inelAlice[i] + statInelAlice[i];
39 relStatInelAlice[i] = statInelAlice[i] / inelAlice[i];
40 lowSystInelAlice[i] = inelAlice[i] - systInelAlice[i];
41 highSystInelAlice[i] = inelAlice[i] + systInelAlice[i];
42 relSystInelAlice[i] = systInelAlice[i] / inelAlice[i];
43 errInelAlice[i] = systInelAlice[i] + statInelAlice[i];
44 lowErrInelAlice[i] = inelAlice[i] - errInelAlice[i];
45 highErrInelAlice[i] = inelAlice[i] + errInelAlice[i];
46 relErrInelAlice[i] = errInelAlice[i] / inelAlice[i];
47 err2InelAlice[i] = sqrt(systInelAlice[i]*systInelAlice[i] + statInelAlice[i]*statInelAlice[i]);
48 lowErr2InelAlice[i] = inelAlice[i] - err2InelAlice[i];
49 highErr2InelAlice[i] = inelAlice[i] + errInelAlice[i];
50 relErr2InelAlice[i] = err2InelAlice[i] / inelAlice[i];
54 ptInel2PiPtAlice[i] = ptInelAlice[i];
55 centerPtInel2PiPtAlice[i] = centerPtInelAlice[i];
56 widthPtInel2PiPtAlice[i] = widthPtInelAlice[i];
57 errPtInel2PiPtAlice[i] = errPtInelAlice[i];
58 lowPtInel2PiPtAlice[i] = lowPtInelAlice[i];
59 highPtInel2PiPtAlice[i] = highPtInelAlice[i];
61 inel2PiPtAlice[i] = inelAlice[i]*centerPtInelAlice[i]*2*M_PI;
62 statInel2PiPtAlice[i] = statInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
63 lowStatInel2PiPtAlice[i] = lowStatInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
64 highStatInel2PiPtAlice[i] = highStatInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
65 relStatInel2PiPtAlice[i] = relStatInelAlice[i];
66 systInel2PiPtAlice[i] = systInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
67 lowSystInel2PiPtAlice[i] = lowSystInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
68 highSystInel2PiPtAlice[i] = highSystInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
69 relSystInel2PiPtAlice[i] = relSystInelAlice[i];
70 errInel2PiPtAlice[i] = errInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
71 lowErrInel2PiPtAlice[i] = lowErrInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
72 highErrInel2PiPtAlice[i] = highErrInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
73 relErrInel2PiPtAlice[i] = relErrInelAlice[i];
74 err2Inel2PiPtAlice[i] = err2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
75 lowErr2Inel2PiPtAlice[i] = lowErr2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
76 highErr2Inel2PiPtAlice[i] = highErr2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
77 relErr2Inel2PiPtAlice[i] = relErr2InelAlice[i];
81 cout << "ptInelAlice[" << i << "] = " << ptInelAlice[i] <<endl;
82 cout << " centerPtInelAlice[" << i << "] = " << centerPtInelAlice[i] <<endl;
83 cout << " widthPtInelAlice[" << i << "] = " << widthPtInelAlice[i] <<endl;
84 cout << " errPtInelAlice[" << i << "] = " << errPtInelAlice[i] <<endl;
85 cout << " lowPtInelAlice[" << i << "] = " << lowPtInelAlice[i] <<endl;
86 cout << " highPtInelAlice[" << i << "] = " << highPtInelAlice[i] <<endl;
87 cout << "inelAlice[" << i << "] = " << inelAlice[i] <<endl;
88 cout << " statInelAlice[" << i << "] = " << statInelAlice[i] <<endl;
89 cout << " lowStatInelAlice[" << i << "] = " << lowStatInelAlice[i] <<endl;
90 cout << " highStatInelAlice[" << i << "] = " << highStatInelAlice[i] <<endl;
91 cout << " relStatInelAlice[" << i << "] = " << relStatInelAlice[i] <<endl;
92 cout << " systInelAlice[" << i << "] = " << systInelAlice[i] <<endl;
93 cout << " lowSystInelAlice[" << i << "] = " << lowSystInelAlice[i] <<endl;
94 cout << " highSystInelAlice[" << i << "] = " << highSystInelAlice[i] <<endl;
95 cout << " relSystInelAlice[" << i << "] = " << relSystInelAlice[i] <<endl;
96 cout << "errInelAlice[" << i << "] = " << errInelAlice[i] <<endl;
97 cout << " lowErrInelAlice[" << i << "] = " << lowErrInelAlice[i] <<endl;
98 cout << " highErrInelAlice[" << i << "] = " << highErrInelAlice[i] <<endl;
99 cout << " relErrInelAlice[" << i << "] = " << relErrInelAlice[i] <<endl;
100 cout << "err2InelAlice[" << i << "] = " << err2InelAlice[i] <<endl;
101 cout << " lowErr2InelAlice[" << i << "] = " << lowErr2InelAlice[i] <<endl;
102 cout << " highErr2InelAlice[" << i << "] = " << highErr2InelAlice[i] <<endl;
103 cout << " relErr2InelAlice[" << i << "] = " << relErr2InelAlice[i] <<endl;
106 ptInelAliceFit[i] = ptInelAlice[i];
107 InelInelAliceFit[i] = InelInelAlice[i]*ptInelAlice[i];
108 errInelInelAliceFit[i] = errInelInelAlice[i];
109 widthPtInelAliceFit[i] = widthPtInelAlice[i];
110 errPtInelAliceFit[i] = 0; //errPtInelAlice[i];
111 lowPtInelAliceFit[i] = lowPtInelAlice[i];
112 highPtInelAliceFit[i] = highPtInelAlice[i];
116 } // while(!fileInelAlice.eof())
117 fileInelAlice.close();
118 //if (fileInelAlice) { delete fileInelAlice; }
121 cout << "================================" << endl;
123 cout << "Finished reading ALICE INEL data" <<endl;
124 cout << "Number of bins read: " << i <<endl;
126 cout << "================================" << endl;