Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the minimal supersymmetric standard model ( mssm ) is one of the most important candidates for the physics beyond the standard model . the strongly interacting superpartners , gluinos and squarks with mass lighter than 2.5 tev can be discovered at the large hadron collider at cern ( lhc )  @xcite . the physics run of lhc will start in 2008 .    among the supersymmetric ( susy ) particles , the lightest susy particle ( lsp ) plays a key role . in cosmology , it is a candidate for the cold dark matter in the universe . at the lhc , the signals of the supersymmetric particles depends on the nature of the lsp . it may be the lightest neutralino @xmath4 , which escapes from detection and leads to missing @xmath5 in the event . another possibility for the lsp is the gravitino @xmath1 , the superpartner of the graviton . the gravitino coupling to other particles in the mssm sector is extremely small , namely proportional to @xmath6 . the next lightest susy particle ( nlsp ) , whose decay products necessarily include the gravitino , is therefore long - lived . we discuss the scenario where the lsp is the gravitino and the long - lived nlsp carries charge . a natural candidate for a charged nlsp ( cnlsp ) is the lightest scalar tau , @xmath7 , which can be significantly lighter than the other sleptons due to the left - right mixing induced by the off diagonal matrix element @xmath8 . the charged particle leaves tracks in the central detectors ( atlas and cms ) , giving additional information for the susy particle reconstructions . if the nlsp decays in the main detectors , a displaced vertex may be observed as well . the expected lifetime of the cnlsp @xmath9 is unconstrained , because it is proportional to @xmath10 the yet unknown gravitino mass squared . on the other hand , the gravitino mass is proportional to the total susy breaking scale in the hidden sector , therefore the determination of the lifetime is an important physics goal . the lifetime measurement gives direct information on the hidden sector . a particle decays efficiently in the main detectors ( cms and atlas ) if the decay length is sufficiently short , @xmath11 ( 10m)@xmath12 , where @xmath13 is the number of produced susy particles . on the other hand , for typical susy production cross sections , a direct observation of the decay is very difficult for @xmath14 sec . however , it has been pointed out that the cnlsp stopped by the ionization loss in the material may be studied in detail  @xcite . an idea for a stopper based on a water tank is presented in ref.@xcite , where the water can be transported away from the detector site , for concentration and further study . in ref . @xcite , a detector consisting of a tracker and heavy stopping material is proposed , which can measure the arrival time and the location where the cnlsps are stopped , in addition to the energy of the decay products . another possibility , which requires minimal experimental modification , is to study the cnlsp which are stopped in the main detector or surrounding rock  @xcite .    as pointed out in ref . @xcite , the study of the cnlsp decays can probe the underlying supergravity in nature . with the gravitino mass inferred from the kinematics , the additional cnlsp lifetime measurement will test an unequivocal prediction of supergravity . moreover , the study of a rare 3body decay @xmath15 can reveal the peculiar couplings of the gravitino and the gravitino spin @xmath16 .    in this paper , we consider the physics of the cnlsp @xmath9 decays that can be done with the heavy material stopper  detector  @xcite , because only a detector of this type can cover a wide range of lifetime @xmath17 @xmath18 years ) . we find that the mass of the gravitino can be measured if it is sufficiently heavy ( roughly @xmath19 ) . in that case , one can check if the lifetime is consistent with the supergravity interpretation .    as we shall see , the lsp mass resolution is however poor if @xmath20 . in that case , it is very hard to prove that the decay @xmath21 plus an invisible particle @xmath22 is indeed caused by the supergravity interaction involving @xmath23 . we should note that , because of the extremely weak coupling of the gravitino , if there is any unknown supersymmetric particle @xmath22 which couples rather weakly to the mssm particles , the lightest mssm particle may decay into @xmath22 instead of @xmath1 even if the @xmath1 is the lsp . for example , the superpartner of the axion , the axino @xmath2 , can be such a particle . however , we found that it may be possible to discriminate the case of @xmath24 and the case of @xmath25 , by investigating the three body decay @xmath26 as suggested in @xcite , with enough statistics .    this paper is organized as follows . in section [ sec : detector ] , we discuss a possible detector setup . we found an @xmath27 kton detector ( up to 8 kton ) may be placed next to the cms detector without serious modification of the cms experiment itself , but with non - negligible modifications to the cms cavern side walls . in section [ sec : modelpoints ] , we select several model points and estimate the expected number of the stopped particles for @xmath28 fb@xmath29 and 3000 fb@xmath29 . in section [ sec : twobody ] , we discuss the measurement of the two body decay @xmath30 . the three body decay of the cnlsp is studied in section [ sec : three - body ] . section [ sec : discussion ] contains the discussion and comments . 44mcms & 15 m & 12.5kt & 21 m    [ atlascms ]    in this section , we discuss the possibility to install massive stoppers next to the lhc detectors , cms and atlas . it turns out that the cms cavern may allow for an easier installation and more room for a massive stopper , compared to the atlas cavern . the parameters of the two detectors are listed in table . [ atlascms ] . the diameter of the cms detector is smaller , therefore the massive stopper can be placed closer to the interaction point at cms . the weight of the cms detector is about twice as large as that of atlas . because of the large weight , the cavern of cms is designed so that it can sustain a massive object safely , which includes a reinforced floor to spread the detector pressure more equally . a potential massive stopper with a weight of around a few kton can be placed on both sides of the cms detector , but it will need a reorganization of the scaffolding and gallery paths on the cavern walls , to make room for such an additional detector .    the assembly and construction of the two detectors is also very different . most of the cms detector components are assembled on the ground , and about 15 large detector units will be lowered in the cavern for final assembly works . hence the installation is less integrated with the cavern , leaving relatively more freedom and thus changes needed in the cavern to install massive stoppers are somewhat simpler  @xcite . on the other hand atlas detector is assembled mostly in its cavern . the huge magnets toroids and outside muon system fill up the cavern . the cryogenic system in the atlas cavern is also taking space outside the detector  @xcite . is the collision point . right : two stopper  detectors and a circle about the size of cms detector are superimposed on the cross section of cms cavern uxc 55 , drawing taken from ref . @xcite . , title="fig:",width=264 ]    is the collision point . right : two stopper  detectors and a circle about the size of cms detector are superimposed on the cross section of cms cavern uxc 55 , drawing taken from ref . , title="fig:",width=264 ]    we assume two stoppers with the size 3.5 m @xmath31 15 m @xmath31 15 m and the average density @xmath32g/@xmath33 , hence , the total weight of the detector is 8 kton . this is maximum possible rectangular parallelepiped volume that can be placed in the cavern , with the given space to the cavern wall and with its long edge being the same as that of the cms barrel part ( see fig . [ fig : stopper ] ) . the stoppers are thus placed 8.5  m away from the interaction point . we also assume that stopping power of the cms detector is equivalent to @xmath34 g/@xmath35 iron , where @xmath36 is angle between the cnlsp direction and the beam direction . the number comes from the average density of cms detector , 3.37g/@xmath33 , which leads to the weight per @xmath35 for the radial direction of @xmath37 .    as discussed in the previous paper  @xcite , the stopper can be a hadronic and electromagnetic calorimeter simultaneously , if the detector consists of layers of dense stopper and tracking devices . the measurement of the energy of the decay product of the cnlsp is the key ingredient to explore the cnlsp interactions to the @xmath22 particle . in this paper we assume that the cnlsp is the scalar tau lepton @xmath0 , which decays mostly as @xmath38 where @xmath39 or @xmath2 . the @xmath40 decays into @xmath41 , or into @xmath42 and @xmath43 s . we do not consider the decays into @xmath44 , because the muon energy can not be measured unless the stopper contains a magnetic field . the energy of the leptons are much softer than the parent @xmath40 energy anyway , so that they are less useful for the study of the decay kinematics . a large volume detector is advantageous to measure the energy of the @xmath40 decay products , because the detector must contain most of the energy of the showers from the @xmath40 decay products . to fully absorb the hadronic cascade one needs sufficient thickness of the calorimeter . the required thickness for an iron calorimeter ( density is 7.87g/@xmath33 ) is listed in ref . @xcite , and is about 170 ( 120 ) cm for 100  gev single hadron energy for 99% ( 95% ) containment respectively , equivalent to 1337.9g/@xmath35 ( 944.4g/@xmath35 ) . to measure the energy with sufficient accuracy by a stopper with @xmath45 5g/@xmath33 , the distance between the decay position and the end of the detector along the shower path must be at least 190 cm . a simplified monte carlo simulation shows that 63% of the @xmath9 decays satisfy this condition . high energy photons and electrons initiate electromagnetic cascades , which are much easier to contain in a detector . the energy deposition along the axis of the cascade is well described by a gamma distribution , where the maximum occurs at @xmath46 , where @xmath47 is the radiation length . the energy deposition terminates around @xmath48 . ( see fig.27.17 of @xcite . ) for the case of iron ( @xmath49 ) , 200g/@xmath35 thickness of the material is needed to stop the electromagnetic cascade . the difference in the absorption length is used to discriminate isolated photons from hadron , which would be useful to study @xmath50 decay .    in a previous paper  @xcite , we discussed the possibility to re - use of the existing 1 kton detector , such as soudan ii  @xcite as the cnlsp stopper . the soudan ii consists of the layers of @xmath51(m ) long drift tubes and thin iron plates . the physics goal of this detector is the search for proton decays . to be sensitive to the low energies involved in these decay processes , the average density of soudan ii is low , less than 2g/@xmath33 . the size of the detector is probably enough to stop a certain amount of cnlsps and measure the decay rate , and therefore it may well be appropriate for the a first stage of the cnlsp study . on the other hand , it is certainly not enough for a detailed energy measurement of the @xmath0 decay product . most of the hadronic decay cascades are not fully contained for a detector with the geometry allowed by the space in the cms cavern . a high density detector consisting of the layers of drift tubes / scintillators or rpcs to measure the charged particles between iron plates thicker than soudan ii will be more optimized for the cnlsp study . in this paper we assume a conservative energy resolution for hadronic showers , which is around @xmath52 . the value is not unrealistic for a simple massive and affordable detector , if a shower is sufficiently contained in the stopper . an additional complexity may occur for showers which develop parallel to the layers of the tracking devices . if the particles pass mostly through tracking devices , they feel a much lower average density , while the particles going mostly through the iron plates do not give detectable signature efficiently . we assume that the measured energy will be corrected depending on the shower directions . a detector uniform to all directions would be better to measure the energy of decay products of the stopped cnlsp . the calorimetery technologies studied for the ilc may satisfy such conditions ; see @xcite for the cnlsp study at ilc . in this section , we briefly describe supersymmetric models with a long - lived charged next lightest susy particle . we also select some model points , and estimate at each model point the number of cnlsp that can be stopped at the stopper  detector .    in the minimal supergravity ( msugra ) models , the scalar masses and gaugino masses , trilinear couplings are universal at the gut scale @xmath53 which are denoted by @xmath54 , @xmath55 , and @xmath56 respectively . the resulting mass spectrum is obtained by solving renormalization group equations ( rge s ) . when the rge is integrated up to @xmath57 , slepton soft masses and gaugino masses are approximately given by the following convenient formulas @xmath58 where @xmath59 is the gauge couplings at the unification scale.(30% ) if the mass scale is @xmath51(1tev ) . ] in this model , the mass of gravitino @xmath60 is order of @xmath54 or @xmath55 , because both of them are proportional to the @xmath61 with @xmath27 coefficient , where @xmath62 is the fundamental susy breaking scale . depending on the @xmath27 coefficients , the gravitino can be the lightest superparticle in general . also in the gaugino mediation models  @xcite , the gravitino can be the lsp in a large domain of the parameter space  @xcite . in this case the scalar masses are very small at the boundary , i.e. , @xmath63 , and the stau naturally becomes the nlsp .    in the gauge mediation ( gm ) model  @xcite the supersymmetry is broken at lower energy scale , and the susy breaking in the hidden sector is mediated to the mssm sector by gauge interactions . the simplest gm model is described by 6 parameters : @xmath64 , @xmath65 , @xmath66 , @xmath67 , sgn@xmath68 and gravitino mass . here @xmath66 is an effective number of vector - like heavy quarks and leptons in su(5 ) representations , @xmath69 and @xmath70 . ( @xmath71 for @xmath72 quarks and leptons , and @xmath73 for @xmath74 . ) a messenger field @xmath75 couples to the vector - like fields as @xmath76 and develops a vacuum expectation value @xmath77 . the gravitino mass is given by @xmath78 where @xmath79 is the total susy breaking scale of the theory . we take @xmath60 as a free parameter in this model . the masses of mssm particles at @xmath65 are obtained by relatively simple formula . gaugino masses satisfy the gut relation , and @xmath80 squark and slepton masses are given by @xmath81 n_5 \lambda^2 , \cr m^2_{\tilde{u}_r } & = & \left [ { 8\over 3 } \left({\alpha_s\over 4\pi}\right)^2 + { 8\over 15 } \left({\alpha_1\over 4\pi}\right)^2 \right ] n_5 \lambda^2 , \cr m^2_{\tilde{d}_r } & = & \left [ { 8\over 3 } \left({\alpha_s\over 4\pi}\right)^2 + { 2\over 15 } \left({\alpha_1\over 4\pi}\right)^2 \right ] n_5 \lambda^2 , \cr m^2_{\tilde{\ell}_l } & = & \left [ { 3\over 2 } \left({\alpha_2\over 4\pi}\right)^2 + { 3\over 10 } \left({\alpha_1\over 4\pi}\right)^2 \right ] n_5 \lambda^2 , \cr m^2_{\tilde{e}_r } & = & \left [ { 6\over 5 } \left({\alpha_1\over 4\pi}\right)^2 \right ] n_5 \lambda^2 . \label{gmeq}\end{aligned}\ ] ]    all of the above models predict a large mass hierarchy between strongly interacting superpartners and weakly interacting superpartners . heavy gluino and squarks are copiously produced at the lhc , and they decay into the light weakly interacting susy particles . the lightest susy particle in the mssm sector is either the lighter stau @xmath7 or the lightest neutralino @xmath4 . the mass of the @xmath7 is the smaller eigenvalue of the mass matrix , @xmath82 because of the off - diagonal elements of @xmath0 mass matrix , @xmath7 could be significantly lighter than the other sleptons . if the @xmath7 is the nlsp , the stopper  detector is useful to stop it and to study its decay . we therefore consider the phenomenology when @xmath83 in this paper . in the following , we omit the subscript of @xmath7 and denote the nlsp stau as @xmath0 .    in most part of this paper , we will discuss the cnlsp physics as model independent as possible . however , to give some numerical reference , we choose a few model points . for msugra models , we take the points proposed in @xcite . the parameters are listed in table . [ masssugra ] and mass spectrum is given by the isajet ver . 7.69 . in @xcite , the gravitino mass is taken to be @xmath84 in those model points . @xmath85 @xmath55 & 440 & 1000 & 1000 @xmath54 & 20 & 100 & 20 @xmath67 & 15 & 21.5 & 23.7@xmath86[gev ] & 1025 & 2191 & 2190 @xmath87[gev ] & 175 & 417 & 416@xmath88[gev ] & 154.2 & 343.5 & 324.3 @xmath89(susy)[pb ] & @xmath90 & @xmath91 & @xmath92 stopped in the stopper  detector per @xmath93 events & 255 & 250 & 254 stopped for 300 fb@xmath29 & 4636 & 34 & 36     @xmath86[tev ] & 0.93 & 1.13 & 1.34 & 1.54 & 1.74 @xmath87[gev ] & 161.7 & 205.3 & 248.7 & 292.1 & 335.4@xmath88[gev ] & 120.5 & 150.1 & 179.9 & 209.8 & 239.9@xmath89(susy)[pb ] & 5.24 & 1.68 & 0.64 & 0.28 & 0.13 stopped in the stopper  detector per @xmath93 events & 282&274 & 274 & 294 & 302stopped for 300fb@xmath29 & 8830 & 2762 & 1052 & 494 & 236    for the gauge mediation model , we take the model points similar to that for the study in  @xcite . namely , we fix @xmath73 , @xmath94 and @xmath95 , where @xmath96  tev . the mass spectrum and production cross sections are summarized in table . the gravitino mass for these model can be very small , whose minimum value is given by @xmath97 . if the messenger sector couples a fraction of total susy breaking , i.e. @xmath98 , the gravitino mass , which is proportional total susy breaking @xmath62 , can be large . the pattern of the mass spectra are similar for both msugra and gm sample points . the gaugino masses obey the gut relation , and @xmath0 mass is lighter than @xmath99 mass because of the cnlsp assumption . for msugra , this implies @xmath100 , so that @xmath101 mass is about as large as @xmath102 , as it is so in gm models . the relation @xmath103 is realized in our model points .    to estimate the number of stopped cnlsp , the production cross section of the susy particles and the velocity distribution of @xmath0 at lhc must be evaluated . they depend on the @xmath102 and @xmath101 masses , and the mass difference @xmath104 . we estimate the production cross section and the decay distribution by using herwig@xcite , where the mass spectrum and branching ratios are interfaced from isajet@xcite to herwig by using isawig@xcite . we generated @xmath93 susy events for each model point . the flying range @xmath105 of the charged stable massive particle may be calculated by integrating the energy loss equation of heavy ionizing particle ( bethe - bloch equation ) . the result is a function of @xmath106 , with a linear dependence on the mass of the particle @xmath107 . in this paper , the stopping power of the stopper  detector is calculated using the data in @xcite . see also the detailed discussion in ref . @xcite . the maximum length of a particle track through the stopper  detector @xmath108(max ) , which depends on cnlsp direction , is calculated assuming that the track of the cnlsp goes straight from the production point . we regard the particle is stopped in the stopper  detector if the flying range @xmath105 satisfies : @xmath109 where @xmath110g@xmath111@xmath33 is the density of the stopper  detector . the number of stopped cnlsp in a stopper per @xmath93 events , and the number of stopped cnlsp for @xmath112 in the two stoppers , are listed in table [ masssugra ] and [ mass ] . in fig . [ fig : position ] we also show the simulated stopped positions in a stopper  detector . the position distribution is uniform in the detector . tev gm point . here the @xmath113 axis is the beam direction , and the @xmath114 axis is the vertical direction . the origin @xmath115 is the interaction point and we assume the stopper  detector is at 8.5m@xmath11612 m , @xmath117m@xmath1187.5 m , and @xmath117m@xmath1197.5 m ( cf . fig . [ fig : stopper ] ) . the big square is the projection on the @xmath114-@xmath113 plane , the top rectangle is the projection on the @xmath120-@xmath113 plane , and the right rectangle is the projection on the @xmath120-@xmath114 plane . , width=453 ]     distribution of @xmath0 for @xmath121   tev and @xmath122  tev in gm models and for point @xmath123 in msugra models for @xmath93 susy events . light ( dark , gray ) shaded histograms are the number of stopped events in the stopper scaled by factor of 10 , for point @xmath123 and for gm models with @xmath124 tev respectively . , width=264 ]    in fig . [ fig : betagamma ] , we show the @xmath125 distribution of the cnlsps for a few model points . we find that cnlsp tends to be less relativistic when the susy scale is large , because momentum of parent squarks and gluino is reduced . in fig . [ fig : betagamma ] , the peak position of the @xmath125 is at @xmath126 for @xmath127  tev . because of that , the number of cnlsp in the smaller @xmath125 region is increased as the susy scale is increased . on the other hand , as @xmath128 increases , cnlsps with smaller @xmath125 are stopped in the detector while the number of events in the smaller @xmath125 bins are kinetically suppressed . for instance , for @xmath122 and 80  tev , cnlsps in the bins between @xmath129 and @xmath130 are stopped in the stopper  detector , respectively . altogether , the number of stopped @xmath0 for @xmath131 to @xmath132  tev for @xmath93 susy events is roughly constant as we can see in table . [ mass ] . the production cross section reduces when the gluino mass is increased because the parton distribution functions of gluon and quarks are quite small for @xmath133 . for 300fb@xmath29 , the number of events stopped at the assumed two 4 kton stoppers ranges from 8000 events to 30 events in the table . we will see in the next section that accumulation of @xmath134 cnlsps are enough to measure the lifetime with @xmath51(10)% accuracy . we will also estimate the resolution of @xmath40 lepton energy @xmath135 arising from the decay @xmath136 through the end point of the tau jet energy . statistically @xmath137 stopped cnlsps are enough to measure the end point with a few gev error . in the upgrade of lhc ( slhc ) , integrated luminosity of 3000fb@xmath29 is proposed , therefore the number of stopped events ranges from @xmath138 to @xmath139 for the model points presented in table  [ mass ] . in this section we study the two body decay of the cnlsp in the stopper  detector . both in msugra and gm models , the stau can be the cnlsp and decays into the gravitino @xmath1 and the @xmath40lepton . the cnlsp decay width into a gravitino and a lepton is given by  @xcite @xmath140^{3/2}\;.          \nonumber\\          & = &          ( 68~\mathrm{days})^{-1 }          \left(\frac{{m_{\tilde{\tau}}}}{100~\mathrm{gev}}\right)^5          \left(\frac{10~\mathrm{gev}}{{m_{\tilde{g}}}}\right)^2          \times           \nonumber\\          & &          \left(1-          \frac{{m_{\tilde{g}}}^2 + m_\tau^2 }           { { m_{\tilde{\tau}}}^2 }           \right)^4           \left [             1-\frac{4 { m_{\tilde{g}}}^2 m_\tau^2 }             { ( { m_{\tilde{\tau}}}^2-{m_{\tilde{g}}}^2-m_\tau^2)^2 }             \right]^{3/2}\;.           \label{eq : gamma32}\end{aligned}\ ] ] we show the dependence of the stau lifetime on the gravitino mass in fig . [ fig : lifetime ] .     in the case of gravitino lsp , as a function of the gravitino mass , for the stau mass 100 , 200 , and 400  gev ( from top to bottom ) . ]    in general , the two  body decay @xmath141 can be triggered by a single tau jets initiating from the position where @xmath0 is stopped . the tau energy is monochromatic and expressed as @xmath142 here @xmath22 is the invisible particle in the @xmath0 decay , in this case @xmath39 . if one can measure both of the lifetime and the mass of the stau , the gravitino mass can be determined _ assuming _ that eq.([eq : gamma32 ] ) is correct . then the total susy breaking scale @xmath143 is also determined , which is very important to understand the hidden sector physics . the lhc main detectors can determine the mass of @xmath0 through the stau velocity measurement @xmath144 , e.g. , in the muon system of the cms detector . however , the measurement of the lifetime may not be easy at the main detectors if the lifetime is too much longer than the detector size . the cross section is typically @xmath27 pb or less , therefore we have at most @xmath145 event at hand for @xmath146 . thus , only 10 events or less decay inside the detector if @xmath147 km . some of the cnlsp are stopped in the main detector , but measuring the decay precisely in the main detector would be challenging during the beam time . on the other hand a stopper  detector  @xcite can measure the position and the time where a cnlsp is stopped , and its decays without dead - time . the lifetime measurement will be discussed in sec.[subsec : lifetime ] . see @xcite also on the idea to measure the lifetime by triggering muons from the decays of the cnlsp stopped in the surrounding rock . to _ predict _ the cnlsp lifetime , one has to determine the gravitino mass independently . this is possible through the extraction of @xmath135 from the energy distribution of the tau jet from the cnlsp decay , because @xmath135 is a function of @xmath148 and @xmath149 , as can be seen in eq . ( [ etau ] ) . @xmath149 is expressed as a function of @xmath135 as follows ; @xmath150 fig . [ fig : etau ] shows the dependence of the reconstructed lsp mass eq.([eq : mx ] ) on the tau energy @xmath151 for several values of @xmath148 . as can be seen from the figure , it is crucial to measure the @xmath151 as precise as possible , especially for small @xmath149 , in order to determine the mass @xmath149 . hence , the stopper  detector should offer a reliable measurement of @xmath135 . we will discuss the @xmath152 measurement and the @xmath149 reconstruction in sec . [ subsec : etau ] and sec . [ subsec : mx ] , respectively . [ see eq.([eq : mx ] ) ] on the tau energy @xmath151 , for @xmath153 , 150 , 200 , 250 , and 300  gev ( solid lines ) . dashed lines show contours of the lifetime in the case of the gravitino lsp ( @xmath24 ) . ] one can study the supersymmetric version of gravity interaction by studying the consistency between the observation and prediction of the decay rate . if the gravitino is the lsp , the decay width is given by eq.([eq : gamma32 ] ) . now if one can independently determine the gravitino mass by means of kinematics , as described above , eq.([eq : gamma32 ] ) can be used in the other way round , which leads to the measurement of the supergravity planck scale  @xcite @xmath154^{3/2}\;.         \label{eq : mp}\end{aligned}\ ] ] comparison of the obtained @xmath155 with the planck scale of the einstein gravity @xmath156  @xcite would be a crucial test of the supergravity . prospects of the `` planck scale '' measurement will be discussed in sec.[subsec : mpl ] .    it should be noted that the undetectable particle @xmath22 may not be the gravitino @xmath1 . any particle which couples weakly to @xmath0 can be particle @xmath22 . if the decay width into @xmath157 is larger than the decay width into @xmath158 , i.e. @xmath159 , , i.e. , @xmath160 . ] the cnlsp lifetime may be different from the supergravity prediction obtained from measured @xmath148 and tau jet energy distribution . inconsistency between the measured and predicted lifetime immediately means a discovery of a new sector that may not be accessible otherwise . one of well motivated examples of such a non - sugra decay is @xmath161 where @xmath2 is the axino , superpartner of the axion . the cnlsp @xmath0 decay into axino is studied in ref . @xcite for hadronic , or ksvz axion models  @xcite . in this paper we adopt the set  up in ref . @xcite for the axino interaction , which we briefly describe here .    in ksvz axion model , the peccei - quinn ( pq ) mechanism  @xcite is realized in an action with new heavy quarks . when these heavy quarks are integrated out , anomalous terms involving the axion and gauge bosons are generated at low energy effective action . when the axion interaction is supersymmetrized , its fermionic superpartner , the axino @xmath2 , must be introduced . the axino mass @xmath162 can range between the ev and gev scale depending on the model and susy breaking scheme  @xcite , and we regard it as a free parameter in this paper . the coupling of the axino to the bino and the photon/@xmath163-boson at the scale below the peccei - quinn scale @xmath164 is given by the lagrangian @xmath165\tilde{b}(\cos\theta_wf_{\mu\nu } -\sin\theta_wz_{\mu\nu}).\ ] ]    . , width=340 ]    the action does not contain direct @xmath166 coupling and also strongly suppressed by the pq scale @xmath167 . the two body decay @xmath161 is induced by the one loop diagram shown in fig.[fig : twobody ] . the loop integral has a logarithmic divergence . this is because the effective vertex ( [ eq : axinobino ] ) is applicable only if the momentum is smaller than the heavy ( s)quark masses , whereas the loop momentum exceeds that scale . in the full theory , one calculates two  loop diagrams with the heavy ( s)quarks , which leads to a finite result  @xcite . here we regulate the logarithmic divergence with the cut  off @xmath168  @xcite , and the effective @xmath169 coupling is parameterized as @xmath170 here , @xmath171 and we take @xmath172 hereafter . the parameter @xmath173 is an order one parameter to represent the uncertainty coming from the cut  off procedure mentioned above . in this paper , we regard this as a free parameter . the two body decay width is given as  @xcite @xmath174      in this subsection , we estimate the statistical error in the cnlsp lifetime measured by the stopper  detector . the analysis is model  independent and applicable to any long - lived charged particle stopped in the stopper .    for each cnlsp stopped in the stopper  detector , the stopping time @xmath175 and the decaying time @xmath176 will be recorded . the lifetime is then measured by fitting the temporal distribution of the decaying events @xmath177 , where @xmath178 . here , we use a maximum  likelihood fitting and adopt the following procedure : first we define a time @xmath179 such that @xmath180 and @xmath181 , where @xmath182 is total number of stopped event . is defined by @xmath183 , where @xmath184 and @xmath185 is the decay time of the @xmath186-th event : @xmath187 . ] ( for large @xmath182 , this @xmath179 is already a good estimator of the lifetime . ) we then calculate the @xmath188 distribution as a function of lifetime @xmath40 : @xmath189 where @xmath190 is the poisson distribution , @xmath191 is the number of events in the @xmath192-th bin , and @xmath193 is the predicted average number of events in that bin : @xmath194 here , we take 1 bin @xmath195 . an example of the time distribution and the best fit curve for mc distribution are shown in fig . [ fig : lifetime_measure100 ] for @xmath196 .     for @xmath196 . in the left plot , the solid line is the generated events , and dashed line is the best fit . ] the @xmath197 confidence interval can be estimated by the range @xmath198 . [ fig : lifetime_measure100 ] shows that the error of the lifetime is @xmath199 for @xmath196 . we have done the same analysis for @xmath200 and found @xmath201 . so far , we assumed that most of the stopped cnlsps decay within the experimental time scale . we now consider the case where the lifetime is longer . suppose that one observes @xmath200 stopped cnlsps and only 10 events decaying within 1 year . in such a case , the lifetime is estimated from the number of the decaying events . for instance , 95% interval of the mean @xmath202 of poisson variable for @xmath203 is @xmath204 $ ]  @xcite . using @xmath205 , a 95% interval @xmath206 years is obtained . for much longer lifetime or smaller @xmath182 , only a lower bound of the lifetime is obtained . we have assumed that the background is negligible . the background from cosmic neutrino and hard neutrino produced from the main detector interaction point is small ( cf . @xcite ) , however careful study on the accidental background is necessary when statistics is low . this is beyond of the scope of this paper .      in this subsection , we estimate the uncertainty of the tau energy determination . then , in sec.[subsec : mx ] we discuss the kinematical reconstruction of the lsp mass . schematically , the procedure is as follows : @xmath207    when @xmath0 decays into @xmath40 and invisible particle @xmath22 , the tau energy @xmath135 is monochromatic ( see eq.([etau ] ) ) . @xmath151 can be obtained by fitting the @xmath40 jet energy distribution : @xmath208 among the decay products of the @xmath40 lepton we omit the neutrinos and muons . in order to see the prospects of @xmath151 measurement with a finite number of events , we generate the events from @xmath40 decay by using the tauola  @xcite , and we perform a maximum  likelihood fitting of low statistics ( `` experimental '' ) event sets by high statistics ( `` theoretical '' ) distributions .    in practice , the observed jet energy distribution depends not only on the @xmath151 , but also on the detector resolution @xmath209 and the tau polarization @xmath210 : @xmath211 to obtain theoretical predictions for @xmath212 , we have generated tau decay events with high statistics run of tauola for the parameter space @xmath213 30  125  gev and @xmath214 , while fixing @xmath215 @xmath216 @xmath217 . events for each of the parameter sets @xmath218 30 , 31 , 32 , ... 109 , 125  gev and @xmath219 ( i.e. , @xmath220 parameter sets @xmath221 events ) , and interpolated the distribution between those parameter points . ] for each single @xmath40 decay generated by tauola , the jet energy @xmath222 is calculated , and then smeared by a gaussian fluctuation with a variance @xmath223 .    , 70 , 90  gev from left to right . the dashed lines show the spectrum without energy resolution effect , and the solid lines show the spectrum with detector energy resolution @xmath224 @xmath225 . ] fig  [ fig : taujet ] shows examples of the jet energy distribution .    in order to see a realistic situation , we generate a small number of tau lepton decays for a fixed parameter set of ( @xmath151 , @xmath210 , @xmath226 ) , and then fit the result by the `` theoretical '' distribution obtained above . here and hereafter , we assume that the energy resolution will be known in advance . examples are shown in fig . [ fig : lnlfit:1000 ] and fig . [ fig : lnlfit:100 ] for total number of decaying tau leptons @xmath227 and @xmath228 , respectively . the energy of @xmath40 jets which are not contained in the detector may not be measured precisely . the @xmath229 corresponds to the number of well contained events . gev , @xmath230 and @xmath231 ( solid histogram ) , and the best fit distribution ( dotted histogram ) . only the bins between the vertical lines are used for the fit ( see text ) . ( b ) @xmath232 projected onto the @xmath151 axis . ( c ) contour plots of @xmath233 , 4 , 9 projected onto the ( @xmath210 , @xmath151 ) plane . ] [ fig : lnlfit:1000 ]     but with 100 decaying tau leptons . ] [ fig : lnlfit:100 ]    some comments are in order here . ( i ) we take 1 bin = 10  gev . ( ii ) for the maximum  likelihood fitting we take only the bins with number of events @xmath234 , and the bins above the peak energy , which in the examples of fig . [ fig : lnlfit:1000](a ) and [ fig : lnlfit:100](a ) correspond to the bins between the two vertical lines . ( iii ) we then calculate the @xmath188 distribution in the parameter space of ( @xmath151 , @xmath210 ) as follows @xmath235 where @xmath190 is the poisson distribution , @xmath236 is the number of events in the @xmath192th bin for the low statistics run , and @xmath237 is the predicted number of events in the @xmath192th bin as a function of the parameter set @xmath238 , normalized by the total number of decaying tau leptons . here , we assume that the total number of stopped @xmath0 is known , i.e. , the total number of events is _ not _ taken as a free parameter for the fit .    the @xmath197 confidence interval can be estimated by the range @xmath198 projected onto the @xmath151 axis [ see figs . [ fig : lnlfit:1000](b ) and [ fig : lnlfit:100](b ) ] . one can see that the primary tau energy can be determined within an error of a few  gev . from fig . [ fig : lnlfit:1000](c ) one can also see that the polarization is hardly determined by the energy distribution analysis even with @xmath239 . this is because the sensitivity to the polarization becomes very weak once the finite detector resolution is taken into account ( cf . [ fig : taujet ] ) .    to estimate the statistical error in @xmath135 measurement we have generated the event sets with same statistics @xmath240 @xmath241 and repeat the fit on @xmath135 and @xmath242 to obtain the best fit value of the @xmath135 , @xmath243 . in figure [ fig : deletau ] , we show the distribution of @xmath243 for @xmath244  gev , @xmath245 , and @xmath246 .     for @xmath247 , 70 and 90 gev ( from left to the right ) , for @xmath239 ( left ) and @xmath248 ( right ) . ]    from the variances of these distributions , we estimate the @xmath249 uncertainty of @xmath151 as @xmath250      now we can estimate the uncertainty of the lsp mass . in fig .  [ fig : mgmg ] we plot the range of reconstructed lsp mass @xmath251 where @xmath252 and @xmath253 . uncertainty of @xmath151 , for @xmath239 and 100 , and @xmath254 and 300  gev . solid lines represent the case without an error of @xmath148 , and dashed lines include 1% uncertainty of @xmath148 . vertical dotted lines represent the stau lifetime of 1 month , 1 year , and 10 years in the case of gravitino lsp ( @xmath39 ) . ]    note that the reconstructed lsp mass depends on not only @xmath151 but also the measured cnlsp mass @xmath255 . in fig . [ fig : mgmg ] we show the range of lsp mass for @xmath256 ( solid lines ) and @xmath257 ( dashed lines ) , the latter corresponding to 1% uncertainty of the cnlsp mass . the stau mass determination from time of flight was discussed in @xcite . for cms detector the mass resolution is estimated as 1020% in each event and less than 1% for @xmath258 events .    as can be seen in fig . [ fig : mgmg ] , the kinematical reconstruction of the lsp mass is possible if @xmath149 is sufficiently large , @xmath259 for @xmath260 and @xmath261 for @xmath262 . otherwise one can get only an upper bound on the mass @xmath149 .      finally , if the lsp is the gravitino ( @xmath39 ) , the uncertainty of the reconstructed gravitino mass @xmath263 translates into an uncertainty of the supergravity planck scale , which is obtained by substituting @xmath263 for @xmath60 in eq.([eq : mp ] ) . to take into account the error of @xmath148 , we also substitute @xmath255 for @xmath148 . eq.([eq : mp ] ) then becomes @xmath264 which is shown in fig . [ fig : mgmpl ] for @xmath256 ( solid lines ) and @xmath265 ( dashed lines ) . here we do not include the uncertainty of the lifetime measurement , which simply affects the measured planck scale as @xmath266 . neglecting the @xmath40lepton mass , eq . ( [ calmpl ] ) is simplified as @xmath267     uncertainty of @xmath151 , for @xmath239 and 100 , and @xmath268 , 200 , and 300  gev . solid lines represent the case without an error of @xmath148 measurement , and dashed lines include 1% uncertainty of @xmath148 . vertical dotted lines represent the stau lifetime of 1 month , 1 year , and 10 years . ] as discussed in the previous section , the kinematical reconstruction of the gravitino mass is possible only if @xmath60 is sufficiently large , @xmath269 for @xmath260 and @xmath270 for @xmath262 . for smaller values of the gravitino mass , one can get only a lower bound on the planck scale @xmath271 . one can see whether the determined ` planck scale ' is inconsistent with the planck scale of the einstein gravity , @xmath272  gev . in other words , one can test the assumption of the decay @xmath273 by comparing the observed lifetime with the predicted lifetime . for example , if the nlsp dominantly decays into axino , the lifetime would be far shorter compared with gravitino assumption for a fixed mass @xmath149 , leading to a smaller value of fitted @xmath271 . suppose for instance one measures @xmath274  gev and @xmath275  gev . when @xmath25 , the lifetime becomes @xmath276 for @xmath277 and @xmath278 [ cf . ( [ eq : axinogamma ] ) ] . if one uses eq . ( [ calmpl ] ) , a `` planck scale '' @xmath279 would be obtained , thereby falsifying the gravitino assumption . we now discuss msugra model points discussed in section [ sec : modelpoints ] , together with cosmological constraints . in the early universe , the @xmath0 cnlsp has been in thermal equilibrium until its decoupling , @xmath280 . if such a particle has a very long lifetime , as discussed in this paper , its decay during or after the big - bang nucleosynthesis ( bbn ) , @xmath281  mev , may spoil the successful prediction of bbn  @xcite . in the models with @xmath0 nlsp and @xmath1 lsp , this leads to severe constraints on the parameter space of ( @xmath148 , @xmath60 ) , in particular to upper bounds on the gravitino mass for a given stau mass  @xcite . furthermore , it has been recently pointed out that a heavy charged particle can form a bound state with light elements during the bbn , which can lead to new effects and/or severer constraints  @xcite . in this paper we do not discuss these effects because it is difficult to evaluate the net effect quantitatively and it still awaits detailed analysis . we should also mention that those bbn constraints may disappear if there is entropy production between the stau decoupling ( @xmath282 ) and the bbn ( @xmath283  mev ) , because the stau abundance is diluted before its decay  @xcite .    keeping in mind the possibilities of severer bounds and also a possible loophole , let us discuss the cases of @xmath268  gev and @xmath284  gev , corresponding to the msugra model points @xmath285 and @xmath123 . according to the latest analyses  @xcite including the effects of the hadronic decay  @xcite , the bounds on the gravitino are @xmath286@xmath287 for @xmath288 and @xmath289@xmath290 for @xmath291.he . if one adopts it , the constraints become severer ( cf . @xcite ) . ] the ranges of upper bounds correspond to the uncertainties of various bounds from primordial light elements .    for @xmath274 gev ( model point @xmath285 ) , as can be seen in table  [ masssugra ] , one could collect more than 1000 staus for 300 fb@xmath29 . the bound @xmath286@xmath292  gev then suggests that the measurement of the gravitino mass and the `` planck scale '' may become possible if one assumes conservative bbn bounds and if the gravitino mass is sufficiently large ( cf . [ fig : mgmg ] and [ fig : mgmpl ] ) . for @xmath293  gev ( model point @xmath123 ) , the measurement would become easier if one could collect the same number of staus ( cf . figs . [ fig : mgmg ] and [ fig : mgmpl ] ) . however , from table  [ masssugra ] we find that expected number of stopped cnlsp is around 30 for 300 fb@xmath29 . this is because gluino mass is above 2 tev for this point and production cross section is small . one needs slhc ( @xmath294fb@xmath29 ) to collect @xmath134 events .    in the case of axino lsp , the bbn bound is much weaker because the lifetime of the cnlsp stau becomes much shorter [ see eq.([eq : axinogamma ] ) ] . hence , the axino mass measurement is plausible for sufficiently large @xmath295 . when @xmath296 , it is difficult to determine the gravitino mass from @xmath40 energy measurement at stopper  detector . axino is a motivated candidate which couples weakly to the mssm particle with comparable strength to the gravitino . in this section we therefore compare the decay @xmath297 with the decay @xmath298 . for simplicity , we will assume the nlsp is pure ` right - handed ' stau , @xmath299 throughout this section . extension to the case with mixing with @xmath300 is straightforward , and the mixing angle dependence is expected to be small . the gravitino is a spin-3/2 particle . however , in the limit where @xmath301 , the effective interaction to mssm particles would be reduced to that of spin 1/2 particle , goldstino @xmath302 . the effective action relevant to the @xmath303 decay is given as follows , @xmath304\tilde{b }   ( \cos\theta_wf_{\mu\nu } -\sin\theta_wz_{\mu\nu } ) . \label{eq : actiongv}\ ] ] the action is similar to that of axino given in eqs.([eq : axinobino ] ) and ( [ eq : axinostau ] ) except the coupling coefficients . the relative weight of the two terms in eq.([eq : actiongv ] ) are fixed by the supergravity , while for the axino the coefficient in eq.([eq : axinostau ] ) is induced from eq.([eq : axinobino ] ) by the radiative corrections . note that the term proportional to @xmath305\tilde{b}f_{\mu\nu}$ ] is a non renormalizable coupling of the photon to gravitino or axino and induces significantly different @xmath306 , @xmath40 distribution . the axino three body decay @xmath307 proceeds through the diagrams shown in fig.[fig : threebody ] , where the hatched triangle express the effective vertex shown in eq . ( [ eq : axinostau ] ) . on the other hand , the relevant diagrams for the three body decay into goldstino @xmath308 , are given in fig.[fig : threebody_gra ] . the diagram corresponding to the top right of fig . [ fig : threebody ] does not exist for the goldstino case . the difference of the actions and the relevant diagrams will appear as the deviation of the decay distributions . in the appendix , we list the three body decay differential width into gravitino / axino in the limit where the gravitino / axino mass can be neglected compared to @xmath148 . the formula for the massive gravitino and axino are given in ref . @xcite .    . , width=377 ]     in the goldstino limit . , width=377 ]      the three body decay @xmath309 should be visible in the stopper  detector if it has an ability to measure charged tracks , and also segmented into small units . the position where @xmath0 decays in the detector is the position where the charged track by the @xmath310 , @xmath44 and @xmath311 from @xmath40 is initiated . for hadronic tau decays , a @xmath312 is always in the decay products , sometimes with photons coming from @xmath43 decays . photons are converted into electron after passing @xmath313 one radiation length @xmath47 . therefore for the iron based detector with @xmath314g/@xmath33 , the photon shower starts @xmath315 cm from the decay point . in summary , the three body decay of @xmath316 is identified as a charged track ( which might be associated with collinear em showers ) + an isolated hard em shower pointing back to the point where charged track is started . if the segmentation is not good enough , the efficiency to discriminate em showers from @xmath312 would be reduced .    as can be seen in the appendix , the three body decay amplitude can be written as a function of the angle between photon and @xmath40 , @xmath36 , and @xmath317 . to be conservative we assume the energy resolution for the isolated photon shower is @xmath318 and ignore the angular resolution of the photon momentum . ( in the following @xmath319 denotes the photon energy after taking into account of this finite resolution effect . ) note that the shortest length of the detector is 3.5 m , which corresponds to 1750g/@xmath35 for @xmath320g/@xmath321 . the em showers are likely contained in the stopper because we only need 200g/@xmath35 to fully absorb them . we need to require several cuts for the accepted events . * experimentally , the angle between @xmath40 and @xmath306 must be large enough to avoid the overlap between @xmath40 decay products and prompt @xmath306 . we only use the events where @xmath322 . we do not lose sensitivity to the differece between the two decays @xmath323 by cutting these events . in the collinear region , the amplitude is dominated by the contributions from qed vertex , which is common for both of the decays . * the three body decay amplitude suffers soft and collinear singularity . because we only adopt simple leading order calculation , we require @xmath324 10  gev and @xmath325  gev . we define the ratio @xmath326 the dependence of @xmath105 on the mssm parameter is quite different between @xmath327 and @xmath39 . while @xmath328 is one loop process controlled by the parameter @xmath173 , the three body decay contains a tree level contribution which depends on non - renormalizable axino-@xmath99-gauge coupling ( top left of fig . [ fig : threebody ] ) . when @xmath173 is small , the tree level contribution plays a dominant role in the three body decay into axino . this can be seen in fig . [ fig : br ] , where the ratio @xmath329 is plotted as a function of the mass difference @xmath330 . we also fix @xmath331  gev but @xmath332 is insensitive to @xmath333 . @xmath329 is enhanced when the mass difference between @xmath334 and @xmath148 is small relative to the @xmath0 mass . this is because the two body decay of axino is suppressed by @xmath335 . the branching ratio ranges from 8% to 0.5% for the model parameters given in the figure .     as a function of @xmath330 . the solid , dashed , dot - dashed , short - dashed curves correspond to the case where @xmath336  gev respectively . lines which increase toward smaller mass difference corresponds to @xmath327 . the @xmath39 lines do not show significant bino mass dependence . , width=264 ]    on the other hand , the branching ratio @xmath337 is well below 1% for the parameter given in fig.[fig : br ] . @xmath338 @xmath339 for @xmath340  gev respectively . when the three body decay branching ratio turns out be above 5% , 200 stopped @xmath0 is enough to see the 3 @xmath89 deviation from the gravitino assumption . we estimate the number of events @xmath341 stopped in the detector , which is required to find 3@xmath89 deviation from the gravitino scenario as follows ; @xmath342 the @xmath343 as a function of @xmath330 is given in fig.[fig : total ] . each curve increases as @xmath330 is increased up to the value where @xmath344 coincides with @xmath345 . when @xmath346  gev , @xmath137 ( @xmath347 ) stopped @xmath0 are enough to see the deviation from the measurement of @xmath105 for @xmath348 respectively .     deviates more than 3 sigma from the prediction for @xmath349 . the lines corresponds to different stau mass as in fig.[fig : br ] . , width=264 ]     in some @xmath36 intervals for @xmath350 decay . left panel : @xmath327 @xmath351 , right panel : @xmath39 . @xmath153  gev , and @xmath352  gev . , title="fig:",width=188 ]   in some @xmath36 intervals for @xmath350 decay . left panel : @xmath327 @xmath351 , right panel : @xmath39 . @xmath153  gev , and @xmath352  gev . , title="fig:",width=188 ]    gev and right for @xmath353  gev , title="fig:",width=257 ]  gev and right for @xmath353  gev , title="fig:",width=264 ]    not only the branching ratio , but also the decay distribution contain the information of the invisible particle . the axino decay distribution can be enhanced at the region where @xmath354 is large and @xmath355 is small , namely hard photon and @xmath40 is back to back  @xcite . this occurs when @xmath173 is small and relative importance of the direct axino - bino - gauge boson coupling is enhanced .    in fig . [ fig : egamdist ] , left panel shows the @xmath317 distribution for different @xmath355 intervals . here we fix @xmath356 , @xmath352   gev , and @xmath153  gev . the distribution has significantly hard component for @xmath357 when compared with that of gravitino in the right panel . the enhancement of back - to - back events is a signature of an axino , and clear @xmath40 and @xmath306 separation is not required to distinguish them . we now estimate the the number of cnlsp decays in the stopper to see the 3 @xmath89 deviation between @xmath358 and @xmath359 from the decay distribution . we define @xmath360 like function from the differential decay width ; @xmath361 here , @xmath362 is the number of event in a @xmath192-th bin for @xmath363 when the number of stopped @xmath0 is @xmath364 . @xmath365 and @xmath366 is the bin size in @xmath317 and @xmath36 . we divide the @xmath317 and @xmath36 into 8 bins and 3 bins respectively , for the range @xmath367 and @xmath368 ; @xmath369 the photon energy resolution of the detector is assumed @xmath370 . we apply the cut @xmath371 , @xmath372  gev and @xmath373  gev . on the other hand , @xmath374 is the number of event in a @xmath192-th bin for @xmath375 with @xmath376 normalized so that the total number is same to that of axino three body decay . namely , we do not use the information for @xmath377 in our fit .    from @xmath378 , we define @xmath379 , the number of stopped @xmath0 required to see the 3@xmath89 deviation between @xmath363 and @xmath375 as follows , @xmath380 we show the @xmath381 as a function of @xmath330 in fig . [ fig : dist ] . the sensitivity is significantly increased from the estimate using the branching ratio only . for @xmath382  gev , the deviation may be visible for @xmath137 events for @xmath383 ( @xmath347 events for @xmath356 ) .    finally , we estimate sensitivity at our model points in section 3 . for simplicity we assume @xmath384 . in general @xmath303 is mixed with @xmath300 . the mixing angle is defined as @xmath385 for model points discussed in table [ mass ] , the angle is @xmath386 . the effect of the mixing angle in the axino decay is small because the amplitude of @xmath387 is suppressed by both by the small @xmath388 factor and smaller hypercharge , and can be safely ignored .    in fig.[fig : lambda ] the expected sensitivity at the stopper  detector is shown . here , long - dashed ( long - short - dashed , dashed ) line corresponds to the required statistics for @xmath383 ( 0.75 , 1 ) for different @xmath389 , while the upper and lower solid lines correspond to the number of stopped cnlsp for 300 fb@xmath29 and 3000 fb@xmath29 . one can address the difference between axino and gravitino for @xmath390tev or less for @xmath391 fb@xmath29 and @xmath356 ( 0.75 ) . note that the expected reach in @xmath389 is essentially determined by the susy production cross section , as they decrease steeply with increasing @xmath389 . we therefore show the gluino masses scale on the top of the figure . luminosity for gm points with @xmath392tev . dashed lines are number of required stopped events to see 3-@xmath89 deviation in the @xmath393 distribution of @xmath394 from those expected for the decay @xmath395 . , width=264 ] in this paper we investigate the physics of the long - lived charged next lightest susy particle ( cnlsp ) , which may be explored at a massive stopper  detector placed next to the cms detector at the lhc . we assume the cnlsp is the lighter scalar tau , @xmath9 , which decay into @xmath396 where @xmath22 is a invisible particle . a natural candidate of the particle @xmath22 is a gravitino but we also consider the case where @xmath22 is axino @xmath2 . in this paper , we assume very large stopper  detector next to the cms detector , with total mass of 8 kton . the stopper must have a capability to measure the position where the cnlsp stopped , and also the energy of the @xmath40 decay products . if the detector can be highly segmented , it is also possible to identify the @xmath40 decay products separately . the number of stopped nlsp ranges from @xmath397 to @xmath134 for gluino and squark with mass below 2 tev . if the size of the detector should be smaller , the number of the events must be scaled down linearly . we estimate the statistical error for @xmath135 , which can be determined from the end point of the tau jet energy . we assume that the energy resolution of the detector is 150%/@xmath398 . @xmath399 ( @xmath400 ) can be achieved for the statistics of the order of @xmath401 . when lsp mass is above @xmath402 , one can constrain the gravitino mass both from above and below . from the lifetime and gravitino mass measurement , the supergravity planck scale @xmath271 can be measured . if many cnlsp can be accumulated , one can even study the three body decay of the cnlp @xmath50 . when the dominant decay mode of the @xmath9 is @xmath403 , we may be able to see the deviation of the decay branching ratio and the distribution from the @xmath39 case . finally we comment on the strategy to proceed this experiment . susy particles will be found in the early stage of the lhc experiment if susy scale is @xmath404 . if lsp is gravitino or axino and the nlsp is charged and long - lived , it would also be recognized easily . the detector proposed in this paper may be placed after the existence of long - lived cnlsp is observed , roughly at the same time to the high luminosity run of the lhc , or the proposed super lhc run . lhc experiment has a great potential to explore new physics in tev regions . it is important to explore new possibilities that can be done with lhc . in this paper we have proposed a large additional detector in the cms cavern when the long - lived cnlsp is found . this requires a significant modification of the cms experiment . the reward is low systematics study of the cnlsp decay which primarily serve for the study of the gravitino sector . the determination of gravitino mass either from the lifetime or ( independently ) from @xmath152 measurement would give us a direct information of the the total susy breaking scale . together with high precision determination of the mssm sector expected with the cnslp momentum information , the nature of the interaction of the mssm sector and hidden sector can be studied in detail . we hope this paper is useful for further , and more realistic studies . we thank j.  ellis , h.  u.  martyn and f.  d.  steffen for discussions . this work is supported in part by the grant - in - aid for science research , ministry of education , science and culture , japan ( no.16081207 , 18340060 ) for mmn . in this appendix we list the differential decay width of @xmath303 decay into @xmath405 where @xmath22 is either gravitino or axino in limit where @xmath406 . formulas in this appendix are obtained by taking massless limits @xmath407 of the formulas in ref . @xcite . the @xmath303 decay width to @xmath408 is given as follows ; @xmath409 ^ 2 } \sum_{\rm spins } \vert { \cal m}(\tilde{\tau}_r\rightarrow \tau\gamma { \tilde{g}})\vert^2\ ] ] where @xmath410 and @xmath411 \cr & & + \frac{1-{x_{\gamma}}}{2-{x_{\gamma}}(1-{\cos\theta } ) } \left ( -4 + 2x-\frac{4}{{x_{\gamma } } } \right ) + \frac{4(1-{x_{\gamma}})^2}{\{2-{x_{\gamma}}(1-{\cos\theta})\}^2 } \cr & & + \frac{1}{\ { { x_{\gamma}(1+\cos\theta)-2a_{\tilde{b}}+a_{\tilde{b}}x_{\gamma}(1-\cos\theta ) } \}^2}\biggl [ \cr & & ~~ 2{a_{\tilde{b}}}^2(1-{x_{\gamma } } ) \left\{{x_{\gamma}}^2 - 2{x_{\gamma}}+\frac{4}{{x_{\gamma}}}-({x_{\gamma}}^2 - 2{x_{\gamma}}-2){\cos\theta}\right\ } \cr & & ~~ + { a_{\tilde{b}}}(1+{\cos\theta})(4 - 3{x_{\gamma}}+{x_{\gamma}}{\cos\theta})({x_{\gamma}}^2-{x_{\gamma}}-2 )   \cr & & ~~ -(1+{\cos\theta})^2{x_{\gamma}}({x_{\gamma}}^2-{x_{\gamma}}-2 ) \biggr ] .\end{aligned}\ ] ] here @xmath412 and @xmath36 is the angle between @xmath40 and @xmath306 . k.  hamaguchi , y.  kuno , t.  nakaya and m.  m.  nojiri , phys . d * 70 * ( 2004 ) 115007 [ arxiv : hep - ph/0409248 ] . j.  l.  feng and b.  t.  smith , phys . d * 71 * ( 2005 ) 015004 [ erratum - ibid . d * 71 * ( 2005 ) 0109904 ] [ arxiv : hep - ph/0409278 ] . a.  de roeck , j.  r.  ellis , f.  gianotti , f.  moortgat , k.  a.  olive and l.  pape , arxiv : hep - ph/0508198 . w.  buchmuller , k.  hamaguchi , m.  ratz and t.  yanagida , phys . b * 588 * ( 2004 ) 90 [ arxiv : hep - ph/0402179 ] ; arxiv : hep - ph/0403203 . a.  brandenburg , l.  covi , k.  hamaguchi , l.  roszkowski and f.  d.  steffen , phys . b * 617 * ( 2005 ) 99 [ arxiv : hep - ph/0501287 ] . cms technical proposal , cern / lhcc/94 - 38 ( 1994 ) . atlas technical proposal cern / lhcc/94 - 43 ( 1994 ) ; + atlas liquid argon tdr , cern / lhcc/96 - 41 ( 1996 ) . s.grohmann , cern , st - cv , `` concepts of refrigeration circuit design for cms tracker cooling , '' at the homepage of ts / cv / dc section of cern ; + ` http://st-support-cooling-electronics.web.cern.ch/st-support-cooling-electronics/ `    s.  eidelman _ et al . _ [ particle data group ] , phys . b * 592 * ( 2004 ) 1 . r.  k.  bock , t.  hansl - kozanecka and t.  p.  shah , nucl . instrum . * 186 * ( 1981 ) 533 . w.  w.  m.  allison _ et al . _ [ soudan-2 collaboration ] , nucl . .  meth . a * 376 * ( 1996 ) 36 . ` http://hepunx.rl.ac.uk/soudan2/ ` h.  u.  martyn , eur . j.  c * 48 * ( 2006 ) 15 [ arxiv : hep - ph/0605257 ] . d.  e.  kaplan , g.  d.  kribs and m.  schmaltz , phys . d * 62 * ( 2000 ) 035010 [ arxiv : hep - ph/9911293 ] ; + z.  chacko , m.  a.  luty , a.  e.  nelson and e.  ponton , jhep * 0001 * ( 2000 ) 003 [ arxiv : hep - ph/9911323 ] . w.  buchmuller , k.  hamaguchi and j.  kersten , phys . b * 632 * ( 2006 ) 366 [ arxiv : hep - ph/0506105 ] . g.  f.  giudice and r.  rattazzi , phys . * 322 * ( 1999 ) 419 [ arxiv : hep - ph/9801271 ] .            t.  goto and m.  yamaguchi , phys .  lett . b * 276 * ( 1992 ) 103 ; + e.  j.  chun , j.  e.  kim , and h.  p.  nilles , phys . b * 287 * ( 1992 ) 123 [ arxiv : hep - ph/9205229 ] ; + e.  j.  chun and a.  lukas , phys . b * 357 * ( 1995 ) 43 [ arxiv : hep - ph/9503233 ] ; + and references therein . s.  jadach , j.  h.  kuhn and z.  was , comput . commun . * 64 * ( 1990 ) 275 ; + m.  jezabek , z.  was , s.  jadach and j.  h.  kuhn , comput . commun . * 70 * ( 1992 ) 69 . + s.  jadach , z.  was , r.  decker and j.  h.  kuhn , comput . commun . * 76 * ( 1993 ) 361 . i.  hinchliffe and f.  e.  paige , phys . d * 60 * ( 1999 ) 095002 [ arxiv : hep - ph/9812233 ] . i.  v.  falomkin , g.  b.  pontecorvo , m.  g.  sapozhnikov , m.  y.  khlopov , f.  balestra and g.  piragino , nuovo cim . a * 79 * ( 1984 ) 193 [ yad . * 39 * ( 1984 ) 990 ] ; + m.  y.  khlopov and a.  d.  linde , phys .  lett . b * 138 * ( 1984 ) 265 ; + j.  r.  ellis , j.  e.  kim and d.  v.  nanopoulos , phys . b * 145 * ( 1984 ) 181 . m.  pospelov , arxiv : hep - ph/0605215 ; + k.  kohri and f.  takayama , arxiv : hep - ph/0605243 ; + m.  kaplinghat and a.  rajaraman , phys . d * 74 * ( 2006 ) 103004 [ arxiv : astro - ph/0606209 ] . + r.  h.  cyburt , j.  ellis , b.  d.  fields , k.  a.  olive and v.  c.  spanos , arxiv : astro - ph/0608562 .
Abstract:

Answer: if the scalar tau @xmath0 is the next lightest supersymmetric particle and decays into a gravitino(@xmath1 ) being the lightest supersymmetric particle , it will have generally a very long lifetime . in this paper , we investigate the possibility to study the decay of such a long lived scalar tau at the lhc . if we can add to the present lhc experiments additional detectors which are able to stop the stau particles and measure the produced decay products , the decay characteristics can be studied precisely at the lhc . we identify a maximum `` stopper detector '' that could be added in the cms cavern , and estimate the sensitivity to the lifetime of the stau and to the mass of gravitino with this detector . the decay of the scalar tau may be significantly modified if the decay channel to the axino @xmath2 is open . we study the possibility to distinguish such decays from decays into gravitinos by measuring the process @xmath3 using the stopper detector .