|
1 | 1 |
|
| 2 | +using JetReconstruction |
| 3 | +using JetReconstruction.C_JetReconstruction: jetreconstruction_PseudoJet_init, |
| 4 | + C_ClusterSequence, |
| 5 | + jetreconstruction_ClusterSequence_free_members_, |
| 6 | + C_JetsResult, |
| 7 | + jetreconstruction_JetsResult_free_members_, |
| 8 | + jetreconstruction_jet_reconstruct, |
| 9 | + jetreconstruction_inclusive_jets, |
| 10 | + jetreconstruction_exclusive_jets_njets, |
| 11 | + jetreconstruction_exclusive_jets_dcut |
| 12 | + |
| 13 | +pseudoJets_len = Csize_t(2) |
| 14 | +pseudoJets_ptr = Ptr{PseudoJet}(Libc.malloc(pseudoJets_len * sizeof(PseudoJet))) |
| 15 | + |
| 16 | +jetreconstruction_PseudoJet_init(pseudoJets_ptr, 0.0, 1.0, 2.0, 3.0) |
| 17 | +jetreconstruction_PseudoJet_init(pseudoJets_ptr + sizeof(PseudoJet), 1.0, 2.0, 3.0, 4.0) |
| 18 | + |
| 19 | +strategy = RecoStrategy.Best |
| 20 | +algorithm = JetAlgorithm.CA |
| 21 | +R = 2.0 |
| 22 | + |
| 23 | +clustersequence_ptr = Ptr{C_ClusterSequence{PseudoJet}}(Libc.malloc(sizeof(C_ClusterSequence{PseudoJet}))) |
| 24 | +jetreconstruction_jet_reconstruct(pseudoJets_ptr, pseudoJets_len, algorithm, R, |
| 25 | + RecoStrategy.Best, |
| 26 | + clustersequence_ptr) |
| 27 | + |
| 28 | +results = C_JetsResult{PseudoJet}(C_NULL, 0) |
| 29 | +results_ptr = Base.unsafe_convert(Ptr{C_JetsResult{PseudoJet}}, Ref(results)) |
| 30 | +jetreconstruction_exclusive_jets_njets(clustersequence_ptr, Csize_t(2), results_ptr) |
| 31 | +jetreconstruction_JetsResult_free_members_(results_ptr) |
| 32 | + |
| 33 | +jetreconstruction_exclusive_jets_dcut(clustersequence_ptr, 1.0, results_ptr) |
| 34 | +jetreconstruction_JetsResult_free_members_(results_ptr) |
| 35 | + |
| 36 | +jetreconstruction_inclusive_jets(clustersequence_ptr, 0.0, results_ptr) |
| 37 | +jetreconstruction_JetsResult_free_members_(results_ptr) |
| 38 | + |
| 39 | +jetreconstruction_ClusterSequence_free_members_(clustersequence_ptr) |
| 40 | +Libc.free(pseudoJets_ptr) |
| 41 | +Libc.free(clustersequence_ptr) |
0 commit comments