ShwFit Class Reference

#include <ShwFit.h>

List of all members.

Public Member Functions
	ShwFit ()
	~ShwFit ()
void	Reset ()
void	Insert (double ph, double z)
void	Fit3d (int psf=0)
void	Insert3d (double ph, double u, double v, double z)
void	SetDetector (int mydet)
Double_t	GetMaximumX (TF1 *efit, Double_t xmin=0, Double_t xmax=0)
void	BuildPlaneMap ()
int	GetPlaneFromZ (double z)
Public Attributes
double	par_a
double	par_b
double	par_e0
double	par_a_err
double	par_b_err
double	par_e0_err
double	chisq
double	ndf
double	prob
double	shwmax
double	shwmaxplane
double	conv
double	pp_chisq
int	pp_ndf
int	pp_igood
double	pp_p
double	cmp_chisq
int	cmp_ndf
double	peakdiff
double	pred_e_a
double	pred_g_a
double	pred_b
double	pred_e0
double	pred_e_chisq
double	pred_e_ndf
double	pred_g_chisq
double	pred_g_ndf
double	zshift
double	pre_over
double	pre_under
double	post_over
double	post_under
std::vector< double >	v_ph
std::vector< double >	v_u
std::vector< double >	v_v
std::vector< double >	v_z
int	detector
int	debug
Static Public Attributes
TH1F *	lenepl = 0
TF1 *	efit = 0
TH1F *	lenepl3d = 0
std::map< double, int > *	planemap = 0

---

Constructor & Destructor Documentation

ShwFit::ShwFit (   )

