BDSBOptrMultiParticleChangeCrossSection.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2022.
 
This file is part of BDSIM.
 
BDSIM 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 version 3 of the License.
 
BDSIM 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 BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSBOPTRMULTIPARTICLECHANGECROSSSECTION_H
#define BDSBOPTRMULTIPARTICLECHANGECROSSSECTION_H
 
#include "G4Version.hh"
#if G4VERSION_NUMBER > 1009
 
#include "G4VBiasingOperator.hh"
class BDSBOptrChangeCrossSection;
class G4ParticleDefinition; 
 
#include <map>
 
class BDSBOptrMultiParticleChangeCrossSection: public G4VBiasingOperator
{
public:
  BDSBOptrMultiParticleChangeCrossSection();
  virtual ~BDSBOptrMultiParticleChangeCrossSection();
  
  void AddParticle(const G4String& particleName);
  void SetBias(const G4String& biasObjectName,
               const G4String& particleName,
               const G4String& process,
               G4double        dBias,
               G4int           iPrimary);
  void StartTracking(const G4Track* track) override;
 
private: 
  virtual G4VBiasingOperation* ProposeOccurenceBiasingOperation(const G4Track*                   track,
                                                                const G4BiasingProcessInterface* callingProcess) override;
  // -- Methods not used:
  virtual G4VBiasingOperation* ProposeFinalStateBiasingOperation(const G4Track*, const G4BiasingProcessInterface*) override
  {return 0;}
  virtual G4VBiasingOperation* ProposeNonPhysicsBiasingOperation(const G4Track*, const G4BiasingProcessInterface*) override
  {return 0;}
  virtual void OperationApplied(const G4BiasingProcessInterface* callingProcess,
                                G4BiasingAppliedCase             biasingCase,
                                G4VBiasingOperation*             occurenceOperationApplied,
                                G4double                         weightForOccurenceInteraction,
                                G4VBiasingOperation*             finalStateOperationApplied, 
                                const G4VParticleChange*         particleChangeProduced) override;
  // prevent compiler warning (since second G4VBiasingOperator::OperationApplied is hidden)
  using G4VBiasingOperator::OperationApplied;
  
  std::map<const G4ParticleDefinition*, BDSBOptrChangeCrossSection*> fBOptrForParticle;
  std::vector<const G4ParticleDefinition*>                           fParticlesToBias;
  BDSBOptrChangeCrossSection*                                        fCurrentOperator;
 
  G4int  fnInteractions; 
  G4bool debug;
};
 
#endif
 
#endif