| Rev | Line | |
|---|
| [1350] | 1 | MACRO NTUPLE
|
|---|
| 2 |
|
|---|
| 3 | * MERGES THE ASCII DATA (hbk/file.txt) INTO A NTUPLE (hbk/file.hbk)
|
|---|
| 4 |
|
|---|
| 5 | * FILE : ASCII INPUT FILE ( taken from the shell environment )
|
|---|
| 6 | * HBK : FILE CONTAINING THE RESULTING NTUPLE
|
|---|
| 7 | * NEVENT : NUMBER OF EVENTS TO READ, DEFAULT IN PAW IS 1E6
|
|---|
| 8 |
|
|---|
| 9 | * NTUPLE VARIABLES:
|
|---|
| 10 | * u,v : surface parameters (eg. angle phi, z-position. See PrimaryGeneratorAction)
|
|---|
| 11 | * x,y,z: generated (true) surface point in cartesian coordinates
|
|---|
| 12 | * px,py,pz: vertex position in cartesian coordinates
|
|---|
| 13 | * dist: distance between vertex and intersection
|
|---|
| 14 | * vx,vy,vz: momentum of particle at vertex
|
|---|
| 15 | * rx,ry,rz: reconstructed intersection
|
|---|
| 16 | * delta : distance between reconstructed and true intersection
|
|---|
| 17 | * theta: angle between momentum and surface normal at intersection
|
|---|
| 18 |
|
|---|
| 19 | FILE = $ENV(arg1)
|
|---|
| 20 | HBK = $ENV(arg2)
|
|---|
| 21 | NEVENT = $ENV(arg3)
|
|---|
| 22 |
|
|---|
| 23 | *file = hbk/surf.txt
|
|---|
| 24 | *hbk = hbk/surf.hbk
|
|---|
| 25 | *nevent = 100000
|
|---|
| 26 |
|
|---|
| 27 | MESS PAW: Process [NEVENT] events.
|
|---|
| 28 |
|
|---|
| 29 | HI/FILE 1 [HBK] 1024 N
|
|---|
| 30 | NT/CRE 1 'merge' 17 //lun1 ! u v x y z px py pz dist vx vy vz rx ry rz delta theta
|
|---|
| 31 | NT/READ 1 [FILE] ! ! [NEVENT]
|
|---|
| 32 |
|
|---|
| 33 | HROUT 0
|
|---|
| 34 | CLOSE 1
|
|---|
| 35 |
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
