GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » PndKinVtxFitter Resoultion
PndKinVtxFitter Resoultion [message #15030] Mon, 29 July 2013 14:57 Go to previous message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *ext.kfa-juelich.de
Hello!

I tried to do a resolution plot for my PndKinVtxFitter Vertex Resolution.
If I understand correctly the distance is in cm. If this is true I have a vertex resolution in the order of mm Confused which is in my opinion really bad.
Maybe I compare the wrong values with eatch other. Can someone have a look at my code?

Best regards

Simone


What I did:

	lambda2.Cleanup();
	lambda2.Combine(protonplus, pionminus); // Daugther 0 proton, Daughter 1 pion
	lambda2.SetType(3122);

	// **********************
	// PndKinVtxFitter
	// **********************
	double smallestchi2 =100000000000;
	int indexofsmallestchi2 =-1;
	RhoCandidate* lambdaPndKinVtxFitterfit;

	for(int k=0;k < lambda2.GetLength();k++)
	{
		PndKinVtxFitter kinvtxfitter(lambda2[k]);
		kinvtxfitter.AddMassConstraint(1.115);
		kinvtxfitter.Fit();
		if(kinvtxfitter.GetChi2()< smallestchi2 && kinvtxfitter.GetChi2() > 0)
		{
			smallestchi2= kinvtxfitter.GetChi2();
			indexofsmallestchi2=k;
			lambdaPndKinVtxFitterfit=lambda2[k]->GetFit();
			std::cout << "smallest chi2 " << smallestchi2<< std::endl;
		}
	}

	if(indexofsmallestchi2!=-1)
	{
fLambdaMassPndKinVtxFitter->Fill((lambdaPndKinVtxFitterfit->Daughter(0)->P4()+lambdaPndKinVtxFitterfit->Daughter(1)->P4()).M());
		fLambdaMassPndKinVtxFitterChi2->Fill(smallestchi2);
	    fPndKinVtxFitter_VertexRes_X->Fill(lambdaPndKinVtxFitterfit->Daughter(0)->GetPosition().X()-lambda2[indexofsmallestchi2]->Daughter(0)->GetMcTruth()->GetPosition().X());
	    fPndKinVtxFitter_VertexRes_Y->Fill(lambdaPndKinVtxFitterfit->Daughter(0)->GetPosition().Y()-lambda2[indexofsmallestchi2]->Daughter(0)->GetMcTruth()->GetPosition().Y());
	    fPndKinVtxFitter_VertexRes_Z->Fill(lambdaPndKinVtxFitterfit->Daughter(0)->GetPosition().Z()-lambda2[indexofsmallestchi2]->Daughter(0)->GetMcTruth()->GetPosition().Z());
}


index.php?t=getfile&id=7469&private=0
index.php?t=getfile&id=7470&private=0
index.php?t=getfile&id=7471&private=0
  • Attachment: resx.png
    (Size: 8.76KB, Downloaded 502 times)
  • Attachment: resy.png
    (Size: 8.70KB, Downloaded 590 times)
  • Attachment: resz.png
    (Size: 8.31KB, Downloaded 606 times)
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: [FIXED] why is total energyloss always too small?
Next Topic: Documentation for Rho classes
Goto Forum:
  


Current Time: Fri Mar 29 09:22:41 CET 2024

Total time taken to generate the page: 0.01121 seconds