31 #include "TCompIsoReg.h"
39 fCorrectedMass = Mass;
45 fAngulo.resize(fCorrectedMass.size());
46 fModulo.resize(fCorrectedMass.size());
49 for(
int i = fCorrectedMass.size() - 1; i >= 0; i--) {
50 Convert2Polar(fCorrectedMass[i], fCompRatio[i] - 1.);
51 fAngulo[k] = Convert2Polar.Ang;
53 if(fAngulo[k] < atan((fCompRatio[i] - fCompRatio[i - 1]) / (fCorrectedMass[i] - fCorrectedMass[i - 1]))) {
54 std::cout <<
"error interpolacion polar" << std::endl;
57 fModulo[k] = Convert2Polar.Mod;
60 fAngMax = MaxComponent(fAngulo);
61 fAngMin = MaxComponent(fAngulo);
67 TCompIsoReg::~TCompIsoReg() {
73 double TCompIsoReg::CompRatio(
double Mass) {
74 return CurvaRC->interp(Mass);
77 double TCompIsoReg::Efficiency(
double Mass) {
78 return CurvaEf->interp(Mass);
81 double TCompIsoReg::CompRatioPol(
double Ang) {
82 return CurvaPol->interp(Ang);
85 double TCompIsoReg::AnguloMaximo() {
89 double TCompIsoReg::AnguloMinimo() {
93 #pragma package(smart_init)