GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Vertex and kinematic fitter for full simulation.
Vertex and kinematic fitter for full simulation. [message #16255] Tue, 08 April 2014 13:52 Go to previous message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *dip0.t-ipconnect.de
Hi Ralf,

Now I am trying to fit for full simulation data.
PndKinFitter with mass constraint is properly working, mass distribution is clearly improved but
PndKinVtxFitter (vertex fitter) doesn't work! All entries in mass histogram are existing but content is empty.

Could you check below code where I made mistake during the fit?
I am talking about full simulation, all spectrum without fit are surely correct.
I used mar14 and test again in scrut14. results are same for both.

Quote:

//Pnd Kin Vertex fitter
PndKinVtxFitter vtxfitter(dp[j]); // instantiate a vertex fitter
vtxfitter.Fit();
double chi2_vtx = vtxfitter.GetChi2(); // access chi2 of fit
double prob_vtx = vtxfitter.GetProb(); // access probability of fit
RhoCandidate *dpfit_vtx = dp[j]->GetFit(); // access the fitted cand
TVector3 Vtx_kin=dpfit_vtx->Pos(); // and the decay vertex position

hdp_fit_mass_vtx_all->Fill( dpfit_vtx->M() );
if (fAnalysis->McTruthMatch(dp[j]))
{
hdp_fit_mass_vtx_ftm->Fill( dpfit_vtx->M() );
}

//Pnd Kin Mass Fitter
PndKinFitter mfitter(dp[j]); // instantiate the kin fitter in psi(2S)
mfitter.AddMassConstraint(DpmMass); // set 4 constraint
mfitter.Fit(); // do fit
double chi2_m = mfitter.GetChi2(); // get chi2 of fit
double prob_m = mfitter.GetProb(); // access probability of fit
RhoCandidate *dpfit_mass = dp[j]->GetFit(); // get fitted J/psi
hdp_fit_mass_mass_all->Fill( dpfit_mass->M() );
if (fAnalysis->McTruthMatch(dp[j]))
{
hdp_fit_mass_mass_ftm->Fill( dpfit_mass->M() );
}


Thank you in advance...
Yours,
donghee
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Fast Simulation- neutral particle
Next Topic: (SOLVED)Obtain Raw Info from Calorimeter (old GetMicroCandidate from TCandList class)
Goto Forum:
  


Current Time: Thu Mar 28 17:50:47 CET 2024

Total time taken to generate the page: 0.00711 seconds