- // Common "deformation" tendencies
- //
- fstring+="(sin(atan2(gy.fElements,gx.fElements)))++";
- fstring+="(cos(atan2(gy.fElements,gx.fElements)))++";
- //
- fstring+="(sin(atan2(gy.fElements,gx.fElements)*2))++";
- fstring+="(cos(atan2(gy.fElements,gx.fElements)*2))++";
- fstring+="(sin(atan2(gy.fElements,gx.fElements)*3))++";
- fstring+="(cos(atan2(gy.fElements,gx.fElements)*3))++";
- //
- fstring+="(sin(atan2(gy.fElements,gx.fElements)*2))*lxr++";
- fstring+="(cos(atan2(gy.fElements,gx.fElements)*2))*lxr++";
- fstring+="(sin(atan2(gy.fElements,gx.fElements)*3))*lxr++";
- fstring+="(cos(atan2(gy.fElements,gx.fElements)*3))*lxr++";
- //
-
- TString *fit1 =stat.FitPlane(tree,"dt",fstring.Data(),"cutF&&cutCE",chi2,npoints,vec1,mat);
- tree->SetAlias("f1",fit1->Data());
- //
- // Central electrode "deformation"
- //
- fstring+="(side*sin(atan2(gy.fElements,gx.fElements)))++";
- fstring+="(side*cos(atan2(gy.fElements,gx.fElements)))++";
- //
- fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*2))++";
- fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*2))++";
- fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*3))++";
- fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*3))++";
- //
- fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*2))*lxr++";
- fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*2))*lxr++";
- fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*3))*lxr++";
- fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*3))*lxr++";
-
- TString *fit2 =stat.FitPlane(tree,"dt",fstring.Data(),"cutF&&abs(dt-f0)<0.7&&cutCE",chi2,npoints,vec2,mat);
- tree->SetAlias("f2",fit2->Data());