detectors/DIRC/DIRC.cc

0001 // Copyright 2022,2023 Christopher Dilks, Nilanga Wickramaarachchi
0002 // Subject to the terms in the LICENSE file found in the top-level directory.
0003 //
0004 //
0005
0006 // SPDX-License-Identifier: LGPL-3.0-or-later
0007 // Copyright (C) 2024, Dmitry Kalinkin
0008
0009 #include <DD4hep/Detector.h>
0010 #include <JANA/JApplication.h>
0011 #include <stddef.h>
0012 #include <algorithm>
0013 #include <gsl/pointers>
0014 #include <memory>
0015 #include <stdexcept>
0016 #include <utility>
0017 #include <vector>
0018
0019 // algorithm configurations
0020 #include "algorithms/digi/PhotoMultiplierHitDigiConfig.h"
0021 #include "algorithms/pid_lut/PIDLookupConfig.h"
0022 #include "extensions/jana/JOmniFactoryGeneratorT.h"
0023 // factories
0024 #include "global/digi/PhotoMultiplierHitDigi_factory.h"
0025 #include "global/pid_lut/PIDLookup_factory.h"
0026 #include "services/geometry/dd4hep/DD4hep_service.h"
0027
0028
0029 extern "C" {
0030   void InitPlugin(JApplication *app) {
0031     InitJANAPlugin(app);
0032
0033     using namespace eicrecon;
0034
0035     // configuration parameters ///////////////////////////////////////////////
0036
0037     // digitization
0038     PhotoMultiplierHitDigiConfig digi_cfg;
0039     digi_cfg.seed            = 5; // FIXME: set to 0 for a 'unique' seed, but
0040                                   // that seems to delay the RNG from actually randomizing
0041     digi_cfg.hitTimeWindow   = 20.0; // [ns]
0042     digi_cfg.timeResolution  = 1/16.0; // [ns]
0043     digi_cfg.speMean         = 80.0;
0044     digi_cfg.speError        = 16.0;
0045     digi_cfg.pedMean         = 200.0;
0046     digi_cfg.pedError        = 3.0;
0047     digi_cfg.enablePixelGaps = false;
0048     digi_cfg.safetyFactor    = 1.0;
0049     digi_cfg.quantumEfficiency.clear();
0050
0051     std::vector<double> sensor_qe{
0052       0,    0,    14.0, 14.8, 14.5, 14.9, 14.4, 14.2, 13.9, 14.6, 15.2, 15.7, 16.4, 16.9, 17.5,
0053       17.7, 18.1, 18.8, 19.3, 19.8, 20.6, 21.4, 22.4, 23.1, 23.6, 24.1, 24.2, 24.6, 24.8, 25.2,
0054       25.7, 26.5, 27.1, 28.2, 29.0, 29.9, 30.8, 31.1, 31.7, 31.8, 31.6, 31.5, 31.5, 31.3, 31.0,
0055       30.8, 30.8, 30.4, 30.2, 30.3, 30.2, 30.1, 30.1, 30.1, 29.8, 29.9, 29.8, 29.7, 29.7, 29.7,
0056       29.8, 29.8, 29.9, 29.9, 29.8, 29.9, 29.8, 29.9, 29.8, 29.7, 29.8, 29.7, 29.8, 29.6, 29.5,
0057       29.7, 29.7, 29.8, 30.1, 30.4, 31.0, 31.3, 31.5, 31.8, 31.8, 31.9, 32.0, 32.0, 32.0, 32.0,
0058       32.2, 32.2, 32.1, 31.8, 31.8, 31.8, 31.7, 31.6, 31.6, 31.7, 31.5, 31.5, 31.4, 31.3, 31.3,
0059       31.2, 30.8, 30.7, 30.5, 30.3, 29.9, 29.5, 29.3, 29.2, 28.6, 28.2, 27.9, 27.8, 27.3, 27.0,
0060       26.6, 26.1, 25.9, 25.5, 25.0, 24.6, 24.2, 23.8, 23.4, 23.0, 22.7, 22.4, 21.9, 21.4, 21.2,
0061       20.7, 20.3, 19.8, 19.6, 19.3, 18.9, 18.7, 18.3, 17.9, 17.8, 17.8, 16.7, 16.5, 16.4, 16.0,
0062       15.6, 15.6, 15.2, 14.9, 14.6, 14.4, 14.1, 13.8, 13.6, 13.3, 13.0, 12.8, 12.6, 12.3, 12.0,
0063       11.9, 11.7, 11.5, 11.2, 11.1, 10.9, 10.7, 10.4, 10.3, 9.9,  9.8,  9.6,  9.3,  9.1,  9.0,
0064       8.8,  8.5,  8.3,  8.3,  8.2,  7.9,  7.8,  7.7,  7.5,  7.3,  7.1,  6.9,  6.7,  6.6,  6.3,
0065       6.2,  6.0,  5.8,  5.7,  5.6,  5.4,  5.2,  5.1,  4.9,  4.8,  4.6,  4.5,  4.4,  4.2,  4.1,
0066       4.0,  3.8,  3.7,  3.5,  3.3,  3.2,  3.1,  3.0,  2.9,  2.5,  2.4,  2.4,  2.3,  2.3,  2.1,
0067       1.8,  1.6,  1.5,  1.5,  1.6,  1.8,  1.9,  1.4,  0.8,  0.9,  0.8,  0.7,  0.6,  0.3,  0.3,
0068       0.5,  0.3,  0.4,  0.3,  0.1,  0.2,  0.1,  0.2,  0.3,  0.0};
0069
0070     for(size_t i=0; i < sensor_qe.size(); i++)
0071       {
0072         double wavelength = 180 + i * 2; // wavelength units are [nm]
0073         digi_cfg.quantumEfficiency.push_back({wavelength, sensor_qe[i]*0.01});
0074       }
0075
0076
0077     // digitization
0078     app->Add(new JOmniFactoryGeneratorT<PhotoMultiplierHitDigi_factory>(
0079           "DIRCRawHits",
0080           {"DIRCBarHits"},
0081           {"DIRCRawHits", "DIRCRawHitsAssociations"},
0082           digi_cfg,
0083           app
0084           ));
0085
0086     int BarrelDIRC_ID = 0;
0087     try {
0088         auto detector = app->GetService<DD4hep_service>()->detector();
0089         BarrelDIRC_ID = detector->constant<int>("BarrelDIRC_ID");
0090     } catch(const std::runtime_error&) {
0091         // Nothing
0092     }
0093     PIDLookupConfig pid_cfg {
0094       .filename="calibrations/hpdirc.lut.gz",
0095       .system=BarrelDIRC_ID,
0096       .pdg_values={11, 211, 321, 2212},
0097       .charge_values={-1, 1},
0098       .momentum_edges={0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0, 6.2, 6.4, 6.6, 6.8, 7.0, 7.2, 7.4, 7.6, 7.8, 8.0, 8.2, 8.4, 8.6, 8.8, 9.0, 9.2, 9.4, 9.6, 9.8, 10.0, 10.2},
0099       .polar_edges={25.0, 26.0, 27.0, 28.0, 29.0, 30.0, 31.0, 32.0, 33.0, 34.0, 35.0, 36.0, 37.0, 38.0, 39.0, 40.0, 41.0, 42.0, 43.0, 44.0, 45.0, 46.0, 47.0, 48.0, 49.0, 50.0, 51.0, 52.0, 53.0, 54.0, 55.0, 56.0, 57.0, 58.0, 59.0, 60.0, 61.0, 62.0, 63.0, 64.0, 65.0, 66.0, 67.0, 68.0, 69.0, 70.0, 71.0, 72.0, 73.0, 74.0, 75.0, 76.0, 77.0, 78.0, 79.0, 80.0, 81.0, 82.0, 83.0, 84.0, 85.0, 86.0, 87.0, 88.0, 89.0, 90.0, 91.0, 92.0, 93.0, 94.0, 95.0, 96.0, 97.0, 98.0, 99.0, 100.0, 101.0, 102.0, 103.0, 104.0, 105.0, 106.0, 107.0, 108.0, 109.0, 110.0, 111.0, 112.0, 113.0, 114.0, 115.0, 116.0, 117.0, 118.0, 119.0, 120.0, 121.0, 122.0, 123.0, 124.0, 125.0, 126.0, 127.0, 128.0, 129.0, 130.0, 131.0, 132.0, 133.0, 134.0, 135.0, 136.0, 137.0, 138.0, 139.0, 140.0, 141.0, 142.0, 143.0, 144.0, 145.0, 146.0, 147.0, 148.0, 149.0, 150.0, 151.0, 152.0, 153.0, 154.0, 155.0, 156.0, 157.0, 158.0, 159.0, 160.0},
0100       .azimuthal_binning={0.0, 30.5, 0.5}, // lower, upper, step
0101     };
0102
0103     app->Add(new JOmniFactoryGeneratorT<PIDLookup_factory>(
0104           "DIRCTruthSeededLUTPID",
0105           {
0106           "ReconstructedTruthSeededChargedWithPFRICHTOFPIDParticles",
0107           "ReconstructedTruthSeededChargedWithPFRICHTOFPIDParticleAssociations",
0108           },
0109           {
0110           "ReconstructedTruthSeededChargedWithPFRICHTOFDIRCPIDParticles",
0111           "ReconstructedTruthSeededChargedWithPFRICHTOFDIRCPIDParticleAssociations",
0112           "DIRCTruthSeededParticleIDs",
0113           },
0114           pid_cfg,
0115           app
0116           ));
0117
0118     app->Add(new JOmniFactoryGeneratorT<PIDLookup_factory>(
0119           "DIRCLUTPID",
0120           {
0121           "ReconstructedChargedWithPFRICHTOFPIDParticles",
0122           "ReconstructedChargedWithPFRICHTOFPIDParticleAssociations",
0123           },
0124           {
0125           "ReconstructedChargedWithPFRICHTOFDIRCPIDParticles",
0126           "ReconstructedChargedWithPFRICHTOFDIRCPIDParticleAssociations",
0127           "DIRCParticleIDs",
0128           },
0129           pid_cfg,
0130           app
0131           ));
0132   }
0133 }