packages/CaloReco/BEmcRecCEMC.h

0001 #ifndef CALORECO_BEMCRECCEMC_H
0002 #define CALORECO_BEMCRECCEMC_H
0003
0004 #include "BEmcRec.h"
0005
0006 #include <string>
0007 #include <vector>
0008 #include <limits>  // for std::numeric_limits<float>::quiet_NaN()
0009
0010 // class BEmcProfile;
0011 class EmcModule;
0012
0013 class BEmcRecCEMC : public BEmcRec
0014 {
0015  public:
0016   BEmcRecCEMC();
0017   ~BEmcRecCEMC() override = default;
0018   void CorrectEnergy(float energy, float x, float y, float &ecorr) override;
0019   void CorrectECore(float ecore, float x, float y, float &ecorecorr) override;
0020   void CorrectPosition(float energy, float x, float y, float &xcorr, float &ycorr) override;
0021   void CorrectShowerDepth(int ix, int iy, float energy, float x, float y, float z, float &xc, float &yc, float &zc) override;
0022
0023   void LoadProfile(const std::string &fname) override;
0024   //  float GetProb(std::vector<EmcModule> HitList, float e, float xg, float yg, float zg, float &chi2, int &ndf) override;
0025   void GetImpactThetaPhi(float xg, float yg, float zg, float &theta, float &phi) override;
0026
0027   bool calculateIncidence(float x, float y,
0028                           float& a_phi_sgn, float& a_eta_sgn);
0029   // optional accessors for the cached signed angles
0030   float lastSignedAlphaPhi() const { return m_lastAlphaPhiSigned; }
0031   float lastSignedAlphaEta()  const { return m_lastAlphaEtaSigned; }
0032
0033  private:
0034   //  BEmcProfile *_emcprof;
0035   // cache of the last signed mechanical incidence angles (in radians)
0036   float m_lastAlphaPhiSigned { std::numeric_limits<float>::quiet_NaN() };
0037   float m_lastAlphaEtaSigned { std::numeric_limits<float>::quiet_NaN() };
0038   // Average tower angle with respect to the transverse plane for each bin in eta
0039   // Only used when the detailed RawTowerGeom objects are not available
0040   double angles[96] = {2.382132067280038,2.3708136440088796,2.355383841518272,2.3434466029558525,2.32455480538339,2.3135158793606005,
0041                        2.2999027706894153,2.289397309387001,2.2727994372634956,2.259380682790801,2.2392521417328277,2.226335700118477,
0042                        2.2095791811782735,2.1935361258633024,2.1754991989110875,2.160931303742618,2.142623616095435,2.1259546238847804,
0043                        2.1054808114318946,2.0900705081463244,2.0670824189635053,2.050851798834259,2.026500612990249,2.005259135338716,
0044                        1.9868346921870874,1.966226760986188,1.9455349050274326,1.924706496569077,1.9017832210973762,1.8805537731028699,
0045                        1.8568682630187272,1.8337068107715466,1.8112862676694754,1.7884337950298892,1.76612732335177,1.7439109866141163,
0046                        1.7137229887165837,1.6887545076082588,1.6601011126718555,1.633546284539203,1.6126575096324751,1.594318479859694,
0047                        1.5699679790263024,1.5683722241543814,1.570109507416994,1.5686069317665148,1.5706670797900908,1.5655401209498148,
0048                        1.5760525326399784,1.5709255737997023,1.5729857218232786,1.5714831461727994,1.573220429435412,1.5716246745634908,
0049                        1.5472741737300992,1.528935143957318,1.50804636905059,1.4814915409179379,1.4528381459815345,1.4278696648732097,
0050                        1.3976816669756769,1.375465330238023,1.353158858559904,1.3303063859203177,1.3078858428182465,1.284724390571066,
0051                        1.2610388804869233,1.2398094324924172,1.2168861570207163,1.1960577485623607,1.1753658926036052,1.154757961402706,
0052                        1.1363335182510774,1.1150920405995441,1.0907408547555344,1.074510234626288,1.051522145443469,1.0361118421578988,
0053                        1.015638029705013,0.9989690374943581,0.9806613498471752,0.9660934546787058,0.948056527726491,0.9320134724115198,
0054                        0.9152569534713161,0.9023405118569655,0.8822119707989924,0.8687932163262975,0.852195344202792,0.8416898829003777,
0055                        0.8280767742291925,0.817037848206403,0.798146050633941,0.7862088120715213,0.7707790095809137,0.7594605863097555};
0056
0057   double slopes_[96] = {0, 0.00023905, 4.90448e-05, 0.00027207, 4.62457e-05, 0.000252483, 4.07806e-05, 0.000248309, 3.76606e-05, 0.000257749, 3.96662e-05, 0.000268067, 3.79489e-05, 0.000258554, 4.64794e-05, 0.000252937, 4.00813e-05, 0.000256408, 5.76557e-05, 0.000254708, 6.78359e-05, 0.00025745, 9.10421e-05, 0.000249737, 9.98244e-05, 0.000248255, 0.000129006, 0.000248438, 0.000137128, 0.000252377, 0.000162495, 0.000255381, 0.000178851, 0.000255713, 0.000196615, 0.000257513, 0.000218168, 0.000260783, 0.00023512, 0.000266362, 0.000253771, 0.000258489, 0.000261129, 0.000237658, 0.000229224, 0.000219019, 0.0002103555, 0.000194438, 0.000191487, 0.000169857, 0.000163719, 0.000145249, 0.000140684, 0.00011482, 0.000114154, 0.000100124, 0.00011904, 7.79723e-05, 0.000107485, 5.65278e-05, 9.40123e-05, 3.13618e-05, 8.63381e-05, 6.76016e-06, 7.83596e-05, -1.41566e-05, 7.57772e-05, -3.53334e-05, 6.71255e-05, -5.46557e-05, 6.50159e-05, -7.18416e-05, 6.88217e-05, -8.98952e-05, 6.77945e-05, -0.000104923, 6.87795e-05, -0.000117168, 7.67749e-05, -0.000130278, 8.43921e-05, -0.000137269, 8.63373e-05, -0.000142283, 8.78225e-05, -0.000145142, 9.22318e-05, -0.000121506, 0.000104699, -0.000114551, 0.00010633, -0.0001077, 9.84409e-05, -0.000109789, 8.14774e-05, -8.70686e-05};
0058
0059   double factor_[8] = {-0.07, 0.03, -0.07, 0.20, 0.15, -0.05, 0.00, -0.07}; // Correction factors for the phi bias within a sector of 8 towers (phi direction)
0060
0061 };
0062
0063
0064
0065 #endif