-Webster's Ninth New College Dictionary
iversion to characterize the
data with different reconstruction version, iversion
can be 5, 6, 7 or 8 (the highest version now is 8).
n_part, n_res and rc
before doing the kinematic fit.
n_part, n_res and rc
before doing the kinematic fit with a resonance.
iversion.le.7
Integer cd, code, ichg, ifit(2), irc, itk1, itk2, itrack Integer i1, i2, i3, i4, i5, i6, i7, i8, i9, i10 Integer i_squid, mem(*), mfit, mzfit, n_part, n_res Integer p_mem(*), p_spec(4,*), q1, q2, rc Integer spec1, spec2, spec3, spec4, type_i Logical newsquid, squidjr Real bc, bchisq, chisq, chi2xy, chi2z, crxy(2), ct Real err(6), error(15), ers1(6), ers2(6), er1(15), er2(15) Real fitcos, fitpgg(4), fitpg1(4), fitpg2(4), ggmass, lxy Real m_res(*), m_res_i, mass(*), mass1, mass2, mass_i, mpar Real p_fit(4,*), pm(5), pull(*), pxyz, p4par(4), rawcos Real rawmgg, rm, seclm(4), sqrt_s, tkp(3), tkpar(6) Real tks1(3), tks2(3), tk1(6), tk2(6), xi0, xmas, zswum(2) Real*8 b(5,6)
This is the source of > 90% of all Telesis problems!
TRY THIS FIRST!
Telesis routines have lots of arguments, so the chances of making a mistake with an argument list are increased.
rc = 0 Successful fit.
1 Neutral track or missing track specified with
inappropriate routine. (Sould only Use pstrkcp or
psextern for these kinds of tracks.)
2 Charged track did not have MFIT=2 at BES.
3 Beam fit failed.
6 Newsquid or Squidjr failed because the error matrix
passed to Leqt2p (IMSL routine) was not positive-
definite. It is normal for this to occur
occasionally. (However the two error messages are not
printed at BES for save the space of log file)
7 Newsquid or Squidjr failed because chi-squared > 50.
This should be the most common reason for failure.
NOTE: rc can not equal 4 or 5 for this version of telesis at BES now.
If you still can't get your kinematic fits to work, the routine
dmpsqdl may be useful for dumping the input and output info,
and some of the contents of the Newsquid Common block.
Performing a kinematic fit with Telesis requires a sequence of subroutine calls in which all of the particles and resonances are specified. The sequence of calls FOR EACH FIT must have the following elements:
List of examples in this file:
XMCHTY COMMON has been included in your
program, so the particle masses xmpich, xmk0, etc. are defined.
(Naturally you may input the necessary masses in any way.)For Newsquid fits it is assumed that the center-of-mass energy has been looked up (and corrected if desired), and is contained in the real variable 'energy'.
It is also assumed that the user is looping over some set of candidate track combinations, or that one particular combination has been chosen. The pions are Trklst track numbers ip1, ip2, ...; the kaons are ik1, ik2, ...; the gammas are ig1, ig2, ..., etc.
Call psinit(energy,npart,nres,rc)
Call psdotrkc(1,xmkch,1,1,1,ik1,npart,mass,pspec,rc,iversion)
Call psdotrkc(2,xmkch,1,1,1,ik2,npart,mass,pspec,rc,iversion)
Call psdotrkc(3,xmkch,1,1,1,ik3,npart,mass,pspec,rc,iversion)
Call psdotrkc(4,xmkch,1,1,1,ik4,npart,mass,pspec,rc,iversion)
If (.not. newsquid(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
Call psinit(energy,npart,nres,rc)
Call pstrkcp(1,0.,1,1,1,ig1,npart,mass,pspec,rc,iversion)
Call psdobmft(2,xmpich,1,1,1,ip1,npart,mass,pspec,rc,bchi1,xi1,
*iversion)
Call psdobmft(3,xmpich,1,1,1,ip2,npart,mass,pspec,rc,bchi2,xi2,
*iversion)
Call pstrkcp(4,0.,1,1,1,ig2,npart,mass,pspec,rc,iversion)
Call pstrkcp(5,0.,1,1,1,ig3,npart,mass,pspec,rc,iversion)
Call pstrkcp(6,0.,1,1,1,ig4,npart,mass,pspec,rc,iversion)
Call pstrkcp(7,0.,1,1,1,ig5,npart,mass,pspec,rc,iversion)
Call psres2(xmpi0,4,5,nres,mres,pmem,mem,rc)
Call psres4(xmeta,2,3,4,5,nres,mres,pmem,mem,rc)
Call psres2(xmeta,6,7,nres,mres,pmem,mem,rc)
If (.not. newsquid(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
C These additional declarations and dimensions are needed:
Integer irc
Real tks1(3), tks2(3), ers1(6), ers2(6)
Real crxy(2), zswum(2), p4k0(4), mpipi, pxyz
Real chi2xy, chi2z, lxy, ct
.
.
C Selection of K0 candidates:
clhb The two pion couple or pion-proton couple have been identified
c before you call klamcorr to find the k0 candidate.
Call Klamcorr (ip1,ip2,xmpich,xmpich,
*irc,tks1,tks2,ers1,ers2,crxy,zswum,p4k0,mpipi,pxyz,chi2xy,
*chi2z,lxy,ct,chisqfit,fitflag,iversion)
clhb-----if fitflag=.true.,then do the kinematical fitting for the
clhb-----found k0(lambda0) vertex by Call FITVEE,the chi square of
clhb-----this fit is chisqfit,otherwise,not do the kinematic fitting
clhb-----and the chisqfit=999. by default.
If (irc .ne. 0) go to 100
C Cut on swum pi pi mass, momentum alignment, etc. if desired:
If (abs(mpipi-xmk0) .gt. 0.100) go to 100
probxy = prob(chi2xy,1)
clhb-------the following condition can vary depending on the different
things you are doing.
If (probxy .lt. 0.01) go to 100
C Kinematic fit:
Call psinit(energy,npart,nres,rc)
Call psextern(1,xmpich,1,1,1,ip1,1,tks1,ers1,
* npart,mass,pspec,rc)
Call psextern(2,xmpich,1,1,1,ip2,1,tks2,ers2,
* npart,mass,pspec,rc)
Call psdobmft(3,xmkch,1,1,1,ik1,npart,mass,pspec,rc,bchik,xik,
*iversion)
Call psdobmft(4,xmmu,1,1,1,imu,npart,mass,pspec,rc,bchimu,
*ximu,iversion)
Call pstrkcp(5,0.,0,0,0,0,npart,mass,pspec,rc,iversion)
if(.not.fitflag)then
Call psres2(xmk0,1,2,nres,mres,pmem,mem,rc,iversion)
endif
If (.not. newsquid(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
Call psinit(energy,npart,nres,rc)
Call psdotrkc(1,xmkch,1,1,1,ik1,npart,mass,pspec,rc,iversion)
Call psdotrkc(2,xmpich,1,1,1,ip1,npart,mass,pspec,rc,iversion)
Call psdotrkc(3,xmkch,1,1,1,ik2,npart,mass,pspec,rc,iversion)
Call psdotrkc(4,xmpich,1,1,1,ip2,npart,mass,pspec,rc,iversion)
Call psres2(-1.,1,2,nres,mres,pmem,mem,rc)
Call psres2(-1.,3,4,nres,mres,pmem,mem,rc)
If (.not. newsquid(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
C These arrays are defined to save the track parameters and
C error matrices.
Real tkpi(3,50), errpi(6,50), tkk(3,50), errk(6,50)
Integer irc, ichg
.
.
C The DEDX and MCS corrections are done once and for all.
C (ichg is just discarded.)
Do 50 i=1,ncharged
If (i is a pion candidate)
* Call Sundae (i,xmpich,irc,ichg,tkpi(1,i),errpi(1,i),iversion)
okpi(i) = (irc .eq. 0)
If (i is a kaon candidate)
* Call Sundae (i,xmkch,irc,ichg,tkk(1,i),errk(1,i),iversion)
okk(i) = (irc .eq. 0)
50 Continue
.
.
C Now we are looping over all pi and K candidates which
C had okpi or okk=true.
Call psinit(energy,npart,nres,rc)
Call psextern(1,xmkch,1,1,1,ik1,1,tkk(1,ik1),errk(1,ik1),
* npart,mass,pspec,rc)
Call psextern(2,xmkch,1,1,1,ik2,1,tkpi(1,ik2),errk(1,ik2),
* npart,mass,pspec,rc)
Call psextern(3,xmpich,1,1,1,ip1,1,tkpi(1,ip1),errpi(1,ip1),
* npart,mass,pspec,rc)
Call pstrkcp(4,-1.,0,0,0,0,npart,mass,pspec,rc,iversion)
Call psres3(-1.,1,2,3,nres,mres,pmem,mem,rc)
If (.not. newsquid(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
Call psinitjr(npart,nres,rc)
Call pstrkcp(1,0.,1,1,1,ig1,npart,mass,pspec,rc,iversion)
Call pstrkcp(2,0.,1,1,1,ig2,npart,mass,pspec,rc,iversion)
Call psres2(xmpi0,1,2,nres,mres,pmem,mem,rc)
IF (.not. squidjr(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
C Alternative: Use the routine Ggfit, which does the
C kinematic fit and returns raw mass, decay angles, etc.
Integer irc
Real rawmgg, rawcos
Real chisq, fitcos, fitpg1(4), fitpg2(4), fitpgg(4)
.
.
Call Ggfit (ig1,ig2,xmpi0,irc,rawmgg,rawcos,
* chisq,fitcos,fitpg1,fitpg2,fitpgg)
C These additional declarations and dimensions are needed:
Integer irc
Real tks1(3), tks2(3), ers1(6), ers2(6)
Real crxy(2), zswum(2), p4k0(4), mpipi, pxyz
Real chi2xy, chi2z, lxy, ct
.
.
C Swim to crossing point, calculate various quantities:
Call Klamcorr (ip1,ip2,xmpich,xmpich,irc,tks1,tks2,ers1,ers2,
*crxy,zswum,p4k0,mpipi,pxyz,chi2xy,chi2z,lxy,ct,chisqfit,
*fitflag,iversion)
Call psinit(energy,npart,nres,rc)
Call pstrkcp(1,0.,1,1,1,ig1,npart,mass,pspec,rc,iversion)
Call pstrkcp(2,0.,1,1,1,ig2,npart,mass,pspec,rc,iversion)
Call psdobmft(3,xmpich,1,1,1,ip1,npart,mass,pspec,rc,bchi1,
*xi1,iversion)
Call psdobmft(4,xmpich,1,1,1,ip2,npart,mass,pspec,rc,bchi2,
*xi2,iversion)
Call psdobmft(5,xmpich,1,1,1,ip3,npart,mass,pspec,rc,bchi3,
*xi3,iversion)
Call psres2(xmeta,1,2,nres,mres,pmem,mem,rc)
Call psres4(0.95757,1,2,3,4,nres,mres,pmem,mem,rc)
Call psres5(1.9714,1,2,3,4,5,nres,mres,pmem,mem,rc)
IF (.not. squidjr(npart,mass,pspec,nres,
* mres,pmem,mem,rc,chisq,pfit)) go to 100
A fit usually involves several iterations, as follows:
Suppose there are Npart particles in our fit. The outputs of the fit are the fit four-momenta of these particles or 4*Npart numbers. (Of course the chi-squared is also output, but that doesn't count here.)
The inputs to the fit are: the measured track momenta (3*Npart numbers), the assumed masses of the particles (Npart numbers), and the four-momentum of the entire event (four numbers; we must assume e+e- annihilation, of course). In addition, we may impose resonance mass constraints which give one input quantity per resonance.
We therefore obtain
Nc = [3*Npart + Npart + 4 + Nres] - [4*Npart],
or
Nc = 4 + Nres,
where Nres is the number of resonance constraints.You may just think that for an e+e- events, there is only one energy conservation law exists (so we have 4 constraints), if we have Nres independent resonances among these final state, then we come to
Nc = 4 + Nres,
The following complications must also be considered:
(0,0,0,sqrt_s). The complete event topology hypothesis must
be specified, allowing at most one missing particle.Squidjr is used to impose resonance (mass) constraints on some subset of the tracks in an event.
Performing a kinematic fit with Telesis requires a sequence of subroutine calls in which all of the particles and resonances are specified. The sequence of calls FOR EACH FIT must have the following elements:
sqrt_s. The other arrays and variables are
being initialized.
Routines to specify a Particle:
iversion.le.7)
Also use this routine for missing particles.
Subroutine pstrkcp (i_squid,mass_i,spec1,spec2,spec3,spec4,
n_part,mass,p_spec,rc,iversion)
Subroutine psdotrkc (i_squid,mass_i,spec1,spec2,spec3,spec4,
n_part,mass,p_spec,rc,iversion)
Subroutine psdobmft (i_squid,mass_i,spec1,spec2,spec3,spec4,
n_part,mass,p_spec,rc,bchisq,xi0,iversion)
Subroutine psextern (i_squid,mass_i,spec1,spec2,spec3,spec4,
type_i,tkp,err,n_part,mass,p_spec,rc)
Notice that the value of spec4 (the Trklst track number)
is NOT used by
psextern. The value is loaded as usual
into p_spec but does not affect the fit in any way.
Particles may be specified in any order, and you can use any
combination of calls from the above list in the same fit. For
each of these routines, the first parameter (i_squid) gives
the SQUID particle number (as opposed to the Trklst track
number). The first particle to be specified must have
i_squid=1, the second i_squid=2,
and so on. [If any resonances are to
be specified, the particle number assignments can then be read
off instantly.]
Routines to specify a Resonance:
For each resonance, choose the routine which has the
correct number of daughter particles:
Subroutine psres2 (m_res_i,i1,i2,
n_res,m_res,p_mem,mem,rc)
Subroutine psres3 (m_res_i,i1,i2,i3,
n_res,m_res,p_mem,mem,rc)
Subroutine psres4 (m_res_i,i1,i2,i3,i4,
n_res,m_res,p_mem,mem,rc)
Subroutine psres5 (m_res_i,i1,i2,i3,i4,i5,
n_res,m_res,p_mem,mem,rc)
Subroutine psres6 (m_res_i,i1,i2,i3,i4,i5,i6,
n_res,m_res,p_mem,mem,rc)
Subroutine psres7 (m_res_i,i1,i2,i3,i4,i5,i6,i7,
n_res,m_res,p_mem,mem,rc)
Subroutine psres8 (m_res_i,i1,i2,i3,i4,i5,i6,i7,i8,
n_res,m_res,p_mem,mem,rc)
Subroutine psres9 (m_res_i,i1,i2,i3,i4,i5,i6,i7,i8,i9,
n_res,m_res,p_mem,mem,rc)
Subroutine psres10(m_res_i,i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,
n_res,m_res,p_mem,mem,rc)
Resonance and particle calls may be interspersed at will.
You can even define resonances using particles that are defined
later.The routines mentioned so far are sometimes known as the 'pre-squid' routines (hence the initials 'Ps').
Function Newsquid (n_part,mass,p_spec,n_res,
m_res,p_mem,mem,rc,chisq,p_fit)
Function Squidjr (n_part,mass,p_spec,n_res,
m_res,p_mem,mem,rc,chisq,p_fit)
by referring to Newsquid or Squidjr in a logical expression
such as
If (.not. newsquid(...)) go to 100
or
goodfit = newsquid(...) .
Newsquid or Squidjr = .TRUE. if the fit was successful. The
other outputs are rc, chisq, and p_fit. From the value of rc
one can determine the reasons for fit failure.
Subroutine squidpul (p_spec,rc,pull)
The first 42 words of Trklst track subtype 3 contain the UNCORRECTED track momentum and an error matrix which is computed assuming the particle is a pion. The first six elements of this matrix are reasonable under this assumption, but the last nine elements are not correctly calculated. Since the kinematic fitters in BES only use the first six elements, good results are usually obtained by simply correcting the magnitude of the track momentum for DEDX before performing fits. The chi-squared distributions are incorrect when any charged tracks other than pions are involved, but the fit momenta come out pretty good. [Notice that track subtype 8 vertex-constrained fits are performed assuming all tracks are pions, and using the incorrect track error matrix elements which pertain to xi and eta.]
delta(P) = A * (P/Pxy) / beta**u, (beta<0.94)
delta(P) = A * (P/Pxy) / 0.94**u, (beta>=0.94)
where A is determined from the amount of matter traversed (a
radial path from the beampipe through Layer 5 was considered) and
u=2.74, as determined from a plot of the more complicated
phenomenological function. As mentioned in the memo, slow tracks
could be more accurately treated by dividing the material into
several parts and accounting for the changing beta of the
particle as the material is traversed.
The Telesis routine superice recovers(approximately)the geometrical error matrix from the trklst track subtype 3 error matrix by removing the 'corrections' applied by SCATR and PQSCAT (at BES7).Because some intermediate results are not available when reading back the reconstructed event, only an approximate inversion is possible. (It is good enough, though.)
In Sundae the mean momentum loss is calculated by integrating the Hitlin formula and referring to a look-up table, yielding good accuracy and considerable speed. (For beta>=0.94 the Hitlin formula is retained.) Let X (running from 0 to 1) denote the fractional amount of material traversed. We have (for beta<0.94)
dP = dX * A * (P/Pxy) / beta**u,
or
beta**u
------------------ * d(beta) = dX * A * (P/Pxy) / M,
(1 - beta**2)**1.5
where M is the hypothesized mass of the particle. The left-hand
side of this equation is numerically integrated and saved in the
array F1. The index of this array is 1000*beta at the upper
limit of integration. (The lower limit is zero.) The material in
the detector is divided into four parts (to accomodate the
multiple scattering calculations described below), so the A is
replaced with an Ai for each part. To determine the momentum
correction for a track passing through one of the scatterers, the
quantity Ai*(P/Pxy)/M is calculated and the value of the integral
corresponding to the exit value of beta is interpolated
(linearly) from the array F1. The sum of these quantities is the
value of the integral corres- ponding to the entrance value of
beta. This beta is determined from F1 and the new momentum is
computed from beta.the track phi correction is
d(phi) = CHG * PATHXY * 0.03 * B * (1/PxyNEW - 1/PxyOLD),
sigma(k) = 2*alpha*t*sec(lambda) / beta**3 / P**2,
with sigma(phi), and the error matrix contributions as before.
(14.1 MeV)**2
d(DELTH2) = ------------- * RADLi * (P/Pxy) * dX,
(beta*P)**2
where RADLi is the number of radiation lengths of material in a
radial path through the scatterer. Inserting d(beta)/dX from
above, we obtain
(14.1 MeV)**2 * RADLi d(beta)
d(DELTH2) = --------------------- * -----------------------------
M * Ai beta**(4-u) * SQRT(1-beta**2)
The integral of the second factor on the right-hand side is
evaluated up to various values of beta and tabulated in the array
F2. Since the integral diverges when the interval includes beta
= 0, I "renormalize" by setting the integrand to zero at this
point. (Adding a constant to the values in F2 does not affect the
subsequent calculations.) To determine DELTH2 for the passage of
a track through scatterer i, we look up the values of F2
corresponding to the entrance and exit betas determined from the
DEDX section. (Again the array index is 1000*beta.) For
beta>0.94 we are safe to treat beta and P as constants in the
first equation of this section. From DELTH2 the error matrix
adjustments are found as before.
Important beam fit outputs:
These quantities provide useful track selection criteria. To select tracks which are consistent with coming from the interaction region, a cut of Prob(chisq,1) > 0.001 or 0.01 is very effective. You may also wish to apply a very loose cut on xi0 (< 0.025 meters, for example) to eliminate tracks which have small bchisq because they have enormous errors on impact distance. Naturally it is a mistake to make such cuts or to use beam fit momenta if you are looking for the pions coming from K0 decays. Beam fits may be useful to reject some K0 daughters, however.
Beam fit track information is easy to use in kinematic fits under Telesis.
To perform beam fits for use outside of kinematic fits, the routine banaspi is usually most convenient. The user specifies only the Trklst track number and the desired particle mass. The outputs include the fit track parameters, xi0, and bchisq.
____________________________________________________________________________
NEWSQUID | SQUID
A.In general
A.0.NEWSQUID and SQUID are the same in |
logical, but NEWSQUID has a more variable |
'rc' to charactrize the return mode: |
OKSQ=NEWSQUID(npart,mass,pspec,nres,mres | OKSQ=SQUID(npart,mass,pspec,nres,
,pmem,mem,rc,sqchi,pfit) | mres,pmem,mem,sqchi,pfit)
|
A.1 One more equal mass constraints | Not available in the old code
can be treated: |
Assuming there is four(real or miss |
ing) particles in the final state,we |
can impose |
M(X1X2)=M(X3X4) or M(X1)=M(X2X3X4) |
A.2.It can deal with vary kinds of track | It can only treat the track from
with DEDX,MCS and BMFIT correction | Track list,and do a mean DEDX
or Track from K0 or Lambda directly. | correction for the chrged track
It does not do any DEDX correction | in subroutine kfinit for the
in subroutine kfinitl!!! | amendment.
__________________________________________|_____________________________________
B. Method of INPUT and OUTPUT
B.1.INPUT
Assuming there are Npart(real and missing) particles in the final state.
the first two particles are pion and kaon,and the final particle is kaon.
NEWSQUID SQUID
o Energy of Center Mass System(can vary) | Etot=Certain value specified in
Call Psinit(Etot,npart,nres,rc) | Kfinit
*:you just give the value to Etot(4.03, |
or 3.097 or 3.77),the other variable |
s are output |
|
o Array Pspec(4,npart) and Mass(npart) | DO i=1,npart
Call pstrkcp(1,0.1396,1,1,1,itrk1, | Do J=1,3
npart,mass,pspec,rc,iversion) | Pspec(j,i)=1
Call pstrkcp(2,0.49367,1,1,1,itrk2, | enddo
npart,mass,pspec,rc,iversion) | enddo
... | Pspec(4,1)=itrk1
... | pspec(4,2)=itrk2
Call pstrkcp(npart,0.49367,1,1,1, | ...
itrknpart,npart,mass,pspec,rc,iversion) | ...
*:where itrk1,itrk2 and itrknpart are | Pspec(4,npart)=itrknpart
the track number of the first,second and |
the last track. | Mass(1)=0.1396
*:iversion is the reconstruction version | Mass(2)=0.49367
it can be 5,6,7 or 8(now the highest | ...
version) | ...
*:the other variables are output | Mass(npart)=0.49367
*:pstrkcp just take the track paramete |
rs from track list,it can be substituted |
by: |
1.psdotrkc(...,iversion):Doing DEDX |
and MCS correction for the tracks. |
2.psdobmft(...,bchk1,xik1,iversion): |
Doing DEDX,MCS and BMFIT for the track. |
the '...' are the same as that in |
pstrkcp. |
the bchk1 characterize the BMFIT chi |
square,and Xik1 is just the closet |
distance of the track to the I.R. point |
3.psextern(1,0.1396,1,1,1,itrk1,ktype, |
tkp,err,npart,mass,pspec,rc): |
Any tracks that are modified or swimmed |
by any means. |
*:ktype has the same meannig as that in |
SQUID(NEWSQUID),the origin of the track |
MDC,ESC or WSC |
__________________________________________|___________________________________
o Resonance array Mres(Nres),Pmem(...), Mem(pmem(...)+...),assuming we
consider the decay e+e- -----> gamma eta(gammagamma)eta(gammagamma):
|
NEWSQUID | SQUID
Call psres2(0.54745,1,2,Nres,Mres, | Nres=2
pmem,mem,rc) | Mres(1)=0.54745
Call psres2(0.54745,3,4,Nres,Mres, | Mres(2)=0.54745
pmem,mem,rc) | Pmem(1)=1
*:the Nres,mres,pmem,mem and rc are |
output | Mem(pmem(1)+0)=2
*:psresN means there are N partic- | Mem(pmem(1)+1)=1
les in this resonance,the telesis can | Mem(pmem(1)+2)=2
deal with N up to 10. | Pmem(2)=4
| mem(pmem(2)+0)=2
| mem(pmem(2)+1)=3
| mem(pmem(2)+2)=4
__________________________________________|______________________________________
o equal mass constraints.
NEWSQUID | SQUID
1.for the above example,we can require |
M(gamma1gamma2)=M(gamma3gammma4): |
Call psres2(-1.,1,2,nres,mres,pmem | It could not deal with
,mem,rc) | this kind of constraint.
Call psres2(-1.,3,4,nres,mres,pmem |
,mem,rc) |
2.for the decay e+e----kkpiDs we can |
require(assuming the missing Ds is the |
fourth particle) |
M(KKPI)=M(Ds) |
Call pstrkcp(4,-1.,0,0,0,npart,mass, |
pspec,rc) |
Call psres3(-1.,1,2,3,nres,mres, |
pmem,mem,rc) |
__________________________________________|_______________________________________
B.2.OUTPUT
o Return code
NEWSQUID | SQUID
In every steps,there is the rc parameter | You can only know whether
,so it can tell you whether that step is | the fitting is ok by seeing the
ok. | chi square in the final step.
__________________________________________|_______________________________________
B.3.CHECKING
NEWSQUID SQUID
You can | Not available
Call squidpul(pspec,rc,pull) |
to check the goodness of the fit |
*:a reasonable fitting results should |
give the pull distribution of the momen- |
tum not far from a Gaussian distribution |
with the 0 mean value and 1 width |
__________________________________________|_______________________________________