AngClusterFitAna Class Reference

#include <AngClusterFitAna.h>

Inheritance diagram for AngClusterFitAna:

List of all members.

Public Types
typedef std::deque< Int_t >	DeqInt_t
typedef DeqInt_t::iterator	IterDeqInt_t
typedef DeqInt_t::size_type	SizeDeqInt_t
typedef std::deque< Float_t >	DeqFloat_t
typedef DeqFloat_t::iterator	IterDeqFloat_t
typedef std::vector< Int_t >	VecInt_t
typedef VecInt_t::iterator	IterVecInt_t
typedef std::vector< Float_t >	VecFloat_t
typedef VecFloat_t::iterator	IterVecFloat_t
typedef std::deque< std::deque< Int_t > >	DeqDeqInt_t
typedef DeqDeqInt_t::iterator	IterDeqDeqInt_t
typedef DeqDeqInt_t::size_type	SizeDeqDeqInt_t
typedef std::deque< TMatrixD >	DeqTMatrixD
typedef DeqTMatrixD::iterator	IterDeqTMatrixD
typedef std::deque< TPolyMarker3D * >	DeqTPoly
typedef DeqTPoly::iterator	IterDeqTPoly
Public Member Functions
	AngClusterFitAna (AngClusterFit &acf)
	~AngClusterFitAna ()
void	Set3DHit (DeqFloat_t &x, DeqFloat_t &y, DeqFloat_t &z, DeqFloat_t &e)
void	SetAngCluster (Int_t &primShow, DeqDeqInt_t &clusterMap, TVector3 &primDir)
void	Analyze (int evtn, RecRecordImp< RecCandHeader > *srobj)
void	FitShower (Float_t totalHitEnergy)
void	TransVarCalc (TH1F *transHisto)
void	Draw (TPad *pad)
void	Reset ()
Public Attributes
TH1F *	fLongHitHisto
TH1F *	fLongHitHistoRad
TH1F *	fTransHitHisto
TH1F *	fTransHitHistoRad
TH3F *	fScatterAxis
TF1 *	fLongFit
TF1 *	fTransFit
Private Member Functions
Double_t	GetMaximum (TF1 *efit, Double_t xmin=0, Double_t xmax=0)
Private Attributes
AngClusterFit &	fAngClusterFit
DeqFloat_t	fX
DeqFloat_t	fY
DeqFloat_t	fZ
DeqFloat_t	fE
Int_t	fPrimShow
DeqDeqInt_t	fClusterMap
TVector3	fPrimDir

---

Member Typedef Documentation

typedef std::deque< std::deque <Int_t> > AngClusterFitAna::DeqDeqInt_t

Definition at line 53 of file AngClusterFitAna.h. Referenced by Analyze(), and Draw().

typedef std::deque<Float_t> AngClusterFitAna::DeqFloat_t

Definition at line 44 of file AngClusterFitAna.h.

typedef std::deque<Int_t> AngClusterFitAna::DeqInt_t

Definition at line 40 of file AngClusterFitAna.h.

typedef std::deque<TMatrixD> AngClusterFitAna::DeqTMatrixD

Definition at line 57 of file AngClusterFitAna.h.

typedef std::deque<TPolyMarker3D *> AngClusterFitAna::DeqTPoly

Definition at line 60 of file AngClusterFitAna.h. Referenced by Draw().

typedef DeqDeqInt_t::iterator AngClusterFitAna::IterDeqDeqInt_t

Definition at line 54 of file AngClusterFitAna.h. Referenced by Analyze(), and Draw().

typedef DeqFloat_t::iterator AngClusterFitAna::IterDeqFloat_t

Definition at line 45 of file AngClusterFitAna.h.

typedef DeqInt_t::iterator AngClusterFitAna::IterDeqInt_t

Definition at line 41 of file AngClusterFitAna.h. Referenced by Analyze(), and Draw().

