//Macro for data. //Needs to be compiled with .L ....C++ in root //Author: Tiger #include "TCanvas.h" #include "TFile.h" #include "TH2F.h" #include "TH3F.h" #include "TClonesArray.h" #include "TGeoManager.h" #include "TTree.h" #include "TLine.h" #include "TProfile.h" #include "TObjArray.h" #include "TPolyMarker.h" #include "TPolyMarker3D.h" #include "TPolyLine.h" #include "TVector3.h" #include "TBox.h" #include "TSystem.h" #include "TROOT.h" #include "TMinuit.h" #include "TView3D.h" #include "TApplication.h" #include "TMath.h" #include "TStyle.h" #include "TString.h" #include "TColor.h" #include "THStack.h" #include "TVector3.h" #include "tpc/par/TpcAlignmentManager.h" #include "tpc/base/TpcCluster.h" #include "tpc/base/TpcDigi.h" #include "tpc/base/TpcSample.h" #include #include #include #include #include #include #include Double_t t2kconv(Double_t * x , Double_t*par); Double_t t2k(Double_t * x , Double_t*par); const bool DEBUG = false; void TestIntPSA(TString filename, unsigned int EvMax=0) { while(gROOT->GetListOfCanvases()->GetSize()) delete gROOT->GetListOfCanvases()->At(0); gStyle->SetOptStat(10); gStyle->SetPalette(1); int a = 0; TFile* recofile = new TFile(filename); if(recofile->IsZombie()) { std::cerr<<"Reco file not existing! Aborting."<Get("cbmsim"); if(!recotree) { std::cerr<<"Reco Tree not existing! Aborting."< Tbranches; //I CAN HAZ ARRAYS :3 TClonesArray* tpcArr=new TClonesArray("TpcCluster"); recotree->SetBranchAddress("TpcCluster", &tpcArr); Tbranches["TpcCluster"] = tpcArr; if(tpcArr==0) {std::cout<< "Init TpcCluster Array not found!"<< std::endl;return;} //Histograms diffinitions: TH2F * hsigt0max = new TH2F("hsigt0max","hsigt0max",40,-20,20,200,0,1.2); TH2F * hsigt0fit = new TH2F("hsigt0fit","hsigt0fit",200,-20,20,200,0,1.2); TH2F * hsigt0diff = new TH2F("hsigt0diff","hsigt0diff",200,-20,20,200,0,1.2); TH2F * hsigt0 = new TH2F("hsigt0","hsigt0",40,-20,20,200,0,1.2); TH2F * hsigt0cf = new TH2F("hsigt0cf","hsigt0cf",200,-20,20,200,0,1.2); TH1F * ht0fitt0max = new TH1F("ht0fitt0max","ht0fitt0max",200,-5,5); TH1F * ht0fitt0diff = new TH1F("ht0fitt0diff","ht0fitt0diff",200,-5,5); TH1F * hsigma = new TH1F ("hsigma","hsigma",500,0,500); TH1F * ht0max = new TH1F ("ht0max","ht0max",2000,0,511); TH1F * ht0cf = new TH1F ("ht0cf","ht0cf",2000,0,511); TH1F * ht0diff = new TH1F ("ht0diff","ht0diff",200,0,511); TH1F * ht0fit = new TH1F ("ht0fit","ht0fit",2000,0,511); unsigned int MAX =recotree->GetEntries(); std::cout << "Number of events: " << MAX <GetEvent(k); bool paded[11000]; for(int i =0;i<11000;++i) paded[i] = false; int ntpccl = tpcArr->GetEntries(); for(int c=0; cAt(c); //if(cl->amp()<12000) continue; // std::cout << "Number of digi: " << cl->nDigi() <nDigi();di++) { const TpcDigi * curdi = cl->getDigi(di); //if(curdi->nSample()<5) continue; //if(curdi->t()<100) continue; TH1F * hsig = new TH1F("hsig","hsig",511,-0.5,510.5); TF1 * f2 = new TF1("f1t2k",t2k,0,511,3); f2->SetLineColor(kBlue); f2->SetParNames("A","t0","sigmag"); f2->SetParLimits(0,0,2000); f2->SetParLimits(1,0,511); f2->SetParLimits(2,0,500); f2->SetNpx(5000); double t0max = curdi->t(); double maxamp = 0;//curdi->maxamp(); double firstsa = 511;double ampt1 = 0.; double diffmax = 0.; double sample1 = -2.; double sample2 = -2.; double ci_1 =-10000.; double t0=-1.;double t0cf=-1.;double t0diff = -1; int DELAY = 2; //for CF // std::cout << "Number of sampel" << curdi->nSample() << std::endl; for (int sa = 0; sanSample();sa++) { const TpcSample * cursa = curdi->getSample(sa); if(cursa->t()t(); if(t0<0) t0 = cursa->t(); hsig->Fill(cursa->t(),cursa->amp()); if(cursa->amp()>maxamp) maxamp = cursa->amp(); if(ampt1 && diffmaxamp()-ampt1 && cursa->t()amp()-ampt1 << " at (" << cursa->t() << "," << cursa->amp() <<") t0= " <t()-1.0*(cursa->amp()/(cursa->amp()-ampt1)) << std::endl; diffmax = cursa->amp()-ampt1; t0diff = cursa->t()-(cursa->amp()/diffmax);} ampt1=cursa->amp(); //CF if(sa>=DELAY) { const TpcSample * prevsa = curdi->getSample(sa-DELAY); sample1=prevsa->amp(); double ci = cursa->amp()-2.0*sample1; if(sample1amp())//rising edge if (ci_1>0) if (ci<0) { t0cf=ci_1/(ci_1-ci)+cursa->t()-1; } ci_1=ci; } } // if(t0diff<0) //|| paded[curdi->padId()]) // {delete(hsig);delete(f2);continue;} //t0diff = t0max; paded[curdi->padId()] = true; if(DEBUG)std::cout << "T0DIFF == "<< t0diff<< " \n ---------------------------" <SetParameters(maxamp/0.835065,firstsa,30.); TFitResultPtr r; int nfit = 0; do{ r = hsig->Fit(f2, "SQ"); nfit++; }while(((int)(r)&&nfit<2)); if(r==0)//Fit OK { double sigma = f2->GetParameter(2);//f2->GetParameter(0); hsigma->Fill(sigma); double amp = curdi->amp(); double ampc; } double t0fit = f2->GetParameter(1); bool plot=false; if(t0diff>0 && t0max>0 && t0fit>0 && t0cf>0)//time ok if(maxamp) for (int sa = 0; sanSample();sa++) { const TpcSample * cursa = curdi->getSample(sa); hsigt0fit->Fill(cursa->t()-t0fit,cursa->amp()/maxamp); hsigt0->Fill(cursa->t()-t0,cursa->amp()/maxamp); if(t0diff!=t0max) hsigt0diff->Fill(cursa->t()-t0diff,cursa->amp()/maxamp); hsigt0max->Fill(cursa->t()-t0max,cursa->amp()/maxamp); // if (t0diff == 3.5) plot=true; if(t0cf>0) hsigt0cf->Fill(cursa->t()-t0cf,cursa->amp()/maxamp); } ht0fitt0max->Fill(t0max-t0fit); ht0fitt0diff->Fill(t0max-t0diff); ht0max->Fill(t0max); ht0fit->Fill(t0fit); ht0cf->Fill(t0cf); ht0diff->Fill(t0diff); if(1)//DISPLAY (plot) { TCanvas * cx = new TCanvas("cx","cx",1200,600); hsig->Draw(); TLine * lt0max = new TLine(t0max,0,t0max,maxamp); lt0max->SetLineWidth(2.); lt0max->SetLineColor(kBlack);lt0max->Draw("same"); TLine * lt0fit = new TLine(t0fit,0,t0fit,maxamp); lt0fit->SetLineWidth(2.); lt0fit->SetLineColor(kBlue);lt0fit->Draw("same"); TLine * lt0diff = new TLine(t0diff,0,t0diff,maxamp); lt0diff->SetLineWidth(2.); lt0diff->SetLineColor(kRed);lt0diff->Draw("same"); TLine * lt0cf = new TLine(t0cf,0,t0cf,maxamp); lt0cf->SetLineWidth(2.); lt0cf->SetLineColor(kCyan);lt0cf->Draw("same"); std::cout << "f2->SetParameters("<GetChisquare()<< "/" << f2->GetNDF() << "=" << f2->GetChisquare()/f2->GetNDF() << " sigma = " << f2->GetParameter(2) << "ns" << std::endl; else std::cout << "Fit Failed" << std::endl; cx->WaitPrimitive(); delete(cx); delete(lt0fit); delete(lt0max); } delete(f2); delete(hsig); }//digi l00p }//cluster l00p }//end event loop std::cout << "DRAWING " << std::endl; gROOT->SetStyle("Plain"); gStyle->SetPalette(1); const Int_t NCont = 250; const Int_t NRGBs = 5; Double_t stops[NRGBs] = { 0.00, 0.34, 0.61, 0.84, 1.00 }; Double_t red[NRGBs] = { 0.00, 0.00, 0.87, 1.00, 0.51 }; Double_t green[NRGBs] = { 0.00, 0.81, 1.00, 0.20, 0.00 }; Double_t blue[NRGBs] = { 0.51, 1.00, 0.12, 0.00, 0.00 }; TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont); TCanvas * canvconv= new TCanvas("canvconv","canvconv", 1000,1000); canvconv->Divide(2,2); canvconv->cd(1); hsigt0fit->Draw("colz"); canvconv->cd(2); hsigt0max->Draw("colz"); canvconv->cd(3); hsigt0diff->Draw("colz"); canvconv->cd(4); // hsigt0->Draw("colz"); hsigt0cf->Draw("colz"); // ht0fitt0max->Draw(); // ht0fitt0diff->SetLineColor(kRed); // ht0fitt0diff->Draw("same"); TCanvas * canv= new TCanvas("canv","canv", 1000,800); hsigma->Draw(); //ht0diff->SetLineColor(kRed); // ht0diff->Draw(); // ht0max->Draw("same"); // ht0fit->SetLineColor(kBlue); // ht0fit->Draw("same"); //canvconv->SaveAs("all.root"); //canvconv->SaveAs("all.eps"); return; }//end of macro Double_t t2k(Double_t * x , Double_t*par){ double y = (x[0]-par[1])*50.; double A =par[0]; double t =par[2];//time constant in ns if(y<0) return 0.; else return A*pow((y/t),3)*sin(y/(3*t))*exp(-1.*y/t); } Double_t t2kconv(Double_t * x , Double_t*par) { double y = (x[0]-par[1])*50.;//in ns double A = par[0]; double t = 2.05*20; double b = par[2]*50.;//mean gaussian in ns double s = par[3]*50.;//mean sigma in ns double E = exp(1); double Pi = TMath::Pi(); if(y<-50||y>1000) return 0.; else { double sss = (pow(E,pow(s,2)/(2.*pow(t,2)) + b/t - y/t)* ((105*sqrt(Pi/2.))/pow(pow(s,-2),4.5) + (pow(b,8)*sqrt(Pi/2.))/sqrt(pow(s,-2)) + (sqrt(Pi/2.)*pow(s,16))/(sqrt(pow(s,-2))*pow(t,8)) + (4*b*sqrt(2*Pi)*pow(s,14))/(sqrt(pow(s,-2))*pow(t,7)) + (14*pow(b,2)*sqrt(2*Pi)*pow(s,12))/ (sqrt(pow(s,-2))*pow(t,6)) + (28*pow(b,3)*sqrt(2*Pi)*pow(s,10))/ (sqrt(pow(s,-2))*pow(t,5)) + (35*pow(b,4)*sqrt(2*Pi)*pow(s,8))/ (sqrt(pow(s,-2))*pow(t,4)) + (28*pow(b,5)*sqrt(2*Pi)*pow(s,6))/ (sqrt(pow(s,-2))*pow(t,3)) + (14*pow(b,6)*sqrt(2*Pi)*pow(s,4))/ (sqrt(pow(s,-2))*pow(t,2)) + (4*pow(b,7)*sqrt(2*Pi)*pow(s,2))/(sqrt(pow(s,-2))*t) - (384*pow(s,8)*(pow(s,2) + b*t))/t + (210*sqrt(2*Pi)*pow(pow(s,2) + b*t,2))/ (pow(pow(s,-2),3.5)*pow(t,2)) - (448*pow(s,6)*pow(pow(s,2) + b*t,3))/pow(t,3) + (105*sqrt(2*Pi)*pow(pow(s,2) + b*t,4))/ (pow(pow(s,-2),2.5)*pow(t,4)) - (112*pow(s,4)*pow(pow(s,2) + b*t,5))/pow(t,5) + (14*sqrt(2*Pi)*pow(pow(s,2) + b*t,6))/ (pow(pow(s,-2),1.5)*pow(t,6)) - (8*pow(s,2)*pow(pow(s,2) + b*t,7))/pow(t,7) + (8*(-1 + pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*pow(s,2)* pow(pow(s,2) + b*t,7))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (112*pow(s,4)*pow(pow(s,2) + b*t,5)* (pow(s,2)/2. + b*t - (-1 + pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))) - pow(b,2)/(2.*pow(s,2)))* pow(t,2)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (224*pow(s,6)*pow(pow(s,2) + b*t,3)* (pow(s,4)/4. + b*pow(s,2)*t + (1 + (3*pow(b,2))/(2.*pow(s,2)))*pow(s,2)*pow(t,2) + 2*b*(1 + pow(b,2)/(2.*pow(s,2)))*pow(t,3) - (-2 + 2*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))) - pow(b,4)/(4.*pow(s,4)) - pow(b,2)/pow(s,2))* pow(t,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (64*pow(s,8)*(pow(s,2) + b*t)* (pow(s,6)/8. + (3*b*pow(s,4)*t)/4. + (3*(1 + (5*pow(b,2))/(2.*pow(s,2)))*pow(s,4)*pow(t,2))/ 4. + b*(3 + (5*pow(b,2))/(2.*pow(s,2)))*pow(s,2)* pow(t,3) + (3*(2 + (5*pow(b,4))/(4.*pow(s,4)) + (3*pow(b,2))/pow(s,2))*pow(s,2)*pow(t,4))/2. + 3*b*(2 + pow(b,4)/(4.*pow(s,4)) + pow(b,2)/pow(s,2))* pow(t,5) - (-6 + 6*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))) - pow(b,6)/(8.*pow(s,6)) - (3*pow(b,4))/(4.*pow(s,4)) - (3*pow(b,2))/pow(s,2))* pow(t,6)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (14*sqrt(2)*pow(b + pow(s,2)/t,7)* pow(pow(s,2) + b*t,2)* (15*sqrt(Pi) - (sqrt(2)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* (15*pow(t,4) + (5*pow(t,2)*pow(pow(s,2) + b*t,2))/pow(s,2) + pow(pow(s,2) + b*t,4)/pow(s,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,4)) + 15*sqrt(Pi)* (-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (pow(pow(b + pow(s,2)/t,2)/pow(s,2),3.5)*pow(t,2)) - (pow(b + pow(s,2)/t,9)* (105*sqrt(Pi) - (sqrt(2)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* (105*pow(s,6)*pow(t,6) + 35*pow(s,4)*pow(t,4)*pow(pow(s,2) + b*t,2) + 7*pow(s,2)*pow(t,2)*pow(pow(s,2) + b*t,4) + pow(pow(s,2) + b*t,6)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(s,6)*pow(t,6)) + 105*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))))/ (sqrt(2)*pow(pow(b + pow(s,2)/t,2)/pow(s,2),4.5)) - (pow(b,8)*sqrt(Pi/2.)*(b + pow(s,2)/t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/sqrt(pow(b + pow(s,2)/t,2)/pow(s,2)) - (sqrt(Pi/2.)*pow(s,16)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,9)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (4*b*sqrt(2*Pi)*pow(s,14)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,8)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (14*pow(b,2)*sqrt(2*Pi)*pow(s,12)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,7)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (28*pow(b,3)*sqrt(2*Pi)*pow(s,10)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,6)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (35*pow(b,4)*sqrt(2*Pi)*pow(s,8)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,5)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (28*pow(b,5)*sqrt(2*Pi)*pow(s,6)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,4)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (14*pow(b,6)*sqrt(2*Pi)*pow(s,4)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,3)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) - (4*pow(b,7)*sqrt(2*Pi)*pow(s,2)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,2)*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) + (14*sqrt(2)*pow(pow(s,2) + b*t,9)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* sqrt(Pi)*TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,9)*pow(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)),1.5)) + (70*pow(s,6)*pow(pow(s,2) + b*t,3)* ((pow(pow(s,2) + b*t,2)* (pow(s,2) + 2*b*t + (3 + pow(b,2)/pow(s,2))*pow(t,2)) )/pow(s,2) - 3*pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* sqrt(Pi/2.)*pow(t,4)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7))) + (4860* pow(E,(pow(s,2) + 2*b*t - 2*t*y)/(2.*pow(t,2)))*sqrt(2*Pi)* pow(t,3)*pow(y,4)*(t* sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2))) - sqrt(pow(s,-2))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (sqrt(pow(s,-2))*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (19440*pow(t,2)*pow(y,3)* (-(t*(-2*sqrt(pow(s,-2))*pow(s,2)*t + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*(pow(s,2) + b*t))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2))))\ + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,2) + b*t)*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (540*pow(y,5)*(-(t*(-2*sqrt(pow(s,-2))*pow(s,2)*t + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*(pow(s,2) + b*t))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2))))\ + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,2) + b*t)*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (4*pow(y,7)*(-(t*(-2*sqrt(pow(s,-2))*pow(s,2)*t + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*(pow(s,2) + b*t))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2))))\ + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,2) + b*t)*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,2)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (19440*y*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t)) ) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,4) + pow(s,2)*t*(2*b + 2*t + y) + pow(t,2)*(pow(b,2) + b*y + pow(y,2)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (1800*pow(y,3)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t)) ) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,4) + pow(s,2)*t*(2*b + 2*t + y) + pow(t,2)*(pow(b,2) + b*y + pow(y,2)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,2)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (28*pow(y,5)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t)) ) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,4) + pow(s,2)*t*(2*b + 2*t + y) + pow(t,2)*(pow(b,2) + b*y + pow(y,2)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,6) + pow(b,3)*pow(t,3) + 3*pow(s,4)*t*(b + t) + 3*b*pow(s,2)*pow(t,2)*(b + t))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,4)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (29160*t*pow(y,2)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))\ - 2*sqrt(pow(s,-2))*pow(s,2)*t*(pow(s,2) + t*(b + y))) - pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (1350*pow(y,4)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))\ - 2*sqrt(pow(s,-2))*pow(s,2)*t*(pow(s,2) + t*(b + y))) - pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*t* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (14*pow(y,6)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))\ - 2*sqrt(pow(s,-2))*pow(s,2)*t*(pow(s,2) + t*(b + y))) - pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + t))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,3)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (4860*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)) )) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,6) + pow(s,4)*t*(3*b + 5*t + y) + pow(s,2)*pow(t,2)* (3*pow(b,2) + 5*b*t + 2*b*y + 3*t*y + pow(y,2)) + pow(t,3)*(pow(b,3) + pow(b,2)*y + b*pow(y,2) + pow(y,3)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*t* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (1350*pow(y,2)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)) )) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,6) + pow(s,4)*t*(3*b + 5*t + y) + pow(s,2)*pow(t,2)* (3*pow(b,2) + 5*b*t + 2*b*y + 3*t*y + pow(y,2)) + pow(t,3)*(pow(b,3) + pow(b,2)*y + b*pow(y,2) + pow(y,3)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,3)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (35*pow(y,4)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)) )) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,6) + pow(s,4)*t*(3*b + 5*t + y) + pow(s,2)*pow(t,2)* (3*pow(b,2) + 5*b*t + 2*b*y + 3*t*y + pow(y,2)) + pow(t,3)*(pow(b,3) + pow(b,2)*y + b*pow(y,2) + pow(y,3)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,8) + pow(b,4)*pow(t,4) + 2*pow(s,6)*t*(2*b + 3*t) + 2*pow(b,2)*pow(s,2)*pow(t,3)*(2*b + 3*t) + 3*pow(s,4)*pow(t,2)*(2*pow(b,2) + 4*b*t + pow(t,2)))* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,5)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (540*y*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,10) + pow(b,5)*pow(t,5) + 5*pow(s,8)*t*(b + 2*t) + 5*pow(b,3)*pow(s,2)*pow(t,4)*(b + 2*t) + 5*pow(s,6)*pow(t,2)* (2*pow(b,2) + 6*b*t + 3*pow(t,2)) + 5*b*pow(s,4)*pow(t,3)* (2*pow(b,2) + 6*b*t + 3*pow(t,2)))) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,8) + pow(s,6)*t*(4*b + 9*t + y) + pow(t,4)*(pow(b,4) + pow(b,3)*y + pow(b,2)*pow(y,2) + b*pow(y,3) + pow(y,4)) + pow(s,4)*pow(t,2)* (6*pow(b,2) + 8*pow(t,2) + 7*t*y + pow(y,2) + 3*b*(6*t + y)) + pow(s,2)*pow(t,3)* (4*pow(b,3) + 3*pow(b,2)*(3*t + y) + pow(y,2)*(4*t + y) + b*y*(7*t + 2*y)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,10) + pow(b,5)*pow(t,5) + 5*pow(s,8)*t*(b + 2*t) + 5*pow(b,3)*pow(s,2)*pow(t,4)*(b + 2*t) + 5*pow(s,6)*pow(t,2)*(2*pow(b,2) + 6*b*t + 3*pow(t,2)) + 5*b*pow(s,4)*pow(t,3)*(2*pow(b,2) + 6*b*t + 3*pow(t,2))) *(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,4)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (28*pow(y,3)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,10) + pow(b,5)*pow(t,5) + 5*pow(s,8)*t*(b + 2*t) + 5*pow(b,3)*pow(s,2)*pow(t,4)*(b + 2*t) + 5*pow(s,6)*pow(t,2)* (2*pow(b,2) + 6*b*t + 3*pow(t,2)) + 5*b*pow(s,4)*pow(t,3)* (2*pow(b,2) + 6*b*t + 3*pow(t,2)))) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,8) + pow(s,6)*t*(4*b + 9*t + y) + pow(t,4)*(pow(b,4) + pow(b,3)*y + pow(b,2)*pow(y,2) + b*pow(y,3) + pow(y,4)) + pow(s,4)*pow(t,2)* (6*pow(b,2) + 8*pow(t,2) + 7*t*y + pow(y,2) + 3*b*(6*t + y)) + pow(s,2)*pow(t,3)* (4*pow(b,3) + 3*pow(b,2)*(3*t + y) + pow(y,2)*(4*t + y) + b*y*(7*t + 2*y)))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,10) + pow(b,5)*pow(t,5) + 5*pow(s,8)*t*(b + 2*t) + 5*pow(b,3)*pow(s,2)*pow(t,4)*(b + 2*t) + 5*pow(s,6)*pow(t,2)*(2*pow(b,2) + 6*b*t + 3*pow(t,2)) + 5*b*pow(s,4)*pow(t,3)*(2*pow(b,2) + 6*b*t + 3*pow(t,2))) *(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,6)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) + (90*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,12) + pow(b,6)*pow(t,6) + 3*pow(s,10)*t*(2*b + 5*t) + 3*pow(b,4)*pow(s,2)*pow(t,5)*(2*b + 5*t) + 15*pow(s,8)*pow(t,2)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 15*pow(b,2)*pow(s,4)*pow(t,4)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 5*pow(s,6)*pow(t,3)* (4*pow(b,3) + 18*pow(b,2)*t + 18*b*pow(t,2) + 3*pow(t,3)))) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,10) + pow(s,8)*t*(5*b + 14*t + y) + pow(s,6)*pow(t,2)* (10*pow(b,2) + 42*b*t + 33*pow(t,2) + 4*b*y + 12*t*y + pow(y,2)) + pow(t,5)*(pow(b,5) + pow(b,4)*y + pow(b,3)*pow(y,2) + pow(b,2)*pow(y,3) + b*pow(y,4) + pow(y,5)) + pow(s,2)*pow(t,4)* (5*pow(b,4) + 3*pow(b,2)*y*(4*t + y) + pow(y,3)*(5*t + y) + 2*pow(b,3)*(7*t + 2*y) + b*pow(y,2)*(9*t + 2*y)) + pow(s,4)*pow(t,3)* (10*pow(b,3) + 6*pow(b,2)*(7*t + y) + 3*b*(11*pow(t,2) + 8*t*y + pow(y,2)) + y*(15*pow(t,2) + 9*t*y + pow(y,2))))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,12) + pow(b,6)*pow(t,6) + 3*pow(s,10)*t*(2*b + 5*t) + 3*pow(b,4)*pow(s,2)*pow(t,5)*(2*b + 5*t) + 15*pow(s,8)*pow(t,2)*(pow(b,2) + 4*b*t + 3*pow(t,2)) + 15*pow(b,2)*pow(s,4)*pow(t,4)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 5*pow(s,6)*pow(t,3)* (4*pow(b,3) + 18*pow(b,2)*t + 18*b*pow(t,2) + 3*pow(t,3)))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,5)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (14*pow(y,2)*(t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (-(pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)* (pow(s,12) + pow(b,6)*pow(t,6) + 3*pow(s,10)*t*(2*b + 5*t) + 3*pow(b,4)*pow(s,2)*pow(t,5)*(2*b + 5*t) + 15*pow(s,8)*pow(t,2)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 15*pow(b,2)*pow(s,4)*pow(t,4)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 5*pow(s,6)*pow(t,3)* (4*pow(b,3) + 18*pow(b,2)*t + 18*b*pow(t,2) + 3*pow(t,3)))) + 2*sqrt(pow(s,-2))*pow(s,2)*t* (pow(s,10) + pow(s,8)*t*(5*b + 14*t + y) + pow(s,6)*pow(t,2)* (10*pow(b,2) + 42*b*t + 33*pow(t,2) + 4*b*y + 12*t*y + pow(y,2)) + pow(t,5)*(pow(b,5) + pow(b,4)*y + pow(b,3)*pow(y,2) + pow(b,2)*pow(y,3) + b*pow(y,4) + pow(y,5)) + pow(s,2)*pow(t,4)* (5*pow(b,4) + 3*pow(b,2)*y*(4*t + y) + pow(y,3)*(5*t + y) + 2*pow(b,3)*(7*t + 2*y) + b*pow(y,2)*(9*t + 2*y)) + pow(s,4)*pow(t,3)* (10*pow(b,3) + 6*pow(b,2)*(7*t + y) + 3*b*(11*pow(t,2) + 8*t*y + pow(y,2)) + y*(15*pow(t,2) + 9*t*y + pow(y,2))))) + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(2*Pi)*sqrt(pow(s,-2))* (pow(s,12) + pow(b,6)*pow(t,6) + 3*pow(s,10)*t*(2*b + 5*t) + 3*pow(b,4)*pow(s,2)*pow(t,5)*(2*b + 5*t) + 15*pow(s,8)*pow(t,2)*(pow(b,2) + 4*b*t + 3*pow(t,2)) + 15*pow(b,2)*pow(s,4)*pow(t,4)* (pow(b,2) + 4*b*t + 3*pow(t,2)) + 5*pow(s,6)*pow(t,3)* (4*pow(b,3) + 18*pow(b,2)*t + 18*b*pow(t,2) + 3*pow(t,3)))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(b - y,2)/(2.*pow(s,2)))*sqrt(pow(s,-2))*pow(t,7)* sqrt(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)))) - (sqrt(2)*pow(E,pow(s,2)/(2.*pow(t,2)) + b/t - y/t)* ((4*sqrt(2)*(-1 + pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2))))* sqrt(pow(s,-2))*pow(s,2)*pow(pow(s,2) + b*t,7))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (4*sqrt(2)* (-1 + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2))))*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,7))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,7)) - (28*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)*pow(pow(s,2) + b*t,5)* (pow(s,4) + pow(b,2)*pow(t,2) + 2*pow(s,2)*t*(b + t - pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*t)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (28*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,3)* (pow(s,8) + pow(b,4)*pow(t,4) + 4*pow(s,6)*t*(b + t) + 4*pow(b,2)*pow(s,2)*pow(t,3)*(b + t) + 2*pow(s,4)*pow(t,2)* (3*pow(b,2) + 4*b*t - 4*(-1 + pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2))))* pow(t,2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) - (4*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)* (pow(s,2) + b*t)* (pow(s,12) + pow(b,6)*pow(t,6) + 6*pow(s,10)*t*(b + t) + 6*pow(b,4)*pow(s,2)*pow(t,5)*(b + t) + 3*pow(s,8)*pow(t,2)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 3*pow(b,2)*pow(s,4)*pow(t,4)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 4*pow(s,6)*pow(t,3)* (5*pow(b,3) + 9*pow(b,2)*t + 12*b*pow(t,2) - 12*(-1 + pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2))))* pow(t,3))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) + (28*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,5)* (pow(s,4) + pow(t,2)*pow(b - y,2) + 2*pow(s,2)*t*(b + t - pow(E,pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2)))*t - y)))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,7)) + (28*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)*pow(pow(s,2) + b*t,3)* (pow(s,8) + pow(t,4)*pow(b - y,4) + 4*pow(s,6)*t*(b + t - y) + 4*pow(s,2)*pow(t,3)*pow(b - y,2)*(b + t - y) + 2*pow(s,4)*pow(t,2)* (3*pow(b,2) - 4*(-1 + pow(E, pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,2) + b*(4*t - 6*y) - 4*t*y + 3*pow(y,2))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,7)) + (4*sqrt(2)*sqrt(pow(s,-2))*pow(s,2)*(pow(s,2) + b*t)* (pow(s,12) + pow(t,6)*pow(b - y,6) + 6*pow(s,10)*t*(b + t - y) + 6*pow(s,2)*pow(t,5)*pow(b - y,4)*(b + t - y) + 3*pow(s,8)*pow(t,2)* (5*pow(b,2) + 8*pow(t,2) + 2*b*(4*t - 5*y) - 8*t*y + 5*pow(y,2)) + 3*pow(s,4)*pow(t,4)*pow(b - y,2)* (5*pow(b,2) + 8*pow(t,2) + 2*b*(4*t - 5*y) - 8*t*y + 5*pow(y,2)) + 4*pow(s,6)*pow(t,3)* (5*pow(b,3) - 12*(-1 + pow(E, pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,3) + 3*pow(b,2)*(3*t - 5*y) - 12*pow(t,2)*y + 9*t*pow(y,2) - 5*pow(y,3) + 3*b*(4*pow(t,2) - 6*t*y + 5*pow(y,2)))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,7)) - (14*sqrt(pow(s,-2))*pow(pow(s,2) + b*t,9)* (15*sqrt(Pi) - (sqrt(2)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* (15*pow(s,4)*pow(t,4) + 5*pow(s,2)*pow(t,2)*pow(pow(s,2) + b*t,2) + pow(pow(s,2) + b*t,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(s,4)*pow(t,4)) + 15*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (pow(t,9)*pow(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)),3.5)) - (sqrt(pow(s,-2))*pow(pow(s,2) + b*t,9)* (105*sqrt(Pi) - (sqrt(2)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* (105*pow(s,6)*pow(t,6) + 35*pow(s,4)*pow(t,4)*pow(pow(s,2) + b*t,2) + 7*pow(s,2)*pow(t,2)*pow(pow(s,2) + b*t,4) + pow(pow(s,2) + b*t,6)))/ (pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(s,6)*pow(t,6)) + 105*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (2.*pow(t,9)*pow(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)),4.5)) - (sqrt(Pi)*sqrt(pow(s,-2))*pow(s,16)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,9)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (4*b*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,14)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,8)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (14*pow(b,2)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,12)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,7)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (28*pow(b,3)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,10)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,6)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (35*pow(b,4)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,8)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,5)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (28*pow(b,5)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,6)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,4)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (14*pow(b,6)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,4)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,3)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (4*pow(b,7)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,2)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (pow(b,8)*sqrt(Pi)*sqrt(pow(s,-2))*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*t*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))) + (14*pow(pow(s,-2),1.5)*pow(s,2)*pow(pow(s,2) + b*t,9)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* sqrt(Pi)*TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,9)*pow(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)),1.5)) + (35*sqrt(pow(s,-2))*pow(s,2)*pow(pow(s,2) + b*t,3)* (sqrt(2)*pow(pow(s,2) + b*t,2)* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + 3*t)) - 3*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)*pow(s,4)*pow(t,4)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,7)) + (14*pow(pow(s,-2),3.5)*pow(s,6)* pow(pow(s,2) + b*t,2)*pow(pow(s,2) + t*(b - y),7)* (15*sqrt(Pi) - (sqrt(2)* sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (pow(s,8) + pow(s,6)*t*(4*b + 5*t - 4*y) + pow(s,2)*pow(t,3)*(4*b + 5*t - 4*y)*pow(b - y,2) + pow(t,4)*pow(b - y,4) + pow(s,4)*pow(t,2)* (6*pow(b,2) + 15*pow(t,2) - 10*t*y + 6*pow(y,2) - 2*b*(-5*t + 6*y))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(s,4)*pow(t,4)) + 15*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (pow(t,9)*pow(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)),3.5)) + (pow(pow(s,-2),4.5)*pow(s,8)*pow(pow(s,2) + t*(b - y),9)* (105*sqrt(Pi) - (sqrt(2)* (105*pow(t,6) + (35*pow(t,4)*pow(pow(s,2) + t*(b - y),2))/ pow(s,2) + (7*pow(t,2)*pow(pow(s,2) + t*(b - y),4))/ pow(s,4) + pow(pow(s,2) + t*(b - y),6)/pow(s,6)) *sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,6)) + 105*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (2.*pow(t,9)*pow(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)),4.5)) + (sqrt(Pi)*sqrt(pow(s,-2))*pow(s,16)*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,9)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (4*b*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,14)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,8)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (14*pow(b,2)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,12)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,7)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (28*pow(b,3)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,10)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,6)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (35*pow(b,4)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,8)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,5)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (28*pow(b,5)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,6)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,4)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (14*pow(b,6)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,4)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,3)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (4*pow(b,7)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,2)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (pow(t,2)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (pow(b,8)*sqrt(Pi)*sqrt(pow(s,-2))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (14*pow(pow(s,-2),1.5)*pow(s,2)*pow(pow(s,2) + b*t,6)* pow(pow(s,2) + t*(b - y),3)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,9)* pow(pow(pow(s,2) + t*(b - y),2)/(pow(s,2)*pow(t,2)), 1.5)) - (35*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,4)* (sqrt(2)*pow(pow(s,2) + t*(b - y),2)* (pow(s,4) + pow(s,2)*t*(2*b + 3*t - 2*y) + pow(t,2)*pow(b - y,2)) - 3*pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)*pow(s,4)*pow(t,4)* sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,7)* (pow(s,2) + t*(b - y)))))/sqrt(pow(s,-2)) - 8*y*(-(pow(E,pow(s,2)/(2.*pow(t,2)) + b/t - y/t)* (-((pow(b,7)*sqrt(Pi/2.))/sqrt(pow(s,-2))) + 48*pow(s,8) - (sqrt(Pi/2.)*pow(s,14))/(sqrt(pow(s,-2))*pow(t,7)) - (7*b*sqrt(Pi/2.)*pow(s,12))/(sqrt(pow(s,-2))*pow(t,6)) - (21*pow(b,2)*sqrt(Pi/2.)*pow(s,10))/ (sqrt(pow(s,-2))*pow(t,5)) - (35*pow(b,3)*sqrt(Pi/2.)*pow(s,8))/ (sqrt(pow(s,-2))*pow(t,4)) - (35*pow(b,4)*sqrt(Pi/2.)*pow(s,6))/ (sqrt(pow(s,-2))*pow(t,3)) - (21*pow(b,5)*sqrt(Pi/2.)*pow(s,4))/ (sqrt(pow(s,-2))*pow(t,2)) - (7*pow(b,6)*sqrt(Pi/2.)*pow(s,2))/(sqrt(pow(s,-2))*t) - (105*sqrt(Pi/2.)*(pow(s,2) + b*t))/ (pow(pow(s,-2),3.5)*t) + (168*pow(s,6)*pow(pow(s,2) + b*t,2))/pow(t,2) - (105*sqrt(Pi/2.)*pow(pow(s,2) + b*t,3))/ (pow(pow(s,-2),2.5)*pow(t,3)) + (70*pow(s,4)*pow(pow(s,2) + b*t,4))/pow(t,4) - (21*sqrt(Pi/2.)*pow(pow(s,2) + b*t,5))/ (pow(pow(s,-2),1.5)*pow(t,5)) + (7*pow(s,2)*pow(pow(s,2) + b*t,6))/pow(t,6) - (7*(-1 + pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*pow(s,2)* pow(pow(s,2) + b*t,6))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (70*pow(s,4)*pow(pow(s,2) + b*t,4)* (pow(s,2)/2. + b*t - (-1 + pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))\ - pow(b,2)/(2.*pow(s,2)))*pow(t,2)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (84*pow(s,6)*pow(pow(s,2) + b*t,2)* (pow(s,4)/4. + b*pow(s,2)*t + (1 + (3*pow(b,2))/(2.*pow(s,2)))*pow(s,2)* pow(t,2) + 2*b*(1 + pow(b,2)/(2.*pow(s,2)))*pow(t,3) - (-2 + 2*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))\ - pow(b,4)/(4.*pow(s,4)) - pow(b,2)/pow(s,2))* pow(t,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (pow(s,2)* (pow(s,12) + pow(b,6)*pow(t,6) + 6*pow(s,10)*t*(b + t) + 6*pow(b,4)*pow(s,2)*pow(t,5)*(b + t) + 3*pow(s,8)*pow(t,2)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 3*pow(b,2)*pow(s,4)*pow(t,4)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 4*pow(s,6)*pow(t,3)* (5*pow(b,3) + 9*pow(b,2)*t + 12*b*pow(t,2) - 12*(-1 + pow(E, pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,3))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (7*pow(b + pow(s,2)/t,7)* (pow(s,2) + b*t)* (15*sqrt(Pi) - (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))* (15*pow(t,4) + (5*pow(t,2)*pow(pow(s,2) + b*t,2))/pow(s,2) + pow(pow(s,2) + b*t,4)/pow(s,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,4)) + 15*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (sqrt(2)*pow(pow(b + pow(s,2)/t,2)/pow(s,2),3.5)*t) + (pow(b,7)*sqrt(Pi/2.)*(b + pow(s,2)/t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/sqrt(pow(b + pow(s,2)/t,2)/pow(s,2)) + (sqrt(Pi/2.)*pow(s,14)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,8)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (7*b*sqrt(Pi/2.)*pow(s,12)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,7)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (21*pow(b,2)*sqrt(Pi/2.)*pow(s,10)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,6)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (35*pow(b,3)*sqrt(Pi/2.)*pow(s,8)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,5)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (35*pow(b,4)*sqrt(Pi/2.)*pow(s,6)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,4)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (21*pow(b,5)*sqrt(Pi/2.)*pow(s,4)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,3)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (7*pow(b,6)*sqrt(Pi/2.)*pow(s,2)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (pow(t,2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) - (21*pow(pow(s,2) + b*t,8)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (sqrt(2)*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,8)*pow(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)),1.5)) - (35*pow(s,6)*pow(pow(s,2) + b*t,2)* ((pow(pow(s,2) + b*t,2)* (pow(s,2) + 2*b*t + (3 + pow(b,2)/pow(s,2))*pow(t,2)))/pow(s,2) - 3*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi/2.)*pow(t,4)* sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (pow(E,pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)))) + (sqrt(2)*pow(E,pow(s,2)/(2.*pow(t,2)) + b/t - y/t)* ((-7*(-1 + pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*sqrt(pow(s,-2))* pow(s,2)*pow(pow(s,2) + b*t,6))/ (sqrt(2)*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (7* (-1 + pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2))))*sqrt(pow(s,-2))* pow(s,2)*pow(pow(s,2) + b*t,6))/ (sqrt(2)*pow(E, pow(pow(s,2) + t*(b - y),2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (35*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,4)* (pow(s,4) + pow(b,2)*pow(t,2) + 2*pow(s,2)*t* (b + t - pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* t)))/ (sqrt(2)*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (21*sqrt(pow(s,-2))* pow(pow(s,3) + b*s*t,2)* (pow(s,8) + pow(b,4)*pow(t,4) + 4*pow(s,6)*t*(b + t) + 4*pow(b,2)*pow(s,2)*pow(t,3)*(b + t) + 2*pow(s,4)*pow(t,2)* (3*pow(b,2) + 4*b*t - 4*(-1 + pow(E, pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,2))))/ (sqrt(2)*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (sqrt(pow(s,-2))*pow(s,2)* (pow(s,12) + pow(b,6)*pow(t,6) + 6*pow(s,10)*t*(b + t) + 6*pow(b,4)*pow(s,2)*pow(t,5)*(b + t) + 3*pow(s,8)*pow(t,2)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 3*pow(b,2)*pow(s,4)*pow(t,4)* (5*pow(b,2) + 8*b*t + 8*pow(t,2)) + 4*pow(s,6)*pow(t,3)* (5*pow(b,3) + 9*pow(b,2)*t + 12*b*pow(t,2) - 12*(-1 + pow(E, pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,3))))/ (sqrt(2)*pow(E, pow(pow(s,2) + b*t,2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) - (35*sqrt(pow(s,-2))*pow(s,2)* pow(pow(s,2) + b*t,4)* (pow(s,4) + pow(t,2)*pow(b - y,2) + 2*pow(s,2)*t* (b + t - pow(E, pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2)))*t - y)))/ (sqrt(2)*pow(E, pow(pow(s,2) + t*(b - y),2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) - (21*sqrt(pow(s,-2))* pow(pow(s,3) + b*s*t,2)* (pow(s,8) + pow(t,4)*pow(b - y,4) + 4*pow(s,6)*t*(b + t - y) + 4*pow(s,2)*pow(t,3)*pow(b - y,2)*(b + t - y) + 2*pow(s,4)*pow(t,2)* (3*pow(b,2) - 4*(-1 + pow(E, pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,2) + b*(4*t - 6*y) - 4*t*y + 3*pow(y,2))))/ (sqrt(2)*pow(E, pow(pow(s,2) + t*(b - y),2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) - (sqrt(pow(s,-2))*pow(s,2)* (pow(s,12) + pow(t,6)*pow(b - y,6) + 6*pow(s,10)*t*(b + t - y) + 6*pow(s,2)*pow(t,5)*pow(b - y,4)*(b + t - y) + 3*pow(s,8)*pow(t,2)* (5*pow(b,2) + 8*pow(t,2) + 2*b*(4*t - 5*y) - 8*t*y + 5*pow(y,2)) + 3*pow(s,4)*pow(t,4)*pow(b - y,2)* (5*pow(b,2) + 8*pow(t,2) + 2*b*(4*t - 5*y) - 8*t*y + 5*pow(y,2)) + 4*pow(s,6)*pow(t,3)* (5*pow(b,3) - 12*(-1 + pow(E, pow(pow(s,2) + b*t - t*y,2)/ (2.*pow(s,2)*pow(t,2))))*pow(t,3) + 3*pow(b,2)*(3*t - 5*y) - 12*pow(t,2)*y + 9*t*pow(y,2) - 5*pow(y,3) + 3*b*(4*pow(t,2) - 6*t*y + 5*pow(y,2)))))/ (sqrt(2)*pow(E, pow(pow(s,2) + t*(b - y),2)/(2.*pow(s,2)*pow(t,2)))* pow(t,6)) + (7*sqrt(pow(s,-2))*pow(pow(s,2) + b*t,7)* (pow(s,8) + b*pow(s,6)*t)* (15*sqrt(Pi) - (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))* (15*pow(s,4)*pow(t,4) + 5*pow(s,2)*pow(t,2)*pow(pow(s,2) + b*t,2) + pow(pow(s,2) + b*t,4)))/ (pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(s,4)*pow(t,4)) + 15*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (2.*pow(s,6)*pow(t,8)* pow(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)),3.5)) + (sqrt(Pi)*sqrt(pow(s,-2))*pow(s,14)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,8)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (7*b*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,12)*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,7)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (21*pow(b,2)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,10)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,6)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (35*pow(b,3)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,8)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,5)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (35*pow(b,4)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,6)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,4)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (21*pow(b,5)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,4)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,3)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (7*pow(b,6)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,2)* (pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*pow(t,2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))) + (pow(b,7)*sqrt(Pi)*sqrt(pow(s,-2))*(pow(s,2) + b*t)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)))/ sqrt(2)))/ (2.*t*sqrt(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2))))\ - (21*pow(pow(s,-2),1.5)*pow(s,2)* pow(pow(s,2) + b*t,8)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (2.*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,8)* pow(pow(pow(s,2) + b*t,2)/(pow(s,2)*pow(t,2)),1.5)) - (35*sqrt(pow(s,-2))*pow(pow(s,3) + b*s*t,2)* (sqrt(2)*pow(pow(s,2) + b*t,2)* (pow(s,4) + pow(b,2)*pow(t,2) + pow(s,2)*t*(2*b + 3*t)) - 3*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)*pow(s,4)* pow(t,4)*sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (2.*pow(E,pow(pow(s,2) + b*t,2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,6)) - (7*pow(pow(s,-2),3.5)*(pow(s,8) + b*pow(s,6)*t)* pow(pow(s,2) + t*(b - y),7)* (15*sqrt(Pi) - (sqrt(2)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* (pow(s,8) + pow(s,6)*t*(4*b + 5*t - 4*y) + pow(s,2)*pow(t,3)*(4*b + 5*t - 4*y)* pow(b - y,2) + pow(t,4)*pow(b - y,4) + pow(s,4)*pow(t,2)* (6*pow(b,2) + 15*pow(t,2) - 10*t*y + 6*pow(y,2) - 2*b*(-5*t + 6*y))))/ (pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(s,4)*pow(t,4)) + 15*sqrt(Pi)*(-1 + TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))))/ (2.*pow(t,8)*pow(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)),3.5)) - (sqrt(Pi)*sqrt(pow(s,-2))*pow(s,14)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,8)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (7*b*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,12)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,7)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (21*pow(b,2)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,10)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,6)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (35*pow(b,3)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,8)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,5)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (35*pow(b,4)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,6)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,4)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (21*pow(b,5)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,4)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,3)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (7*pow(b,6)*sqrt(Pi)*sqrt(pow(s,-2))*pow(s,2)* (pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*pow(t,2)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) - (pow(b,7)*sqrt(Pi)*sqrt(pow(s,-2))*(pow(s,2) + t*(b - y))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2)))/ (2.*t*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))) + (21*pow(pow(s,-2),1.5)*pow(s,2)*pow(pow(s,2) + b*t,5)* pow(pow(s,2) + t*(b - y),3)* (sqrt(2)*sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2))) - pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (2.*pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,8)* pow(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)),1.5)) + (35*sqrt(pow(s,-2))*pow(s,2)*pow(pow(s,2) + b*t,3)* (sqrt(2)*pow(pow(s,2) + t*(b - y),2)* (pow(s,4) + pow(s,2)*t*(2*b + 3*t - 2*y) + pow(t,2)*pow(b - y,2)) - 3*pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*sqrt(Pi)*pow(s,4)* pow(t,4)*sqrt(pow(pow(s,2) + t*(b - y),2)/ (pow(s,2)*pow(t,2)))* TMath::Erf(sqrt(pow(pow(s,2) + b*t - t*y,2)/ (pow(s,2)*pow(t,2)))/sqrt(2))))/ (2.*pow(E,pow(pow(s,2) + t*(b - y),2)/ (2.*pow(s,2)*pow(t,2)))*pow(t,6)* (pow(s,2) + t*(b - y)))))/sqrt(pow(s,-2))) - 90*pow(E,(pow(s,2) + 2*b*t - 2*t*y)/(2.*pow(t,2)))*sqrt(2*Pi)*s* pow(t,2)*pow(y,6)*(-1 + sqrt(pow(s,-2))*s + TMath::Erfc((pow(s,2) + t*(b - y))/(sqrt(2)*s*t))) + pow(E,(pow(s,2) + 2*b*t - 2*t*y)/(2.*pow(t,2)))*sqrt(Pi/2.)*s* pow(y,8)*(-1 + sqrt(pow(s,-2))*s + TMath::Erfc((pow(s,2) + t*(b - y))/(sqrt(2)*s*t))))/ (29160.*sqrt(2*Pi)*s*pow(t,8)); return A*sss;} }