include/EvtGenBase/EvtRelBreitWignerBarrierFact.hh

0001
0002 /***********************************************************************
0003 * Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
0004 *                                                                      *
0005 * This file is part of EvtGen.                                         *
0006 *                                                                      *
0007 * EvtGen is free software: you can redistribute it and/or modify       *
0008 * it under the terms of the GNU General Public License as published by *
0009 * the Free Software Foundation, either version 3 of the License, or    *
0010 * (at your option) any later version.                                  *
0011 *                                                                      *
0012 * EvtGen is distributed in the hope that it will be useful,            *
0013 * but WITHOUT ANY WARRANTY; without even the implied warranty of       *
0014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
0015 * GNU General Public License for more details.                         *
0016 *                                                                      *
0017 * You should have received a copy of the GNU General Public License    *
0018 * along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
0019 ***********************************************************************/
0020
0021 #ifndef EVTRELBREITWIGNERBARRIERFACT_HH
0022 #define EVTRELBREITWIGNERBARRIERFACT_HH
0023
0024 #include "EvtGenBase/EvtAbsLineShape.hh"
0025
0026 class EvtRelBreitWignerBarrierFact : public EvtAbsLineShape {
0027   public:
0028     EvtRelBreitWignerBarrierFact() = default;
0029     EvtRelBreitWignerBarrierFact( double mass, double width, double maxRange,
0030                                   EvtSpinType::spintype sp );
0031     //figure the m1 and l on the fly
0032     //                 double mDaug1, double mDaug2, int l);
0033     EvtRelBreitWignerBarrierFact& operator=( const EvtRelBreitWignerBarrierFact& x );
0034     EvtRelBreitWignerBarrierFact( const EvtRelBreitWignerBarrierFact& x );
0035
0036     EvtAbsLineShape* clone() override;
0037
0038     double getMassProb( double mass, double massPar, int nDaug,
0039                         double* massDau ) override;
0040     // othDaugId is the other daughter of the parent in the case of a two body decay (only!)
0041     // ie B->rho K rho->pipi, othDaugId = K
0042     double getRandMass( EvtId* parId, int nDaug, EvtId* dauId, EvtId* othDaugId,
0043                         double maxMass, double* dauMasses ) override;
0044
0045     void reSetBlatt( double blatt ) override { _blattDecay = blatt; }
0046     void reSetBlattBirth( double blatt ) override { _blattBirth = blatt; }
0047
0048   protected:
0049     double _blattDecay;
0050     double _blattBirth;
0051     bool _errorCond;
0052 };
0053
0054 #endif