Extracting the pole and Breit-Wigner properties of nucleon and $\Delta$ resonances from the $\gamma N\to K\Sigma$ photoproduction

We have developed a covariant isobar model to phenomenologically explain the four possible isospin channels of $K\Sigma$ photoproduction. To obtain a consistent reaction amplitude, which is free from the lower-spin background problem, we used the consistent electromagnetic and hadronic interactions described in our previous report. We used all available experimental data, including the recent MAMI A2 2018 data obtained for the $K^0\Sigma^+$ and $K^0\Sigma^0$ isospin channels. By fitting the calculated observables to these data we extract the resonance properties, including their Breit-Wigner and pole parameters. Comparison of the extracted parameters with those listed by the Particle Data Group yields a nice agreement. An extensive comparison of the calculated observables with experimental data is also presented. By using the model we investigated the effects of three different form factors used in the hadronic vertex of each diagram. A brief discussion on the form factors is given.


I. INTRODUCTION
For more than 50 years kaon photoproduction off a nucleon has gained a special interest in the hadronic physics community. Early theoretical work was reported in 1957 [1], but a more comprehensive analysis with fitting to experimental data was just started in 1966 [2]. Since then many efforts have been devoted to explain this reaction, ranging from quark to hadronic coupled-channel models, as briefly mentioned in the introduction part of our previous report [3].
In the beginning, the main motivations to study kaon photoproduction off the nucleon were merely to obtain theoretical explanation of the reaction process. However, it has been soon realized that an accurate theoretical model describing this elementary process is also useful in many branches of hadronic and nuclear physics. In hadronic physics the model is indispensable in the investigation of the missing resonances that have considerably large branching ratio to the strangeness channel [4], the narrow resonance which is also predicted to have a large branching ratio to the KΛ channel [5], and the resonance hadronic coupling constant that measures the strength of the interaction between kaon-hyperon final state and the resonance [6]. In the nuclear physics this elementary model is required in calculating the cross section of hypernuclear photoproduction [7], which is the main observable in the investigation of hypernuclear spectroscopy.
Recently, the nnΛ electroproduction on a tritium target, i.e, the 3 H(e, e ′ K + )nnΛ process, has been performed at Jefferson Lab Hall A (JLab E12-17-003) and the result is currently being analyzed [8]. There have been intensive discussions on whether the nnΛ system could be bound or would lead to a resonant state. Thus, this experiment is expected to shed light on the nnΛ puzzle. Furthermore, an accurate measurement of hypertriton electroproduction has been proposed and conditionally approved as the JLab C12-19-002 experiment [8]. The experiment is expected to further elucidate the Λ binding energy in the hypertriton, since the result of previous measurement indicated a stronger binding energy [9]. Therefore, an accurate elementary model, describing the photo-and electroproduction of kaon on the nucleon target is timely and urgently required.
Our previous model used to this end was Kaon-Maid [10], which includes kaon photo-and electroproduction off the nucleon in six isospin channels. However, due to tremendous increase in the number of experimental data, especially in the case of photoproduction, Kaon-Maid started to show its deficiency. To get rid of this problem we have started to improve the model in the photoproduction sector, for which experimental data dominate our present database.
In the previous works we have developed a new and modern covariant isobar model for kaon photoproduction off the nucleon γ + p → K + + Λ and γ + n → K 0 + Λ [11]. The model fits nearly 9000 experimental data points and employs the consistent hadronic and electromagnetic interactions that eliminate contributions of high-spin resonance background [12,13]. The latter is widely known as an intrinsic problem that plagues the formalism of high-spin propagators used to describe the contribution of nucleon and delta resonances. In this paper we extend the model, based on the covariant effective Lagrangian method, to include the other four isospin channels in the KΣ photoproduction.
We have organized this paper as follows. In Sec. II we discuss the formalism used in our model. In principle, we use the same interaction Lagrangians as described in our previous paper for the KΛ photoproduction, γ +p → K + +Λ [13]. In addition, we also briefly discuss the formalism used to extract resonance masses and widths at their pole positions in this section. In Sec. III we present the result of our analysis and compare the calculated observables with the available experimental data.
To describe the accuracy of our model in details, an extensive comparison of polarization observables is given in this section. In Sec. IV we summarize our analysis and conclude the important findings. The extracted Breit-Wigner masses and widths of the resonances used in the model are listed in Table VII of Appendix A. II. FORMALISM

