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.