#include #include #include #include #include using namespace std; //// for non-standard W bins (W_start is usually 0.0) int GetWIndex(double W, double W_start, double W_width) { double aa = W/W_width ; cout<<"aa="<