A. The Model
In the present work we consider the photoproduction process of KΣ on a nucleon, i.e., γ(k) + N (p) → K(q) + Σ(p Σ ) . (1) Based on the isospin and strangeness conservations, Eq. (1) implies four different photoproduction processes given in Table I.   TABLE I. Four possible isospin channels for KΣ photoproduction off the nucleon. The corresponding threshold energies are also listed in terms of the photon laboratory energy k thr 0,lab and total c.m. energy W thr .
No. Channel The scattering amplitude of these reactions is calculated from the first-order Feynman diagrams shown in Fig. 1. According to their intermediate states the diagrams can be grouped into three main channels, i.e., the s-, t-and u-channel, with the corresponding Mandelstam variables are defined as (2) Note that the notation of the momenta written in Eq. (2) is given explicitly in Eq. (1). The corresponding vertex factors can be obtained from the effective Lagrangian approach, specifically by using the prescription given in Refs. [14][15][16]. In our previous study of the K + Λ photoproduction [13] the Lagrangians were constructed according to the method proposed in Ref. [12] in order to be consistent with the formulation of high spin (J > 3/2) propagators. In the case of KΣ photoproduction, the hadronic and electromagnetic Lagrangians of the s-channel spin-(n + 1/2) particle with positive parity reads respectively, where Ψ is the field of the Σ particle, F µν is the antisymmetric tensor of photon field, andΨ µ1···µn is the modified RS-field of spin-(n + 1/2) particles constructed to make the interaction consistent as proposed in Ref. [12]. The modified RS-field reads whereψ is the original RS-field for spin-(n+1/2) particles and the interaction operator O is defined by where P (ν) and P (λ) indicate the permutations of all possible ν and λ indices, respectively, and The propagator used for calculating the scattering amplitude is obtained from the completeness relation of the RS-fields. This propagator, however, bears unphysical lower spin projection operator and eventually would yield unphysical contribution to the scattering amplitude if it was not properly handled. In this work, by using the interaction Lagrangians given by Eqs. (3) and (4), the remaining lower spin terms are automatically removed from the amplitude, leaving only the pure spin-(n + 1/2) contribution that comes from the projection operator with the projection operator for spin-n particle where k max is equal to n/2 if n is even and to (n − 1)/2 if n is odd, P µν (p) = (−g µν + p µ p ν /p 2 ), and the coefficient A n k is defined by 3 .  As mentioned above the constructed scattering amplitude is free from the unphysical lower-spin contribution that originates from the RS-fields. In the compact form the amplitude can be written as

