BaryonResonanceDecayer.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::BaryonResonanceDecayer
 
\brief    Baryon resonance decayer module.
 
          A simple resonance decay simulation built using resonance branching
          fraction data and an N-body phase space generator.
          Since the resonance can be produced off-the-mass-shell, decay
          channels with total-mass > W are suppressed.
 
          Is a concerete implementation of the EventRecordVisitorI interface.
 
\author   Costas Andreopoulos <c.andreopoulos \at cern.ch>
          University of Liverpool
 
\created  November 27, 2004
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
*/
//____________________________________________________________________________
 
#ifndef _BARYON_RESONANCE_DECAYER_H_
#define _BARYON_RESONANCE_DECAYER_H_
 
#include <TGenPhaseSpace.h>
#include <TLorentzVector.h>
 
#include "Physics/Decay/Decayer.h"
 
namespace genie {
 
class GHepParticle;
class BaryonResonanceDecayer : protected Decayer {
 
public:
  BaryonResonanceDecayer();
  BaryonResonanceDecayer(string config);
  virtual ~BaryonResonanceDecayer();
 
  // Implement the EventRecordVisitorI interface
  void ProcessEventRecord(GHepRecord * event) const;
 
  virtual void LoadConfig    (void);
 
private:
 
  void           Initialize        (void) const;
  bool           IsHandled         (int pdgc) const;
  void           InhibitDecay      (int pdgc, TDecayChannel * ch=0) const;
  void           UnInhibitDecay    (int pdgc, TDecayChannel * ch=0) const;
  double         Weight            (void) const;
  bool           Decay             (int dec_part_id, GHepRecord * event) const;
  TDecayChannel* SelectDecayChannel(int dec_part_id, GHepRecord * event, bool & to_be_deleted ) const;
  // the flag to_be_deleted is referred to the returned decay channel 
  bool           DecayExclusive    (int dec_part_id, GHepRecord * event, TDecayChannel * ch) const;
 
  // Methods specific for Delta decay
  TObjArray *    EvolveDeltaBR        (int dec_part_pdgc, TObjArray * decay_list, double W) const;
  double         EvolveDeltaDecayWidth(int dec_part_pdgc, TDecayChannel * ch, double W) const;
  bool           AcceptPionDecay( TLorentzVector lab_pion, int dec_part_id, const GHepRecord * event ) const ;
 
  double         FinalStateMass    ( TDecayChannel * ch ) const;
  bool           IsPiNDecayChannel ( TDecayChannel * ch ) const;
 
  static bool IsDelta( int dec_part_pdgc ) ; 
  static bool HasEvolvedBRs( int dec_part_pdgc ) ; 
 
  // utilities for pion angular distribution with phi dependency 
 
  static double PionAngularDist( const double * x, const double * par ) ;
  static double MinusPionAngularDist( const double * x, const double * par ) {  // this is used to find the maxima of the previous function
    return -1. * BaryonResonanceDecayer::PionAngularDist( x, par ) ; 
  }
 
  double FindDistributionExtrema( unsigned int i /*q2_bin_index*/, 
                                  bool find_maximum = false  ) const ;
 
 
  mutable TGenPhaseSpace fPhaseSpaceGenerator;
  mutable double         fWeight;
 
  bool   fDeltaThetaOnly ;
 
  double fMaxTolerance ;
 
  std::vector<double> fR33, fR31, fR3m1 ;
  std::vector<double*> fRParams ;  // this contains same parameter as fR33, fR31 and fR3m1 alingned by Q2 bin
 
  std::vector<double> fQ2Thresholds ;
 
  std::vector<double> fW_max ;
 
  double fFFScaling ;  // Scaling factor of the form factor of the Delta wrt to Q2
 
};
 
}         // genie namespace
#endif    // _BARYON_RESONANCE_DECAYER_H_