| 1016 |
Double_t E22x = ptt->xtr_tof[3];//tr->x[3]; |
Double_t E22x = ptt->xtr_tof[3];//tr->x[3]; |
| 1017 |
Double_t E22y = ptt->ytr_tof[3];//tr->y[3]; |
Double_t E22y = ptt->ytr_tof[3];//tr->y[3]; |
| 1018 |
Double_t E22z = zin[3]; |
Double_t E22z = zin[3]; |
| 1019 |
if ( E11x < 100. && E11y < 100. && E22x < 100. && E22y < 100. ){ |
if ( (E11x < 100. && E11y < 100. && E22x < 100. && E22y < 100.) || ptt->trkseqno != -1 ){ |
| 1020 |
Double_t norm = sqrt(pow(E22x-E11x,2)+pow(E22y-E11y,2)+pow(E22z-E11z,2)); |
Double_t norm = sqrt(pow(E22x-E11x,2)+pow(E22y-E11y,2)+pow(E22z-E11z,2)); |
| 1021 |
Double_t MyAzim = TMath::RadToDeg()*atan(TMath::Abs(E22y-E11y)/TMath::Abs(E22x-E11x)); |
Double_t MyAzim = TMath::RadToDeg()*atan(TMath::Abs(E22y-E11y)/TMath::Abs(E22x-E11x)); |
| 1022 |
if(E22x-E11x>=0 && E22y-E11y <0) MyAzim = 360. - MyAzim; |
if(E22x-E11x>=0 && E22y-E11y <0) MyAzim = 360. - MyAzim; |