1. 1
2. 2
3. 3
4. 4
5. 5
6. 6
7. 7
8. 8
9. 9
10. 10
11. 11
12. 12
13. 13
14. 14
15. 15
16. 16
17. 17
18. 18
19. 19
20. 20
21. 21
22. 22
23. 23
24. 24
25. 25
26. 26
27. 27
28. 28
29. 29
30. 30
31. 31
32. 32
33. 33
34. 34
35. 35
36. 36
37. 37
38. 38
39. 39
40. 40

bool EvtDecay::isCharm(EvtId xid){
EvtId d0    = EvtPDL::getId(std::string("D0"));
EvtId d0bar = EvtPDL::getId(std::string("anti-D0"));
EvtId dp    = EvtPDL::getId(std::string("D+"));
EvtId dm    = EvtPDL::getId(std::string("D-"));
EvtId d0h   = EvtPDL::getId(std::string("D0H"));
EvtId d0l    = EvtPDL::getId(std::string("D0L"));
EvtId dstp   = EvtPDL::getId(std::string("D*+"));
EvtId dstm   = EvtPDL::getId(std::string("D*-"));
EvtId ds0    = EvtPDL::getId(std::string("D*0"));
EvtId ds0bar = EvtPDL::getId(std::string("anti-D*0"));
EvtId dsp    = EvtPDL::getId(std::string("D_s+"));
EvtId dsm    = EvtPDL::getId(std::string("D_s-"));
EvtId dsstp  = EvtPDL::getId(std::string("D_s*+"));
EvtId dsstm  = EvtPDL::getId(std::string("D_s*-"));
EvtId ds0stp = EvtPDL::getId(std::string("D_s0*+"));
EvtId ds0stm = EvtPDL::getId(std::string("D_s0*-"));

 std::vector<EvtId> Vid; Vid.clear();
 Vid.push_back(d0);
 Vid.push_back(d0bar);
 Vid.push_back(dp);
 Vid.push_back(dm);
 Vid.push_back(d0h);
 Vid.push_back(d0l);
 Vid.push_back(dstp);
 Vid.push_back(dstm);
 Vid.push_back(ds0);
 Vid.push_back(ds0bar );
 Vid.push_back(dsp );
 Vid.push_back(dsm );
 Vid.push_back(dsstp );
 Vid.push_back(dsstm );
 Vid.push_back(ds0stp );
 Vid.push_back(ds0stm );

 bool flag=true;
 for(int i=0;i<Vid.size();i++){ if(xid == Vid[i]) return flag;}
 return false;
}