X-ray emission spectrum modelling of a tungsten target anode at 40 kVp I have an X-ray tube which has the voltage potential from $4kV$ to $60kV$ and the target material is Tungsten $(W, Z=74)$. The emission spectrum provided by the manufacturer is given below:

The spectrum is measured directly out the end of the window using a SiPIN detector, energized to $40 kV /16.2 µA$, and integrated over $80$ seconds live time.  The detector is approximately $8$ cm away from the tube, with a $1.6$ mm collimator immediately in front of the tube window, and a $0.1$ mm aperture immediately in front of the detector window. The space in between is open air.
It can be seen from the spectrum graph that there are Characteristic radiations involved in the form of sudden peaks besides the continuum that is due to the Bremsstrahlung radiation. The manufacturer has not provided the details regarding the spectrum behavior and I am trying to simulate the spectrum under same test conditions but my simulated spectrum shows no such peaks. Looking at the energies of those peaks (from the manufacturer's graph), it seems to be the L shell emissions $(L\alpha,$ $ L\beta)$ since the energies are in those range for Tungsten.
The simulation algorithm (TASMIP) that I have seems to model only the K shell emissions $(K\alpha,$ $ K\beta)$ which are beyond the peak energy of the test, i.e. $40keV$, and therefore absent. My simulations only have the Bremsstrahlung continuum, how do I simulate the L shell emissions so that it approximately matches the manufacturer's spectrum?
 A: I think the answer can be found by looking at the original paper "An accurate method for computer-generating tungsten anode x-ray spectra from 30 to 140 kV" by Boone and Seibert. In the abstract they state

X-ray output measurements (mR/mAs measured at 1 m) were made on a calibrated constant potential generator in our laboratory from 50 kV to 124 kV, and with 0–5 mm added aluminum filtration. 

The actual algorithm consists of interpolating the different measurements that were done - so only things that were observed in their data will be reflected in the spectra. This in contrast with Monte Carlo methods (e.g. MNCP) which include "all the physics". 
Next, upon reading Sisniega, Desgo and Vaquera (2013) I notice that the authors observed that

Inherent filtration for the radiology spectra from TASMIP was found to be equivalent to 1.68 mm Al, as compared to spectra obtained from the mammography model.

Aha! Since the original data was only measured with Xray energies above 50 keV, and with 1.68 mm Al equivalent filtration, they may never have observed these L shell emissions. You need a different model (perhaps the Sisniega paper I linked above will be helpful).
