本文整理汇总了C++中Flavour::Charge方法的典型用法代码示例。如果您正苦于以下问题：C++ Flavour::Charge方法的具体用法？C++ Flavour::Charge怎么用？C++ Flavour::Charge使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Flavour的用法示例。
在下文中一共展示了Flavour::Charge方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: AnalyseEvent
void AnalyseEvent(Blob_List* blobs)
{
#ifdef USING__ROOT
// int outgoing = 1;
// int incoming = -1;
// Particle_List outparts = blobs->ExtractParticles(part_status::active, outgoing);
///////////////////////////////////////////////////////////////////////////////////
// analyse primary decay blob, ignore subsequent decays //
///////////////////////////////////////////////////////////////////////////////////
Blob * primarydecayblob = blobs->FindFirst(btp::Hadron_Decay);
// msg_Out()<<"primary decay blob:"<<endl<<*primarydecayblob<<endl;
// photon multiplicity and decay frame radiated energy (total)
unsigned int photmult = 0;
double photener = 0.;
for (int i=0; i<primarydecayblob->NOutP(); i++) {
if ((primarydecayblob->OutParticle(i)->Flav().IsPhoton() == true) &&
(primarydecayblob->OutParticle(i)->Info() == 'S')) {
photmult++;
photener = photener + primarydecayblob->OutParticle(i)->Momentum()[0];
}
}
photonmultiplicity->Fill(photmult);
if (photener != 0.) decayframeenergy->Fill(photener);
// multipole rest frame angles
Vec4D multipolesum = Vec4D(0.,0.,0.,0.);
Vec4D axis = Vec4D(0.,0.,0.,1.);
std::list<Vec4D> multipole;
std::list<Vec4D> newphot;
for (int i=0; i<primarydecayblob->NOutP(); i++) {
if (primarydecayblob->OutParticle(i)->Flav().Charge() != 0.) {
multipolesum = multipolesum + primarydecayblob->OutParticle(i)->Momentum();
multipole.push_back(primarydecayblob->OutParticle(i)->Momentum());
}
}
if (primarydecayblob->InParticle(0)->Flav().Charge() != 0) {
multipolesum = multipolesum + primarydecayblob->InParticle(0)->Momentum();
multipole.push_front(primarydecayblob->InParticle(0)->Momentum());
}
Poincare boost(multipolesum);
Poincare rotate;
// charged initial state: rotate such that initial state at theta = 0
if (mother_flav.Charge() != 0.) {
Vec4D inmom = *multipole.begin();
boost.Boost(inmom);
rotate = Poincare(inmom,axis);
}
// neutral initial state: rotate such that heaviest charged final state at theta = 0
else {
std::list<Vec4D>::iterator heaviest = multipole.begin();
for (std::list<Vec4D>::iterator iter=multipole.begin(); iter!=multipole.end(); iter++) {
if (abs((iter->Abs2() - heaviest->Abs2())/(iter->Abs2() + heaviest->Abs2())) > 1E-6) {
heaviest = iter;
}
}
boost.Boost(*heaviest);
rotate = Poincare(*heaviest,axis);
}
for (int i=0; i<primarydecayblob->NOutP(); i++) {
if (primarydecayblob->OutParticle(i)->Flav().IsPhoton() == true) {
Vec4D mom = primarydecayblob->OutParticle(i)->Momentum();
boost.Boost(mom);
rotate.Rotate(mom);
double theta = acos((Vec3D(mom)*Vec3D(axis))/(Vec3D(mom).Abs()*Vec3D(axis).Abs()));
multipoleframeangles->Fill(theta);
}
}
///////////////////////////////////////////////////////////////////////////////////
// inclusive analysis of whole decay chain //
///////////////////////////////////////////////////////////////////////////////////
// to be done ..
#endif
}