In high energy scattering, multiproduction process is unique in its relevance to the total cross section and in its global property such as rapidity and other kinematic distributions. If there is hard interaction, the jet rate and structure is a good arena of perturbative chromodynamics. However, once any hadron is specified, its property and structure must make sense while the global and/or perturbative chromodynamics mechanism still can put important constrains. The relation of the hadron property and structure with its production cross section, distribution etc. can be much more complex than its decay width. In one hand there are many difficulties and challenges in calculation; on the other hand, production process provides unique way to study the details of property and structure of the hadron, which is beyond the approach of its decay process. Here I review our works on such topic in recent years, mainly on the multiquark state production in multiproduction process and the Bethe-Salpeter wave function in exclusive process. For the former, I emphasize on the unitarity of hadronization process and relevant models, so that Almost No multiquark state is observed in multiproduction process; I also address how to calculate hadron molecule production in multiproduction process; the most recent observed T + cc is also stated with its relevance to colour and baryon number fluctuation of the preconfinement clusters. For the latter, I emphasize the Dirac structure of the hadron-quark coupling vertex. Contents 1 Introduction 2 2 Unitarity, colour and baryon fluctuation 2 3 Bethe-Salpeter (BS) wave function in production 3 4 Conclusion 3 References 4