GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » [FIXED] PndAnalysis::McTruthMatch always returns false
icon9.gif  [FIXED] PndAnalysis::McTruthMatch always returns false [message #14960] Mon, 15 July 2013 14:21 Go to next message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *vpn.kfa-juelich.de
Hi everyone!

at the moment I try to reconstruct Lambda LambdaBar events with the new Rho Classes and PndAnalysis.
There is output which I dont understand


I have done the following:

theAnalysis->FillList(mctrack,"McTruth","McTruth");
theAnalysis->FillList(protonplus,"ProtonLoosePlus","PidAlgoIdealCharged");
theAnalysis->FillList(pionminus,"PionLooseMinus","PidAlgoIdealCharged");

lambda.Cleanup();
lambda.Combine(protonplus, pionminus); 
lambda.SetType(3122);
for (int r=0; r< lambda.GetLength();r++)
	{
		fLambdaMass->Fill( lambda[r]->M() );
		std::cout << "ndaughters " << lambda[r]->NDaughters() << std::endl;

		lambda[r]->Daughter(0)->PrintOn(std::cout);
		std::cout << std::endl;
		lambda[r]->Daughter(0)->GetMcTruth()->PrintOn(std::cout);
		std::cout << std::endl;
		std::cout << std::endl;
                theAnalysis->McTruthMatch(lambda[r],0,true)
	}



This gives me the following output

ndaughters 2
uid:2 trk:1 (-0.22592,0.0670787,0.84685;1.28571) 1e 0.938272 GeV/c2  (0,0,0,4) daughters: 0() pdg: 2212 PID:0,0,0,0,1,  mc truth pointer: 0x7fb13f203680
uid:5 trk:-1 (-0.225709,0.0602535,0.849309;1.28696) 1e 0.938272 GeV/c2  (0,0,0,32) daughters: 0() pdg: 2212 PID:-999999,-999999,-999999,-999999,-999999,  mc truth pointer: 0

uid:3 trk:2 (0.0527879,0.0460279,0.232573;0.280134) -1e 0.13957 GeV/c2  (0,0,0,8) daughters: 0() pdg: -211 PID:0,0,1,0,0,  mc truth pointer: 0x7fb13f2038b0
uid:6 trk:-1 (0.056507,0.047805,0.224331;0.274377) -1e 0.13957 GeV/c2  (0,0,0,64) daughters: 0() pdg: -211 PID:-999999,-999999,-999999,-999999,-999999,  mc truth pointer: 0

Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=2212)
Info in <PndMcTruthMatch::MctMatch>: rejected composite (pdg=3122) by non-matching daughter: idau=0



My questions are now:
Why is the treknumber trk from my MCTruth Candidates always -1? Should it not be the same number as above ( This is the position in the TCA, right? )

Why are ALL my MctMatch false? Somehow this dosnt work or I am doing it wrong

External Packages: Apr13
Pandaroot: REV 20654

Best Regards

Simone

[Updated on: Fri, 19 July 2013 09:42] by Moderator

Report message to a moderator

Re: PndAnalysis::McTruthMatch always returns false [message #14966 is a reply to message #14960] Wed, 17 July 2013 12:08 Go to previous messageGo to next message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *pool.mediaWays.net
I think there are two "return" statements missing. That is the reason why my MCTruthMatch always returns false.

In PndAnalysis.h in function
Bool_t MctMatch(RhoCandidate& c, RhoCandList& mct, Int_t level = 2, bool verbose=false)
(Line 55)

and in PndAnalysis.cxx in function
Bool_t PndAnalysis::McTruthMatch(RhoCandidate* cand, Int_t level, bool verbose)
(Line 764)

Can someone check this and add this to svn?

Thanks and Best Regards

Simone
Re: PndAnalysis::McTruthMatch always returns false [message #14968 is a reply to message #14966] Wed, 17 July 2013 22:00 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *dip0.t-ipconnect.de
Hi Simone,


thanks for that hint! I think you are right, and I guess Ralf will correct it asap.

Did you modify it in your own software to check whether this indeed fixes the problem?


Best,
Klaus
Re: PndAnalysis::McTruthMatch always returns false [message #14969 is a reply to message #14968] Thu, 18 July 2013 09:23 Go to previous messageGo to next message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *ext.kfa-juelich.de
Hi Klaus,

yes I modified my own version.
Before McTruthMatch found no single event which reached matchlevel 2, now there are ~40% of events with matchlevel 2 for Lambdas.

Is there any other way to control if the McTruthMatch works properly? I expect to have more events with Lambdas.


Best Regards

Simone
Re: PndAnalysis::McTruthMatch always returns false [message #14972 is a reply to message #14969] Thu, 18 July 2013 14:49 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *dip0.t-ipconnect.de
Hi Simone,


did you make sure, that all your Lambdas decay to p pi- in your decay file?

What you can do is to make the truth match 'by hand'. Since you have the McTruth list, you should be able to find the proper indices of the p and the pi- from the true Lambda and check, whether there are particles in your reco list with corresponding McIndex.

You can then also check the properties of the true candidates, which don't get a reco match. Maybe they have very small transverse momenta or whatsoever, so that they don't get reconstructed properly.


Best regards,
Klaus
Re: PndAnalysis::McTruthMatch always returns false [message #14975 is a reply to message #14960] Thu, 18 July 2013 15:45 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Hi Simone.

I put up the fix.

Ralf

Edit: MC candidates now also know the mc track number.

[Updated on: Thu, 18 July 2013 15:51]

Report message to a moderator

Re: PndAnalysis::McTruthMatch always returns false [message #14976 is a reply to message #14975] Thu, 18 July 2013 19:09 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
Good,
can Simone give the final OK, so that I can put [fixed] in the thread?
icon7.gif  Re: PndAnalysis::McTruthMatch always returns false [message #14978 is a reply to message #14976] Fri, 19 July 2013 06:24 Go to previous messageGo to next message
Simone Esch is currently offline  Simone Esch
Messages: 37
Registered: August 2010
Location: Forschungszentrum Jülich...
continuous participant

From: *pool.mediaWays.net
Everthing fine for me now Smile

Thank you all!


Simone
Re: PndAnalysis::McTruthMatch always returns false [message #14979 is a reply to message #14978] Fri, 19 July 2013 09:41 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Good,
one problem less Smile
Previous Topic: RHO: Major Modifications
Next Topic: Vertex fitters' problems
Goto Forum:
  


Current Time: Mon Nov 11 06:23:46 CET 2024

Total time taken to generate the page: 0.00941 seconds