EvtOpenCharm.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See CVS repository
//      Copyright (A) 2011      Ping Rong-Gang
//
// Module: EvtOpenCharm.cc
//
// Description: Routine to decay charmonium-> DD, DDpi according the
// cross section measurement by CLEO  PRD 80, 072001.
//
// Modification history:
//
//    Ping R.-G.    December, 2011       Module created
//
//------------------------------------------------------------------------
//
 
 
#ifndef EVTOPENCHARM_HH
#define EVTOPENCHARM_HH
 
#include "EvtGenBase/EvtDecayIncoherent.hh"
 
 
class EvtParticle;
 
typedef EvtDecayBase* EvtDecayBasePtr;
 
#include <iosfwd>
 
class EvtOpenCharm:public  EvtDecayIncoherent  {
 
public:
 
  EvtOpenCharm();
  virtual ~EvtOpenCharm();
 
  void getName(std::string& name);
  EvtDecayBase* clone();
  void decay(EvtParticle *p); 
 
  std::string commandName();
  void command(std::string cmd);
 
  void init();
 
  void initProbMax();
  int getTotalEvt(){return nevt;}
  static void OpencrmInit(int f);
 
  static int myiter;
  
  bool isbelong(EvtId myid);
  int which_mode(EvtId myid);
 
private:
 
  void store(EvtDecayBase* jsdecay);
  void fixPolarizations(EvtParticle* p);
 
  static int nOpencharmdecays;
  static EvtDecayBasePtr* Opencharmdecays;
  static int ntable;
 
  static int ncommand;
  static int lcommand;
  static std::string* commands;
  static int nevt;
 
  int theMode;
  std::vector<EvtId> theId;
  int ndaugjs;
  static std::vector<EvtId> mypar;
  static std::vector<int> vmode;
  static std::vector<double> Vcms;
  int _index;
};
 
#endif