]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1996/03/08 17:27:17 mclareni | |
6 | * Isajet714 | |
7 | * | |
8 | * | |
9 | ||
10 | ISAJET ZEBRA INTERFACE | |
11 | F.E. Paige and S.D. Protopopescu | |
12 | Physics Department | |
13 | Brookhaven National Laboratory | |
14 | Upton, NY 11973 | |
15 | March 30, 1991 | |
16 | ||
17 | ||
18 | [1] Overview | |
19 | ============ | |
20 | ||
21 | The ZEBRA version of ISAJET contains subroutines which replace the | |
22 | usual output package for ISAJET with one using the CERN package ZEBRA. | |
23 | This version provides some features which are not part of the standard | |
24 | ISAJET output. Momenta and masses of the partons are rescaled from the | |
25 | final particles so that they match exactly (this is not so in ISAJET). | |
26 | Resonances decaying strongly are not part of the particle list while | |
27 | particles with heavy quarks decaying weakly are part of a vertex list with | |
28 | the corresponding vertex position and appropriate links to the particles | |
29 | belonging to that vertex. Thus there is no need to generate secondary | |
30 | vertices for charmed or other heavy quark particles and the event structure | |
31 | makes it easy to find which particles come from the decay of a heavy quark. | |
32 | The association of particles with parent partons is also made easy by | |
33 | reference links. For simple calorimeter studies banks are also provided | |
34 | with energy sums over calorimeter cells and banks containing only the | |
35 | leptons (see below). If only those are written instead of particle banks | |
36 | the resulting files are more compact and require considerably less CPU time | |
37 | to analyze (albeit the information is rather limited). In addition jets are | |
38 | calculated by using a cone algorithm on the partons, these can be useful | |
39 | for studying efficiencies of jet finding algorithms and comparing results | |
40 | to theoretical jets. The results of parton jets are in Zebra banks PJET. | |
41 | ||
42 | NOTE: The present version can only handle one run at a time, so you | |
43 | cannot supply a command file with many runs for now. This problem will be | |
44 | fixed in a future release. | |
45 | ||
46 | [2] Creating and running a job | |
47 | ============================== | |
48 | ||
49 | There are command files in BNLCL6::$2$DUA14:[ISAJET.ISALIBRARY] to make | |
50 | EXE files for event generation and analysis. | |
51 | ||
52 | To make an EXE generating ISAJET events which are then written to | |
53 | a file in ZEBRA format users should run command file: | |
54 | ||
55 | $ @isa$area:MAKE_ISAJET_ZEBRA | |
56 | ||
57 | this starts from ISAJET.PAM and makes ISAJET_ZEBRA.FOR, ISAJET_ZEBRA.OLB | |
58 | and ISAJET_ZEBRA.EXE. It requires that ISA$AREA de defined as the area | |
59 | where ISAJET.PAM and command files reside. It also requires that CERN$AREA | |
60 | be defined as the area for the CERN libraries. The above command file | |
61 | assumes symbol YPATCHY is defined. If running on a machine other than a VAX | |
62 | looking at the command file should give a pretty good idea of what needs to | |
63 | be done. | |
64 | ||
65 | If you already have ISAJET_ZEBRA.OLB you can run command file | |
66 | ||
67 | $ @isa$area:LINK_ISAJET_ZEBRA | |
68 | ||
69 | to make ISAJET_ZEBRA.EXE, this program runs with an interactive | |
70 | dialog to generate ISAJET instructions. You can supply your own | |
71 | ISARUN main program if you wish to bypass the dialog, see command | |
72 | file MAKE_ISAJET.COM for the simplest example. | |
73 | ||
74 | User can provide any of 3 logical functions to reduce the | |
75 | size of output files: | |
76 | 1) EDIT(I) = if false event will not be written out | |
77 | I is event number | |
78 | Zebra banks have not yet been generated, so | |
79 | user must work from standard ISAJET common blocks | |
80 | 2) ISEDIT() = same as EDIT but Zebra banks have been filled | |
81 | 3) PEDIT(ID,P) = if false particle will not be included in | |
82 | the list of particles written out | |
83 | ID= particle ID, P(4)= particle momentum vector | |
84 | ||
85 | ||
86 | To generate a user program that reads an ISAJET output file in ZEBRA | |
87 | format users should run | |
88 | ||
89 | $ @isa$area:MAKE_ISZRUN | |
90 | this makes a sample user analysis subroutine ISZUSR and an ISZRUN.EXE | |
91 | program. ISZUSR gives an example of how to use the utility subroutines | |
92 | mentioned in next section for analysis. | |
93 | ||
94 | If user already has an ISZUSR subroutine he can generate ISZRUN.EXE by | |
95 | ||
96 | $ @isa$area:LINK_ISZRUN | |
97 | ||
98 | ||
99 | [3] Utilities | |
100 | ============= | |
101 | ||
102 | A number of utilities are provided for the user to make it easier for | |
103 | him to use the information in ZEBRA form without him having to learn much | |
104 | about ZEBRA. It is nevertheless recommended that the user familiarize | |
105 | himself with ZEBRA by reading the first chapter of the ZEBRA manual. Also | |
106 | to make full use of the power of ZEBRA users may eventually want to | |
107 | familiarize themselves with some of the utilities provided with it such as | |
108 | MZDROP and ZSORT. | |
109 | ||
110 | o For every ZEBRA bank of name ISxx (most ISAJET ZEBRA banks start with | |
111 | IS) there exist: | |
112 | 1) GZISxx : an integer function returning pointer to first bank in | |
113 | the linear structure | |
114 | 2) GTISxx : a subroutine returning the contents of the bank | |
115 | 3) PRISxx : a subroutine generating a printout of the banks with name ISxx | |
116 | ||
117 | o To generate a full dump: | |
118 | CALL PRTEVZ(unit #) | |
119 | This subroutine calls all the PRISxx subroutines. | |
120 | ||
121 | o Users can provide a primary vertex for each event by: | |
122 | CALL SETVTX(XYZ) | |
123 | where XYZ is a 3-vector with values x,y,z | |
124 | To supply a function for the vertex user should supply subroutine | |
125 | INTVTX(XYZ) which returns an XYZ vector with each call. The | |
126 | default returns 0,0,0. | |
127 | ||
128 | o The parameters for generating PJET banks can be controlled by | |
129 | CALL PJPSET(N,ALG,ETCUT,DRCUT,MAXIT,IR,MUON) | |
130 | where | |
131 | Default values | |
132 | N = NUMBER OF PJET ALGORITHMS 1 | |
133 | ALG = ALGORITHM TPYE 'CONE' OR 'ANGLE' 'CONE' | |
134 | ETCUT = ET CUT FR PJET JETS 5.0 | |
135 | DRCUT = CONE CUT (IN R OR ANGLE) 0.45 | |
136 | MAXIT = MAXIMUM NUMBER OF ITERATIONS 1 | |
137 | IR = INIT RADIATION FLAG 0=USE IR, 1=NO IR 0 | |
138 | MUON = MUON FLAG 0=NO MUONS, 1=use MUONS in PJET 0 | |
139 | ||
140 | ||
141 | Other subroutines of interest are: | |
142 | o ISZJET find jets using ISCL banks | |
143 | o ISZUSR example of a user subroutine for analysis | |
144 | ||
145 | ||
146 | ||
147 | [4] Other Subroutines | |
148 | ===================== | |
149 | ||
150 | We list here the subroutines that are part of ISAJET_ZEBRA: | |
151 | ||
152 | GENVTX generate secondary vertices (for heavy quark decays) | |
153 | ISABFL fill ISAB bank (begin of run bank) | |
154 | ISACFL fill ISAC and ISCL banks (calorimeter description and | |
155 | calorimeter cells with non-zero energy deposited). | |
156 | Entry point ISACIN calculates all numbers needed for | |
157 | calorimeter description. | |
158 | ISAEFL fill ISAE bank (main header) | |
159 | ISAFFL fill ISAF bank (end of run bank) | |
160 | ISAJFL fill ISAJ banks (primary jets) | |
161 | ISALFL fill ISAL banks (leptons) | |
162 | ISAPFL fill ISV1 and ISP1 banks (vertices and particles) | |
163 | ISAQFL fill ISAQ banks (secondary partons, both initial and final) | |
164 | ISCMFL fill ISCM (copy of command instructions) | |
165 | ISMEAR generates smeared calorimeter data (modifies ISCL banks) | |
166 | PJETFL fill PJET banks (jets obtained from partons) | |
167 | QRECAL to recalculate momenta and masses of partons starting | |
168 | from the final particles | |
169 | ||
170 | ||
171 | ||
172 | ||
173 | ||
174 | [5] Brief description of Zebra banks | |
175 | ==================================== | |
176 | ||
177 | The beginning and end of run records have Zebra banks named ISAB and | |
178 | ISAF while the whole event hangs from Zebra bank ISAE. A complete | |
179 | description of the banks is given in section [6]. | |
180 | ||
181 | ISAB data describing event generation (in beginning record) | |
182 | ISCM copy of command file used for event generation( " ) | |
183 | ISAF cumulative results from run (in end record) | |
184 | ISAE general event information such as weight, type, etc. | |
185 | ISAJ primary jet momenta and masses | |
186 | ISAQ final and initial partons (momenta and masses) | |
187 | ISV1 primary and short lived vertices (heavy quarks), includes | |
188 | id,mass and 4-momenta of decaying particle | |
189 | ISP1 particles associated with ISV1 (connected via a structural link) | |
190 | id,mass and 4-momenta | |
191 | ISAC simple calorimeter description (no depth) | |
192 | ISCL energy deposition and cell description in calorimeter cell | |
193 | ISMR description of smearing parameters | |
194 | ISAL stable leptons (including neutrinos):id,mass and 4-momenta | |
195 | with reference links to calorimeter cell and ISP1 if they | |
196 | exist | |
197 | ISJT jets found by subroutine ISZJET. | |
198 | PJHD header for PJET banks | |
199 | PJET jets found adding over partons | |
200 | PJPT pointers to parton banks (ISAQ) contributing to a jet in PJET | |
201 | ||
202 | Whether banks ISV1,ISP1,ISAC,ISCL and ISAL are written out is | |
203 | optional. The ISARUN job asks user to chose an option ISAP, ISAC or ISAL. | |
204 | If ISAP is selected then ISV1 and ISP1 are written out, if ISAC then ISAC | |
205 | and ISCL are written out and if ISAL then ISAL is written out. Any | |
206 | combination can be selected, the instruction ISAPISACISAL causes all banks | |
207 | to be written out. | |
208 | ||
209 | The banks ISJT are created by calling ISZJET, the banks ISAC and ISCL | |
210 | must be available. One can create the banks ISAC,ISCL and/or ISAL starting | |
211 | from a file containing ISV1 and ISP1 by calling ISACFL and/or ISALFL for | |
212 | each event. | |
213 | ||
214 | ||
215 | ||
216 | [6] Zebra banks documentation | |
217 | ============================= | |
218 | ||
219 | ||
220 | +------------------+ | |
221 | | ISAE | | |
222 | | event descriptor | | |
223 | | | | |
224 | +------------------+ | |
225 | struc. link ||| | | | |
226 | ||| | +------------------------------+ | |
227 | ||| +-----------------------------+ | | |
228 | +-------+------------+| | | | |
229 | | | | | | | |
230 | ______ ______ ______ ______ | | | |
231 | | ISAJ \ | ISAQ \ | ISV1 | | ISV1 \ | | | |
232 | | -1 > | -2 > | -3 |--| (-3) > | | | |
233 | |______/ |______/ |______| |______/ | | | |
234 | : ^ | :...<..: | | | | |
235 | : : | [-2] | | | | |
236 | [PJET -1].<.. : | | | | | |
237 | : [-2] ______ ______ | | | |
238 | :.....| ISP1 \ | ISP1 \ | | | |
239 | [PJET -3] ..<........| -1 > | -1 > | | | |
240 | |______/ |______/ | | | |
241 | : |[-1] : : : | | | |
242 | [-4]: v <....: : :[-4] | | | |
243 | : free <......: : | | | |
244 | v v | | | |
245 | ISV2 ISV2 | | | |
246 | (end vertex) (end vertex) | | | |
247 | | | | |
248 | +------------+---------------------+ | | |
249 | | | | | |
250 | | | | | |
251 | _________ ______ ______ | |
252 | | ISAC | | ISAL \ | PJHD \ | |
253 | | -6 | | -7 > | -8 > | |
254 | |_________| |______/ |______/ | |
255 | | | | |
256 | +-------+---------+--+ ______ | |
257 | | | | [ISAQ -1]..>.....| PJET \ | |
258 | ______ ______ ______ | -1 > | |
259 | | ISCL \ | ISMR \ | ISJT \ |______/ | |
260 | | -1 > | -2 > | -3 > | | |
261 | |______/ |______/ |______/ ______ | |
262 | | PJPT | | |
263 | [ISAQ]...<...| -1 | | |
264 | |______| | |
265 | ||
266 | +-----------+ +-----------+ | |
267 | | ISAB | | ISAF | | |
268 | | begin run | | end run | | |
269 | | | | | | |
270 | +-----------+ +-----------+ | |
271 | | | |
272 | +----------+ | |
273 | | ISCM | | |
274 | | -1 | | |
275 | +----------+ | |
276 | ||
277 | ||
278 | ||
279 | C======================================================================= | |
280 | C | |
281 | C Bank name: ISAC | |
282 | C Author : SDP | |
283 | C Date : June 14,1989 | |
284 | C | |
285 | C Calorimeter descriptor | |
286 | C | |
287 | C LQ Q/IQ | |
288 | C ___________________________________________________________________________ | |
289 | C -3 struct. to ISJT (jets found using ISZJET) | |
290 | C -2 struct. to smearing parameters (ISMR) | |
291 | C -1 struct. to calorimeter cells (ISCL) | |
292 | C 0 next to next ISAC (only one at present) | |
293 | C +1 up to ISAE | |
294 | C +2 origin to ISAE for first and previous ISAC for others | |
295 | C ............................................................................ | |
296 | C I-5 calorimeter number | |
297 | C -4 bank name, 'ISAC' | |
298 | C -3 NL=2 | |
299 | C -2 NS=0 | |
300 | C -1 ND=6 | |
301 | C 0 STATUS | |
302 | C +1 I number of phi cells | |
303 | C 2 I " eta " | |
304 | C 3 F phi cell size | |
305 | C 4 F eta " " | |
306 | C 5 F minimum eta | |
307 | C 6 F maximum eta | |
308 | C======================================================================= | |
309 | ||
310 | ||
311 | ||
312 | C======================================================================= | |
313 | C | |
314 | C Bank name: ISAE | |
315 | C Author : SDP | |
316 | C Date : June 19,1986 DH add ISAM | |
317 | C | |
318 | C Top level bank and event descriptor for ISAZEB | |
319 | C | |
320 | C LQ Q/IQ | |
321 | C ___________________________________________________________________________ | |
322 | C -8 struct. to PJHD (parton jet structure header) | |
323 | C -7 struct. to lepton bank (ISAL) | |
324 | C -6 struct. to pseudo cal. (ISAC) | |
325 | C -5 struct. to unassociated particles (ISP3) | |
326 | C -4 struct. to long lived vertices (ISV2) | |
327 | C -3 struct. to short " " (ISV1) | |
328 | C -2 struct. to final/initial bank (ISAQ) | |
329 | C -1 struct. to primary jet (ISAJ) | |
330 | C 0 next to next ISAE | |
331 | C +1 up | |
332 | C +2 origin | |
333 | C ............................................................................ | |
334 | C I-5 bank number | |
335 | C -4 bank name, 'ISAE' | |
336 | C -3 NL=7 | |
337 | C -2 NS=7 | |
338 | C -1 ND=18 | |
339 | C 0 STATUS | |
340 | C +1 I event id(1) | |
341 | C 2 I event id(2) | |
342 | C 3 I event number | |
343 | C 4 I reaction type | |
344 | C 5 I number of primary jet banks | |
345 | C 6 I " of stable parton banks (final+initial) | |
346 | C 7 I " of PJET banks | |
347 | C 8 I " of particle banks | |
348 | C 9 I " of vertex banks | |
349 | C 10 I " of lepton banks | |
350 | C 11 F cross section in microbarns | |
351 | C 12 F weight | |
352 | C 13 F effective q**2 | |
353 | C 14 F hard scattering invariant s | |
354 | C 15 F " " " t | |
355 | C 16 F " " " u | |
356 | C 17 D Seed (part 1) | |
357 | C 18 D Seed (part 2) | |
358 | C======================================================================= | |
359 | ||
360 | ||
361 | ||
362 | C======================================================================= | |
363 | C | |
364 | C Bank name: ISAJ | |
365 | C Author : SDP | |
366 | C Date : June 19,1986 | |
367 | C | |
368 | C Primary jet information | |
369 | C | |
370 | C LQ Q/IQ | |
371 | C ___________________________________________________________________________ | |
372 | C L-1 ref. to parent primary (ISAJ) (0 if original) | |
373 | C +0 next to next jet (ISAJ) | |
374 | C +1 up to ISAE | |
375 | C +2 origin to ISAE | |
376 | C ............................................................................ | |
377 | C I-5 Primary jet number | |
378 | C -4 bank name, 'ISAJ' | |
379 | C -3 NL=0 | |
380 | C -2 NS=0 | |
381 | C -1 ND=9 | |
382 | C 0 STATUS | |
383 | C +1 I jet ID | |
384 | C 2 F px | |
385 | C 3 F py | |
386 | C 4 F pz | |
387 | C 5 F p | |
388 | C 6 F mass | |
389 | C 7 F phi | |
390 | C 8 F theta | |
391 | C 9 F eta | |
392 | C======================================================================= | |
393 | ||
394 | ||
395 | ||
396 | C======================================================================= | |
397 | C | |
398 | C Bank name: ISAL | |
399 | C Author : SDP | |
400 | C Date : June 19,1986 | |
401 | C | |
402 | C Lepton bank (from primary and short lived vertices) | |
403 | C | |
404 | C LQ Q/IQ | |
405 | C ___________________________________________________________________________ | |
406 | C L-4 ref. to particle bank (ISP1) | |
407 | C -3 ref. to primary jet (ISAJ) | |
408 | C -2 ref. to initial parton (ISAQ) | |
409 | C -1 struct. (free) | |
410 | C 0 next to next ISAL | |
411 | C +1 up to ISAE | |
412 | C +2 origin to ISAE for first and previous ISAL for others | |
413 | C ............................................................................ | |
414 | C I-5 particle number | |
415 | C -4 bank name, 'ISAL' | |
416 | C -3 NL=4 | |
417 | C -2 NS=1 | |
418 | C -1 ND=9 | |
419 | C 0 STATUS | |
420 | C +1 I ISAJET particle ID | |
421 | C 2 F px | |
422 | C 3 F py | |
423 | C 4 F pz | |
424 | C 5 F p | |
425 | C 6 F mass | |
426 | C 7 F phi | |
427 | C 8 F theta | |
428 | C 9 F eta | |
429 | C======================================================================= | |
430 | ||
431 | ||
432 | ||
433 | C======================================================================= | |
434 | C | |
435 | C Bank name: ISAQ | |
436 | C Author : SDP | |
437 | C Date : June 19,1986 | |
438 | C | |
439 | C Final and initial partons | |
440 | C | |
441 | C LQ Q/IQ | |
442 | C ___________________________________________________________________________ | |
443 | C L-2 ref. to PJET (parton structure jets) | |
444 | C -1 ref. to primary jet (ISAJ) | |
445 | C 0 for initial parton | |
446 | C 0 next to next ISAQ | |
447 | C +1 up to ISAE | |
448 | C +2 origin to ISAE for first and previous ISAQ for others | |
449 | C ............................................................................ | |
450 | C I-5 final/initial parton number | |
451 | C -4 bank name, 'ISAQ' | |
452 | C -3 NL=2 | |
453 | C -2 NS=0 | |
454 | C -1 ND=9 | |
455 | C 0 STATUS | |
456 | C +1 I final/initial parton type | |
457 | C 2 F px | |
458 | C 3 F py | |
459 | C 4 F pz | |
460 | C 5 F E | |
461 | C 6 F mass | |
462 | C 7 F phi | |
463 | C 8 F theta | |
464 | C 9 F eta | |
465 | C======================================================================= | |
466 | ||
467 | ||
468 | ||
469 | C======================================================================= | |
470 | C | |
471 | C Bank name: ISCL | |
472 | C Author : SDP | |
473 | C Date : June 19,1986 | |
474 | C | |
475 | C Calorimeter cell bank | |
476 | C | |
477 | C LQ Q/IQ | |
478 | C ___________________________________________________________________________ | |
479 | C L-2 ref. to reconstructed jet (ISJT) | |
480 | C -1 ref. to parent parton (ISAQ) | |
481 | C 0 next to next ISCL | |
482 | C +1 up to ISAC | |
483 | C +2 origin to ISAC for first and previous ISAL for others | |
484 | C ............................................................................ | |
485 | C I-5 cell number (ordinal) | |
486 | C -4 bank name, 'ISCL' | |
487 | C -3 NL=2 | |
488 | C -2 NS=0 | |
489 | C -1 ND=10 | |
490 | C 0 STATUS | |
491 | C +1 I 10000*(phi index) + y index | |
492 | C 2 I 10000*(# of charged trks) + # of gammas | |
493 | C 3 F e.m. energy deposited | |
494 | C 4 F had. " " | |
495 | C 5 F sin(theta) at center of cell | |
496 | C 6 F cos(theta) " | |
497 | C 7 F sin(phi) " | |
498 | C 8 F cos(phi) " | |
499 | C 9 F phi " | |
500 | C 10 F eta " | |
501 | C=========================================================================== | |
502 | ||
503 | ||
504 | ||
505 | ||
506 | C======================================================================= | |
507 | C | |
508 | C Bank name: ISJT | |
509 | C Author : SDP | |
510 | C Date : June 19,1986 | |
511 | C | |
512 | C Jets found by ISZJET | |
513 | C | |
514 | C LQ Q/IQ | |
515 | C ___________________________________________________________________________ | |
516 | C -1 to possible parent parton (ISAQ) | |
517 | C +0 next to next jet (ISJT) | |
518 | C +1 up to ISJT | |
519 | C +2 origin to ISAC | |
520 | C ............................................................................ | |
521 | C I-5 jet number | |
522 | C -4 bank name, 'ISJT' | |
523 | C -3 NL=0 | |
524 | C -2 NS=0 | |
525 | C -1 ND=6 | |
526 | C 0 STATUS | |
527 | C +1 F ET (transverse energy) | |
528 | C 2 F px | |
529 | C 3 F py | |
530 | C 4 F pz | |
531 | C 5 F E | |
532 | C 6 F mass | |
533 | C 7 F phi | |
534 | C 8 F theta | |
535 | C 9 F eta | |
536 | C======================================================================= | |
537 | ||
538 | ||
539 | ||
540 | ||
541 | C======================================================================= | |
542 | C | |
543 | C Bank name: ISP1 | |
544 | C Author : SDP | |
545 | C Date : June 19,1986 | |
546 | C | |
547 | C Particle bank (from primary and short lived vertices) | |
548 | C | |
549 | C LQ Q/IQ | |
550 | C ___________________________________________________________________________ | |
551 | C L-5 ref. to GCAH | |
552 | C L-4 ref. to secondary vertex (ISV2) | |
553 | C -3 ref. to primary jet (ISAJ) | |
554 | C -2 ref. to initial parton (ISAQ) | |
555 | C -1 struct. (free) | |
556 | C 0 next to next ISP1 | |
557 | C +1 up to ISV1 | |
558 | C +2 origin to ISV1 for first and previous ISP1 for others | |
559 | C ............................................................................ | |
560 | C I-5 particle number | |
561 | C -4 bank name, 'ISP1' | |
562 | C -3 NL=4 | |
563 | C -2 NS=1 | |
564 | C -1 ND=9 | |
565 | C 0 STATUS | |
566 | C +1 I ISAJET particle ID | |
567 | C 2 F px | |
568 | C 3 F py | |
569 | C 4 F pz | |
570 | C 5 F E | |
571 | C 6 F mass | |
572 | C 7 F phi | |
573 | C 8 F theta | |
574 | C 9 F eta | |
575 | C======================================================================= | |
576 | ||
577 | ||
578 | ||
579 | C======================================================================= | |
580 | C | |
581 | C Bank Name : PJET | |
582 | C Author : Chip Stewart | |
583 | C Date : 7-NOV-1989 18:10:09.84 | |
584 | C Version : 2.0 | |
585 | C | |
586 | C Bank description : PARTON JET BANK | |
587 | C PJET is constructed by applying a jet algorithm (CONE or | |
588 | C OPENING ANGLE) to either Partons given in the ISAQ bank | |
589 | C or Paricles given in the ISP1 bank. | |
590 | C | |
591 | C LQ Q/IQ | |
592 | C----------------------------------------------------------------------- | |
593 | C +2 Down link to SPARE | |
594 | C -1 Down link to PJPT | |
595 | C 0 Next link to | |
596 | C +1 Up link to PJHD | |
597 | C +2 Origin link to PJHD | |
598 | C....................................................................... | |
599 | C -5 Bank number | |
600 | C -4 Bank name, 'PJET' | |
601 | C -3 NL = 2 | |
602 | C -2 NS = 2 | |
603 | C -1 ND = 12 | |
604 | C 0 Status | |
605 | C +1 I NV = 2 Version Number | |
606 | C 2 F Et ( =Pt in version NV=1 ) | |
607 | C 3 F Px | |
608 | C 4 F Py | |
609 | C 5 F Pz | |
610 | C 6 F E | |
611 | C 7 F mass | |
612 | C 8 F phi | |
613 | C 9 F theta | |
614 | C 10 F eta | |
615 | C 11 I CLASS - (eg. for mapping PJET with RECO ) | |
616 | C 12 I NEXT | |
617 | C======================================================================= | |
618 | C | |
619 | C Definitions: | |
620 | C ----------- | |
621 | C | |
622 | C E - Sum(Ei) over all the Partons/Particles included in the jet. | |
623 | C | |
624 | C Px - Sum(Pxi). | |
625 | C Py - Sum(Pyi). | |
626 | C Pz - Sum(Pzi). | |
627 | C | |
628 | C Et - Sum[Ei*sin(Thetai)] | |
629 | C | |
630 | C Comment: In version NV=1 word 2 was Pt = SQRT ( Px**2 + Py**2 ) | |
631 | C | |
632 | C Mass - SQRT [ E**2 - ( Px**2 + Py**2 + Pz**2 )] | |
633 | C | |
634 | C For Theta, Phi and Eta the following code was used | |
635 | C to calculate them from Px,Py,Pz | |
636 | C | |
637 | C PARAMETER( SMALL = 1.0E-5 ) | |
638 | C Phi=ATAN2(Py,Px+SMALL) ; IF(Phi.LT.0) Phi=Phi+TWOPI | |
639 | C Theta=ACOS((Pz+SMALL)/(SQRT(Px*Px+Py*Py+Pz*Pz)+SMALL)) | |
640 | C Eta=-ALOG(TAN(Theta/2.)+SMALL) | |
641 | C | |
642 | C======================================================================= | |
643 | ||
644 | ||
645 | ||
646 | C======================================================================= | |
647 | C | |
648 | C Bank Name : PJHD | |
649 | C Author : Chip Stewart | |
650 | C Date : 7-NOV-1989 17:57:58.00 | |
651 | C | |
652 | C Bank description : Header for PJET structure | |
653 | C | |
654 | C LQ Q/IQ | |
655 | C----------------------------------------------------------------------- | |
656 | C -2 DOWN spare | |
657 | C -1 DOWN LINK TO PJET | |
658 | C 0 Next link to PJHD | |
659 | C +1 Up link to ISAE | |
660 | C +2 Origin link to ISAE | |
661 | C....................................................................... | |
662 | C -5 Bank number | |
663 | C -4 Bank name, 'PJHD' | |
664 | C -3 NL = 2 | |
665 | C -2 NS = 2 | |
666 | C -1 ND = 8 | |
667 | C 0 Status | |
668 | C 1 I NV=2 Version number | |
669 | C 2 I N_ALG - 1 for DR CONE , 2 for CMS OPEN ANGLE | |
670 | C 3 I NJET - NUMBER OF PARTON JETS | |
671 | C 4 F algorithm parameters DR_CONE_CUT/OPENING_ANGLE_CUT | |
672 | C 5 F algorithm parameters JET_ET_CUT | |
673 | C 6 I algorithm parameters MAX ITERATIONS | |
674 | C 7 I Initial Radiation switch 0-use IR 1-no IR | |
675 | C 8 I MUON switch (IF NV=1 MUON NOT USED) | |
676 | C 0-no MUONS 1-use MUONS | |
677 | C======================================================================= | |
678 | ||
679 | ||
680 | ||
681 | C======================================================================= | |
682 | C | |
683 | C Bank Name : PJPT | |
684 | C Author : Chip Stewart | |
685 | C Date : 15-Dec-1989 | |
686 | C | |
687 | C Bank description : pointers to Partons (ISAQ) in a given jet | |
688 | C | |
689 | C LQ Q/IQ | |
690 | C----------------------------------------------------------------------- | |
691 | C -N+1 Ref Link to parton bank ISAQ for Nth parton in this PJET | |
692 | C -2 Ref Link to parton bank ISAQ for 1st " " " " | |
693 | C -1 not used | |
694 | C 0 Next link to none | |
695 | C +1 Up link to PJET | |
696 | C +2 Origin link to PJET | |
697 | C....................................................................... | |
698 | C -5 Bank number | |
699 | C -4 Bank name, 'PJPT' | |
700 | C -3 NL = N+1 | |
701 | C -2 NS = 0 | |
702 | C -1 ND = 1 - number of data words in PJPT bank = 1 | |
703 | C 0 Status | |
704 | C +1 I Bank version (=1) | |
705 | C======================================================================= | |
706 |