waveSetGenerator.h

Go to the documentation of this file.

//
//    Copyright 2010
//
//    This file is part of rootpwa
//
//    rootpwa is free software: you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation, either version 3 of the License, or
//    (at your option) any later version.
//
//    rootpwa is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//
//    You should have received a copy of the GNU General Public License
//    along with rootpwa. If not, see <http://www.gnu.org/licenses/>.
//
//-------------------------------------------------------------------------
// File and Version Information:
// $Rev::                             $: revision of last commit
// $Author::                          $: author of last commit
// $Date::                            $: date of last commit
//
// Description:
//      class that generates set of isobar decay topologies from a
//      template and according to user defined criteria
//
//
// Author List:
//      Boris Grube          TUM            (original author)
//
//
//-------------------------------------------------------------------------


#ifndef WAVESETGENERATOR_H
#define WAVESETGENERATOR_H


#include <boost/tuple/tuple.hpp>

#include "isobarDecayTopology.h"


namespace rpwa {


        class waveSetGenerator {

        public:

                // some typedefs for convenience
                typedef isobarDecayTopology::nodeDesc    nodeDesc;
                typedef isobarDecayTopology::adjIterator adjIterator;


                waveSetGenerator();
                virtual ~waveSetGenerator();

                bool setWaveSetParameters(const std::string& templateKeyFileName);  

                // wave set parameter accessors
                // !note! isospin and angular momentum quantum numbers are in units of hbar / 2
                void setIsospinRange         (const int maxI = 2,
                                              const int minI = 0) { _isospinRange = boost::tuples::make_tuple(minI, maxI); }
                void setJRange               (const int maxJ = 6,
                                              const int minJ = 0) { _JRange       = boost::tuples::make_tuple(minJ, maxJ); }
                void setLRange               (const int maxL = 6,
                                              const int minL = 0) { _LRange       = boost::tuples::make_tuple(minL, maxL); }
                void setSRange               (const int maxS = 6,
                                              const int minS = 0) { _SRange       = boost::tuples::make_tuple(minS, maxS); }
                void setIsobarBlackList      (const std::vector<std::string>& isobarList)
                { _isobarBlackList = isobarList; }
                void setIsobarWhiteList      (const std::vector<std::string>& isobarList)
                { _isobarWhiteList = isobarList; }
                void setAllowSpinExotics     (const bool   flag     ) { _allowSpinExotics      = flag;  }
                void setRequireMinIsobarMass (const bool   flag     ) { _requireMinIsobarMass  = flag;  }
                void setForceDecayCheck      (const bool   flag     ) { _forceDecayCheck       = flag;  }
                void setIsobarMassWindowSigma(const double sigma = 1) { _isobarMassWindowSigma = sigma; }

                std::size_t generateWaveSet();  

                std::vector<isobarDecayTopology>&       waveSet()       { return _waveSet; }  
                const std::vector<isobarDecayTopology>& waveSet() const { return _waveSet; }  

                bool writeKeyFiles(const std::string& dirName                  = "",
                                   const bool         newKeyFileNameConvention = false);  

                virtual void reset();  

                virtual std::ostream& print(std::ostream& out) const;  

                static bool debug() { return _debug; }                             
                static void setDebug(const bool debug = true) { _debug = debug; }  


        private:

                const isobarDecayTopologyPtr createNewDecayTopology(const isobarDecayTopology&  parentDecay,
                                                                    const isobarDecayVertexPtr& parentVertex,
                                                                    const int                   L,
                                                                    const int                   S,
                                                                    const particleProperties&   isobar,
                                                                    const int                   parentCharge);

                std::set<std::size_t> findBoseSymDecays() const;  

                boost::tuples::tuple<int, int> _isospinRange;           
                boost::tuples::tuple<int, int> _JRange;                 
                boost::tuples::tuple<int, int> _spinProjRange;          
                int                            _reflectivity;           
                bool                           _useReflectivity;        
                bool                           _allowSpinExotics;       
                boost::tuples::tuple<int, int> _LRange;                 
                boost::tuples::tuple<int, int> _SRange;                 
                std::vector<std::string>       _isobarBlackList;        
                std::vector<std::string>       _isobarWhiteList;        
                bool                           _requireMinIsobarMass;   
                double                         _isobarMassWindowSigma;  
                bool                           _forceDecayCheck;        

                isobarDecayTopologyPtr _templateTopo;  

                std::vector<isobarDecayTopology> _waveSet;  

                static bool _debug;  

        };


        inline
        std::ostream&
        operator <<(std::ostream&           out,
                    const waveSetGenerator& gen)
        {
                return gen.print(out);
        }


} // namespace rpwa


#endif  // WAVESETGENERATOR_H