C++ Flavour::Strong方法代码示例

void SM_U1_B::FillSpectrum(const PDF::ISR_Handler_Map& isr) {
  p_dataread->RereadInFile();
  p_constants->insert(make_pair(string("M_Z'"),    
				p_dataread->GetValue<double>("M_Z'",-1.)));
  p_constants->insert(make_pair(string("g'_1"),    
				p_dataread->GetValue<double>("g'_1",0.)));

  double MZprime(ScalarConstant(string("M_Z'"))),GZprime(0.),inc;
  double g1prime(sqr(ScalarConstant(std::string("g'_1"))));
  FixMix();
  msg_Out()<<"Calculate width for Z' with mass = "<<MZprime<<" GeV.\n";

  for (short int i=1;i<=6;i++) {
    Flavour quark = Flavour((kf_code)(i));
    double  massq = quark.HadMass(); 
    if (!quark.IsOn() || !quark.Strong()) continue;
    if (i%2) {
      if (2.*massq>MZprime) continue;
      GZprime += inc = 
	3.*MZprime*g1prime/(24.*M_PI)*pow(1.-sqr(2.*massq/MZprime),3./2.);
      msg_Out()<<"   Add Z' --> "<<quark<<"+"<<quark.Bar()<<" ("<<i<<"), "
	       <<" add "<<inc<<" GeV.\n";
    }
    else {
      double M2 = MZprime*MZprime, mq2 = massq*massq;
      for (short int j=2;j<=7;j+=2) {
	Flavour anti = Flavour((kf_code)(j)).Bar();
	if (!anti.IsOn() || !anti.Strong()) continue;
	double massa  = anti.HadMass(), ma2 = massa*massa; 
	if (massq+massa>MZprime) continue;
	Complex mix   = ComplexMatrixElement(std::string("UpMix"),i/2-1,j/2-1);
	double absmix = ATOOLS::sqr(mix.real())+ATOOLS::sqr(mix.imag());
	if (ATOOLS::IsZero(absmix)) continue;
	GZprime += inc = 
	  3.*g1prime*absmix/(24.*M_PI*M2)*
	  ((2.*M2-mq2-ma2)/2.-3.*massq*massa-ATOOLS::sqr(ma2-mq2)/(2.*M2)) *
	  sqrt(ATOOLS::sqr(M2-mq2-ma2)-4.*ma2*mq2)/(2.*MZprime);
	msg_Out()<<"   Add Z' --> "<<quark<<"+"<<anti<<" ("<<i<<"), "
		 <<" add "<<inc<<" GeV.\n";
      }
    }
  }

  Flavour flav;
  flav = Flavour(kf_Z0_2);
  flav.SetMass(ScalarConstant(string("M_Z'")));
  flav.SetHadMass(MZprime);
  flav.SetMassOn(true);
  flav.SetWidth(GZprime);

  msg_Out()<<METHOD<<" initializes Z' boson with \n"
	   <<"    mass = "<<MZprime<<" GeV and width = "<<GZprime<<" GeV\n"
	   <<"    for g'_1 = "<<ScalarConstant(std::string("g'_1"))<<".\n";
}

void Interaction_Model_QCD::c_FFV(std::vector<Single_Vertex>& vertex,int & vanz)
{
  Kabbala kcpl0 = -g3*M_I;
  Kabbala kcpl1 = kcpl0;

  for (short int i=1;i<=(m_extension==2?8:6);i++) {
    Flavour flav = Flavour((kf_code)(i));
    if (flav.Strong() && flav.IsOn() && Flavour(kf_gluon).IsOn()) { 
      vertex[vanz].in[0]         = flav;
      vertex[vanz].in[1]         = Flavour(kf_gluon);
      vertex[vanz].in[2]         = flav;

      vertex[vanz].cpl[0]        = kcpl0;
      vertex[vanz].cpl[1]        = kcpl1;
      vertex[vanz].Str           = (kcpl0*PR+kcpl1*PL).String();      

      
      vertex[vanz].Color.push_back(Color_Function(cf::T,1,2,0,'1','2','0'));     

      vertex[vanz].Lorentz.push_back(LF_Getter::GetObject("Gamma",LF_Key()));     
      vertex[vanz].Lorentz.back()->SetParticleArg(1);     
                  
      vertex[vanz].on            = 1;
      vertex[vanz].oqcd          = 1;
      vertex[vanz].oew           = 0;
      vertex.push_back(Single_Vertex());vanz++;
    } 
  }
}

 void Interaction_Model_SM_Zprime::c_FFV(vector<Single_Vertex>& vertex,int& vanz)
 {
   // create the vertices for the standard model
   Interaction_Model_SM::c_FFV(vertex,vanz);
   
   
   // create FFV vertices with Z' if it's on
   kf_code kfZp=32;
   Flavour flZprime(kfZp);
   if (flZprime.IsOn()) {
     
     // parse through all fermions than couple to Z' and create vertices
     int PossibleFermions[12] = {1,2,3,4,5,6,11,12,13,14,15,16};
     for (int i=0; i<12; i++) {
       
       // initialize the currently parsed fermion
       int FermionNumber = PossibleFermions[i];
       Flavour flFermion = Flavour((kf_code)(FermionNumber));
       Kabbala B = Kabbala(string("B_{")+flFermion.TexName()+string("}"),
                           flFermion.BaryonNumber());
       Kabbala L = Kabbala(string("L_{")+ flFermion.TexName()+string("}"),
                           flFermion.LeptonNumber());
       Kabbala Y3R = Kabbala(string("YR_{")+flFermion.TexName()+string("}"),
                             flFermion.IsoWeak());
       
       if (flFermion.IsOn()) {
         // create the vertex for that particular fermion and a Z'.
         // Right-handed neutrinos will not take part in any interaction.
         Kabbala kcpl0;
         if ((FermionNumber==12)||(FermionNumber==14)||(FermionNumber==16))
         {kcpl0 = Kabbala("0.0", 0.);}
         else {kcpl0 = -M_I * g2 * (Y3R * alphaLR + (L-B)/(alphaLR*2));};
         Kabbala kcpl1 = -M_I * g2 * (L-B) / (alphaLR*2);
         
         // set couplings and particle info for current vertex
         vertex[vanz].in[0] = flFermion;
         vertex[vanz].in[1] = flZprime;
         vertex[vanz].in[2] = Flavour((kf_code)(FermionNumber));
         vertex[vanz].cpl[0] = kcpl0;
         vertex[vanz].cpl[1] = kcpl1;
         vertex[vanz].Str = (kcpl0*PR+kcpl1*PL).String(); 
         
         // Color Function for vertex
         
         if (flFermion.Strong()) {
           vertex[vanz].Color.push_back(Color_Function(cf::D));;
           vertex[vanz].Color.back().SetParticleArg(0,2);
           vertex[vanz].Color.back().SetStringArg('0','2');
         } 
         else 
           vertex[vanz].Color.push_back(Color_Function(cf::None));;
         
         // Lorenz function for vertex
         
         vertex[vanz].Lorentz.push_back(LF_Getter::GetObject("Gamma",LF_Key()));
         vertex[vanz].Lorentz.back()->SetParticleArg(1);
         
         vertex[vanz].on     = 1;
         vertex.push_back(Single_Vertex());vanz++; 
       }; 
     };
   };
 }