- if (fGRF!=0) {
- fGRF->Delete();
- fGRF=0;
- }
- if (strncmp(fType,"User",3)==0){
- fGRF= new TF2;
- R__b>>fGRF;
- }
- if (strncmp(fType,"Gauss",3)==0)
- fGRF = new TF2("fun",funGauss2D,-5.,5.,-5.,5.,4);
- if (strncmp(fType,"Cosh",3)==0)
- fGRF = new TF2("fun",funCosh2D,-5.,5.,-5.,5.,4);
- if (strncmp(fType,"Gati",3)==0)
- fGRF = new TF2("fun",funGati2D,-5.,5.,-5.,5.,5);
- //read interpolation parameters
- R__b >>fY1;
- R__b >>fY2;
- R__b >>fNYdiv;
- R__b >>fDStep;
- R__b >>fNPRF;
- R__b >>fNChargeArray;
- if (fChargeArray!=0) delete [] fChargeArray;
- if (fNChargeArray>0) {
- fChargeArray = new Float_t[fNChargeArray];
- R__b.ReadFastArray(fChargeArray,fNChargeArray);
+ if (strncmp(fType,"User",3)!=0){
+ delete fGRF;
+ if (strncmp(fType,"Gauss",3)==0)
+ fGRF = new TF2("funGauss2D",funGauss2D,-5.,5.,-5.,5.,4);
+ if (strncmp(fType,"Cosh",3)==0)
+ fGRF = new TF2("funCosh2D",funCosh2D,-5.,5.,-5.,5.,4);
+ if (strncmp(fType,"Gati",3)==0)
+ fGRF = new TF2("funGati2D",funGati2D,-5.,5.,-5.,5.,5);
+ if (fGRF!=0) fGRF->SetParameters(funParam);