HydjetHadronizer.h

Go to the documentation of this file.

#ifndef HydjetHadronizer_h
#define HydjetHadronizer_h


#define PYCOMP pycomp_

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "GeneratorInterface/Core/interface/BaseHadronizer.h"

#include <map>
#include <string>
#include <vector>
#include <math.h>

namespace CLHEP {
  class HepRandomEngine;
}

namespace HepMC {
  class GenEvent;
  class GenParticle;
  class GenVertex;
}

namespace gen
{
   class Pythia6Service;

  class HydjetHadronizer : public BaseHadronizer {
  public:
    HydjetHadronizer(const edm::ParameterSet &);
    virtual ~HydjetHadronizer();

    bool generatePartonsAndHadronize();
    bool hadronize();
    bool decay();
    bool residualDecay();
    bool readSettings( int );
    bool initializeForExternalPartons();
    bool initializeForInternalPartons();
    bool declareStableParticles( const std::vector<int>& );
    bool declareSpecialSettings( const std::vector<std::string>& ) { return true; }

    void finalizeEvent();
    void statistics();
    const char* classname() const;

  private:

    virtual void doSetRandomEngine(CLHEP::HepRandomEngine* v) override;
    virtual std::vector<std::string> const& doSharedResources() const override { return theSharedResources; }

    static const std::vector<std::string> theSharedResources;

    void                    add_heavy_ion_rec(HepMC::GenEvent *evt);
    HepMC::GenParticle*                         build_hyjet( int index, int barcode );  
    HepMC::GenVertex*                           build_hyjet_vertex(int i, int id);
    bool                    get_particles(HepMC::GenEvent* evt);
    bool                                        call_hyinit(double energy, double a, int ifb, double bmin,
                                double bmax,double bfix,int nh);
    bool                    hydjet_init(const edm::ParameterSet &pset);
    inline double                   nuclear_radius() const;
    void                                        rotateEvtPlane();

    HepMC::GenEvent   *evt;
    edm::ParameterSet  pset_;
    double             abeamtarget_;           // beam/target atomic mass number 
    int            angularspecselector_;   // angular emitted gluon spectrum selection
    double             bfixed_;                // fixed impact param (fm); valid only if cflag_=0
    double             bmax_;                  // max impact param; 
                                               // units of nucl radius
    double             bmin_;                  // min impact param; 
                                               // units of nucl radius
    int                cflag_;                 // centrality flag 
                                               // =  0 fixed impact param, 
                                               // <> 0 between bmin and bmax
    bool             embedding_;               // Switch for embedding mode
    double             comenergy;              // collision energy   
    bool               doradiativeenloss_;     
    bool               docollisionalenloss_;   
    double             fracsoftmult_;          // fraction of soft hydro induced hadronic multiplicity
                                               // proportional to no of nucleon participants
                                               // (1-fracsoftmult_)--- fraction of soft 
                                               // multiplicity proportional to the numebr 
                                               // of nucleon-nucleon binary collisions
                                               // DEFAULT=1., allowed range [0.01,1]
    double             hadfreeztemp_;          // hadron freez-out temperature
                                               // DEFAULT=0.14MeV, allowed ranges [0.08,0.2]MeV
    std::string        hymode_;                // Hydjet running mode
    unsigned int       maxEventsToPrint_;      // Events to print if verbosity  
    double             maxlongy_;              // max longitudinal collective rapidity: 
                                               // controls width of eta-spectra
                                               // DEFAULT=4, allowed range [0.01,7.0]
    double             maxtrany_;              // max transverse collective rapidity: 
                                               // controls slope of low-pt spectra
                                               // DEFAULT=1.5, allowed range [0.01,3.0]
    int                nsub_;                  // number of sub-events
    int                nhard_;                 // multiplicity of PYTHIA(+PYQUEN)-induced particles in event
    int                nmultiplicity_;         // mean soft multiplicity in central PbPb
                                               // automatically calculated for other centralitie and beams         
    int                nsoft_;                 // multiplicity of HYDRO-induced particles in event 
    unsigned int       nquarkflavor_;          
    unsigned int       pythiaPylistVerbosity_; // pythia verbosity; def=1 
    double             qgpt0_;                 // initial temperature of QGP
                                               // DEFAULT = 1GeV; allowed range [0.2,2.0]GeV; 
    double             qgptau0_;               // proper time of QGP formation
                                               // DEFAULT = 0.1 fm/c; allowed range [0.01,10.0]fm/

    double             phi0_;                  // Event plane angle
    double             sinphi0_;
    double             cosphi0_;
    bool               rotate_;                // Switch to rotate event plane

    unsigned int       shadowingswitch_;       // shadowing switcher
                                               // 1-ON, 0-OFF
    double             signn_;                 // inelastic nucleon nucleon cross section [mb]
                                               // DEFAULT= 58 mb

    Pythia6Service* pythia6Service_;
    edm::InputTag   src_;
  };

  double HydjetHadronizer::nuclear_radius() const
  {
    // Return the nuclear radius derived from the 
    // beam/target atomic mass number.

    return 1.15 * pow((double)abeamtarget_, 1./3.);
  }
} /*end namespace*/
#endif