| 191 |
// |
// |
| 192 |
// |
// |
| 193 |
fM = 2. + 0.096; // real position from cbar |
fM = 2. + 0.096; // real position from cbar |
| 194 |
fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc. |
// fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc. |
| 195 |
|
fM1 = 2. - 0.122 + 0.096; // due to calculation of xe1 etc. |
| 196 |
if ( fM1 < 0. ) fM1 = 0.; |
if ( fM1 < 0. ) fM1 = 0.; |
| 197 |
// |
// |
| 198 |
CaloStrip *cs = new CaloStrip(fSimu); |
CaloStrip *cs = new CaloStrip(fSimu); |
| 379 |
printf(" energyxo:.............. %f \n",energyxo); |
printf(" energyxo:.............. %f \n",energyxo); |
| 380 |
printf(" energyye:.............. %f \n",energyye); |
printf(" energyye:.............. %f \n",energyye); |
| 381 |
printf(" energyyo:.............. %f \n",energyyo); |
printf(" energyyo:.............. %f \n",energyyo); |
| 382 |
|
printf(" fXEen_maxplane:........ %f \n",fXEen_maxplane); |
| 383 |
|
printf(" fXOen_maxplane:........ %f \n",fXOen_maxplane); |
| 384 |
|
printf(" fYEen_maxplane:........ %f \n",fYEen_maxplane); |
| 385 |
|
printf(" fYOen_maxplane:........ %f \n",fYOen_maxplane); |
| 386 |
printf(" x0max :.............. %f \n",x0max); |
printf(" x0max :.............. %f \n",x0max); |
| 387 |
printf(" debug :.............. %i \n",debug); |
printf(" debug :.............. %i \n",debug); |
| 388 |
|
|
| 727 |
// event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section |
// event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section |
| 728 |
// |
// |
| 729 |
if ( fXosel ){ |
if ( fXosel ){ |
| 730 |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)) ;iplm++) fXOen_maxplane += en_xop[iplm]; |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)) ;iplm++){ |
| 731 |
|
fXOen_maxplane += en_xop[iplm]; |
| 732 |
|
if ( debug ) printf(" XO iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXOen_maxplane,en_xop[iplm]); |
| 733 |
|
}; |
| 734 |
fEnergyxo = fXOen_maxplane/fConv_rxo; |
fEnergyxo = fXOen_maxplane/fConv_rxo; |
| 735 |
}; |
}; |
| 736 |
}; |
}; |
| 779 |
fXesel = true; |
fXesel = true; |
| 780 |
}; |
}; |
| 781 |
if ( fXesel ){ |
if ( fXesel ){ |
| 782 |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexe+nPl)) ;iplm++) fXEen_maxplane += en_xep[iplm]; |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexe+nPl)) ;iplm++){ |
| 783 |
|
fXEen_maxplane += en_xep[iplm]; |
| 784 |
|
if ( debug ) printf(" XE iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXEen_maxplane,en_xep[iplm]); |
| 785 |
|
}; |
| 786 |
fEnergyxe = fXEen_maxplane/fConv_rxe; |
fEnergyxe = fXEen_maxplane/fConv_rxe; |
| 787 |
}; |
}; |
| 788 |
}; |
}; |
| 997 |
} else { |
} else { |
| 998 |
fXOen_maxplane = clong->Get_E0(); |
fXOen_maxplane = clong->Get_E0(); |
| 999 |
}; |
}; |
| 1000 |
|
fMax_plane = clong->Get_tmax(); |
| 1001 |
fYOen_maxplane = 0.; |
fYOen_maxplane = 0.; |
| 1002 |
fYEen_maxplane = 0.; |
fYEen_maxplane = 0.; |
| 1003 |
fXEen_maxplane = 0.; |
fXEen_maxplane = 0.; |