* * *
where p R is the four-momentum of resonance particle and the vertex factors Γ had µ1···µn and Γ em ν1···νn are derived directly from the interaction Lagrangians given by Eqs. (3) and (4), respectively. For the purpose of numerical computation of observables we need to calculate the total scattering amplitude, which is obtained by adding the back-ground M back fi and resonance M R fi contributions, i.e., where the summation on the right hand side is performed over all nucleon resonances considered in the present work. The total amplitude is decomposed into the gauge and Lorentz invariant matrices M i where M i denotes the gauge and Lorentz invariant matrices given by [13,31] with P = 1 2 (p + p Σ ) and ǫ µνρσ is the four dimensional Levi-Civita tensor. The required cross section and polarization observables can be calculated from the functions A i given by Eq. (13), which depend on the Mandelstam variables given in Eq. (2).
Note that as seen from Eq. (11) the present formalism introduces the high momentum dependence p 2n R that might lead to a non-resonance behavior of scattering amplitude at high energies. To alleviate this problem we need a stronger hadronic form factor that can sufficiently suppresses the divergence of the amplitude at high energies. The widely use form factor is the dipole one, i.e., with Λ the cutoff parameter and M the resonance mass. However, we found that such a form factor does not have a sufficiently strong suppression for this purpose. Therefore, in this study we propose the use of the multi-dipole form factor as well as the Gaussian one and investigate the effects of these form factors on the constructed model. For the multi-dipole form factor we will present the result with n = 3, for which we obtained the best agreement with experimental data, and denote the model with HFF-P3. The models that use the form factors given by Eqs. (18) and (20) are denoted with HFF-P1 and HFF-G, respectively. Equations (3) and (4) indicate that for each resonance with J ≥ 3/2 there are four unknown coupling constants which might be considered as free parameters. These coupling constants can be extracted from fitting the calculated observables to the corresponding experimental data. In the present work the fitting process was performed by using the CERN-MINUIT code [32] to minimize the value of where N data and N par indicate the numbers of experimental data and free parameters used in the fit, respectively, σ exp i and σ th i are the i-th values of experimental and theoretical observables, and ∆σ exp i is the corresponding experimental error bar.
Experimental data used in the present analysis are obtained from a number of experimental collaborations as listed in Table II. Note that the KΣ photoproduction offers 4 possible isospin channels. Among them the γp → K + Σ 0 channel has the largest number of experimental data, as can be seen in Table II. Furthermore, experimental data of this channel are available for different types of observables, which can be expected to complete our understanding of this reaction. Nevertheless, although with a limited number and observable types, the existence of experimental data in the other three isospin channels is very important to constrain the extracted coupling constants as well as the predicted observables [33].

B. Resonances Properties
In this study, we extract a number of important resonance properties, i.e., their masses and total widths evaluated at pole, their partial widths, and their individual contribution to the process. The evaluation of mass and width at pole starts with a complex root equation of the denominator of the scattering amplitude, which reads where m R and Γ R are the Breit-Wigner mass and width, respectively. The variable s p is defined as where m p and Γ p are the mass and width evaluated at the pole, respectively. We can clearly see that the solutions of m p and Γ p were actually simple. However, in the present work, we use an energy-dependent width Γ(s) that is directly proportional to the total width Γ R [34]. As a consequence, we cannot analytically solve the root equation given in Eq. (22). Therefore, in the present analysis we solve this equation numerically. Furthermore, we can also compute the partial decay widths Γ γp Γ KΣ /Γ tot to conveniently compare the values of extracted coupling constants in this study. To this end, we start with the interaction Lagrangians for spin J = n + 1/2 resonances and for each J we obtain the decay width formula. Note that in the present analysis we can only obtain the product of hadronic and electromagnetic decay widths, because in the single channel analysis only the product of the electromagnetic and hadronic coupling constants can be extracted. The formulation of these decay widths can be found in our previous study [13].
The significance of each resonance can be also evaluated by excluding the specific resonance in the fitting process. Mathematically, the significance of an N * resonance is defined through where χ 2 all is obtained from fitting the data by using all resonances, while χ 2 all−N * is obtained by using all but a specific N * resonance. In the present work Eq. (24) is also used for investigating the significance of ∆ resonances. The significance of a resonance is not only very useful for investigation of the role of each resonance in the reaction, but also for a practical guidance to simplify the model if minimizing the number of resonances used in the model is important. A simple covariant isobar model is very important for the application in few-body nuclear physics, for which numerical computation and accuracy are very demanding.

