8550415 />and respectively) in p + p collisions at = 7 TeV at LHC would be,
For the calculation of the rapidity distribution from the Equation (4) we can make use of a standard relation as given below:
For p + p collisions, the calculated rapidity distribution equation at RHIC-energy = 200 GeV is
In a similar fashion, the SCM-based rapidity distribution equation at LHC-energy = 7 TeV in p + p collisions has been given hereunder
4.2. Invariant Yields in d + Au, Cu + Cu and Au + Au Collisions at RHIC-Energy
= 200 GeV
From the expression (4), we arrive at the invariant yields for the J/Ψ production in reactions for mid and forward-rapidities.
For the case of Cu + Cu most central collisions (0% - 20%) at RHIC, the SCM-based calculated theoretical invariant yields for the rapidities and are given by the following equations respectively;
Similarly, for Au + Au collisions at = 200 GeV at RHIC, the equations of transverse momenta spectra for 0% - 20% centrality regions are given by the undernoted relations.
For calculating the values of NR, in general, we have used the values of from Refs.  and .
4.3. The Nuclear Modification Factor RAB
There is yet another very important observable called nuclear modification factor (NMF), denoted here by which for the production of J/Ψ is defined by 
The SCM-based results on NMFs for Cu + Cu and Au + Au collisions for forward rapidities are deduced on the basis of Equation (6), Equation (15), Equation (17) and Equation (19) and they are given by the undernoted relations
herein the value of to be used is  for Cu + Cu collisions and for Au + Au collisions it is taken as .
4.4. Analysis of the Figures
In Figures 1(a) and (b), we have drawn the total inclusive cross-sections for the production of J/Ψ-mesons in p + p collisions in different rapidities at RHIC and LHC energies = 200 GeV and 7 TeV respectively. The solid lines in the figures are depicting the SCM modelbased results with the help of the Equations (6)-(9) while the experimental measurements are taken from Ref. [19, 20] respectively.
In Figure 2 we have plotted the rapidity distributions for J/Ψ-production in p + p collisions at = 200 GeV. Data in the figure are taken from Ref.  and the line shows the SCM-based output.
Similarly, in Figure 3, we have drawn the rapidity distribution for J/Ψ production in p + p collisions at = 7 TeV as function of y. The solid lines depict the SCM-based results [Equation (12)] and the points indicate the experimental measurements .
Figure 1. Plot of the invariant cross-section for J/Ψ production in proton-proton collisions at (a) = 200 GeV and (b) = 7 TeV as function of. The data points are from  for (a) and from  for (b). The solid curves show the SCM-based results.
Figure 2. Plot of the rapidity distribution for J/Ψ production in proton-proton collisions at = 200 GeV as function of y. The data points are from [15,21]. The solid curves show the SCM-based results (Equation (11)).
Figure 3. Plot of the rapidity distribution for J/Ψ production in p + p collisions at = 7 TeV as function of y. The data points are from . The solid curves show the SCM-based results (Equations (15) and (16)).
In Figure 4, we have drawn the solid lines depicting the SCM-based results for invariant yields for J/Ψ production in d + Au collisions at = 200 GeV with the help of two Equations (13) and (14) against the experimental measurements .
The experimental results for the invariant yields of J/Ψ production as a function of transverse momenta for Cu + Cu collisions at = 200 GeV are taken from Ref.  at centrality 0% - 20% and are plotted in Figure 5. The solid lines in the figure show the SCM-induced results.
Figure 4. Plot of the invariant yields for J/Ψ production in d + Au collisions at = 200 GeV as function of. The data points are from . The solid curves show the SCMbased results.
In the Figure 6, the solid lines are the plots of SCMbased invariant yields vs. as described by Equations (17) and (18) at forward and mid-rapidities for Au + Au collisions at = 200 GeV, while the dotted curve in the Fig. shows results of coalescence model . The experimental data points in the Figure 6 for the invariant yields of J/Ψ production as a function of transverse momenta at centrality values 0% - 20% and at the rapidities and respectively are taken from the PHENIX Collaboration .
Figure 5. Transverse momenta spectra at and for J/Ψ production in Cu + Cu central collisions at = 200 GeV. The data are taken from . The solid curves depict the SCM-based results.
Figure 6. Plot of the invariant yields for J/Ψ production in Au + Au collisions at = 200 GeV as function of. The data points are from . The solid curve shows the SCM-based results while the dotted curve depicts the Coalescence Model .
In Figure 7(a), we plot vs. for 0% - 20% central region in Cu + Cu and Au + Au collisions at = 200 GeV. The solid lines in the figure show the SCM-based results (Equations (20) and (21)) against the experimentally measured results [15,16]. The dotted lines in the figure. represent the double Color Filtering approach .
And in Figure 7(b), we plot vs. for 0% - 20% central region in Au + Au collisions at =
200 GeV and for Pb + Pb collisions at = 2.76 TeV. The solid lines in the figure show the SCM-based calculations for different
5. Summary and Outlook
Let us first concentrate on what we have achieved here: 1) The features related to pT—spectra for J/Ψ production in some particle-particle and nuclear collisions at various high energies have been reproduced quite successfully; 2) The characteristics of rapidity spectra in p + p collisions at TeV energies have been brought out somewhat modestly satisfactorily; 3) The features of nuclear modification factors in Cu + Cu, Au + Au and Pb + Pb reactions at RHIC and LHC energies have also been brought out with the help of the applied model. Besides, some of our model-based results have also been compared with the performances on the same observables by some competing models of “standard” variety. And these comparisons with data and the results obtained by some other models reveal that SCM-based results describe the features of the data, at least, not inferior to the performance by the other approaches grounded on the “Standard” model ilk. In the past such were the recurrent observations made by us valid for many other obsevables measured in the various high statistics high energy particle and nuclear experiments.
Thus, summing up our past experiences and consider-
Figure 7. (a) Plot of the for J/Ψ production in Cu + Cu and Au + Au collisions at = 200 GeV as function of. The data points are from [15,16]. The solid and dotted curves show respectively the SCM and the Double Color Filter-oriented  results. (b) Plot of the vs. for J/Ψ production in Au + Au collisions at = 200 GeV and Pb + Pb collisions at = 2.76 TeV. The red circles depict Au + Au collisions  while the black squares represent Pb + Pb collisions . The solid curves show the SCM-based results while the dotted curve depicts the Gluon Saturation Approach .
ing the weightage of the results reported here, we are forced to comment finally that this work essentially represents a case of paradigm shift in the domain of particle theory, as we have eschewed the conventional views of approach to J/Ψ production in the “standard” framework. And this is just the reflection of our radical views about the particle structure and the nature of particle collisions. Obviously we obtain the fair agreement with data on some observables without inductions of 1) any QGP concept; 2) any prognosis of suppression or enhancement of J/Ψ-production. The production of J/Ψ- particles resembles all other hadrons.
The authors would like to express their thankful gratitude to the learned Referee for his/her valuable remarks and constructive suggestions in improving the earlier draft of the manuscript.