Projecting dipole into eeg electrode

You can create a scout, define some activity for it (=simulate a signal), and use both to simulate recordings. This would use the orientation defined in your source model (= the normal to the cortex):

From the forward model (matrix Gain), if you know how to obtain a matrix [Nchannels x 1], you can use the function view_topography() to display it, by setting it as the input variable "F" (read the help at the beginning of the function):
[hFig, iDS, iFig] = view_topography(DataFile, Modality, TopoType, F)

changing the orientation of that dipole and again project it to the electrodes

This is not that easy to do, it would involve some scripting and deep understanding of the Brainstorm data structures involves.
A research group in Oldenburg has recently done a simulation study with Brainstorm, they were observing the effect of dipoles orientation and depth on the scalp EEG. The article will be published soon, and maybe they have some example scripts they could already share.
@ArndMeiser @martinbleichner @stedeben

where in BST we do not have "theta"

As mentioned by @SBeumer, you can use [theta,phi,r] = cart2sph(x,y,z).

OK thank you, can we choose 1 dipole and project it to the electrode?

To get the index of one dipole: either use the panel "Get coordinates" or create a scout with only one vertex and get the corresponding index from its structure in the cortex surface.
https://neuroimage.usc.edu/brainstorm/Tutorials/ExploreAnatomy#Get_coordinates
https://neuroimage.usc.edu/brainstorm/Tutorials/ExploreAnatomy#On_the_hard_drive:_Surface
https://neuroimage.usc.edu/brainstorm/Tutorials/SourceEstimation#On_the_hard_drive
https://neuroimage.usc.edu/brainstorm/Tutorials/HeadModel#On_the_hard_drive