1. Calo jets
In your current code "jet_particles_list_calo" never contains CASTOR data, since it is only filled in the loop over CaloTowersp4
btw. in the code version from 21.9.2020 you commented this code in lines 187 to 191 -> thus currently calo jets will NOT work.
After you fixed the above, ALSO add CASTOR-p4 to "jet_particles_list_calo" before you reconstruct calo jets.
2. Jet matching
The jet-matching procedure we need is explained on pages 85 and 86 of:
https://publikationen.bibliothek.kit.edu/1000074295
It works in the following steps
- reject all calo jets with pt<1GeV
- reject all gen jets with pt<1GeV
- select only calo jets which include at least one CASTOR tower (plus additional HF towers)
- select only gen jets between "-5.2 < eta<-6.6+R"
- select the single, highest-pT calo jet only
- find a gen jet within delta-R=0.2 or calo jet
- if found: this gen jet and the calo jet are matched. This event is signal.
- if not found: event is background.