Definition at line 91 of file ShwFit.cxx. References BuildPlaneMap(), debug, MsgService::Instance(), MsgService::IsActive(), MSG, and Reset(). { debug=0; if(MsgService::Instance()->IsActive("ShwFit",Msg::kDebug))debug=1; Reset(); //const Int_t LHISTBINS=100; // if(lenepl)delete lenepl; // lenepl = new TH1F("lenepl","longitudinal energy by plane",LHISTBINS,0.0,LHISTBINS); //lenepl->SetDirectory(0); if(!planemap)BuildPlaneMap(); if(debug)MSG("ShwFit",Msg::kError)<<"ParticleFinder ShwFit is in debug mode... there will be a purposeful memory leak\n"; }

ShwFit::~ShwFit (   )

Definition at line 184 of file ShwFit.cxx. References efit, lenepl, and lenepl3d. { if(!debug) { if(lenepl){delete lenepl;lenepl=0;} if(lenepl3d){delete lenepl3d;lenepl3d=0;} if(efit){delete efit;efit=0;} } }

---

Member Function Documentation

void ShwFit::BuildPlaneMap (   )

Definition at line 113 of file ShwFit.cxx. References UgliGeomHandle::GetSteelPlnHandleVector(), infile, MSG, outfile, and planemap. Referenced by ShwFit(). { //try to open plane location file ... for faster running! ifstream infile("planemap.txt"); if(infile.is_open()) { MSG("ShwFit",Msg::kDebug)<<"loading shwfit planemap from file!.... "; planemap=new std::map<double,int>; string line; int cnt=0; while(!infile.eof()) { getline(infile,line); int planei=0; double planez=0; sscanf(line.c_str(),"%lf %d",&planez, &planei); planemap->insert(make_pair(planez,planei)); cnt++; } MSG("ShwFit",Msg::kDebug)<<"loaded "<<cnt<<" planes...."; infile.close(); }else{ VldContext vldc(Detector::kFar , SimFlag::kData , VldTimeStamp()); UgliGeomHandle geo(vldc); planemap=new std::map<double,int>; //iterate over the planes, and insert the upper z and plane number into the map std::vector<UgliSteelPlnHandle> plns = geo.GetSteelPlnHandleVector(); ofstream outfile("planemap.txt"); for(unsigned int i=0;i<plns.size();i++) { planemap->insert(make_pair(plns[i].GetZ0(),i)); if(outfile.is_open()) { outfile << (double)plns[i].GetZ0() << " "<<i<<"\n"; } } outfile.close(); } }

void ShwFit::Fit3d (  int  psf = 0  )

see if we need a z shift! fast version //N do not draw detail version detail version Definition at line 365 of file ShwFit.cxx. References chisq, cmp_chisq, cmp_ndf, conv, debug, UgliGeomHandle::Draw(), efit, GetMaximumX(), GetPlaneFromZ(), lenepl3d, MSG, ndf, par_a, par_a_err, par_b, par_b_err, par_e0, par_e0_err, peakdiff, post_over, post_under, pp_chisq, pp_igood, pp_ndf, pp_p, pre_over, pre_under, pred_b, pred_e0, pred_e_a, pred_e_chisq, pred_e_ndf, pred_g_a, pred_g_chisq, pred_g_ndf, prob, shwfunc3d(), shwmax, shwmaxplane, v_ph, v_u, v_v, and v_z. Referenced by Finder::FinalizeParticles3D(), and PrimaryShowerFinder::FindPrimaryShower(). { double Rsteel = 1.445; // for steel double Rscint=0.0241; //for scint only double R = Rscint+Rsteel; double offset=Rsteel / 2; //offset from front of scint plane to count as 0 radiation lengths... /* //build a vector measuring distance in radiation lengths from first hit VldContext vldc(Detector::kFar , SimFlag::kData , VldTimeStamp()); UgliGeomHandle geo(vldc); */ std::vector<int> planes; for(unsigned int i=0;i<v_z.size();i++) { /* PlexPlaneId p; p = geo.GetPlaneIdFromZ(v_z[i]); planes.push_back(p.GetPlane()); */ planes.push_back(GetPlaneFromZ(v_z[i])); } std::vector<double> rdlen; double rddist=offset; rdlen.push_back(rddist); int sm1=0; int sm2=0; if(v_z[0]<15.5)sm1=1; if(v_z[0]>=15.5)sm2=1; for(unsigned int i=1;i<v_z.size();i++) { double dz = v_z[i]-v_z[i-1]; double du = v_u[i]-v_u[i-1]; double dv = v_v[i]-v_v[i-1]; double tdist = sqrt(dz*dz +du*du +dv*dv); if(planes[i]==planes[i-1]) //in the same scint plane { rddist+=tdist * Rscint; }else{ int dp = planes[i]-planes[i-1]; double dt = sqrt(du*du+dv*dv); rddist += sqrt( R*R*dp*dp + dt*dt); } rdlen.push_back(rddist); if(v_z[i]<15.5)sm1=1; if(v_z[i]>=15.5)sm2=1; } if(sm1 && sm2) { MSG("ShwFit",Msg::kError)<<"Attempting to fit between super modules! exiting\n"; return; } //for(unsigned int i=0;i<v_z.size();i++) //std::cout << "z " << v_z[i] << " plane " <<planes[i] <<" rdlen "<<rdlen[i] <<endl; //decide what size histo to make, and make it int nbins = v_z.size()*2; double rlmax = rdlen[rdlen.size()-1]*2; //don't delete the histo if we want to be viewing it later //this will leak... so just don't run this module in debug mode //except when debugging! if(!debug) if (lenepl3d){delete lenepl3d;lenepl3d=0;} if(debug && psf) lenepl3d = new TH1F("lenepl3dPSF","longitudinal energy by plane",nbins,0.0,rlmax); else lenepl3d = new TH1F("lenepl3d","longitudinal energy by plane",nbins,0.0,rlmax); //hide the histo unless we are debugging the primary shower if(!debug) lenepl3d->SetDirectory(0); //try a smear fill... for(unsigned int i=0;i<rdlen.size();i++) { double ph = v_ph[i]; ph/=R; //cout<<"adding to hist "<<rdlen[i]<<" e "<<ph<<"\n"; if(i==0) { lenepl3d->Fill(rdlen[i],ph);//*0.8); //lenepl3d->Fill(rdlen[i]+lenepl3d->GetBinWidth(1),ph*0.2); //lenepl3d->Fill(rdlen[i]+lenepl3d->GetBinWidth(1)*2,ph*0.1); }else{ //lenepl3d->Fill(rdlen[i]+lenepl3d->GetBinWidth(1),ph*0.2); lenepl3d->Fill(rdlen[i],ph);//*0.60); //lenepl3d->Fill(rdlen[i]-lenepl3d->GetBinWidth(1),ph*0.2); //lenepl3d->Fill(rdlen[i]+lenepl3d->GetBinWidth(1),ph*0.05); //lenepl3d->Fill(rdlen[i]-lenepl3d->GetBinWidth(1),ph*0.05); } } //calculate expected values here.... double chisq_elec=-1; double chisq_gamma=-1; //what should a be for this energy? //a = 1 +b(lnE/Ec -0.5) (or+0.5 for gamma) b=0.5 //Ec is about 36.64 MeV assuming steel has z 26 and scint is about z=7 //calculate total e double etot=0; for(int i=0;i<(int)v_ph.size();i++) etot+=v_ph[i]; //make it MeV with roughly 25 mip/gev etot=etot/23*1000;// *2; if(etot<1e-5)return; pred_b=0.55; pred_e0=etot; double a=1+0.5*(log(etot/36.64)-0.5); TF1* predElec = 0; if(psf)predElec = new TF1("predElecPSF",shwfunc3d,0.1,rlmax,3); else predElec = new TF1("predElec",shwfunc3d,0.1,rlmax,3); predElec->SetParNames("a","b","e0"); predElec->SetParameter(0,a); predElec->SetParameter(1,0.5); predElec->SetParameter(2,etot); predElec->SetLineColor(2); pred_e_a=a; MSG("ShwFit",Msg::kDebug)<<"elec a b e0 "<<a<<" "<<0.5<<" "<<etot<<"\n"; double maxX = GetMaximumX(predElec,0,rlmax); int maxbin = lenepl3d->GetMaximumBin(); double maxH = lenepl3d->GetBinCenter(maxbin); MSG("ShwFit",Msg::kDebug)<<"pred max "<<maxX<<" hist max "<<maxH<<"\n"; /* //do we need to reset the offset of the histogram? if(fabs(maxH-maxX)>0.05) { zshift = maxX-maxH; //lenepl3d->Reset(); if(lenepl3d)delete lenepl3d;lenepl3d=0; lenepl3d = new TH1F("lenepl3d","longitudinal energy by plane",nbins,0.0,rlmax); lenepl3d->SetDirectory(0); //try a smear fill... for(unsigned int i=0;i<rdlen.size();i++) { double ph = v_ph[i]; ph/=R; if(i==0) { lenepl3d->Fill(rdlen[i]+zshift,ph);// *0.8); //lenepl3d->Fill(rdlen[i]+zshift+lenepl3d->GetBinWidth(1),ph*0.2); //lenepl3d->Fill(rdlen[i]+zshift+lenepl3d->GetBinWidth(1)*2,ph*0.1); }else{ // lenepl3d->Fill(rdlen[i]+zshift+lenepl3d->GetBinWidth(1),ph*0.2); lenepl3d->Fill(rdlen[i]+zshift,ph);// *0.60); // lenepl3d->Fill(rdlen[i]+zshift-lenepl3d->GetBinWidth(1),ph*0.2); //lenepl3d->Fill(rdlen[i]+zshift+lenepl3d->GetBinWidth(1),ph*0.05); //lenepl3d->Fill(rdlen[i]+zshift-lenepl3d->GetBinWidth(1),ph*0.05); } } //store zshift in terms of actual distance, not radiation length zshift/=R; //to get number of planes zshift*=0.0354;//average plane width int maxbin = lenepl3d->GetMaximumBin(); double maxH = lenepl3d->GetBinCenter(maxbin); printf("new hist max %f\n",maxH); } */ a=1+0.5*(log(etot/36.64)+0.5); //gamma TF1* predGamma = 0; if(psf) predGamma = new TF1("predGammaPSF",shwfunc3d,0.1,rlmax,3); else predGamma = new TF1("predGamma",shwfunc3d,0.1,rlmax,3); predGamma->SetParNames("a","b","e0"); predGamma->SetParameter(0,a); predGamma->SetParameter(1,0.5); predGamma->SetParameter(2,etot); predGamma->SetLineColor(3); pred_g_a=a; MSG("ShwFit",Msg::kDebug)<<"gamma a b e0 "<<a<<" "<<0.5<<" "<<etot<<"\n"; //calculate chi sq's .... chisq_elec=0; chisq_gamma=0; int dc_e=0; int dc_g=0; for(int i=1;i<lenepl3d->GetNbinsX()+1;i++) { double val = lenepl3d->GetBinContent(i); double x = lenepl3d->GetBinCenter(i); MSG("ShwFit",Msg::kDebug)<<"v "<<val<<" x "<<x<<" eval "<<predElec->Eval(x)<<" "<<predGamma->Eval(x)<<"\n"; // double predval = predElec->Integral( lenepl3d->GetBinLowEdge(i),lenepl3d->GetBinLowEdge(i+1),(double *)0,1e-3); double bincenter = lenepl3d->GetBinCenter(i); double predval = predElec->Eval(bincenter); // if(predval>1e-4 && val>0) if(predval)chisq_elec+=(val-predval)*(val-predval)/(predval); //else dc_e++; //predval = predGamma->Integral( lenepl3d->GetBinLowEdge(i),lenepl3d->GetBinLowEdge(i+1),(double *)0,1e-3); predval = predGamma->Eval(bincenter); //if(predval>1e-4 && val>0) if(predval) chisq_gamma+=(val-predval)*(val-predval)/(predval); //else dc_g++; } pred_e_chisq=chisq_elec; pred_e_ndf=lenepl3d->GetNbinsX()-dc_e; pred_g_chisq=chisq_gamma; pred_g_ndf=lenepl3d->GetNbinsX()-dc_g; chisq_elec/=lenepl3d->GetNbinsX()-dc_e; chisq_gamma/=lenepl3d->GetNbinsX()-dc_g; post_over=0; post_under=0; pre_over=0; pre_under=0; int m = lenepl3d->GetMaximumBin(); //do the fit int hmin=0; int npare=3; double pulseheight = lenepl3d->Integral(); if(!debug) if(efit){delete efit;efit=0;} efit = new TF1("efit",shwfunc3d,hmin+0.001,rlmax,npare); efit->SetParNames("a","b","e0","zshift"); //efit->SetParLimits(0,hmin+0.001,rlmax); //efit->SetParLimits(1,0.001,20000); //efit->SetParLimits(2,0+0.001,1000000); efit->SetParameters(3,0.5,300); efit->SetParLimits(0,1.5,5); efit->SetParLimits(1,0.01,1.5); efit->SetParLimits(2,0,1000000); //efit->SetParameters(3,-5,5); //efit->SetParameters(3.,0.5,pulseheight *60*25,0); //gues e as sum of ph in projection.... TCanvas * c=0; if(debug){ MSG("ShwFit",Msg::kDebug)<<(psf!=0)<<" "<<(c!=0)<<"\n"; if(psf){ c = new TCanvas("psfcan","psfcan"); c->cd(); lenepl3d->Draw(); } MSG("ShwFit",Msg::kDebug)<<(psf!=0)<<" "<<(c!=0)<<"\n"; } if(!debug){ lenepl3d->Fit(efit,"NWLLQ"); }else{ if(psf) lenepl3d->Fit(efit,"WLLV+"); else lenepl3d->Fit(efit,"NWLLV"); } TH1F * predhist= 0; if(psf){ predhist = new TH1F("predhistPSF","longitudinal energy by plane",nbins,0.0,rlmax); // predhist->SetDirectory(0); predhist->SetLineColor(2); }else{ predhist = new TH1F("predhist","longitudinal energy by plane",nbins,0.0,rlmax); predhist->SetDirectory(0); predhist->SetLineColor(2); } for(int i=1;i<lenepl3d->GetNbinsX()+1;i++) { double x = lenepl3d->GetBinCenter(i); // double pred = predElec->Integral(lenepl3d->GetBinLowEdge(i),lenepl3d->GetBinLowEdge(i+1),(double *)0,1e-3); double pred = predElec->Eval(lenepl3d->GetBinCenter(i)); double val = lenepl3d->GetBinContent(i); if(predhist)predhist->Fill(x,pred); if(i>m) { if(pred<val) post_over+=val-pred; else post_under+=pred-val; }else{ if(pred<val) pre_over+=val-pred; else pre_under+=pred-val; } } if(debug&&psf&&predhist)predhist->Draw("same"); //draw the predicted fit here... if(debug){ MSG("ShwFit",Msg::kDebug)<<(psf!=0)<<" "<<(c!=0)<<"\n"; if(psf && c){ c->cd(); predElec->Draw("same"); predGamma->Draw("same"); MSG("ShwFit",Msg::kDebug)<<"ASDFASDFASDF\n"; } } pp_chisq=0; pp_ndf=0; pp_igood=0; pp_p=0; if(psf){ //require nonempty histograms to prevent endless loop! if(predhist->Integral()>0.01 && lenepl3d->Integral()>0.01) { MSG("ShwFit",Msg::kDebug)<<"int pred "<<(predhist->Integral()>0)<<" shw "<< (lenepl3d->Integral()>0)<<"\n"; //pp_p=lenepl3d->Chi2TestX(predhist,pp_chisq,pp_ndf,pp_igood,"WW",(double*)0); } } //lets directly compare the predicted histogram with the original one cmp_chisq=0; cmp_ndf = 0; peakdiff=0; if(predhist) { cmp_ndf=predhist->GetNbinsX(); for(int k=1;k<cmp_ndf+1;k++) { double oc = lenepl3d->GetBinContent(k); double pc = predhist->GetBinContent(k); if(pc)cmp_chisq+=(oc-pc)*(oc-pc)/pc; else if(oc)cmp_ndf--; } if(cmp_chisq>0)cmp_chisq=sqrt(cmp_chisq); int maxbin = predhist->GetMaximumBin(); if(maxbin>1) { peakdiff+= lenepl3d->GetBinContent(maxbin-1)-predhist->GetBinContent(maxbin-1); peakdiff+= lenepl3d->GetBinContent(maxbin)-predhist->GetBinContent(maxbin); peakdiff+= lenepl3d->GetBinContent(maxbin+1)-predhist->GetBinContent(maxbin+1); }else if(maxbin==1) { peakdiff+= lenepl3d->GetBinContent(maxbin)-predhist->GetBinContent(maxbin); peakdiff+= lenepl3d->GetBinContent(maxbin+2)-predhist->GetBinContent(maxbin+2); peakdiff+= lenepl3d->GetBinContent(maxbin+1)-predhist->GetBinContent(maxbin+1); } // printf("-- %f %d %f\n",cmp_chisq,cmp_ndf, peakdiff); } MSG("ShwFit",Msg::kDebug)<<" STATUS: "<<gMinuit->fCstatu <<" "<<efit->GetParameter(0) <<" "<<efit->GetNDF()<<" "<<pulseheight<<endl; string fitstatus = (string)(gMinuit->fCstatu); MSG("ShwFit",Msg::kDebug)<<" STATUS: "<<fitstatus <<", "<<efit->GetParameter(0) <<", "<<efit->GetNDF()<<" "<<pulseheight<<endl; if(fitstatus=="CONVERGED "&&efit->GetParameter(0)<29.9&& efit->GetNDF()>0&&pulseheight>0){ MSG("ShwFit",Msg::kDebug)<<" filling vars"<<endl; par_a=efit->GetParameter(0); par_b=efit->GetParameter(1); par_e0=efit->GetParameter(2); par_a_err=efit->GetParError(0); par_b_err=efit->GetParError(1); par_e0_err=efit->GetParError(2); chisq=efit->GetChisquare(); ndf = efit->GetNDF(); //ndf should be equal to the number of bins... because we are giving 0 bins weight 1 ndf = nbins; if(debug && psf){ c->cd(); TLatex l; l.SetTextAlign(12); l.SetNDC(1); char txt[100]; sprintf(txt,"\\chi^2=%f",efit->GetChisquare()); l.DrawLatex(0.5,0.8,txt); sprintf(txt,"NDF=%f",ndf); l.DrawLatex(0.5,0.75,txt); sprintf(txt,"\\chi^2/NDF=%f",efit->GetChisquare()/ndf); l.DrawLatex(0.5,0.7,txt); sprintf(txt,"par a %f \\pm %f",par_a,par_a_err); l.DrawLatex(0.5,0.65,txt); sprintf(txt,"par b %f \\pm %f",par_b,par_b_err); l.DrawLatex(0.5,0.6,txt); sprintf(txt,"e0 %f \\pm %f",par_e0,par_e0_err); l.DrawLatex(0.5,0.55,txt); sprintf(txt,"elec \\chi^2/NDF = %f ",chisq_elec); l.DrawLatex(0.5,0.50,txt); sprintf(txt,"gamma \\chi^2/NDF = %f ",chisq_gamma); l.DrawLatex(0.5,0.45,txt); sprintf(txt,"pre over %f under %f",pre_over,pre_under); l.DrawLatex(0.5,0.40,txt); sprintf(txt,"post over %f under %f",post_over,post_under); l.DrawLatex(0.5,0.35,txt); sprintf(txt,"%f --- %f %d %f", pp_p ,pp_chisq,pp_ndf,pp_chisq/pp_ndf); l.DrawLatex(0.5,0.3,txt); if(cmp_ndf) { sprintf(txt,"cmp c2 %f ndf %d c2ndf %f ", cmp_chisq,cmp_ndf,cmp_chisq/(double)cmp_ndf); l.DrawLatex(0.5,0.25,txt); } sprintf(txt,"pd %f ",peakdiff); l.DrawLatex(0.5,0.2,txt); } shwmax=1.*GetMaximumX(efit); prob=efit->GetProb(); shwmaxplane=(int)(lenepl3d-> GetBinContent(lenepl3d->FindBin(shwmax))); conv=1; } if(!psf || !debug) { delete predElec;predElec=0; delete predGamma;predGamma=0; delete predhist;predhist=0; } if(!psf) { delete efit;efit=0; } }

Double_t ShwFit::GetMaximumX (  TF1 *  efit, Double_t  xmin = 0, Double_t  xmax = 0 )

Definition at line 263 of file ShwFit.cxx. Referenced by Fit3d(). { Double_t fXmin = 0.001; Double_t fXmax = 100; Double_t fNpx = 100; Double_t x,y; if (xmin >= xmax) {xmin = fXmin; xmax = fXmax;} Double_t dx = (xmax-xmin)/fNpx; Double_t xxmax = xmin; Double_t yymax = efit->Eval(xmin+dx); for (Int_t i=0;i<fNpx;i++) { x = xmin + (i+0.5)*dx; y = efit->Eval(x); if (y > yymax) {xxmax = x; yymax = y;} } if (dx < 1.e-9*(fXmax-fXmin)) return TMath::Min(xmax,xxmax); else return GetMaximumX(efit, TMath::Max(xmin,xxmax-dx), TMath::Min(xmax,xxmax+dx)); }

int ShwFit::GetPlaneFromZ (  double  z  )

Definition at line 171 of file ShwFit.cxx. References planemap. Referenced by Fit3d(). { int plane=0; std::map<double,int>::iterator i = planemap->lower_bound(z); if(i!=planemap->end()) plane=i->second; return plane; }

void ShwFit::Insert (  double  ph, double  z )

Definition at line 258 of file ShwFit.cxx. References lenepl. { lenepl->Fill(z+0.5,ph); }

void ShwFit::Insert3d (  double  ph, double  u, double  v, double  z )

Definition at line 355 of file ShwFit.cxx. References v_ph, v_u, v_v, and v_z. Referenced by Finder::FinalizeParticles3D(), and PrimaryShowerFinder::FindPrimaryShower(). { v_ph.push_back(ph); v_u.push_back(u); v_v.push_back(v); v_z.push_back(z); }

void ShwFit::Reset (   )

Definition at line 198 of file ShwFit.cxx. References chisq, cmp_chisq, cmp_ndf, conv, detector, efit, lenepl, lenepl3d, ndf, par_a, par_a_err, par_b, par_b_err, par_e0, par_e0_err, peakdiff, post_over, post_under, pp_chisq, pp_igood, pp_ndf, pp_p, pre_over, pre_under, pred_b, pred_e0, pred_e_a, pred_e_chisq, pred_e_ndf, pred_g_a, pred_g_chisq, pred_g_ndf, prob, shwmax, shwmaxplane, v_ph, v_u, v_v, v_z, and zshift. Referenced by PrimaryShowerFinder::FindPrimaryShower(), and ShwFit(). { if(!debug) { if(lenepl){delete lenepl;lenepl=0;} if(lenepl3d){delete lenepl3d;lenepl3d=0;} if(efit){delete efit;efit=0;} } cmp_chisq=0; cmp_ndf = 0; peakdiff=0; pp_chisq=0; pp_ndf=0; pp_igood=0; pp_p=0; lenepl=0; par_a=0; par_b=0; par_e0=0; par_a_err=0; par_b_err=0; par_e0_err=0; chisq=0; ndf=0; prob=0; shwmax=0; shwmaxplane=0; conv=0; pred_e_a=0; pred_g_a=0; pred_b=0; pred_e0=0; pred_e_chisq=0; pred_e_ndf=0; pred_g_chisq=0; pred_g_ndf=0; detector=0; v_ph.clear(); v_u.clear(); v_v.clear(); v_z.clear(); zshift=0; pre_over=0; pre_under=0; post_over=0; post_under=0; }

void ShwFit::SetDetector (  int  mydet  )

Definition at line 350 of file ShwFit.cxx. References detector. { detector=mydet; }

---

Member Data Documentation

double ShwFit::chisq

Definition at line 27 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::cmp_chisq

Definition at line 41 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

int ShwFit::cmp_ndf

Definition at line 42 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::conv

Definition at line 34 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

int ShwFit::debug

Definition at line 96 of file ShwFit.h. Referenced by Fit3d(), and ShwFit().

int ShwFit::detector

Definition at line 89 of file ShwFit.h. Referenced by Reset(), and SetDetector().

TF1 * ShwFit::efit = 0 [static]

Definition at line 39 of file ShwFit.cxx. Referenced by Fit3d(), Reset(), and ~ShwFit().

TH1F * ShwFit::lenepl = 0 [static]

Definition at line 38 of file ShwFit.cxx. Referenced by Insert(), Reset(), and ~ShwFit().

TH1F * ShwFit::lenepl3d = 0 [static]

Definition at line 40 of file ShwFit.cxx. Referenced by Fit3d(), Reset(), and ~ShwFit().

double ShwFit::ndf

Definition at line 28 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_a

Definition at line 17 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_a_err

Definition at line 22 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_b

Definition at line 18 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_b_err

Definition at line 23 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_e0

Definition at line 19 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::par_e0_err

Definition at line 24 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::peakdiff

Definition at line 43 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

std::map< double, int > * ShwFit::planemap = 0 [static]

Definition at line 42 of file ShwFit.cxx. Referenced by BuildPlaneMap(), and GetPlaneFromZ().

double ShwFit::post_over

Definition at line 72 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::post_under

Definition at line 73 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pp_chisq

Definition at line 36 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

int ShwFit::pp_igood

Definition at line 38 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

int ShwFit::pp_ndf

Definition at line 37 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pp_p

Definition at line 39 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pre_over

Definition at line 70 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pre_under

Definition at line 71 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_b

Definition at line 61 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_e0

Definition at line 62 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_e_a

Definition at line 59 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_e_chisq

Definition at line 63 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_e_ndf

Definition at line 64 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_g_a

Definition at line 60 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_g_chisq

Definition at line 65 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::pred_g_ndf

Definition at line 66 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::prob

Definition at line 30 of file ShwFit.h. Referenced by Finder::FinalizeParticles3D(), PrimaryShowerFinder::FindPrimaryShower(), Fit3d(), and Reset().

double ShwFit::shwmax

Definition at line 32 of file ShwFit.h. Referenced by Fit3d(), and Reset().

double ShwFit::shwmaxplane

Definition at line 33 of file ShwFit.h. Referenced by Fit3d(), and Reset().

std::vector<double> ShwFit::v_ph

Definition at line 84 of file ShwFit.h. Referenced by Fit3d(), Insert3d(), and Reset().

std::vector<double> ShwFit::v_u

Definition at line 85 of file ShwFit.h. Referenced by Fit3d(), Insert3d(), and Reset().

std::vector<double> ShwFit::v_v

Definition at line 86 of file ShwFit.h. Referenced by Fit3d(), Insert3d(), and Reset().

std::vector<double> ShwFit::v_z

Definition at line 87 of file ShwFit.h. Referenced by Fit3d(), Insert3d(), and Reset().

double ShwFit::zshift

Definition at line 68 of file ShwFit.h. Referenced by PrimaryShowerFinder::FindPrimaryShower(), and Reset().

---

The documentation for this class was generated from the following files:

• ShwFit.h
• ShwFit.cxx

---