71 |
fpHits.clear();// we don't need it |
fpHits.clear();// we don't need it |
72 |
|
|
73 |
//From this moment all hits sorted by time and present in tubes |
//From this moment all hits sorted by time and present in tubes |
74 |
|
|
75 |
for(Int_t i = 0; i<NTUBES; i++){ |
for(Int_t i = 0; i<NTUBES; i++){ |
76 |
if( tubeHits[i].size() > 0 ){ |
if( tubeHits[i].size() > 0 ){ |
77 |
for(UInt_t j = 0; j<tubeHits[i].size(); j++){ |
for(UInt_t j = 0; j<tubeHits[i].size(); j++){ |
86 |
Uback += keVtomV*(hit0->GetEREL())*1.e6*Exp((hit0->GetTOF()-/*hit1->GetTOF()*/hit->GetTOF())/Abs(tauHits[i].at(k-1))); |
Uback += keVtomV*(hit0->GetEREL())*1.e6*Exp((hit0->GetTOF()-/*hit1->GetTOF()*/hit->GetTOF())/Abs(tauHits[i].at(k-1))); |
87 |
} |
} |
88 |
|
|
89 |
tauHits[i].push_back(gRandom->Gaus(Tau,SigmaTau)); |
|
90 |
|
tauHits[i].push_back(rnd->Gaus(Tau,SigmaTau)); |
91 |
if (DEBUG){ |
if (DEBUG){ |
92 |
cout<<"++++++HIT+++++++"<<endl; |
cout<<"++++++HIT+++++++"<<endl; |
93 |
cout<<"TubeNo="<<hit->GetPOS()<<endl |
cout<<"TubeNo="<<hit->GetPOS()<<endl |