| 1113 |
Double_t E22z = zin[3]; |
Double_t E22z = zin[3]; |
| 1114 |
if ( (E11x < 100. && E11y < 100. && E22x < 100. && E22y < 100.) || ptt->trkseqno != -1 ){ |
if ( (E11x < 100. && E11y < 100. && E22x < 100. && E22y < 100.) || ptt->trkseqno != -1 ){ |
| 1115 |
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)); |
| 1116 |
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)); |
| 1117 |
if(E22x-E11x>=0 && E22y-E11y <0) MyAzim = 360. - MyAzim; |
// if(E22x-E11x>=0 && E22y-E11y <0) MyAzim = 360. - MyAzim; |
| 1118 |
if(E22x-E11x>=0 && E22y-E11y >=0) MyAzim = MyAzim; |
// if(E22x-E11x>=0 && E22y-E11y >=0) MyAzim = MyAzim; |
| 1119 |
if(E22x-E11x<0 && E22y-E11y >0) MyAzim = 180. - MyAzim; |
// if(E22x-E11x<0 && E22y-E11y >0) MyAzim = 180. - MyAzim; |
| 1120 |
if(E22x-E11x<0 && E22y-E11y <0) MyAzim = 180. + MyAzim; |
// if(E22x-E11x<0 && E22y-E11y <0) MyAzim = 180. + MyAzim; |
| 1121 |
Px = (E22x-E11x)/norm; |
Px = (E22x-E11x)/norm; |
| 1122 |
Py = (E22y-E11y)/norm; |
Py = (E22y-E11y)/norm; |
| 1123 |
Pz = (E22z-E11z)/norm; |
Pz = (E22z-E11z)/norm; |