[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / ChargedHadrons / dNdPt / macros / plots / readUa1.C

Int_t readUa1() 
{
// read UA1 data from file
//

cout << endl;
cout << "================================" << endl;
cout << endl;
cout << "read UA1 YIELD data from file" <<endl;
cout << "Number of bins: " << binsUa1 <<endl;
cout << "Filename:       " << filenameUa1 <<endl;
cout << endl;
cout << "================================" << endl;
cout << endl;

ifstream fileCrossUa1;
fileCrossUa1.open(filenameUa1);

Int_t i = 0;
while(!fileCrossUa1.eof()) {
    if(i == binsUa1) break;
    // textfile content: pt:cross_yield:(stat + syst error added linerarly)
    fileCrossUa1 >> centerPtUa1[i] >> crossUa1[i] >> errCrossUa1[i];
    
    ptUa1[i]          = centerPtUa1[i];
    
    // the width of the pt bins, currently hardwired....
    widthPtUa1[i] = 0.1;
    if (centerPtUa1[i] > 4) widthPtUa1[i] = 0.2;
    if (centerPtUa1[i] > 6) widthPtUa1[i] = 1.0;
    
    errPtUa1[i]       = widthPtUa1[i] / 2.0;              
    lowPtUa1[i]       = centerPtUa1[i] - errPtUa1[i];
    highPtUa1[i]      = centerPtUa1[i] + errPtUa1[i];
    
    /*
    lowStatCrossUa1[i]  = crossUa1[i] - statCrossUa1[i];
    highStatCrossUa1[i] = crossUa1[i] + statCrossUa1[i];
    relStatCrossUa1[i]  = statCrossUa1[i] / crossUa1[i];    
    lowSystCrossUa1[i]  = crossUa1[i] - systCrossUa1[i];
    highSystCrossUa1[i] = crossUa1[i] + systCrossUa1[i];
    relSystCrossUa1[i]  = systCrossUa1[i] / crossUa1[i];    
    */    
    lowErrCrossUa1[i]   = crossUa1[i] - errCrossUa1[i];
    highErrCrossUa1[i]  = crossUa1[i] + errCrossUa1[i];
    relErrCrossUa1[i]   = errCrossUa1[i] / crossUa1[i];
    err2CrossUa1[i]     = errCrossUa1[i];
    lowErr2CrossUa1[i]  = lowErrCrossUa1[i];
    highErr2CrossUa1[i] = highErrCrossUa1[i];
    relErr2CrossUa1[i]  = relErrCrossUa1[i];
    
    yieldUa1[i]         = crossUa1[i] * (avgToHadr / sigmaInelUa1);
    /*
    statYieldUa1[i]     = 
    lowStatYieldUa1[i]  = 
    highStatYieldUa1[i] = 
    relStatYieldUa1[i]  = 
    systYieldUa1[i]     = 
    lowSystYieldUa1[i]  = 
    highSystYieldUa1[i] = 
    relSystYieldUa1[i]  = 
    */
    errYieldUa1[i]      = errCrossUa1[i] * (avgToHadr / sigmaInelUa1);
    lowErrYieldUa1[i]   = yieldUa1[i] - errYieldUa1[i];
    highErrYieldUa1[i]  = yieldUa1[i] + errYieldUa1[i];
    relErrYieldUa1[i]   = errYieldUa1[i] / yieldUa1[i];
    err2YieldUa1[i]     = errYieldUa1[i];
    lowErr2YieldUa1[i]  = lowErrYieldUa1[i];
    highErr2YieldUa1[i] = highErrYieldUa1[i];
    relErr2YieldUa1[i]  = relErrYieldUa1[i];
        
    cout << "ptUa1[" << i << "]               = " << ptUa1[i] <<endl;
    cout << "   centerPtUa1[" << i << "]      = " << centerPtUa1[i] <<endl;
    cout << "   widthPtUa1[" << i << "]       = " << widthPtUa1[i] <<endl;
    cout << "   errPtUa1[" << i << "]         = " << errPtUa1[i] <<endl;
    cout << "   lowPtUa1[" << i << "]         = " << lowPtUa1[i] <<endl;
    cout << "   highPtUa1[" << i << "]        = " << highPtUa1[i] <<endl;
    cout << "crossUa1[" << i << "]            = " << crossUa1[i] <<endl;
    /*
    cout << "   statCrossUa1[" << i << "]     = " << statCrossUa1[i] <<endl;
    cout << "   lowStatCrossUa1[" << i << "]  = " << lowStatCrossUa1[i] <<endl;
    cout << "   highStatCrossUa1[" << i << "] = " << highStatCrossUa1[i] <<endl;
    cout << "   relStatCrossUa1[" << i << "]  = " << relStatCrossUa1[i] <<endl;
    cout << "   systCrossUa1[" << i << "]     = " << systCrossUa1[i] <<endl;
    cout << "   lowSystCrossUa1[" << i << "]  = " << lowSystCrossUa1[i] <<endl;
    cout << "   highSystCrossUa1[" << i << "] = " << highSystCrossUa1[i] <<endl;
    cout << "   relSystCrossUa1[" << i << "]  = " << relSystCrossUa1[i] <<endl;
    */
    cout << "errCrossUa1[" << i << "]         = " << errCrossUa1[i] <<endl;
    cout << "   lowErrCrossUa1[" << i << "]   = " << lowErrCrossUa1[i] <<endl;
    cout << "   highErrCrossUa1[" << i << "]  = " << highErrCrossUa1[i] <<endl;
    cout << "   relErrCrossUa1[" << i << "]   = " << relErrCrossUa1[i] <<endl;
    cout << "err2CrossUa1[" << i << "]        = " << err2CrossUa1[i] <<endl;
    cout << "   lowErr2CrossUa1[" << i << "]  = " << lowErr2CrossUa1[i] <<endl;
    cout << "   highErr2CrossUa1[" << i << "] = " << highErr2CrossUa1[i] <<endl;
    cout << "   relErr2CrossUa1[" << i << "]  = " << relErr2CrossUa1[i] <<endl;

    cout << "yieldUa1[" << i << "]            = " << yieldUa1[i] <<endl;
    /*
    cout << "   statYieldUa1[" << i << "]     = " << statYieldUa1[i] <<endl;
    cout << "   lowStatYieldUa1[" << i << "]  = " << lowStatYieldUa1[i] <<endl;
    cout << "   highStatYieldUa1[" << i << "] = " << highStatYieldUa1[i] <<endl;
    cout << "   relStatYieldUa1[" << i << "]  = " << relStatYieldUa1[i] <<endl;
    cout << "   systYieldUa1[" << i << "]     = " << systYieldUa1[i] <<endl;
    cout << "   lowSystYieldUa1[" << i << "]  = " << lowSystYieldUa1[i] <<endl;
    cout << "   highSystYieldUa1[" << i << "] = " << highSystYieldUa1[i] <<endl;
    cout << "   relSystYieldUa1[" << i << "]  = " << relSystYieldUa1[i] <<endl;
    */
    cout << "errYieldUa1[" << i << "]         = " << errYieldUa1[i] <<endl;
    cout << "   lowErrYieldUa1[" << i << "]   = " << lowErrYieldUa1[i] <<endl;
    cout << "   highErrYieldUa1[" << i << "]  = " << highErrYieldUa1[i] <<endl;
    cout << "   relErrYieldUa1[" << i << "]   = " << relErrYieldUa1[i] <<endl;
    cout << "err2YieldUa1[" << i << "]        = " << err2YieldUa1[i] <<endl;
    cout << "   lowErr2YieldUa1[" << i << "]  = " << lowErr2YieldUa1[i] <<endl;
    cout << "   highErr2YieldUa1[" << i << "] = " << highErr2YieldUa1[i] <<endl;
    cout << "   relErr2YieldUa1[" << i << "]  = " << relErr2YieldUa1[i] <<endl;
    cout << endl;
     
   i++;
} // while(!fileCrossUa1.eof())
fileCrossUa1.close();
//if (fileCrossUa1) { delete fileCrossUa1; }
//fileCrossUa1=0;

cout << "================================" << endl;
cout << endl;
cout << "Finished reading UA1 YIELD data" <<endl;
cout << "Number of bins read: " << i <<endl;
cout << endl;
cout << "================================" << endl;
cout << endl;

return i;
}
Commit	Line	Data
2803ac99	1	Int_t readUa1()
	2	{
	3	// read UA1 data from file
	4	//
	5
	6	cout << endl;
	7	cout << "================================" << endl;
	8	cout << endl;
	9	cout << "read UA1 YIELD data from file" <<endl;
	10	cout << "Number of bins: " << binsUa1 <<endl;
	11	cout << "Filename: " << filenameUa1 <<endl;
	12	cout << endl;
	13	cout << "================================" << endl;
	14	cout << endl;
	15
	16	ifstream fileCrossUa1;
	17	fileCrossUa1.open(filenameUa1);
	18
	19	Int_t i = 0;
	20	while(!fileCrossUa1.eof()) {
	21	if(i == binsUa1) break;
	22	// textfile content: pt:cross_yield:(stat + syst error added linerarly)
	23	fileCrossUa1 >> centerPtUa1[i] >> crossUa1[i] >> errCrossUa1[i];
	24
	25	ptUa1[i] = centerPtUa1[i];
	26
	27	// the width of the pt bins, currently hardwired....
	28	widthPtUa1[i] = 0.1;
	29	if (centerPtUa1[i] > 4) widthPtUa1[i] = 0.2;
	30	if (centerPtUa1[i] > 6) widthPtUa1[i] = 1.0;
	31
	32	errPtUa1[i] = widthPtUa1[i] / 2.0;
	33	lowPtUa1[i] = centerPtUa1[i] - errPtUa1[i];
	34	highPtUa1[i] = centerPtUa1[i] + errPtUa1[i];
	35
	36	/*
	37	lowStatCrossUa1[i] = crossUa1[i] - statCrossUa1[i];
	38	highStatCrossUa1[i] = crossUa1[i] + statCrossUa1[i];
	39	relStatCrossUa1[i] = statCrossUa1[i] / crossUa1[i];
	40	lowSystCrossUa1[i] = crossUa1[i] - systCrossUa1[i];
	41	highSystCrossUa1[i] = crossUa1[i] + systCrossUa1[i];
	42	relSystCrossUa1[i] = systCrossUa1[i] / crossUa1[i];
	43	*/
	44	lowErrCrossUa1[i] = crossUa1[i] - errCrossUa1[i];
	45	highErrCrossUa1[i] = crossUa1[i] + errCrossUa1[i];
	46	relErrCrossUa1[i] = errCrossUa1[i] / crossUa1[i];
	47	err2CrossUa1[i] = errCrossUa1[i];
	48	lowErr2CrossUa1[i] = lowErrCrossUa1[i];
	49	highErr2CrossUa1[i] = highErrCrossUa1[i];
	50	relErr2CrossUa1[i] = relErrCrossUa1[i];
	51
	52	yieldUa1[i] = crossUa1[i] * (avgToHadr / sigmaInelUa1);
	53	/*
	54	statYieldUa1[i] =
	55	lowStatYieldUa1[i] =
	56	highStatYieldUa1[i] =
	57	relStatYieldUa1[i] =
	58	systYieldUa1[i] =
	59	lowSystYieldUa1[i] =
	60	highSystYieldUa1[i] =
	61	relSystYieldUa1[i] =
	62	*/
	63	errYieldUa1[i] = errCrossUa1[i] * (avgToHadr / sigmaInelUa1);
	64	lowErrYieldUa1[i] = yieldUa1[i] - errYieldUa1[i];
65	highErrYieldUa1[i] = yieldUa1[i] + errYieldUa1[i];
66	relErrYieldUa1[i] = errYieldUa1[i] / yieldUa1[i];
67	err2YieldUa1[i] = errYieldUa1[i];
68	lowErr2YieldUa1[i] = lowErrYieldUa1[i];
69	highErr2YieldUa1[i] = highErrYieldUa1[i];
70	relErr2YieldUa1[i] = relErrYieldUa1[i];
71
72	cout << "ptUa1[" << i << "] = " << ptUa1[i] <<endl;
73	cout << " centerPtUa1[" << i << "] = " << centerPtUa1[i] <<endl;
74	cout << " widthPtUa1[" << i << "] = " << widthPtUa1[i] <<endl;
75	cout << " errPtUa1[" << i << "] = " << errPtUa1[i] <<endl;
76	cout << " lowPtUa1[" << i << "] = " << lowPtUa1[i] <<endl;
77	cout << " highPtUa1[" << i << "] = " << highPtUa1[i] <<endl;
78	cout << "crossUa1[" << i << "] = " << crossUa1[i] <<endl;
79	/*
80	cout << " statCrossUa1[" << i << "] = " << statCrossUa1[i] <<endl;
81	cout << " lowStatCrossUa1[" << i << "] = " << lowStatCrossUa1[i] <<endl;
82	cout << " highStatCrossUa1[" << i << "] = " << highStatCrossUa1[i] <<endl;
83	cout << " relStatCrossUa1[" << i << "] = " << relStatCrossUa1[i] <<endl;
84	cout << " systCrossUa1[" << i << "] = " << systCrossUa1[i] <<endl;
85	cout << " lowSystCrossUa1[" << i << "] = " << lowSystCrossUa1[i] <<endl;
86	cout << " highSystCrossUa1[" << i << "] = " << highSystCrossUa1[i] <<endl;
87	cout << " relSystCrossUa1[" << i << "] = " << relSystCrossUa1[i] <<endl;
88	*/
89	cout << "errCrossUa1[" << i << "] = " << errCrossUa1[i] <<endl;
90	cout << " lowErrCrossUa1[" << i << "] = " << lowErrCrossUa1[i] <<endl;
91	cout << " highErrCrossUa1[" << i << "] = " << highErrCrossUa1[i] <<endl;
92	cout << " relErrCrossUa1[" << i << "] = " << relErrCrossUa1[i] <<endl;
93	cout << "err2CrossUa1[" << i << "] = " << err2CrossUa1[i] <<endl;
94	cout << " lowErr2CrossUa1[" << i << "] = " << lowErr2CrossUa1[i] <<endl;
95	cout << " highErr2CrossUa1[" << i << "] = " << highErr2CrossUa1[i] <<endl;
96	cout << " relErr2CrossUa1[" << i << "] = " << relErr2CrossUa1[i] <<endl;
97
98	cout << "yieldUa1[" << i << "] = " << yieldUa1[i] <<endl;
99	/*
100	cout << " statYieldUa1[" << i << "] = " << statYieldUa1[i] <<endl;
101	cout << " lowStatYieldUa1[" << i << "] = " << lowStatYieldUa1[i] <<endl;
102	cout << " highStatYieldUa1[" << i << "] = " << highStatYieldUa1[i] <<endl;
103	cout << " relStatYieldUa1[" << i << "] = " << relStatYieldUa1[i] <<endl;
104	cout << " systYieldUa1[" << i << "] = " << systYieldUa1[i] <<endl;
105	cout << " lowSystYieldUa1[" << i << "] = " << lowSystYieldUa1[i] <<endl;
106	cout << " highSystYieldUa1[" << i << "] = " << highSystYieldUa1[i] <<endl;
107	cout << " relSystYieldUa1[" << i << "] = " << relSystYieldUa1[i] <<endl;
108	*/
109	cout << "errYieldUa1[" << i << "] = " << errYieldUa1[i] <<endl;
110	cout << " lowErrYieldUa1[" << i << "] = " << lowErrYieldUa1[i] <<endl;
111	cout << " highErrYieldUa1[" << i << "] = " << highErrYieldUa1[i] <<endl;
112	cout << " relErrYieldUa1[" << i << "] = " << relErrYieldUa1[i] <<endl;
113	cout << "err2YieldUa1[" << i << "] = " << err2YieldUa1[i] <<endl;
114	cout << " lowErr2YieldUa1[" << i << "] = " << lowErr2YieldUa1[i] <<endl;
115	cout << " highErr2YieldUa1[" << i << "] = " << highErr2YieldUa1[i] <<endl;
116	cout << " relErr2YieldUa1[" << i << "] = " << relErr2YieldUa1[i] <<endl;
117	cout << endl;
118
119	i++;
120	} // while(!fileCrossUa1.eof())
121	fileCrossUa1.close();
122	//if (fileCrossUa1) { delete fileCrossUa1; }
123	//fileCrossUa1=0;
124
125	cout << "================================" << endl;
126	cout << endl;
127	cout << "Finished reading UA1 YIELD data" <<endl;
128	cout << "Number of bins read: " << i <<endl;
129	cout << endl;
130	cout << "================================" << endl;
131	cout << endl;
132
133	return i;
134	}