FlyoDoc_2011 Pisa 2011 by GmP --- 011
|
00001 *** Studio decadimenti K Tirrenia 19.08.01 00002 *** csi = sigma+ e nu 00003 *** file csisp.epc decadimento con solo spazio delle fasi 00004 *** per decadimento con elmento di matrice vedi csimx.epc 00005 *** 00006 **dbginput 00007 illumina 00008 errori 00009 scatter 00010 perdita 00011 **pairprod // 1 pairprod 00012 ***pairbrem // pair prod followed by brems. 00013 **bremsstrhalung 00014 ppp 1.e12 //protoni per puls o burst 00015 tburst 3.0e+9 //lunghezza burst in sec 00016 correla 1 //== 0 no corr., 1 run corr., 2 run and event corr. 00017 run 1 00018 ***wrtdst // Scrive le tracce viste da un rivelatore 00019 ***wrtdstall // Scrive tutto 00020 **wrtnt 00021 wrtntall 00022 maxevnt 300 00023 debgf 0 // debug first 00024 debgl 4 // debug last 00025 ***** 00026 data_base 00027 gam 0. 0. -1. 0 00028 elec 0.00051 0. -1. -1 00029 elep 0.00051 0. -1. 1 00030 nu 0. 0. -1. 0 00031 mum .10565 0. -1. -1 00032 mup .10565 0. -1. 1 00033 n .93957 0. -1. 0 00034 p .93827 0. -1. 1 00035 pip .13956 0. 780.45 1 00036 pim .13956 0. 780.45 -1 00037 pi0 .13497 0. 0. 0 00038 eta .5488 0. 0. 0 00039 Kp .49367 0. 371.3 1 00040 Km .49367 0. 371.3 -1 00041 Kl .49767 0. 1551.0 0 00042 Ks .49767 0. 2.676 0 00043 Ro .789 .150 .0 0 00044 W .783 .008 .0 0 00045 fy 1.0194 .00441 .0 0 00046 sigp 1.18937 0. 2.404 1 00047 csi 1.192642 0. 8.71 0 00048 res .0 .0 .0 0 00049 cms .0 .0 .0 0 00050 dlz .0 .0 .0 0 00051 pic2 .0 .0 .0 0 00052 pi02 .0 .0 .0 0 00053 rsd4 .0 .0 .0 0 00054 ppgg .0 .0 .0 0 00055 nusg .0 .0 .0 0 00056 done 00057 unseen nu 00058 **fisica cms 10 00059 ** 00060 decay csi // rate minus the dalitz contribution 00061 rate 1.0 00062 target 0. 7.19 0.0 0.1 40. 00063 pinc 0. 0. 120. 75. 0.0036 450. 00064 coll 0.0 6.93 480. 0.18 1. 00065 reazione csi = elec sigp nu 00066 csi = elec sigp nu 00067 sigp = p pi0 00068 pi0 = gam gam 00069 fine 00070 end 00071 00072 00073 *********************************************************** 00074 ***************** A p p a r a t o 00075 fiducial 607. 5000. // from defining collimator .. attenzione in 00076 zona 12180. // ultima zeta dell'apparato **** 00077 device corona 00078 **writeout 00079 nome Coll // questo per Kshort 00080 detect mum 00081 detect mup 00082 detect all dump 00083 flags none 00084 centro 0.0 6.887 560. 00085 rin .3 0.0 40. //spessore fittizio 00086 sizeout 50. 50.0 40. // per non intercettare i kl.. 00087 end 00088 ***** 00089 device corona 00090 nome Kvlr 00091 **writeout 00092 detect gam 00093 detect elec 00094 detect elep 00095 detect mup 00096 detect mum 00097 detect char 00098 flags none 00099 centro 0.0 0.0 9568.7 00100 densita 0.8 1.5 00101 radi 0.0034 0.0 00102 rin 7.9 0.0 0.1 00103 rout 140. 0.0 0.1 00104 end 00105 device chamber 00106 nome Cam1 00107 writeout 00108 flags trigger 00109 detect gam noreg 00110 detect elec 00111 detect elep 00112 detect char 00113 centro 0.0 0.0 9708.6 00114 radi 0.0043 0.0 00115 rin 9.6 0.0 15. 00116 rout 135. 0.0 15. 00117 end 00118 *** 00119 device chamber 00120 nome Cam2 00121 writeout 00122 flags trigger 00123 detect gam noreg 00124 detect elec 00125 detect elep 00126 detect char 00127 centro 0.0 0.0 10627.1 00128 radi 0.0043 0.0 00129 rin 9.6 0.0 15. 00130 rout 135. 0.0 15. 00131 end 00132 *** 00133 *** Standard field 3575. 00134 device magnete 00135 nome Magn 00136 **writeout 00137 flags magn 00138 detect char 00139 centro 0.0 0.0 10896.8 00140 densita 0.000125 1.93 *** helium..vedi sotto 00141 radi 0.00 0.00048 00142 field 0.0 3575. 0.0 00143 sizein 120. 120. 120. 00144 sizeout 120. 120. 120. 00145 end 00146 *** 00147 device chamber 00148 nome Cam3 00149 writeout 00150 flags trigger 00151 detec gam noreg 00152 detect elec 00153 detect elep 00154 detect char 00155 centro 0.0 0.0 11166. 00156 radi 0.0043 0.0 00157 rin 9.6 0.0 15. 00158 rout 135. 0.0 15. 00159 end 00160 *** 00161 device chamber 00162 nome Cam4 00163 writeout 00164 flags trigger 00165 detec gam noreg 00166 detect elec 00167 detect elep 00168 detect char 00169 centro 0.0 0.0 11887. 00170 rin 9.6 0.0 15. 00171 rout 135. 0.0 15. 00172 end 00173 *** 00174 device hod48 00175 nome Hodo 00176 writeout 00177 flags trigger 00178 detect elec 00179 detect elep 00180 detect char 00181 centro 0.0 0.0 11988.5 00182 rin 12.8 0.0 1. 00183 rout 121. 0.0 1. 00184 end 00185 *** 00186 device lkr48 00187 **escludi // aggiungere ** all'inizio, per commentare 00188 nome Lkry 00189 writeout 00190 flags trigger 00191 detect gam dump 00192 detect elec dump 00193 detect elep dump 00194 detect mum 00195 detect mup 00196 detect char 00197 centro 0.0 0.0 12171.4 00198 densita 3.7 1.2 ** krypton .... 00199 rin 8.6 0.0 30. 00200 rout 130. 0.0 30. 00201 end 00202 *** 00203 device hac48 00204 **escludi // aggiungere ** all'inizio, per commentare 00205 nome Hadc 00206 writeout 00207 flags trigger 00208 detect mum 00209 detect mup 00210 detect pip dump 00211 detect pim dump 00212 detect Kp dump 00213 detect Km dump 00214 centro 0.0 0.0 12421.3 00215 densita 4.0 1.6 ** ferro e plastica 00216 rin 9.6 0.0 80. 00217 sizeout 135. 135. 80. 00218 end 00219 *** 00220 device rettcili 00221 nome Muvt 00222 writeout 00223 flags trigg 00224 detect char dump 00225 centro 0.0 0.0 12905.5. 00226 rin 10.6 0.0 80. 00227 sizeout 135. 135. 80. 00228 end 00229 ****** 00230 device magnete 00231 nome Blue 00232 **writeout 00233 escludi // va incluso se voglia l'effetto .. asteriscare!! 00234 flags magn 00235 detect char 00236 centro 0.0 0.0 4750. 00237 field .13 0.12 0.0 00238 ***field varibile fornito da User routine 00239 sizein 120. 120. 4250. 00240 sizeout 120. 120. 4250. 00241 end 00242 **** controllare dimensioni 00243 ***** helium gas: Densita 0.178 g/l, x0=94.32 g/cm2 00244 **** L0= 10*94.32/0.178 =5299 mt radi= 2*11.794/5299= 0.0045 00245 **** Min Ener 1.93 MeV/g/cm2 00246 device basedev 00247 nome Base 00248 **writeout 00249 escludi ** va esluso dalla loop dei rivelatori 00250 centro 0.0 0.0 10703.125 00251 densita 0.000125 1.93 00252 radi 0.0045 0.0 00253 sizeout 200. 200. 1171.85 00254 rin 9.6 9.6 1171.85 00255 end 00256 done 00257 *** 00258 end_all