add channelPolarization.cpp