include/EvtGenModels/EvtSVVHelCPMix.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 EVTSVVHELCPMIX_HH
0022 #define EVTSVVHELCPMIX_HH
0023
0024 #include "EvtGenBase/EvtDecayAmp.hh"
0025
0026 //Class to handle decays of the form SCALAR -> VECTOR VECTOR
0027 //according the the helicity amplitudes specified by the
0028 //user.  There are 6 arguements, orders as amplitude then
0029 //phase for H+, H0, and H-, in that order.
0030 //
0031 // Description: Routine to decay scalar -> 2 vectors
0032 //              by specifying the helicity amplitudes, taking appropriate
0033 //      weak phases into account to get mixing and CP violation through
0034 //      interference. Based on EvtSVVHelAmp. Particularly appropriate for
0035 //      Bs->J/Psi+Phi
0036 //
0037 // Model takes the following as user-specified arguments:
0038 //  deltaM, averageM - mass diference and average of light and heavy mass eigenstates (real scalars)
0039 //  gamma, deltagamma - average width and width difference of the l and h eigenstates (real scalars)
0040 //  delta1, delta2 - strong phases (real scalars)
0041 //  direct weak phase (real scalar) (for Bs->JPsiPhi this will be zero)
0042 //  weak mixing phase (real scalar) (this is equal to 2*arg(Vts Vtb) for Bs->JPsiPhi)
0043 //  Magnitudes of helicity amplitudes as in SVV_HELAMP
0044 // See Phys Rev D 34 p1404 - p1417 and chapters 5 and 7 of Physics Reports 370 p537-680 for more details
0045
0046 class EvtAmp;
0047 class EvtParticle;
0048 class EvtId;
0049
0050 class EvtSVVHelCPMix : public EvtDecayAmp {
0051   public:
0052     std::string getName() override;
0053     EvtDecayBase* clone() override;
0054
0055     void init() override;
0056
0057     EvtComplex hp;
0058     EvtComplex h0;
0059     EvtComplex hm;
0060     double averageM;
0061     double deltaM;
0062     double gamma;
0063     double deltagamma;
0064     EvtComplex strongphase1;
0065     EvtComplex strongphase2;
0066     EvtComplex weakmixingphase;
0067     EvtComplex weakdirectphase;
0068
0069     void initProbMax() override;
0070
0071     void decay( EvtParticle* p ) override;
0072
0073     std::string getParamName( int i ) override;
0074     std::string getParamDefault( int i ) override;
0075 };
0076
0077 #endif