| 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 |