typedef DeqTMatrixD::iterator AngClusterFitAna::IterDeqTMatrixD

Definition at line 58 of file AngClusterFitAna.h.

typedef DeqTPoly::iterator AngClusterFitAna::IterDeqTPoly

Definition at line 61 of file AngClusterFitAna.h. Referenced by Draw().

typedef VecFloat_t::iterator AngClusterFitAna::IterVecFloat_t

Definition at line 51 of file AngClusterFitAna.h.

typedef VecInt_t::iterator AngClusterFitAna::IterVecInt_t

Definition at line 48 of file AngClusterFitAna.h.

typedef DeqDeqInt_t::size_type AngClusterFitAna::SizeDeqDeqInt_t

Definition at line 55 of file AngClusterFitAna.h.

typedef DeqInt_t::size_type AngClusterFitAna::SizeDeqInt_t

Definition at line 42 of file AngClusterFitAna.h.

typedef std::vector<Float_t> AngClusterFitAna::VecFloat_t

Definition at line 50 of file AngClusterFitAna.h.

typedef std::vector<Int_t> AngClusterFitAna::VecInt_t

Definition at line 47 of file AngClusterFitAna.h.

---

Constructor & Destructor Documentation

AngClusterFitAna::AngClusterFitAna (  AngClusterFit &  acf  )

Definition at line 51 of file AngClusterFitAna.cxx. : fLongHitHisto(0), fLongHitHistoRad(0), fTransHitHisto(0), fTransHitHistoRad(0), fScatterAxis(0), fLongFit(0), fTransFit(0), fAngClusterFit(acf) { }

AngClusterFitAna::~AngClusterFitAna (   )

