33 {
34 barrelHit1.clear();
35 barrelHit2.clear();
36 ecapHit.clear();
37 wcapHit.clear();
38
39 for ( std::vector<int>::iterator
iter = vtofHit.begin();
iter != vtofHit.end();
iter++ )
40 {
41 int tofId = ( *iter );
42 int part, layer, im;
43 part = int( tofId / 10000. );
44 layer = int( ( tofId - part * 10000 ) / 1000. );
45 im = int( ( tofId - part * 10000 - layer * 1000 ) / 10. );
46 if ( part == 0 )
47 {
48
49 if ( find( ecapHit.begin(), ecapHit.end(), im ) == ecapHit.end() )
50 { ecapHit.push_back( im ); }
51 }
52 else if ( part == 2 )
53 {
54
55 if ( find( wcapHit.begin(), wcapHit.end(), im ) == wcapHit.end() )
56 { wcapHit.push_back( im ); }
57 }
58 else
59 {
60 if ( layer == 0 )
61 {
62
63 if ( find( barrelHit1.begin(), barrelHit1.end(), im ) == barrelHit1.end() )
64 { barrelHit1.push_back( im ); }
65 }
66 if ( layer == 1 )
67 {
68
69 if ( find( barrelHit2.begin(), barrelHit2.end(), im ) == barrelHit2.end() )
70 { barrelHit2.push_back( im ); }
71 }
72 }
73 }
74}