@@ -48,7 +48,7 @@ class GPUTRDRecoParam
4848
4949 GPUd () float getRPhiRes (float snp, float pull = 0 .f, int occupancy = 0 ) const ;
5050 GPUd () float getDyRes (float snp, int occupancy = 0 ) const { return mDyA2 + mDyC2 * (snp - mLorentzAngle ) * (snp - mLorentzAngle ) + mOccDyA * occupancy; } // a^2 + c^2 * (snp - b)^2
51- GPUd () float convertAngleToDy (float snp) const { return 3 .f * snp / CAMath::Sqrt (1 - snp * snp); } // when calibrated, sin(phi) = (dy / xDrift) / sqrt(1+(dy/xDrift)^2) works well
51+ GPUd () float convertAngleToDy (float snp) const { return 3 .f * snp / CAMath::Sqrt (1 - snp * snp); } // when calibrated, sin(phi) = (dy / xDrift) / sqrt(1+(dy/xDrift)^2) works well
5252 GPUd () float getCorrYDy () const { return mCorrYDy ; }
5353 GPUd () float getPileUpProbTracklet (int nBC, bool withChargeInfo, bool Q0 = true , bool Q1 = true ) const ;
5454 GPUd () float getPileUpProbTrack (int nBC, std::array<int , 6 > Q0 , std::array<int , 6 > Q1 ) const ;
@@ -64,9 +64,9 @@ class GPUTRDRecoParam
6464 // tracklet error parameterization depends on the magnetic field
6565 float mLorentzAngle {0 .f };
6666 // rphi
67- float mRPhiA {1 .f }; // /< parameterization for tracklet position resolution
68- float mRPhiATgp {1 .f }; // /< parameterization for tracklet position resolution
69- float mRPhiC2 {0 .f }; // /< parameterization for tracklet position resolution
67+ float mRPhiA {1 .f }; // /< parameterization for tracklet position resolution
68+ float mRPhiATgp {1 .f }; // /< parameterization for tracklet position resolution
69+ float mRPhiC2 {0 .f }; // /< parameterization for tracklet position resolution
7070 // angle
7171 float mDyA2 {1 .225e-3f }; // /< parameterization for tracklet angular resolution
7272 float mDyC2 {0 .f }; // /< parameterization for tracklet angular resolution
@@ -114,7 +114,8 @@ class GPUTRDRecoParam
114114// / \param phi angle of related track
115115// / \return sigma_y^2 of tracklet
116116// / also depend on absolute pull and on chamber occupancy
117- GPUdi () float GPUTRDRecoParam::getRPhiRes (float snp, float pull, int occupancy) const {
117+ GPUdi () float GPUTRDRecoParam::getRPhiRes (float snp, float pull, int occupancy) const
118+ {
118119 // flat uncertainty + radial-alignment uncertainty depending on tan(phi)
119120 float tgp = (CAMath::Abs (snp) < 0 .99999f ) ? CAMath::Abs (snp) / CAMath::Sqrt (1 - snp * snp) : 1e6 ;
120121 float resIdeal = mRPhiA + mRPhiATgp * tgp;
@@ -123,8 +124,8 @@ GPUdi() float GPUTRDRecoParam::getRPhiRes(float snp, float pull, int occupancy)
123124 pull = 10 .f ;
124125 }
125126 float resPull = mPullA * pull * pull + mPullB * pull; // parametrization as pol2 summed in quadrature
126- float resOccupancy = mOccA * occupancy; // parametrization as sqrt() summed in quadrature
127- return (resIdeal * resIdeal + mRPhiC2 * (snp - mLorentzAngle ) * (snp - mLorentzAngle ) + resPull * resPull + resOccupancy);
127+ float resOccupancy = mOccA * occupancy; // parametrization as sqrt() summed in quadrature
128+ return (resIdeal * resIdeal + mRPhiC2 * (snp - mLorentzAngle ) * (snp - mLorentzAngle ) + resPull * resPull + resOccupancy);
128129}
129130
130131GPUdi () float GPUTRDRecoParam::getPileUpProbTracklet (int nBC, bool withChargeInfo, bool Q0 , bool Q1 ) const
0 commit comments