A. General Results
Table III compares the extracted coupling constants of the Born terms obtained for the 3 different hadronic form factor models. By comparing the values of χ 2 , we can clearly see that the model HFF-P1 yields less agreement with experimental data. This indicates that to achieve a better agreement with experimental data we have to use a softer form factor. The result obtained from model HFF-P3 corroborates this. However, since the three form factors used in the present analysis have different forms, see Eqs. (18)- (20), it is difficult to estimate the suppression imposed by the form factors on the amplitude by merely comparing their cutoffs listed in Table III. Therefore, in Fig. 2 we plot the form factors for both Born and resonance terms as a function of the Mandelstam variable s and shows the energy region covered by the experimental data in the fitting database. From the top panel of Fig. 2 we may conclude that the Gaussian form factor is the softest one and, as a consequence, this form factor strongly suppresses the Born amplitude. In contrast to this form factor, the dipole one given by Eq. (18) provides the lightest suppression, whereas the multi-dipole form factor yields the relatively moderate suppression. Nevertheless, even the dipole form factor decreases the contribution of Born terms significantly. Thus, from Fig. 2 we may safely say that all models analyzed in this work are resonance-dominated models because, as shown in the bottom panel of Fig. 2, near the resonance mass the form factors do not suppress the amplitude significantly.
In the case of resonance, within the covered energy the suppression effect of form factor is clearly not symmetrical. This is understandable because the resonance masses used in the present analysis are below 2.3 GeV. However, this asymmetrical suppression is required for the covariant description of a resonance due to the large contribution of a Z-diagram that indicates the existence of a particle and an antiparticle in the intermediate state [35]. This contribution increases quickly as the energy of resonance increases beyond the resonance mass and, therefore, requires an increasing suppression. In our previous study [35] it was shown that the dipole form factor given by Eq. (18) is suitable for this purpose.
In Table III we also show the χ 2 contribution from each channel. It is apparent from this table that the model HFF-P3 shows the best agreement with the experimental data (lowest χ 2 ) from all but the γn → K 0 Σ 0 channel. As we will see later, when we compare the observables, the model deficiency in this channel originates from the discrepancy between the calculated differential cross section and the experimental data at forward angles. Nevertheless, the effect of this discrepancy is less significant compared to those obtained from the other three isospin channels. Furthermore, by analyzing the sources of experimental data we found that the model HFF-P3 yields a nice agreement not only with the SAPHIR [23], but also with the MAMI [30] data.
The calculated total cross sections obtained from the three models are compared with the available experimental data in Fig. 3. Note that although there are no data for the K + Σ − total cross section, both LEPS 2006 [22] and CLAS 2010 [27] collaborations have measured the K + Σ − differential cross section, whereas the LEPS 2006 collaboration has obtained the K + Σ − photon asymmetry [22]. Since these data have been included in our fitting database, we believe that the calculated K + Σ − total cross section shown in Fig. 3 is also accurate. From Fig. 3 we may conclude that all three models can nicely reproduce the total cross section data. It is understandable that this new result is quite different from that of Kaon-Maid, since Kaon-Maid was fitted to the old SAPHIR data [36,37], except for the γn → K 0 Σ 0 reaction, for which Kaon-Maid yields a fairly good prediction to the recent MAMI data. We believe that the latter is pure coincidence and, in fact, the discrepancy between Kaon-Maid prediction and experimental data is clearly seen near the threshold and higher energy regions. It is also important to note that in the K 0 Σ + channel the total cross section indicates two resonance peaks. Interestingly, they originate from the N (1720)P 13 and N (1900)P 13 states that have been found to be important to describe both KΛ and KΣ photoproductions [31,38]. The first peak does not appear in the K 0 Σ + channel,   [4] and different models analyzed in the present work (HFF-P1, HFF-P3, and HFF-G), compared with the presently available experimental data. Notation of the curves is given in the lowest panel, whereas notation of the data is given in the corresponding panels. Note that the experimental data shown in this figure were not included in the fitting process of the present work.  whereas the second one is shifted above 1.9 GeV due to the interference with other resonances, whose extracted masses are heavier than 1.9 GeV.
From Figs. 4 and 5 we can clearly see that all models do not exhibit significant variation in the differential cross sections, except for the extreme kinematics, i.e., in the forward and backward directions and in the higher energy region. In this kinematics, the interference between the resonance and background terms is found stronger than in any other regions. Meanwhile, experimental data in this kinematic are scattered and, in fact, in certain energy intervals there are no data available to constrain the models. Therefore, during the fitting process this condition yields significant variations among the models.
The lack of experimental data in certain energy and angular regions also occurs in the case of polarization observables. It is well known that unlike the cross section, the polarization observables are very sensitive to the ingredient of the reaction amplitude. Therefore, they can severely constrain the flexibility of the model during the   Fig. 4, but for the angular distribution of differential cross section. The corresponding value of the total c.m. energy W in GeV is given in each panel. fitting process. As a result, the calculated χ 2 reported in the present work originates mostly from the polarization data. From their formulations it is easy to understand that the single polarization observables are the simplest polarization observables that depend sensitively on all ingredients of the model, i.e., not only the resonance configuration, but also the background structure. Thus, they are very useful to constrain the models that predict similar trend in differential cross section, but significantly different beam, target, or recoil polarization observables. To this end, we also note that there was a discussion on whether the right model should be resonanceor background-dominated, or both resonance and back-ground are equally contributing [4,39,40]. Hence, the observables can be used to help alleviate the problem. Furthermore, in the certain kinematical region, where differential cross section data are not available, single polarization observables provide an important tool to shape the trend of differential cross section.
From Figs. 6 and 7 we found that in the kinematical region where the differential cross section data are abundant, e.g., at cos θ = −0.70, the three models show a large variation of the calculated recoil polarization P , in contrast to the predicted differential cross sections. On the other hand, in the kinematical region where differential cross section data are not available, e.g., at cos θ = −0.95, we find that the models yield significant variation in both cross section and polarization observables. This phenomenon is also shown by the other single polarization observables, i.e., the photon asymmetry Σ shown in Figs. 8 and 9, as well as the target asymmetry T shown in Figs. 10 and 11. Therefore, we may conclude that single polarization observables are potential to become important tools to determine the correct ingredients of the Born and resonance terms in the KΣ photoproduction channels.
Double-polarization observables are even more sensitive to the constituent of the reaction amplitude compared to the single polarization ones. Therefore, doublepolarization observables can be considered as the main constraint to the models.
For the K + Σ 0 channels we notice that out of 12 possible double polarization observables, experimental data are only available for four types of the beam-recoil observables, i.e., C x , C z , O x , and O z . In the case of the C x and C z , unfortunately, the currently available experimental data have large error bars, especially in the extreme kinematics as shown in Figs. 12-15. As a consequence, all models investigated in the present work yield a large variance in this kinematics. As discussed before, the effect propagates to the calculated cross sections, for which large variation can be observed in the extreme kinematics (see Figs. 4 and 5). Furthermore, we also observe this effect in the recoil polarization as clearly shown in Figs. 6 and 7.
Unlike the double-polarization observables C x and C z , the available data for O x and O z are relatively more accurate, as shown in Figs. [16][17][18][19]. We found that the three models have relatively similar trend in the region where experimental data are sufficiently available. However, a different situation is shown in the extreme kinematics, where experimental data are extremely scarce. Actually, the effects of this polarization is small only at the extreme kinematics, where the cross section data are accidentally scarce. Nevertheless, to obtain a more accurate model, including the background part, these data are urgently required, since contribution of the background part in this kinematic is stronger than in any other kinematical regions.
As stated before, the model HFF-P3 yields the best agreement with experimental data, although the difference is not large compared to the model HFF-G in the case of γp → K 0 Σ + channel. From Figs. 20 and 21, we may conclude that contribution to the first peak in the differential cross section at W ≈ 1.75 GeV (see Fig. 20) stems mostly from the backward angles, since as shown in Fig. 21 the cross section in this energy region is back-ward peaking. The angular distributions of differential cross section obtained from the model HFF-P3 exhibit a clearer picture of the origin of both peaks in the total cross section. By comparing with the results from the other two models we might conclude that the difference between them originates from different ways of reproducing the recoil polarization data shown in Figs It is also obvious that at very forward and backward angles, only the model HFF-P3 shows different recoil polarization. The difference originates from the lack of data in the extreme kinematics from the MAMI collaboration. Thus, we might conclude that the different behavior shown by all models originates from the freedom to predict the missing data in the extreme kinematics.
The other Σ channels that use neutron as target exhibit a similar pattern as shown in Figs. 24-27. As expected,    there are differences in the kinematical regions, where experimental data are lacking. Interestingly, we observe a unique behavior in the prediction of the model HFF-G for the γn → K + Σ − , i.e., there are tiny peaks near the threshold region at backward and forward angles. The peaks, which will be discussed later when we discuss the resonances properties, originate from the contribution of high spin resonances which only appears near threshold in the model HFF-G. Unfortunately, this behavior cannot be further explored due to the lack of experimental data in this kinematics. As shown in Figs. 28 and 29, the same situation also happens in the case of the polarization asymmetry Σ that usually provides a severe constraint to the model. Nevertheless, we will discuss this topic from another perspective in the following subsection.

