#include <ReactorNtuple.hh>
Public Methods | |
ReactorNtuple (int ntId, const char *ntName, char *dirName) | |
ReactorNtuple constructor. More... | |
~ReactorNtuple () | |
ReactorNtuple destructor. | |
ReactorNtuple (const ReactorNtuple &Ntuple) | |
ReactorNtuple copy constructor. | |
ReactorNtuple & | operator= (const ReactorNtuple &rhs) |
ReactorNtuple overloaded = operator. | |
void | Fill (ReactorEvent &Event, ReactorDetector &Detector) |
Definition at line 12 of file ReactorNtuple.hh.
|
ReactorNtuple constructor. Called by ReactorFsim with an ntuple Id, name and directory name. Creates all the ntuple arrays and histograms as set by the user. Definition at line 56 of file ReactorNtuple.cpp.
00058 { 00059 00060 // define ntuple tags here 00061 // to add vars, increment ntVars, and add tag to end of list 00062 // don't change wordsize 00063 // 00064 int wordsize = 8; 00065 int ntVarsMax = 100; 00066 char ntTags[ntVarsMax*wordsize]; 00067 // 00068 int ntVars = 0; 00069 int id = 1; 00070 char Name[100]; 00071 // 00072 // book the general event veriables, nt=1001 00073 // 00074 strcpy(ntTags+ntVars++*wordsize,"nevent "); 00075 strcpy(ntTags+ntVars++*wordsize,"wtc "); 00076 strcpy(ntTags+ntVars++*wordsize,"wtf "); 00077 // 00078 IdGen = 1000*ntId+id++; 00079 sprintf(Name,"%s%s",ntName,", general"); 00080 NvarGen = ntVars; 00081 // 00082 hbookn_(IdGen,Name,NvarGen,dirName,1000,ntTags, 00083 strlen(Name),strlen(dirName),wordsize,NvarGen); 00084 // 00085 // reset 00086 ntVars = 0; 00087 // 00088 // book the nubar event variables, nt=1002 00089 // 00090 // up to 1 nubar per event 00091 strcpy(ntTags+ntVars++*wordsize,"nnu "); 00092 00093 strcpy(ntTags+ntVars++*wordsize,"enu1 "); 00094 strcpy(ntTags+ntVars++*wordsize,"nuproc1 "); 00095 strcpy(ntTags+ntVars++*wordsize,"vnu1r "); 00096 strcpy(ntTags+ntVars++*wordsize,"vnu1cos "); 00097 strcpy(ntTags+ntVars++*wordsize,"vnu1phi "); 00098 strcpy(ntTags+ntVars++*wordsize,"vnu1x "); 00099 strcpy(ntTags+ntVars++*wordsize,"vnu1y "); 00100 strcpy(ntTags+ntVars++*wordsize,"vnu1z "); 00101 // 00102 IdNu = 1000*ntId+id++; 00103 sprintf(Name,"%s%s",ntName,", nubar"); 00104 NvarNu = ntVars; 00105 // 00106 hbookn_(IdNu,Name,NvarNu,dirName,1000,ntTags, 00107 strlen(Name),strlen(dirName),wordsize,NvarNu); 00108 // 00109 // reset 00110 ntVars = 0; 00111 // 00112 // book the electron event variables, nt=1003 00113 // 00114 // up to 1 electron per event 00115 strcpy(ntTags+ntVars++*wordsize,"nele "); 00116 00117 strcpy(ntTags+ntVars++*wordsize,"KEele1 "); 00118 strcpy(ntTags+ntVars++*wordsize,"ele1cos "); 00119 strcpy(ntTags+ntVars++*wordsize,"eleproc1 "); 00120 strcpy(ntTags+ntVars++*wordsize,"gele1r "); 00121 strcpy(ntTags+ntVars++*wordsize,"gele1cos "); 00122 strcpy(ntTags+ntVars++*wordsize,"gele1phi "); 00123 strcpy(ntTags+ntVars++*wordsize,"gele1x "); 00124 strcpy(ntTags+ntVars++*wordsize,"gele1y "); 00125 strcpy(ntTags+ntVars++*wordsize,"gele1z "); 00126 strcpy(ntTags+ntVars++*wordsize,"vele1r "); 00127 strcpy(ntTags+ntVars++*wordsize,"vele1cos "); 00128 strcpy(ntTags+ntVars++*wordsize,"vele1phi "); 00129 strcpy(ntTags+ntVars++*wordsize,"vele1x "); 00130 strcpy(ntTags+ntVars++*wordsize,"vele1y "); 00131 strcpy(ntTags+ntVars++*wordsize,"vele1z "); 00132 strcpy(ntTags+ntVars++*wordsize,"tele1 "); 00133 // 00134 IdEle = 1000*ntId+id++; 00135 sprintf(Name,"%s%s",ntName,", electron"); 00136 NvarEle = ntVars; 00137 // 00138 hbookn_(IdEle,Name,NvarEle,dirName,1000,ntTags, 00139 strlen(Name),strlen(dirName),wordsize,NvarEle); 00140 // 00141 // reset 00142 ntVars = 0; 00143 // 00144 // book the positron event variables, nt=1004 00145 // 00146 // up to 1 positron per event 00147 strcpy(ntTags+ntVars++*wordsize,"npos "); 00148 00149 strcpy(ntTags+ntVars++*wordsize,"KEpos1 "); 00150 strcpy(ntTags+ntVars++*wordsize,"pos1cos "); 00151 strcpy(ntTags+ntVars++*wordsize,"posproc1 "); 00152 strcpy(ntTags+ntVars++*wordsize,"gpos1r "); 00153 strcpy(ntTags+ntVars++*wordsize,"gpos1cos "); 00154 strcpy(ntTags+ntVars++*wordsize,"gpos1phi "); 00155 strcpy(ntTags+ntVars++*wordsize,"gpos1x "); 00156 strcpy(ntTags+ntVars++*wordsize,"gpos1y "); 00157 strcpy(ntTags+ntVars++*wordsize,"gpos1z "); 00158 strcpy(ntTags+ntVars++*wordsize,"vpos1r "); 00159 strcpy(ntTags+ntVars++*wordsize,"vpos1cos "); 00160 strcpy(ntTags+ntVars++*wordsize,"vpos1phi "); 00161 strcpy(ntTags+ntVars++*wordsize,"vpos1x "); 00162 strcpy(ntTags+ntVars++*wordsize,"vpos1y "); 00163 strcpy(ntTags+ntVars++*wordsize,"vpos1z "); 00164 strcpy(ntTags+ntVars++*wordsize,"tpos1 "); 00165 // 00166 IdPos = 1000*ntId+id++; 00167 sprintf(Name,"%s%s",ntName,", positron"); 00168 NvarPos = ntVars; 00169 // 00170 hbookn_(IdPos,Name,NvarPos,dirName,1000,ntTags, 00171 strlen(Name),strlen(dirName),wordsize,NvarPos); 00172 // 00173 // reset 00174 ntVars = 0; 00175 // 00176 // book the neutron event variables, nt=1005 00177 // 00178 // up to 5 neutron per event 00179 strcpy(ntTags+ntVars++*wordsize,"nneu "); 00180 00181 strcpy(ntTags+ntVars++*wordsize,"KEneu1 "); 00182 strcpy(ntTags+ntVars++*wordsize,"neu1cos "); 00183 strcpy(ntTags+ntVars++*wordsize,"neuproc1 "); 00184 strcpy(ntTags+ntVars++*wordsize,"gneu1r "); 00185 strcpy(ntTags+ntVars++*wordsize,"gneu1cos "); 00186 strcpy(ntTags+ntVars++*wordsize,"gneu1phi "); 00187 strcpy(ntTags+ntVars++*wordsize,"gneu1x "); 00188 strcpy(ntTags+ntVars++*wordsize,"gneu1y "); 00189 strcpy(ntTags+ntVars++*wordsize,"gneu1z "); 00190 strcpy(ntTags+ntVars++*wordsize,"vneu1r "); 00191 strcpy(ntTags+ntVars++*wordsize,"vneu1cos "); 00192 strcpy(ntTags+ntVars++*wordsize,"vneu1phi "); 00193 strcpy(ntTags+ntVars++*wordsize,"vneu1x "); 00194 strcpy(ntTags+ntVars++*wordsize,"vneu1y "); 00195 strcpy(ntTags+ntVars++*wordsize,"vneu1z "); 00196 strcpy(ntTags+ntVars++*wordsize,"tneu1 "); 00197 strcpy(ntTags+ntVars++*wordsize,"nsteps1 "); 00198 strcpy(ntTags+ntVars++*wordsize,"absorb1 "); 00199 00200 strcpy(ntTags+ntVars++*wordsize,"KEneu2 "); 00201 strcpy(ntTags+ntVars++*wordsize,"neu2cos "); 00202 strcpy(ntTags+ntVars++*wordsize,"neuproc2 "); 00203 strcpy(ntTags+ntVars++*wordsize,"gneu2r "); 00204 strcpy(ntTags+ntVars++*wordsize,"gneu2cos "); 00205 strcpy(ntTags+ntVars++*wordsize,"gneu2phi "); 00206 strcpy(ntTags+ntVars++*wordsize,"gneu2x "); 00207 strcpy(ntTags+ntVars++*wordsize,"gneu2y "); 00208 strcpy(ntTags+ntVars++*wordsize,"gneu2z "); 00209 strcpy(ntTags+ntVars++*wordsize,"vneu2r "); 00210 strcpy(ntTags+ntVars++*wordsize,"vneu2cos "); 00211 strcpy(ntTags+ntVars++*wordsize,"vneu2phi "); 00212 strcpy(ntTags+ntVars++*wordsize,"vneu2x "); 00213 strcpy(ntTags+ntVars++*wordsize,"vneu2y "); 00214 strcpy(ntTags+ntVars++*wordsize,"vneu2z "); 00215 strcpy(ntTags+ntVars++*wordsize,"tneu2 "); 00216 strcpy(ntTags+ntVars++*wordsize,"nsteps2 "); 00217 strcpy(ntTags+ntVars++*wordsize,"absorb2 "); 00218 00219 strcpy(ntTags+ntVars++*wordsize,"KEneu3 "); 00220 strcpy(ntTags+ntVars++*wordsize,"neu3cos "); 00221 strcpy(ntTags+ntVars++*wordsize,"neuproc3 "); 00222 strcpy(ntTags+ntVars++*wordsize,"gneu3r "); 00223 strcpy(ntTags+ntVars++*wordsize,"gneu3cos "); 00224 strcpy(ntTags+ntVars++*wordsize,"gneu3phi "); 00225 strcpy(ntTags+ntVars++*wordsize,"gneu3x "); 00226 strcpy(ntTags+ntVars++*wordsize,"gneu3y "); 00227 strcpy(ntTags+ntVars++*wordsize,"gneu3z "); 00228 strcpy(ntTags+ntVars++*wordsize,"vneu3r "); 00229 strcpy(ntTags+ntVars++*wordsize,"vneu3cos "); 00230 strcpy(ntTags+ntVars++*wordsize,"vneu3phi "); 00231 strcpy(ntTags+ntVars++*wordsize,"vneu3x "); 00232 strcpy(ntTags+ntVars++*wordsize,"vneu3y "); 00233 strcpy(ntTags+ntVars++*wordsize,"vneu3z "); 00234 strcpy(ntTags+ntVars++*wordsize,"tneu3 "); 00235 strcpy(ntTags+ntVars++*wordsize,"nsteps3 "); 00236 strcpy(ntTags+ntVars++*wordsize,"absorb3 "); 00237 00238 strcpy(ntTags+ntVars++*wordsize,"KEneu4 "); 00239 strcpy(ntTags+ntVars++*wordsize,"neu4cos "); 00240 strcpy(ntTags+ntVars++*wordsize,"neuproc4 "); 00241 strcpy(ntTags+ntVars++*wordsize,"gneu4r "); 00242 strcpy(ntTags+ntVars++*wordsize,"gneu4cos "); 00243 strcpy(ntTags+ntVars++*wordsize,"gneu4phi "); 00244 strcpy(ntTags+ntVars++*wordsize,"gneu4x "); 00245 strcpy(ntTags+ntVars++*wordsize,"gneu4y "); 00246 strcpy(ntTags+ntVars++*wordsize,"gneu4z "); 00247 strcpy(ntTags+ntVars++*wordsize,"vneu4r "); 00248 strcpy(ntTags+ntVars++*wordsize,"vneu4cos "); 00249 strcpy(ntTags+ntVars++*wordsize,"vneu4phi "); 00250 strcpy(ntTags+ntVars++*wordsize,"vneu4x "); 00251 strcpy(ntTags+ntVars++*wordsize,"vneu4y "); 00252 strcpy(ntTags+ntVars++*wordsize,"vneu4z "); 00253 strcpy(ntTags+ntVars++*wordsize,"tneu4 "); 00254 strcpy(ntTags+ntVars++*wordsize,"nsteps4 "); 00255 strcpy(ntTags+ntVars++*wordsize,"absorb4 "); 00256 00257 strcpy(ntTags+ntVars++*wordsize,"KEneu5 "); 00258 strcpy(ntTags+ntVars++*wordsize,"neu5cos "); 00259 strcpy(ntTags+ntVars++*wordsize,"neuproc5 "); 00260 strcpy(ntTags+ntVars++*wordsize,"gneu5r "); 00261 strcpy(ntTags+ntVars++*wordsize,"gneu5cos "); 00262 strcpy(ntTags+ntVars++*wordsize,"gneu5phi "); 00263 strcpy(ntTags+ntVars++*wordsize,"gneu5x "); 00264 strcpy(ntTags+ntVars++*wordsize,"gneu5y "); 00265 strcpy(ntTags+ntVars++*wordsize,"gneu5z "); 00266 strcpy(ntTags+ntVars++*wordsize,"vneu5r "); 00267 strcpy(ntTags+ntVars++*wordsize,"vneu5cos "); 00268 strcpy(ntTags+ntVars++*wordsize,"vneu5phi "); 00269 strcpy(ntTags+ntVars++*wordsize,"vneu5x "); 00270 strcpy(ntTags+ntVars++*wordsize,"vneu5y "); 00271 strcpy(ntTags+ntVars++*wordsize,"vneu5z "); 00272 strcpy(ntTags+ntVars++*wordsize,"tneu5 "); 00273 strcpy(ntTags+ntVars++*wordsize,"nsteps5 "); 00274 strcpy(ntTags+ntVars++*wordsize,"absorb5 "); 00275 // 00276 IdNeu = 1000*ntId+id++; 00277 sprintf(Name,"%s%s",ntName,", neutron"); 00278 NvarNeu = ntVars; 00279 // 00280 hbookn_(IdNeu,Name,NvarNeu,dirName,1000,ntTags, 00281 strlen(Name),strlen(dirName),wordsize,NvarNeu); 00282 // 00283 // reset 00284 ntVars = 0; 00285 // 00286 // book the pmt detector variables, nt=1006 00287 // 00288 strcpy(ntTags+ntVars++*wordsize,"ncosbins "); 00289 strcpy(ntTags+ntVars++*wordsize,"nphibins "); 00290 00291 strcpy(ntTags+ntVars++*wordsize,"phposvis "); 00292 strcpy(ntTags+ntVars++*wordsize,"phposrec "); 00293 strcpy(ntTags+ntVars++*wordsize,"phneuvis "); 00294 strcpy(ntTags+ntVars++*wordsize,"phneurec "); 00295 strcpy(ntTags+ntVars++*wordsize,"phtlvis "); 00296 strcpy(ntTags+ntVars++*wordsize,"phtlrec "); 00297 strcpy(ntTags+ntVars++*wordsize,"phbivis "); 00298 strcpy(ntTags+ntVars++*wordsize,"phbirec "); 00299 strcpy(ntTags+ntVars++*wordsize,"phcfvis "); 00300 strcpy(ntTags+ntVars++*wordsize,"phcfrec "); 00301 strcpy(ntTags+ntVars++*wordsize,"phmuovis "); 00302 strcpy(ntTags+ntVars++*wordsize,"phmuorec "); 00303 strcpy(ntTags+ntVars++*wordsize,"phelevis "); 00304 strcpy(ntTags+ntVars++*wordsize,"phelerec "); 00305 00306 strcpy(ntTags+ntVars++*wordsize,"npepos "); 00307 strcpy(ntTags+ntVars++*wordsize,"attenpos "); 00308 strcpy(ntTags+ntVars++*wordsize,"sfactpos "); 00309 strcpy(ntTags+ntVars++*wordsize,"npeneu "); 00310 strcpy(ntTags+ntVars++*wordsize,"attenneu "); 00311 strcpy(ntTags+ntVars++*wordsize,"sfactneu "); 00312 strcpy(ntTags+ntVars++*wordsize,"npetl "); 00313 strcpy(ntTags+ntVars++*wordsize,"attentl "); 00314 strcpy(ntTags+ntVars++*wordsize,"sfacttl "); 00315 strcpy(ntTags+ntVars++*wordsize,"npebi "); 00316 strcpy(ntTags+ntVars++*wordsize,"attenbi "); 00317 strcpy(ntTags+ntVars++*wordsize,"sfactbi "); 00318 strcpy(ntTags+ntVars++*wordsize,"npecf "); 00319 strcpy(ntTags+ntVars++*wordsize,"attencf "); 00320 strcpy(ntTags+ntVars++*wordsize,"sfactcf "); 00321 strcpy(ntTags+ntVars++*wordsize,"npemuo "); 00322 strcpy(ntTags+ntVars++*wordsize,"attenmuo "); 00323 strcpy(ntTags+ntVars++*wordsize,"sfactmuo "); 00324 strcpy(ntTags+ntVars++*wordsize,"npeele "); 00325 strcpy(ntTags+ntVars++*wordsize,"attenele "); 00326 strcpy(ntTags+ntVars++*wordsize,"sfactele "); 00327 00328 strcpy(ntTags+ntVars++*wordsize,"ngampos "); 00329 strcpy(ntTags+ntVars++*wordsize,"ngamneu "); 00330 strcpy(ntTags+ntVars++*wordsize,"ngamtl "); 00331 strcpy(ntTags+ntVars++*wordsize,"ngambi "); 00332 strcpy(ntTags+ntVars++*wordsize,"ngamcf "); 00333 strcpy(ntTags+ntVars++*wordsize,"ngammuo "); 00334 strcpy(ntTags+ntVars++*wordsize,"ngamele "); 00335 strcpy(ntTags+ntVars++*wordsize,"fpmtpos "); 00336 strcpy(ntTags+ntVars++*wordsize,"fpmtneu "); 00337 strcpy(ntTags+ntVars++*wordsize,"fpmttl "); 00338 strcpy(ntTags+ntVars++*wordsize,"fpmtbi "); 00339 strcpy(ntTags+ntVars++*wordsize,"fpmtcf "); 00340 strcpy(ntTags+ntVars++*wordsize,"fpmtmuo "); 00341 strcpy(ntTags+ntVars++*wordsize,"fpmtele "); 00342 00343 strcpy(ntTags+ntVars++*wordsize,"dposx "); 00344 strcpy(ntTags+ntVars++*wordsize,"dposy "); 00345 strcpy(ntTags+ntVars++*wordsize,"dposz "); 00346 strcpy(ntTags+ntVars++*wordsize,"dneux "); 00347 strcpy(ntTags+ntVars++*wordsize,"dneuy "); 00348 strcpy(ntTags+ntVars++*wordsize,"dneuz "); 00349 // 00350 IdDet = 1000*ntId+id++; 00351 sprintf(Name,"%s%s",ntName,", detector"); 00352 NvarDet = ntVars; 00353 // 00354 hbookn_(IdDet,Name,NvarDet,dirName,1000,ntTags, 00355 strlen(Name),strlen(dirName),wordsize,NvarDet); 00356 // 00357 // define some histograms 00358 // Last argument is need for c++ ==> f strings 00359 // 00360 Offset = 10000*(ntId-1); 00361 00362 hbook1_(Offset+101,"electron PH /PMT/gamma/electron",200,0,0.1,0, 00363 strlen("electron PH /PMT/gamma/electron")); 00364 hbook1_(Offset+102,"positron PH /PMT/gamma/positron",200,0,0.1,0, 00365 strlen("positron PH /PMT/gamma/positron")); 00366 hbook1_(Offset+103,"neutron PH /PMT/gamma/neutron",200,0,0.1,0, 00367 strlen("neutron PH /PMT/gamma/neutron")); 00368 hbook1_(Offset+104,"thallium PH /PMT/gamma/decay",200,0,0.1,0, 00369 strlen("thallium PH /PMT/gamma/decay")); 00370 hbook1_(Offset+105,"bismuth PH /PMT/gamma/decay",200,0,0.1,0, 00371 strlen("bismuth PH /PMT/gamma/decay")); 00372 hbook1_(Offset+106,"cf prompt gamma PH /PMT/gamma/decay",200,0,0.1,0, 00373 strlen("cf prompt gamma PH /PMT/gamma/decay")); 00374 hbook1_(Offset+107,"muon PH /PMT/gamma/muon",200,0,0.1,0, 00375 strlen("muon PH /PMT/gamma/muon")); 00376 // 00377 hbook1_(Offset+201,"electron time /PMT/gamma/electron",250,0,5e8,0, 00378 strlen("electron time /PMT/gamma/electron")); 00379 hbook1_(Offset+202,"positron time /PMT/gamma/positron",250,0,50,0, 00380 strlen("positron time /PMT/gamma/positron")); 00381 hbook1_(Offset+203,"neutron time /PMT/gamma/neutron",250,0,5e5,0, 00382 strlen("neutron time /PMT/gamma/neutron")); 00383 hbook1_(Offset+204,"thallium time /PMT/gamma/decay",250,0,1e6,0, 00384 strlen("thallium time /PMT/gamma/decay")); 00385 hbook1_(Offset+205,"bismuth time /PMT/gamma/decay",250,0,1e6,0, 00386 strlen("bismuth time /PMT/gamma/decay")); 00387 hbook1_(Offset+206,"cf prompt gamma time /PMT/gamma/decay",250,0,50,0, 00388 strlen("cf prompt gamma time /PMT/gamma/decay")); 00389 hbook1_(Offset+207,"muon time /PMT/gamma/muon",250,0,50,0, 00390 strlen("muon time /PMT/gamma/muon")); 00391 // 00392 hbook1_(Offset+301,"number of gamma/electron",10,0,10,0, 00393 strlen("number of gamma/electron")); 00394 hbook1_(Offset+302,"number of gamma/positron",10,0,10,0, 00395 strlen("number of gamma/positron")); 00396 hbook1_(Offset+303,"number of gamma/neutron",20,0,20,0, 00397 strlen("number of gamma/neutron")); 00398 hbook1_(Offset+306,"number of prompt Cf gamma/decay",20,0,20,0, 00399 strlen("number of prompt CF gamma/decay")); 00400 // 00401 hbook1_(Offset+401,"electron gamma energy",100,0,10,0, 00402 strlen("electron gamma energy")); 00403 hbook1_(Offset+402,"positron gamma energy",100,0,10,0, 00404 strlen("positron gamma energy")); 00405 hbook1_(Offset+403,"neutron gamma energy",100,0,10,0, 00406 strlen("neutron gamma energy")); 00407 hbook1_(Offset+406,"prompt Cf gamma energy",100,0,10,0, 00408 strlen("prompt Cf gamma energy")); 00409 // 00410 hbook1_(Offset+501,"electron gamma time",250,0,5e8,0, 00411 strlen("electron gamma time")); 00412 hbook1_(Offset+502,"positron gamma time",250,0,5,0, 00413 strlen("positron gamma time")); 00414 hbook1_(Offset+503,"neutron gamma time",250,0,5e5,0, 00415 strlen("neutron gamma time")); 00416 hbook1_(Offset+506,"prompt Cf gamma time",250,0,5,0, 00417 strlen("prompt Cf gamma time")); 00418 // 00419 hbook1_(Offset+601,"electron gamma radius",300,0,300,0, 00420 strlen("electron gamma radius")); 00421 hbook1_(Offset+602,"positron gamma radius",300,0,300,0, 00422 strlen("positron gamma radius")); 00423 hbook1_(Offset+603,"neutron gamma radius",300,0,300,0, 00424 strlen("neutron gamma radius")); 00425 hbook1_(Offset+606,"prompt Cf gamma radius",300,0,300,0, 00426 strlen("prompt Cf gamma radius")); 00427 // 00428 hbook1_(Offset+701,"electron gamma cos(theta)",20,-1.1,1.1,0, 00429 strlen("electron gamma cos(theta)")); 00430 hbook1_(Offset+702,"positron gamma cos(theta)",20,-1.1,1.1,0, 00431 strlen("positron gamma cos(theta)")); 00432 hbook1_(Offset+703,"neutron gamma cos(theta)",20,-1.1,1.1,0, 00433 strlen("neutron gamma cos(theta)")); 00434 hbook1_(Offset+706,"prompt Cf gamma cos(theta)",20,-1.1,1.1,0, 00435 strlen("prompt Cf gamma cos(theta)"));} |