| 1516 |
double dr = sqrt(dx*dx+dy*dy); |
double dr = sqrt(dx*dx+dy*dy); |
| 1517 |
double theta=atan(dr/76.81); |
double theta=atan(dr/76.81); |
| 1518 |
// |
// |
| 1519 |
|
if ( xtr_tof[1] > 99. || xtr_tof[5] > 99. || ytr_tof[0] > 99. || ytr_tof[4] > 99. ) theta = 0.; |
| 1520 |
|
for (Int_t ii=0; ii<6; ii++){ |
| 1521 |
|
if ( xtr_tof[ii] > 99. ) xtr_tof[ii] = 0.; |
| 1522 |
|
if ( ytr_tof[ii] > 99. ) ytr_tof[ii] = 0.; |
| 1523 |
|
}; |
| 1524 |
|
// |
| 1525 |
|
|
| 1526 |
|
|
| 1527 |
for( int ii=0; ii<48; ii++ ) { |
for( int ii=0; ii<48; ii++ ) { |
| 1528 |
// |
// |
| 1529 |
// printf(" ii %i beta %f atime %u xtr 1 %f ytr 1 %f adc %f \n",ii,betamean,atime,xtr_tof[0],ytr_tof[0],adc[ii]); |
// printf(" ii %i beta %f atime %u xtr 1 %f ytr 1 %f adc %f \n",ii,betamean,atime,xtr_tof[0],ytr_tof[0],adc[ii]); |
| 1619 |
if(adclin<=0) continue; |
if(adclin<=0) continue; |
| 1620 |
// |
// |
| 1621 |
if ( betamean > 99. ){ |
if ( betamean > 99. ){ |
| 1622 |
|
// eDEDXpmt.AddAt((Float_t)adclin,ii); |
| 1623 |
eDEDXpmt[ii]=(Float_t)adclin; |
eDEDXpmt[ii]=(Float_t)adclin; |
| 1624 |
|
// printf(" AAPMT IS %i dedx is %f vector is %f \n",ii,adclin,eDEDXpmt[ii]); |
| 1625 |
continue; |
continue; |
| 1626 |
}; |
}; |
| 1627 |
// |
// |
| 1642 |
if(dEdx<=0) continue; |
if(dEdx<=0) continue; |
| 1643 |
|
|
| 1644 |
eDEDXpmt[ii]=(Float_t)dEdx; |
eDEDXpmt[ii]=(Float_t)dEdx; |
| 1645 |
|
// eDEDXpmt.AddAt((Float_t)dEdx,ii); |
| 1646 |
|
|
| 1647 |
|
// printf(" PMT IS %i dedx is %f vector is %f \n",ii,dEdx,eDEDXpmt[ii]); |
| 1648 |
|
|
| 1649 |
} //end loop on 48 PMT |
} //end loop on 48 PMT |
| 1650 |
|
|