31 {
32 if ( m_part_type == "e" )
33 {
35 trk.setPx( mdcTrk->px() );
36 trk.setPy( mdcTrk->py() );
37 trk.setPz( mdcTrk->pz() );
38 double p3( mdcTrk->p() );
39 trk.setE( sqrt( p3 * p3 +
mmass[0] *
mmass[0] ) );
40 return trk;
41 }
42 if ( m_part_type == "muon" )
43 {
45 trk.setPx( mdcTrk->px() );
46 trk.setPy( mdcTrk->py() );
47 trk.setPz( mdcTrk->pz() );
48 double p3( mdcTrk->p() );
49 trk.setE( sqrt( p3 * p3 +
mmass[1] *
mmass[1] ) );
50 return trk;
51 }
52 if ( m_part_type == "pion" )
53 {
55 trk.setPx( mdcTrk->px() );
56 trk.setPy( mdcTrk->py() );
57 trk.setPz( mdcTrk->pz() );
58 double p3( mdcTrk->p() );
59 trk.setE( sqrt( p3 * p3 +
mmass[2] *
mmass[2] ) );
60 return trk;
61 }
62 if ( m_part_type == "kion" )
63 {
65 trk.setPx( mdcTrk->px() );
66 trk.setPy( mdcTrk->py() );
67 trk.setPz( mdcTrk->pz() );
68 double p3( mdcTrk->p() );
69 trk.setE( sqrt( p3 * p3 +
mmass[3] *
mmass[3] ) );
70 return trk;
71 }
72 if ( m_part_type == "proton" )
73 {
75 trk.setPx( mdcTrk->px() );
76 trk.setPy( mdcTrk->py() );
77 trk.setPz( mdcTrk->pz() );
78 double p3( mdcTrk->p() );
79 trk.setE( sqrt( p3 * p3 +
mmass[4] *
mmass[4] ) );
80 return trk;
81 }
82 else return trk;
83}