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