Definition at line 64 of file AngClusterFitAna.cxx. References fLongFit, fLongHitHisto, fLongHitHistoRad, fScatterAxis, fTransFit, fTransHitHisto, and fTransHitHistoRad. { //Prevent memory leaks if(fLongHitHisto!=0){ fLongHitHisto->Delete(); fLongHitHisto=0; } if(fLongHitHistoRad !=0){ fLongHitHistoRad->Delete(); } if(fTransHitHisto !=0){ fTransHitHisto->Delete(); fTransHitHisto=0; } if(fTransHitHistoRad !=0){ fTransHitHistoRad->Delete(); fTransHitHistoRad=0; } if(fScatterAxis !=0){ fScatterAxis->Delete(); fScatterAxis=0; } if(fLongFit !=0){ fLongFit->Delete(); fLongFit=0; } if(fTransFit !=0){ fTransFit->Delete(); fTransFit=0; } }

---

Member Function Documentation

void AngClusterFitAna::Analyze (  int  evtn, RecRecordImp< RecCandHeader > *  srobj )  [virtual]

Implements NueAnaBase. Definition at line 167 of file AngClusterFitAna.cxx. References DeqDeqInt_t, fAngClusterFit, fClusterMap, fE, FitShower(), fLongHitHisto, fLongHitHistoRad, fPrimDir, fPrimShow, fTransHitHisto, fTransHitHistoRad, fX, fY, fZ, RecRecordImp< T >::GetHeader(), RecPhysicsHeader::GetSnarl(), ANtpDefaultValue::IsDefault(), IterDeqDeqInt_t, IterDeqInt_t, kRadLengthToPlane, MSG, AngClusterFit::Reset(), and TransVarCalc(). Referenced by NueRecordAna::Analyze(), and NueDisplayModule::Analyze(). { if (ANtpDefVal::IsDefault(fPrimShow)){ MSG("AngClusterFitAna",Msg::kInfo)<< "No primary cluster " << " found in event." << endl << endl; fAngClusterFit.Reset(); return; } if(fZ.size() < 3 || fZ.size() > 1000) return; if(fPrimDir.X()==0.0 && fPrimDir.Y()==0.0 && fPrimDir.Z()==0) return; MSG("AngClusterFitAna",Msg::kDebug)<<"In AngClusterFitAna::Analyze"<<endl; MSG("AngClusterFitAna",Msg::kDebug)<<"On Snarl " <<srobj->GetHeader().GetSnarl() <<" event "<<evtn<<endl; if(srobj==0){ return; } if(((dynamic_cast<NtpStRecord *>(srobj))==0)&& ((dynamic_cast<NtpSRRecord *>(srobj))==0)){ return; } Float_t cosZAngle=fPrimDir.CosTheta(); //Angle between //Primary Shower dir //and the Z axis. const Int_t kLHistBin=40; const Int_t kTHistBin=40; const Float_t kTHistXLim=static_cast<Float_t>(kTHistBin)*kStripWidth; char lhh[100]; char thh[100]; char lhhr[100]; char thhr[100]; sprintf(lhh ,"LongHitHisto_%d_%d" ,srobj->GetHeader().GetSnarl(),evtn); sprintf(lhhr,"LongHitHistoRad_%d_%d" ,srobj->GetHeader().GetSnarl(),evtn); sprintf(thh ,"TransHitHisto_%d_%d" ,srobj->GetHeader().GetSnarl(),evtn); sprintf(thhr,"TransHitHistoRad_%d_%d",srobj->GetHeader().GetSnarl(),evtn); fLongHitHisto = new TH1F(lhh,"longitudinal energy hit profile" ,kLHistBin,0.0,kLHistBin); fLongHitHistoRad = new TH1F(lhhr,"longitudinal energy hit profile(X0)" ,kLHistBin ,0.0 ,kLHistBin*kRadLengthToPlane/cosZAngle); fTransHitHisto = new TH1F(thh,"trasverse energy hit profile", kTHistBin,-kTHistXLim,kTHistXLim); fTransHitHistoRad = new TH1F(thhr,"trasverse energy hit profile(Rm)", kTHistBin,0,kTHistBin); TVector3 unitDir; //Unit vector parallel to PrimShow dir unitDir=fPrimDir.Unit(); TRotation r,a; r.SetToIdentity(); a.SetToIdentity(); if (TMath::Abs(fPrimDir.Z()) > 0.0){ //Do not rotate Z axis if shower dir //goes along it. r.SetZAxis(unitDir); //Rotates the Z axis to be parallel to the //primDir vector a=r.Inverse(); //Inverting the rotation means when one //applies the transformation to the hit //vector, one obtains the primDir coordinates } TVector3 hitTemp; DeqDeqInt_t cluTemp; cluTemp.push_back(fClusterMap.at(fPrimShow)); //Get PrimShow cluster Float_t rPrime; //Hit Radius around the shower direction. Float_t totalHitEnergy=0.0; //Loop over primary shower hits and fill Long and Trans histograms for (IterDeqDeqInt_t cluIter = cluTemp.begin() ;cluIter != cluTemp.end() ; ++cluIter ){ for (IterDeqInt_t hitIter = cluIter->begin() ;hitIter < cluIter->end() ; ++hitIter){ hitTemp.SetX(fX.at(*hitIter)); hitTemp.SetY(fY.at(*hitIter)); hitTemp.SetZ(fZ.at(*hitIter)); totalHitEnergy+=(fE.at(*hitIter))*kMeuToGeV; Float_t nPlane=hitTemp.Z()/(kDetectorPitch); Float_t nPlaneDir=nPlane/(cosZAngle); Float_t radDist=nPlaneDir*kRadLengthToPlane; //Longitudinal hit //distance from //vertex in X0 units. //Transform hit coordinate system to Zprime hitTemp.Transform(a); rPrime=TMath::Sqrt(hitTemp.X()*hitTemp.X() +hitTemp.Y()*hitTemp.Y()); Float_t rPrimeRad=rPrime/(kMolRadDet); //Transverse hit distance //from Zprime in Moliere //radius units fLongHitHisto ->Fill(nPlane,fE.at(*hitIter)*kMeuToGeV); fLongHitHistoRad ->Fill(radDist,fE.at(*hitIter)*kMeuToGeV); fTransHitHisto ->Fill(hitTemp.X(),fE.at(*hitIter)*kMeuToGeV); fTransHitHistoRad->Fill(rPrimeRad,fE.at(*hitIter)*kMeuToGeV); }//for (hitIter) }//for(cluIter) //Do the Longitudinal and Transverse fits FitShower(totalHitEnergy); //fill transverse variables TransVarCalc(fTransHitHisto); MSG("AngClusterFitAna",Msg::kDebug)<<"Leaving AngClusterFitAna::Analyze" <<endl; }

void AngClusterFitAna::Draw (  TPad *  pad  )

Definition at line 527 of file AngClusterFitAna.cxx. References DeqDeqInt_t, DeqTPoly, AngClusterFit::Draw(), fAngClusterFit, fClusterMap, fLongFit, fLongHitHistoRad, fPrimDir, fPrimShow, fScatterAxis, fTransFit, fTransHitHisto, fTransHitHistoRad, fX, fY, fZ, ANtpDefaultValue::IsDefault(), IterDeqDeqInt_t, IterDeqInt_t, IterDeqTPoly, MSG, and AngClusterFit::Reset(). Referenced by NueDisplayModule::Analyze(). { if (ANtpDefVal::IsDefault(fPrimShow)){ MSG("AngClusterFitAna",Msg::kInfo)<< "No primary cluster " << " found in event." << endl << endl; fAngClusterFit.Reset(); return; } if(fZ.size() < 3 || fZ.size() > 1000) return; if(fPrimDir.X()==0.0 && fPrimDir.Y()==0.0 && fPrimDir.Z()==0) return; DeqTPoly scatterCluster; TPolyMarker3D *scatterTemp; DeqDeqInt_t scatTemp; Int_t nClu=0; Int_t nHit=0; //Loop over all cluster hits and fill polymarkers for (IterDeqDeqInt_t cluIter = fClusterMap.begin() ;cluIter != fClusterMap.end() ; ++cluIter ){ scatTemp.push_back(fClusterMap.at(nClu)); scatterTemp = new TPolyMarker3D(scatTemp.size(),1); for (IterDeqInt_t hitIter = cluIter->begin() ;hitIter < cluIter->end() ; ++hitIter){ scatterTemp->SetPoint(nHit,fZ.at(*hitIter) ,fX.at(*hitIter) ,fY.at(*hitIter)); nHit++; }//for (hitIter) scatterCluster.push_back(scatterTemp); nHit=0; nClu++; }//for(cluIter) // Find min and max values; Float_t minX=static_cast <Float_t> (*min_element(fX.begin(),fX.end())); Float_t maxX=static_cast <Float_t> (*max_element(fX.begin(),fX.end())); Float_t minY=static_cast <Float_t> (*min_element(fY.begin(),fY.end())); Float_t maxY=static_cast <Float_t> (*max_element(fY.begin(),fY.end())); Float_t minZ=static_cast <Float_t> (*min_element(fZ.begin(),fZ.end())); Float_t maxZ=static_cast <Float_t> (*max_element(fZ.begin(),fZ.end())); if(fScatterAxis) delete fScatterAxis; fScatterAxis = new TH3F("cluHisto","Scatter plot of all hits" , 10,minZ-0.1,maxZ+0.2 , 10,minX-0.2,maxX+0.2 , 10,minY-0.2,maxY+0.2); fScatterAxis->SetXTitle("Z(m)"); fScatterAxis->SetYTitle("X(m)"); fScatterAxis->SetZTitle("Y(m)"); TPolyMarker3D *evtVertex; evtVertex = new TPolyMarker3D(1,8); evtVertex->SetMarkerColor(1); evtVertex->SetPoint(0,0.,0.,0.); TPolyLine3D *lineCluster; lineCluster= new TPolyLine3D(2); lineCluster->SetPoint(0,0.,0.,0.); lineCluster->SetPoint(1,fPrimDir.Z(),fPrimDir.X(),fPrimDir.Y()); pad->cd(1); fScatterAxis->Draw(); Int_t polyColor=1; for (IterDeqTPoly polyIter = scatterCluster.begin() ;polyIter < scatterCluster.end() ; ++polyIter){ (*polyIter)->SetMarkerStyle(7); (*polyIter)->SetMarkerColor(polyColor); (*polyIter)->Draw("sames"); polyColor++; }//for(polyIter) evtVertex->Draw("sames"); lineCluster->SetLineWidth(2); lineCluster->Draw("sames"); pad->cd(3); fLongHitHistoRad->Draw(); fLongFit->Draw("sames"); pad->cd(2); fTransHitHistoRad->Draw(); fTransFit->Draw("sames"); pad->cd(4); fTransHitHisto->Draw(); pad->cd(); pad->Modified(); //Clearing ROOT objects evtVertex->Delete(); scatterTemp->Delete(); lineCluster->Delete(); }

void AngClusterFitAna::FitShower (  Float_t  totalHitEnergy  )

Definition at line 313 of file AngClusterFitAna.cxx. References AngClusterFit::fACluFitE0EnergyRatio, AngClusterFit::fACluFitLongChiSq, AngClusterFit::fACluFitLongConv, AngClusterFit::fACluFitLongNDF, AngClusterFit::fACluFitParA, AngClusterFit::fACluFitParB, AngClusterFit::fACluFitParC12, AngClusterFit::fACluFitParL1, AngClusterFit::fACluFitParL2, AngClusterFit::fACluFitParLongE0, AngClusterFit::fACluFitParTransE0, AngClusterFit::fACluFitShwMax, AngClusterFit::fACluFitTransChiSq, AngClusterFit::fACluFitTransConv, AngClusterFit::fACluFitTransNDF, fAngClusterFit, fLongFit, fLongHitHistoRad, fTransFit, fTransHitHistoRad, GetMaximum(), MSG, shwLongFunc(), and shwTransFunc(). Referenced by Analyze(). { //Configure fit functions Int_t lMin=0; Int_t lMax=100; Int_t tMin=0; Int_t tMax=100; Int_t nParLong=3; Int_t nParTrans=4; char lfn[100]; char tfn[100]; sprintf(lfn,"lfit_%d_%d",lMin,lMax); sprintf(tfn,"tfit_%d_%d",tMin,tMax); fLongFit = new TF1(lfn,shwLongFunc,lMin+0.001,lMax,nParLong); fTransFit = new TF1(tfn,shwTransFunc,tMin+0.001,tMax,nParTrans); fLongFit->SetParNames("a","b","E0"); fLongFit->SetParLimits(0,lMin+0.001,lMax); fLongFit->SetParLimits(1,0.001,20000); fLongFit->SetParLimits(2,0.001,20000); fTransFit->SetParNames("l1,l2,c12,E0"); fTransFit->SetParLimits(0,tMin+0.001,tMax); fTransFit->SetParLimits(1,tMin+0.001,tMax); fTransFit->SetParLimits(2,0.001,50); fTransFit->SetParLimits(3,0.001,20000); //Do the Longitudinal fit fLongFit->SetParameters(3.,0.5,totalHitEnergy); fLongHitHistoRad->Fit(fLongFit,"RLQ0+"); string fitStatus = (string)(gMinuit->fCstatu); MSG("AngClusterFitAna",Msg::kDebug)<<"LongFit status: "<<fitStatus<<endl; if(fitStatus=="CONVERGED " && fLongFit->GetParameter(0)<29.9 && fLongFit->GetNDF()>0 && totalHitEnergy > 0.0){ fAngClusterFit.fACluFitParA =fLongFit->GetParameter(0); fAngClusterFit.fACluFitParB =fLongFit->GetParameter(1); fAngClusterFit.fACluFitParLongE0=fLongFit->GetParameter(2); fAngClusterFit.fACluFitShwMax =GetMaximum(fLongFit); fAngClusterFit.fACluFitLongChiSq=fLongFit->GetChisquare(); fAngClusterFit.fACluFitLongConv =1; fAngClusterFit.fACluFitLongNDF =fAngClusterFit.fACluFitLongChiSq/fLongFit->GetNDF(); fAngClusterFit.fACluFitE0EnergyRatio =fAngClusterFit.fACluFitParLongE0/totalHitEnergy; }//if(fitStatus) //Do the transverse fit fTransFit->SetParameters(1.0,1.0,1.0,totalHitEnergy); fTransHitHistoRad->Fit(fTransFit,"RLQ0+"); fitStatus = (string)(gMinuit->fCstatu); MSG("AngClusterFitAna",Msg::kDebug)<<"TransFit status: "<<fitStatus<<endl; if(fitStatus=="CONVERGED " && fTransFit->GetNDF()>0 ){ fAngClusterFit.fACluFitParL1 =fTransFit->GetParameter(0); fAngClusterFit.fACluFitParL2 =fTransFit->GetParameter(1); fAngClusterFit.fACluFitParC12 =fTransFit->GetParameter(2); fAngClusterFit.fACluFitParTransE0=fTransFit->GetParameter(3); fAngClusterFit.fACluFitTransChiSq=fTransFit->GetChisquare(); fAngClusterFit.fACluFitTransConv =1; fAngClusterFit.fACluFitTransNDF =fAngClusterFit.fACluFitTransChiSq/fTransFit->GetNDF(); } MSG("AngClusterFitAna",Msg::kDebug)<< "total energy " << totalHitEnergy << " from long histo " << fLongHitHistoRad->Integral() << " from trans histo " << fTransHitHistoRad->Integral() << endl; }

Double_t AngClusterFitAna::GetMaximum (  TF1 *  efit, Double_t  xmin = 0, Double_t  xmax = 0 )  [private]

Definition at line 644 of file AngClusterFitAna.cxx. Referenced by FitShower(). { Double_t fXmin = 0.001; Double_t fXmax = 100; Double_t fNpx = 40; 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 GetMaximum(efit, TMath::Max(xmin,xxmax-dx), TMath::Min(xmax,xxmax+dx)); }

void AngClusterFitAna::Reset (   )

Definition at line 638 of file AngClusterFitAna.cxx. { }

void AngClusterFitAna::Set3DHit (  DeqFloat_t &  x, DeqFloat_t &  y, DeqFloat_t &  z, DeqFloat_t &  e )

Definition at line 102 of file AngClusterFitAna.cxx. References fE, fX, fY, and fZ. Referenced by NueRecordAna::Analyze(), and NueDisplayModule::Analyze(). { if(x.size()!=0 && y.size()!=0 && z.size()!=0 && e.size()!=0){ fX=x; fY=y; fZ=z; fE=e; } }

void AngClusterFitAna::SetAngCluster (  Int_t &  primShow, DeqDeqInt_t &  clusterMap, TVector3 &  primDir )

Definition at line 118 of file AngClusterFitAna.cxx. References fClusterMap, fPrimDir, and fPrimShow. Referenced by NueRecordAna::Analyze(), and NueDisplayModule::Analyze(). { if(clusterMap.size()!=0){ fClusterMap=clusterMap; fPrimShow =primShow; fPrimDir =primDir; } else fPrimShow=ANtpDefVal::kInt; //Will trigger stopping of //further computations. }

void AngClusterFitAna::TransVarCalc (  TH1F *  transHisto  )

Definition at line 395 of file AngClusterFitAna.cxx. References AngClusterFit::fACluFitAsymPeak, AngClusterFit::fACluFitAsymVert, AngClusterFit::fACluFitKurt, AngClusterFit::fACluFitMean, AngClusterFit::fACluFitMolRadPeak, AngClusterFit::fACluFitMolRadVert, AngClusterFit::fACluFitRMS, AngClusterFit::fACluFitSkew, fAngClusterFit, MSG, and total(). Referenced by Analyze(). { Int_t kTransHistoBin = static_cast <Int_t> (transHisto->GetNbinsX()); Int_t binMax =transHisto->GetMaximumBin(); Int_t binVert =Int_t(((kTransHistoBin-1)/2)+1); Float_t peakBin=transHisto->Integral(binMax,binMax); Float_t vertBin=transHisto->Integral(binVert,binVert); Float_t total =transHisto->Integral(1,kTransHistoBin); Float_t asymPeak= ANtpDefVal::kFloat; Float_t asymVert= ANtpDefVal::kFloat; MSG("AngClusterFitAna",Msg::kDebug)<<"total: " <<total<<" binMax: "<<binMax <<" binVert "<<binVert <<" peakBin: "<<peakBin <<" kTransHistBin "<<kTransHistoBin <<endl; MSG("AngClusterFitAna",Msg::kDebug)<<"HistoI1: " <<transHisto->Integral(binMax+1 ,kTransHistoBin) <<endl; MSG("AngClusterFitAna",Msg::kDebug)<<"HistoI1: " <<transHisto->Integral(1, binMax-1) <<endl; MSG("AngClusterFitAna",Msg::kDebug)<<"HistoI1: " <<transHisto->Integral(binVert+1 ,kTransHistoBin) <<endl; MSG("AngClusterFitAna",Msg::kDebug)<<"HistoI1: " <<transHisto->Integral(1, binVert-1) <<endl; if(TMath::Abs(total-peakBin) > 1e-5){ MSG("AngClusterFitAna",Msg::kDebug)<<"Evaluating asymPeak: " <<TMath::Abs(transHisto ->Integral(binMax+1 ,kTransHistoBin) -transHisto ->Integral(1,binMax-1)) <<" divided by "<<(total-peakBin) <<endl; asymPeak=(TMath::Abs(transHisto->Integral(binMax+1,kTransHistoBin) -transHisto->Integral(1,binMax-1))) /(total-peakBin); } if(TMath::Abs(total-vertBin) > 1e-5){ MSG("AngClusterFitAna",Msg::kDebug)<<"Evaluating asymVert: "<< TMath::Abs(transHisto->Integral(binVert+1,kTransHistoBin) -transHisto->Integral(1,binVert-1)) <<" divided by "<<(total-peakBin)<<endl; asymVert=(TMath::Abs(transHisto->Integral(binVert+1,kTransHistoBin) -transHisto->Integral(1,binVert-1))) /(total-vertBin); } MSG("AngClusterFitAna",Msg::kDebug)<<"asymPeak: "<<asymPeak <<" asymVert: "<<asymVert<<endl; Float_t ratio,molRadPeak=0,molRadVert=0; if(total){ for(Int_t i=0; i<=binVert-1;i++){ ratio=transHisto->Integral(binMax-i>0?binMax-i :1,binMax+i<kTransHistoBin ?binMax+i:kTransHistoBin)/total; if(ratio>0.90){molRadPeak=i+1; break;} } for(Int_t i=0; i<=binVert-1;i++){ ratio=transHisto->Integral(binVert-i,binVert+i)/total; if(ratio>0.90){molRadVert=i+1; break;} } } Float_t mean=transHisto->GetMean(); Float_t rms =transHisto->GetRMS(); Float_t skew=0,kurt=0; Int_t nCount=0; MSG("AngClusterFitAna",Msg::kDebug)<<"mean: "<<mean<<" rms: "<<rms<<endl; for(Int_t i=1;i<=kTransHistoBin;i++){ if(transHisto->GetBinContent(i)){ skew=skew+(TMath::Power((transHisto->GetBinCenter(i)-mean),3) *transHisto->GetBinContent(i)); kurt=kurt+(TMath::Power((transHisto->GetBinCenter(i)-mean),4) *transHisto->GetBinContent(i)); nCount++; } } if(rms>1e-8 && nCount>1){ skew=skew /(static_cast<Float_t>(nCount-1)*TMath::Power((rms),3.)); kurt=(kurt/(static_cast<Float_t>(nCount-1)*TMath::Power((rms),4.)))-3.; if(TMath::Abs(skew) > 1.e8 ) skew = ANtpDefVal::kFloat; if(TMath::Abs(kurt) > 1.e8 ) kurt = ANtpDefVal::kFloat; } else{skew= ANtpDefVal::kFloat; kurt= ANtpDefVal::kFloat;} //Fill variables fAngClusterFit.fACluFitAsymPeak=asymPeak; fAngClusterFit.fACluFitAsymVert=asymVert; fAngClusterFit.fACluFitMolRadPeak=molRadPeak; fAngClusterFit.fACluFitMolRadVert=molRadVert; fAngClusterFit.fACluFitMean=mean; fAngClusterFit.fACluFitRMS=rms; fAngClusterFit.fACluFitSkew=skew; fAngClusterFit.fACluFitKurt=kurt; MSG("AngClusterFitAna",Msg::kDebug) <<"Leaving AngClusterFitAna::TransCalc" <<endl; }

---

Member Data Documentation

AngClusterFit& AngClusterFitAna::fAngClusterFit [private]

Definition at line 96 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), FitShower(), and TransVarCalc().

DeqDeqInt_t AngClusterFitAna::fClusterMap [private]

Definition at line 110 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and SetAngCluster().

DeqFloat_t AngClusterFitAna::fE [private]

Definition at line 106 of file AngClusterFitAna.h. Referenced by Analyze(), and Set3DHit().

TF1* AngClusterFitAna::fLongFit

Definition at line 75 of file AngClusterFitAna.h. Referenced by Draw(), FitShower(), and ~AngClusterFitAna().

TH1F* AngClusterFitAna::fLongHitHisto

Definition at line 68 of file AngClusterFitAna.h. Referenced by Analyze(), and ~AngClusterFitAna().

TH1F* AngClusterFitAna::fLongHitHistoRad

Definition at line 69 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), FitShower(), and ~AngClusterFitAna().

TVector3 AngClusterFitAna::fPrimDir [private]

Definition at line 111 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and SetAngCluster().

Int_t AngClusterFitAna::fPrimShow [private]

Definition at line 109 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and SetAngCluster().

TH3F* AngClusterFitAna::fScatterAxis

Definition at line 73 of file AngClusterFitAna.h. Referenced by Draw(), and ~AngClusterFitAna().

TF1* AngClusterFitAna::fTransFit

Definition at line 76 of file AngClusterFitAna.h. Referenced by Draw(), FitShower(), and ~AngClusterFitAna().

TH1F* AngClusterFitAna::fTransHitHisto

Definition at line 70 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and ~AngClusterFitAna().

TH1F* AngClusterFitAna::fTransHitHistoRad

Definition at line 71 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), FitShower(), and ~AngClusterFitAna().

DeqFloat_t AngClusterFitAna::fX [private]

Definition at line 103 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and Set3DHit().

DeqFloat_t AngClusterFitAna::fY [private]

Definition at line 104 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and Set3DHit().

DeqFloat_t AngClusterFitAna::fZ [private]

Definition at line 105 of file AngClusterFitAna.h. Referenced by Analyze(), Draw(), and Set3DHit().

---

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

• AngClusterFitAna.h
• AngClusterFitAna.cxx

---