Home » PANDA » PandaRoot » Analysis » Vertex Fitter in Fast simulation
Re: Vertex Fitter in Fast simulation [message #16373 is a reply to message #16371] |
Fri, 18 April 2014 22:26   |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
 From: *ip43.fastwebnet.it
|
|
Hi,
I think I found the problem.
First, the covariance matrix was not filled for neutrals, as you were saying. Second, in the EMC class the error in Energy was set but not the error in momentum. then all the elements with p were set to zero.
I fixed in my private version of the code, and I can see something reasonable. I wait for the Klaus' OK before committing the code.
Meanwhile, if you wanna try, these are the changes:
Index: ../../fsim/PndFsmEmcBarrel.cxx
===================================================================
--- ../../fsim/PndFsmEmcBarrel.cxx (revisione 24552)
+++ ../../fsim/PndFsmEmcBarrel.cxx (copia locale)
@@ -107,12 +107,14 @@
if (wasDetected && fabs(t->charge())<1e-8)
{
result->setdE(dE(t));
+ result->setdp(dE(t));
result->setdphi(dphi(t));
result->setdtheta(dtheta(t));
}
else
{
result->setdE(0.);
+ result->setdp(0.);
result->setdphi(0.);
result->setdtheta(0.);
}
Index: ../../fsim/PndFastSim.cxx
===================================================================
--- ../../fsim/PndFastSim.cxx (revisione 24552)
+++ ../../fsim/PndFastSim.cxx (copia locale)
@@ -465,7 +465,7 @@
RhoCandidate tcand(ft->p4(),ft->charge(),svtx);
- if( (fPropagate||fUseFlatCovMatrix) && fabs(ft->charge())>1e-6)
+ if( (fPropagate||fUseFlatCovMatrix))// && fabs(ft->charge())>1e-6)
{
tcand.SetCov7(ft->Cov7());
pidCand->SetCov7( ft->Cov7() );
Index: ../../fsim/PndFsmEmcBwCap.cxx
===================================================================
--- ../../fsim/PndFsmEmcBwCap.cxx (revisione 24552)
+++ ../../fsim/PndFsmEmcBwCap.cxx (copia locale)
@@ -99,12 +99,14 @@
if (wasDetected && fabs(t->charge())<1e-8)
{
result->setdE(dE(t));
+ result->setdp(dE(t));
result->setdphi(dphi(t));
result->setdtheta(dtheta(t));
}
else
{
result->setdE(0.);
+ result->setdp(0.);
result->setdphi(0.);
result->setdtheta(0.);
}
Index: ../../fsim/PndFsmEmcFwCap.cxx
===================================================================
--- ../../fsim/PndFsmEmcFwCap.cxx (revisione 24552)
+++ ../../fsim/PndFsmEmcFwCap.cxx (copia locale)
@@ -100,12 +100,14 @@
if (wasDetected && fabs(t->charge())<1e-8)
{
result->setdE(dE(t));
+ result->setdp(dE(t));
result->setdphi(dphi(t));
result->setdtheta(dtheta(t));
}
else
{
result->setdE(0.);
+ result->setdp(0.);
result->setdphi(0.);
result->setdtheta(0.);
}
Index: ../../fsim/PndFsmEmcFS.cxx
===================================================================
--- ../../fsim/PndFsmEmcFS.cxx (revisione 24552)
+++ ../../fsim/PndFsmEmcFS.cxx (copia locale)
@@ -100,12 +100,14 @@
if (wasDetected && fabs(t->charge())<1e-8)
{
result->setdE(dE(t));
+ result->setdp(dE(t));
result->setdphi(dphi(t));
result->setdtheta(dtheta(t));
}
else
{
result->setdE(0.);
+ result->setdp(0.);
result->setdphi(0.);
result->setdtheta(0.);
}
|
|
|
Goto Forum:
Current Time: Wed Mar 19 13:53:55 CET 2025
Total time taken to generate the page: 0.00715 seconds
|