// Channel_Generator3V #include "PHASIC++/Channels/Single_Channel.H" #include "ATOOLS/Org/Run_Parameter.H" #include "ATOOLS/Org/MyStrStream.H" #include "ATOOLS/Org/Scoped_Settings.H" #include "PHASIC++/Channels/Channel_Elements.H" #include "PHASIC++/Channels/Vegas.H" using namespace PHASIC; using namespace ATOOLS; namespace PHASIC { class C2_1 : public Single_Channel { double m_thexp; double m_alpha,m_ctmax,m_ctmin; Info_Key m_kTC_0__1__2_3,m_kZS_0; Vegas* p_vegas; public: C2_1(int,int,Flavour*,Integration_Info * const); ~C2_1(); void GenerateWeight(Vec4D *,Cut_Data *); void GeneratePoint(Vec4D *,Cut_Data *,double *); void AddPoint(double); void MPISync() { p_vegas->MPISync(); } void Optimize() { p_vegas->Optimize(); } void EndOptimize() { p_vegas->EndOptimize(); } void WriteOut(std::string pId) { p_vegas->WriteOut(pId); } void ReadIn(std::string pId) { p_vegas->ReadIn(pId); } void ISRInfo(int &,double &,double &); std::string ChID(); }; } extern "C" Single_Channel * Getter_C2_1(int nin,int nout,Flavour* fl,Integration_Info * const info) { return new C2_1(nin,nout,fl,info); } void C2_1::GeneratePoint(Vec4D * p,Cut_Data * cuts,double * _ran) { double *ran = p_vegas->GeneratePoint(_ran); for(int i=0;iGenerateWeight(p_rans); if (wt!=0.) wt = vw/wt/pow(2.*M_PI,2*3.-4.); m_weight = wt; } C2_1::C2_1(int nin,int nout,Flavour* fl,Integration_Info * const info) : Single_Channel(nin,nout,fl) { Settings& s = Settings::GetMainSettings(); m_name = std::string("C2_1"); m_rannum = 2; p_rans = new double[m_rannum]; m_thexp = s["THRESHOLD_EXPONENT"].Get(); m_alpha = s["TCHANNEL_ALPHA"].Get(); m_ctmax = 1.; m_ctmin = -1.; m_kTC_0__1__2_3.Assign(std::string("TC_0__1__2_3"),2,0,info); m_kZS_0.Assign(std::string("ZS_0"),2,0,info); p_vegas = new Vegas(m_rannum,100,m_name); } C2_1::~C2_1() { delete p_vegas; } void C2_1::ISRInfo(int & type,double & mass,double & width) { type = 2; mass = 0; width = 0.; } void C2_1::AddPoint(double Value) { Single_Channel::AddPoint(Value); p_vegas->AddPoint(Value,p_rans); } std::string C2_1::ChID() { return std::string("CG2$TC_0__1__2_3$ZS_0$"); }