Changeset 773 for trunk/MagicDoku
- Timestamp:
- 04/26/01 08:10:07 (24 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicDoku/strategy_mc_ana.tex
r772 r773 3 3 \usepackage{magic-tdas} 4 4 5 \setlength{\unitlength}{1.0cm}6 5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 7 6 %% BEGIN DOCUMENT … … 207 206 208 207 \end{enumerate} 208 209 % &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 210 % &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 211 % &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 209 212 210 213 \section{MC work} … … 227 230 \put (1, 11.5){{\sl Air shower programs}} 228 231 \put (1., 10.){\framebox(3.,1.){MMCS}} 229 \put (2., 10.){\vector(0 .,-1.){.9} }232 \put (2., 10.){\vector(0,-1){.9} } 230 233 \put (1., 8.){\framebox(3.,1.){reflector}} 231 \put (2., 8.){\vector(0 .,-1.){.9}}234 \put (2., 8.){\vector(0,-1){.9}} 232 235 233 236 \put (6, 10.){{\sl star background programs}} 234 \put (6.,8.){\framebox(3.,1.){star fieldadder}}235 \put (6., 8.){\line(0 ., -1.){1.5}}236 \put (10.,8.){\framebox(3 .,1.){starresponse}}237 \put (10., 8.){\line(0 ., -1.){1.5}}238 \put (10., 6.5){\vector(-1 .,0.){6.} }237 \put (6.,8.){\framebox(3.,1.){starresponse}} 238 \put (6., 8.){\line(0, -1){1.5}} 239 \put (10.,8.){\framebox(3,1){starfieldadder}} 240 \put (10., 8.){\line(0, -1){1.5}} 241 \put (10., 6.5){\vector(-1,0){6.} } 239 242 240 243 \put (1., 6.){\framebox(3.,1.){camera}} 241 \put (2., 6.){\vector(3 .,-1.){5.} }244 \put (2., 6.){\vector(3,-1){5.} } 242 245 243 246 244 247 245 248 \put (14, 11.5){{\sl real data programs}} 246 \put (14, 8.){\framebox(3 .,1.){MAGIC DAQ}}247 \put (15, 8.){\vector(0 .,-1.){.9} }249 \put (14, 8.){\framebox(3,1){MAGIC DAQ}} 250 \put (15, 8.){\vector(0,-1){.9} } 248 251 \put (14, 6.){\framebox(3.,1.){MERPP}} 249 \put (15., 6.){\vector(-3 .,-1.){5.} }252 \put (15., 6.){\vector(-3,-1){5.} } 250 253 251 254 \put (8.75, 3.7){\oval(4.,1.)} 252 255 \put (7., 3.5){MAGIC root file} 253 \put (8., 3.2){\vector(0 ., -1.){1.0}}256 \put (8., 3.2){\vector(0, -1){1.0}} 254 257 255 258 \put (7, 1.){\framebox(3.,1.){MARS}} 256 259 257 260 \thicklines 258 \put (5., 11.){\line(0 ., -1.){6.5}}259 \put (13., 12.){\line(0 ., -1.){7.5}}261 \put (5., 11.){\line(0, -1){6.5}} 262 \put (13., 12.){\line(0, -1){7.5}} 260 263 261 264 \end{picture} … … 294 297 \label{sec_exist_progs} 295 298 \subsubsection{MMCS - Magic Monte Carlo Simulation} 296 \subsubsection{reflector } 299 300 This program is based on a CORSIKA simulation. It is used to generate 301 air showers for the MAGIC telecope. At the start one run of the 302 program, one has to define the details of the simulation. 303 One can specify the following parameters of an shower 304 (see also figure \ref{pic_shower}): 305 % 306 \begin{enumerate} 307 \item the type of the particles in one run ($PartID$) 308 \item the energy range of the particles ($E_1, E_2$) 309 \item the slope of the Energy spectra 310 \item the range of the shower core on the ground $r_{core}$. 311 \item the direction of the shower by setting the range of 312 zenith angle ($\Theta_1, \Theta_2$) and 313 azimuth angle ($\phi_1, \phi_2$) 314 \end{enumerate} 315 % 316 \begin{figure}[h] 317 \setlength{\unitlength}{1.5cm} 318 \begin{center} 319 \begin{picture}(9.,6.) 320 \put (0., 0.){\framebox(9.,6.){}} 321 322 \thicklines 323 % telescope 324 \put (5., .5){\oval(.75, .75)[t]} 325 \put (3., 1.){{\sl Telesope position}} 326 \put (4.5, 1.){\vector(1, -1){0.5}} 327 % observation level 328 \put (.5, .5){\line(1, 0){8}} 329 \put (.5, .6){{\sl Observation level}} 330 331 % air shower 332 \put (4. , 5.5 ){\line(2, -3){3.3}} 333 \put (4.5, 5.5 ){{\sl Particle Type ($PartId$)}} 334 \put (4.5, 5.25){{\sl Energy ($E_1 < E < E_2$)}} 335 \put (4.5, 5. ) {$\Theta_1 < \Theta < \Theta_2$} 336 \put (4.5, 4.75) {$\phi_1 < \phi < \phi_2$} 337 \put (7.5, .75){{\sl shower core}} 338 339 \thinlines 340 \put (5., .25){\line(1,0){2.3}} 341 \put (6.1, .25){{\sl $r_{Core}$}} 342 343 \put (5., .5){\line(4,3){1.571}} 344 \put (6., 1.35){{\sl $p$}} 345 346 \end{picture} 347 \end{center} 348 \caption {The parameter of an shower that are possible to define 349 at the begin of an MMCS run.} 350 \label{pic_shower} 351 \end{figure} 352 Other parameters, that will be important in the analysis later, 353 can be calculated. I.e. the impact parameter $p$ is defined by 354 the direction 355 of the shower ($\Theta, \phi$) and the core position 356 ($x_{core}, y_{core}$). 357 358 The program MMCS will track the whole shower development 359 through the atmosphere. All the cerenkov particles that hit a 360 sphere around the telesope (in the figure \ref{pic_shower} 361 drawn as the circle around the telecope position) are stored 362 on disk. It is important to recognize, that up to now no 363 information of the pointing of the telescope was taking into 364 account. 365 This cerenkov photons are the input for the next program, 366 called reflector. 367 368 369 \subsubsection{reflector} 370 371 The aim of the reflector program is the 372 tracking of the cerenkov photons to the camera 373 of the MAGIC telescope. So this 374 is the point where we introduce a specific pointing of 375 the telescope ($\Theta_{MAGIC}, \phi_{MAGIC}$). 376 For all cerenkov photons the program 377 tests if the mirrors are hitten, calculates the 378 probability for the reflection and tracks them to the 379 mirror plane. All the photons that are hitting the 380 camera are written to disk (*.rfl) 381 with their important parameters 382 ($x_{camera}, y_{camera}, \lambda, t_{arrival}$). 383 These parameters are the input from the shower simulation 384 for the next program in the 385 MC simulation chain, the camera program. 386 297 387 \subsubsection{camera} 388 389 The camera program simulates the behaviour of the 390 PMTs and the electronic of the trigger and FAC system. 391 For each photon out of the reflector file (*.rfl) the 392 camera program calculates the probability to generate 393 an photo electron out of the photo cathode. If a photo 394 electrons was ejected, this will create a signal in the 395 trigger and FADC system of the hitted pixel. 396 You have to specify the 397 parameter of the signal shaping 398 (shape, Amplitude, FWHM of signal) 399 at the beginning of the 400 camera, seperatly for the trigger and the FADC system. 401 All signal from all photoelectrons are superimposed for 402 each pixel. As an example you can see the output of 403 the trigger and FADC system in figure \ref{fig_trigger_fadc}. 404 \begin{figure}[h] 405 406 \caption{The response of one shower from the trigger (left) and 407 fadc system (right).} 408 \label{fig_trigger_fadc} 409 \end{figure} 410 411 All these analog signals going into the trigger system are used 412 to check if for a given event a trigger signal was generated or 413 not. But before the start of the camera program on also has to 414 set a few parameters of the trigger system like: 415 \begin{itemize} 416 \item diskriminator threshold 417 \item mulitplicity 418 \item topology 419 \end{itemize} 420 With this set of parameter the camera program will analyse 421 if one event has triggered. For the triggered event all the FADC 422 content will be writen on the file (*.root). In addition all the 423 information about the event ($PartID, E, \Theta$,...) and 424 information of trigger (FirstLevel, SecondLevel, ..) are also 425 be written to the file. 426 427 One of the nice features of the camera program is the possiblity 428 so simulate the NSB, the diffuse and the star light part of it. 429 But before doing this, on has to start other programs 430 (called starresponse and starfieldadder) that are describe 431 below. 432 433 \subsubsection{starresponse} 434 435 This program will simulate the analog response for stars of 436 a given brightness $B$. 437 438 298 439 \subsubsection{starfieldadder} 299 \subsubsection{starresponse} 440 441 442 443 300 444 301 445
Note:
See TracChangeset
for help on using the changeset viewer.