B. Resonances Properties
Having constructed an isobar model that can nicely reproduce all KΣ photoproduction data we are in the position to study the properties of baryon resonances through their electromagnetic and hadronic interactions. The resonance properties that are of interest in the hadronic physics community are the resonance mass and width evaluated at its pole position as well as the extracted resonance partial decay width. Evaluation of the resonance mass and width at pole position has a clear advantage, since the result is model independent. Therefore, the properties of resonances evaluated at pole position in the present study are comparable with those obtained from other model-independent investigations.
The result is shown in Table IV, where we compare the resonance masses and widths evaluated at their poles obtained from the present work with those given by PDG.
For completeness, we present the corresponding Breit-Wigner masses and widths in Table VII of Appendix A. The calculated resonance masses show a good agreement with the PDG values, especially in the case of model HFF-G. On the other hand, the calculated widths show some discrepancies with the PDG ones. There is no obvious pattern observed in these discrepancies and, in fact, it seems to be random for all models. We believe that this phenomenon originates from the use of single channel analysis, where the resonance widths are not unitary defined. In the case of the current best models (HFF-G and HFF-P3), there are two adjacent resonances that show an interesting phenomenon, i.e., the N (1700)D 13 and N (1720)P 13 states. In these models the resonances produce different cross section behavior near the threshold by switching their pole positions. However, in the model HFF-G the pole positions of these resonances are closer to the PDG values. Put in other words, compared to other nucleon resonances the N (1700)D 13 state is very likely to be important in the threshold region. Nevertheless, this finding still needs further investigation once the KΣ photoproduction data near threshold are sufficiently available. It is also interesting to explore how the resonances complement to each other through their significance in the   strong impact on the agreement between model calculation and experimental data. This can be understood because there are no resonances with spin equals to 7/2 and adjacent mass. Another example that shows the importance of high spin resonances is exhibited by the ∆(2400)G 39 state. This resonance has a relatively high significance compared to the other ones. By looking at the PDG estimate in Table IV we obviously find that the extracted masses are significantly heavier. We note that this occurs in both Breit-Wigner and pole position methods. Furthermore, Table IV also indicates that in the model HFF-G the extracted mass is about 500 MeV heavier than the PDG estimate. Therefore, we might conclude that this resonance is responsible for improving the agreement with experimental data in the higher energy region. Nevertheless, for a more conclusive finding, investigation at this energy region is strongly advocated in the future. Since we have observed similar cases in the present work, we might conclude that the high spin resonances are indispensable in the KΣ photoproduction process. The role of spin-7/2 and 9/2 resonances in the K + Λ photoproduction has been thoroughly analyzed in Ref. [13]. For spin-11/2 and 13/2 nucleon and delta resonances the analysis is still ongoing and will be published soon [42]. Figures 30 and 31 also show the less significance resonances, e.g., the N (1680)F 15 and ∆(1930)D 35 , which have ∆χ 2 less than 1%. In the case of the N (1680)F 15 resonance this is understandable, since it contributes to the background part of the model and, furthermore, there is a spin 5/2 resonance with adjacent mass, i.e., the N (1675)D 15 . For the ∆(1930)D 35 resonance, although there is evidence for its branching to the KΣ channels, the present work indicates that this state is less significance. A closer look to the mass position reveals that the mass of this resonance lies among the masses of six ∆ resonances from 1900 to 1950 MeV (see Fig. 31). We believe that the less significance of the ∆(1930)D 35 resonance originates from this phenomenon. As a consequence, if we needed to simplify the model but not its accuracy, this resonance could be excluded.
In Tables V and VI we compare the resonance partial widths obtained in the present analysis and those listed by PDG. Table V shows the comparison for the proton channels, whereas Table VI shows that for the neutron channels. Since the value of partial width for ∆ resonances is the same for both proton and neutron channels, we only show them in Table V. From Tables V and VI, we may conclude that most of the extracted partial widths are in good agreement with those listed by PDG. Nevertheless, we also note that there are large discrepancies in the case of the N (1900)P 13 , N (2060)D 15 , and ∆(1920)P 33 . This would be an interesting phenomenon for future investigation, since for the three resonances the discrepancies between our best model HFF-G and the PDG values are milder.
In our previous study [13] we also observed a similar phenomenon, i.e., the extracted partial widths for a num- In the present work we found that this problem can be overcome by using the Gaussian hadronic form factor given in the HFF-G model, which strongly suppresses the amplitude except near the resonance mass. Comparison among the three models given in Tables V and VI shows that the Gaussian form factor reduces the extracted partial waves to the values closer to the PDG estimates, except for the N (1900)P 13 resonance, for which the extracted value is still much larger than the PDG one. Therefore, future investigation should address this problem, since the N (1900)P 13 is currently a four-star resonance.
A closer look at Tables V and VI reveals that the N (1720)P 13 resonance experiences the strongest reduction if we use the Gaussian form factor (model HFF-G). Compared to the model HFF-P3 the fractional decay width of this resonance is reduced by more than 90%. This is due to the change of the resonance role, between the N (1700)D 13 and N (1720)P 13 states as we mentioned earlier. Since the N (1720)P 13 is closer to the threshold, the branching ratio to this reaction tends to be smaller. In addition to this fact and earlier analysis, we have more confidence to say that the N (1720)P 13 state is more likely to be found near the reaction threshold. FIG. 24. Energy distribution of the γn → K + Σ − differential cross section obtained from all models. Notation of the curves is as in Fig. 3. Experimental data are obtained from LEPS 2006 (open squares [22]) and CLAS 2010 (solid squares [27]) collaborations.

IV. SUMMARY AND CONCLUSION
We have analyzed the KΣ photoproduction data for all four possible isospin channels by using a covariant isobar model and including all appropriate nucleon and   Energy distribution of the γn → K 0 Σ 0 differential cross section obtained from all models. Notation of the curves is as in Fig. 3. Experimental data are obtained from MAMI 2018 (solid squares [30]) collaborations.  delta resonances. We used the consistent Lagrangians for hadronic and electromagnetic interactions to eliminate the problem of lower-spin background contribution. In this analysis, three different form factors have been used in the hadronic vertices, i.e., the dipole, multidipole, and Gaussian ones. The present model yields a nice agreement between calculated observables and experimental data. The best agreement is shown by the models that employ the multidipole and Gaussian form factors. We have also extracted the Breit-Wigner masses and widths of the nucleon and delta resonances as well as