File includes good channels ONLY, so staus == 0. muon_coeffs.open("......vec"); while(muon_coeffs){ muon_coeffs>>lay>>pos>>entr>>r2>>mean>>rms>>gmean>>gsigma>>chi>>coe>>status; coeff[lay][pos] = coe; } entr - number of entries r2 - integral after 6 MIPs / entr Histogram (0.2 MIP - 12.2 MIP) -> mean,rms, from previous iteration: gmean,sigma,chi -> results of gaussian fit at this iteration coe - MPV value to use