GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Fast Simulation- neutral particle
Fast Simulation- neutral particle [message #16126] Mon, 31 March 2014 16:24 Go to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi all,

I try to do fast simulation with this decay chain: pbp->Ds+ Ds-; Ds- ->K+K-pi-; Ds+ ->e+ nv eta; eta->pi+pi-pi0 (gamma gamma). From the examples in macro/scrut, I got some problems for simulating the neutral channel:
1) in simulation macro:
  // generate electro-magnetic / hadronic split offs in the EMC? switch off when running w/o EMC
  Bool_t enableSplitoff   = kFALSE;

  //enable the producting of parametrized neutral (hadronic) split offs
  if (enableSplitoff)
    fastSim->EnableSplitoffs(BaseDir+"/fsim/splitpars.dat");

From the comments there, if I understand correctly, the "enableSplitoff" seems should be switched on for eta->pi+pi-pi0 (gamma gamma). I tried kTRUE, but it crashed.

2) in analysis:
	// *** name of the only PidAlgo TClonesArray in fsim
	TString pidalg = "PidChargedProbability";


Does it mean no neutral algorithm in fast simulation? then, how to fill the neutral list?

Thanks in advance.

Lu
Re: Fast Simulation- neutral particle [message #16127 is a reply to message #16126] Tue, 01 April 2014 08:09 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: *gsi.de
Dear Lu,


actually you don't need the turn on that split-offs switch for your channel. The generation of e.-m. and hadronic split-offs was supposed to somehow remodel neutral secondaries induced by charged tracks. However the code is quite old and was parametrized in the 'old' Babar-like framework. Therefore it should be taken with care. But at least it shouldn't crash, so I'm gonna take a look.


Best,
Klaus
Re: Fast Simulation- neutral particle [message #16130 is a reply to message #16126] Tue, 01 April 2014 09:09 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: *gsi.de
Dear Lu,


I just updated macro/scrut concerning the split-off generation. If you like you can svn up that directory and retry with using split-offs to see the difference in the neutrals multiplicity. Please let me know if it still crashes.


Best,
Klaus
Re: Fast Simulation- neutral particle [message #16136 is a reply to message #16130] Tue, 01 April 2014 11:05 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Dear Klaus,

Thanks. It works now.
Then, coming to the second question, since "PidChargedProbability" is for charged particle, how to fill the neutral list?

Best regards,
Lu
Re: Fast Simulation- neutral particle [message #16145 is a reply to message #16136] Tue, 01 April 2014 11:43 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *dip0.t-ipconnect.de
How do you fill photon list in the full reco? It should be the same.
Re: Fast Simulation- neutral particle [message #16149 is a reply to message #16145] Tue, 01 April 2014 12:46 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
In the full reconstruction, I used
theAnalysis->FillList(gam,"Neutral");

but it makes the ana macro crash in the fast sim with this warning:
Warning in <PndAnaPidCombiner::PndAnaPidCombiner::ReadTCA()>: No "PidAlgoIdealCharged" array found.
 *** Break *** segmentation violation


Best regards,
Lu
Re: Fast Simulation- neutral particle [message #16153 is a reply to message #16145] Tue, 01 April 2014 13:24 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: *gsi.de
Yes, just take the same name 'PidChargedProbability' for neutrals, I think without that parameter the FillList crashes (Ralf is looking at it...).

Best,
Klaus
Re: Fast Simulation- neutral particle [message #16178 is a reply to message #16153] Wed, 02 April 2014 15:40 Go to previous messageGo to next message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ikp.kfa-juelich.de
Dear Klaus and Ralf,

it is still not clear to me how we should fill up the neutral list in the new release.
I am running in the trunk 24275 right now.

For example, in the old release if I type:

theAnalysis->FillList(gammas, "Neutral");

this works ["works"....great word! At least it does not produce a crash]

In the new /macro/scrut for charged particles one has to fill up lists in this way:

theAnalysis->FillList(kplus, "KaonAllPlus", pidalg); //example for kaon+ list


But for neutral particles I think it is wrong to write:

theAnalysis->FillList(gammas, "Neutral", pidalg);


especially if I read in the macro:

// *** name of the only PidAlgo TClonesArray in fsim
TString pidalg = "PidChargedProbability";


Now, to write simply:
theAnalysis->FillList(gammas, "Neutral");

gives a Seg Fault, as the info concerning the pid-algorithm MUST be passed, for how this is defined.


So, my question follows that posed from Lu in the past days: how shall we fill up the list of neutrals in the new release? Without this, here we are stacked, as in all channels we try to reconstruct, at least one photon gets involved.

Thank you for your suggestion,

Elisabetta
Re: Fast Simulation- neutral particle [message #16180 is a reply to message #16178] Wed, 02 April 2014 15:51 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 500
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hi,

I'm on it. You should never have been forced to provide a PidAlgo.

Cheers
Ralf
Re: Fast Simulation- neutral particle [message #16181 is a reply to message #16180] Wed, 02 April 2014 15:59 Go to previous messageGo to next message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ikp.kfa-juelich.de
Hi Ralf,

thank you to take care of this. If I do not write any pid information in the neutral list, I got a Segmentation Fault. So I cannot run that macro at all.

If I define a new string:


TString pidalg0 ="PidNeutralProbability";

and I pass this parameter, actually it runs. Would it be ok, in your opinion?

I see that in PidClassifier/PndPidIdealAssociatorTask.cxx, "PidNeutralProbability"; exists.

cheers, Elisabetta
Re: Fast Simulation- neutral particle [message #16184 is a reply to message #16181] Wed, 02 April 2014 16:18 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 500
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hi Elisabetta,

Yes, it would be OK as long as you don't access PID information. I.e. don't select on the numbers from the charged PID in your neutral candidates.

Cheers
Ralf
Re: Fast Simulation- neutral particle [message #16185 is a reply to message #16180] Wed, 02 April 2014 16:20 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Ralf,

Regarding to the issue of FillList, I just found that scrut/ana_jpsi.C crashes due to this operation:

theAnalysis->FillList(mclist, "McTruth");

It doesn't work any more in the version of 24297, but works with the previous one #24275 (Elisabetta tests).
Is it caused by the on-going modifications of code? and can be fixed in soon?


Thanks in advance.

Best,
Lu
Re: Fast Simulation- neutral particle [message #16186 is a reply to message #16185] Wed, 02 April 2014 16:37 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 500
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hi.

I checked in the fix for not providing a pidalgo.

Regarding the issue with FillList(... "MCTruth") I found nothing wrong in my fresh trunk. Please see with svn diff if you have some accidental changes.

Cheers
Ralf
Re: Fast Simulation- neutral particle [message #16188 is a reply to message #16186] Wed, 02 April 2014 17:33 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi,
I didn't touch the example code...and it still doesn't work even for the newly download. Sad
uh...what can I do anyway?

Best,
Lu
Re: Fast Simulation- neutral particle [message #16189 is a reply to message #16188] Wed, 02 April 2014 17:43 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: *gsi.de
Hi Lu,


I just checked out mar14 and run the example with

root -l -b -q 'simfast.C("psi","psi2s_jpsi2mu_2pi.dec",6.23,1000,"psi(2S)")'
root -l -b 'ana_jpsi.C("psi_fast.root")'


and it runs without any problem.

Did you start a new shell for checking out and compiling and stuff? It can always be that your environment is spoiled, which leads to strange effects.


Best,
Klaus
Re: Fast Simulation- neutral particle [message #16191 is a reply to message #16189] Wed, 02 April 2014 18:09 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Klaus,

Thanks for your checking. My problem sounds pretty strange somehow.. Confused The working environment looks fine, and ana_jpsi.C can work when I comment the line filling the MC list.
By the way, I met a problem in compiling mar14:
[ 93%] Built target LmdTool                                                                                                                  
[ 93%] Building CXX object lmd/CMakeFiles/Model.dir/LmdFit/ModelFramework/fit/ModelFitFacade.cxx.o
/private/pandaroot_mar14/pandaroot/lmd/LmdFit/ModelFramework/fit/ModelFitFacade.cxx: In member function 'ModelFitResult ModelFitFacade::Fit()':
/private/pandaroot_mar14/pandaroot/lmd/LmdFit/ModelFramework/fit/ModelFitFacade.cxx:68:11: error: no match for 'operator==' in '0 == ((ModelFitFacade*)this)->ModelFitFacade::estimator'
/private/pandaroot_mar14/pandaroot/lmd/LmdFit/ModelFramework/fit/ModelFitFacade.cxx:68:11: note: candidate is: operator==(int, int) <built-in>
/private/pandaroot_mar14/pandaroot/lmd/LmdFit/ModelFramework/fit/ModelFitFacade.cxx:81:11: error: no match for 'operator==' in '0 == ((ModelFitFacade*)this)->ModelFitFacade::minimizer'
/private/pandaroot_mar14/pandaroot/lmd/LmdFit/ModelFramework/fit/ModelFitFacade.cxx:81:11: note: candidate is: operator==(int, int) <built-in>
make[2]: *** [lmd/CMakeFiles/Model.dir/LmdFit/ModelFramework/fit/ModelFitFacade.cxx.o] Error 1
make[1]: *** [lmd/CMakeFiles/Model.dir/all] Error 2
make: *** [all] Error 2


Best,
Lu
Re: Fast Simulation- neutral particle [message #16192 is a reply to message #16191] Wed, 02 April 2014 18:17 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: *gsi.de
Hi Lu,


do you work on your own computer? Maybe it's something with the external packages? Because it just worked for me half an hour ago on my GSI machine.


Best,
Klaus
Re: Fast Simulation- neutral particle [message #16193 is a reply to message #16192] Wed, 02 April 2014 18:22 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Klaus,

I'm using own computer with the external package of apr13. There's no problem before today...

Best,
Lu
Re: Fast Simulation- neutral particle [message #16196 is a reply to message #16191] Thu, 03 April 2014 11:10 Go to previous messageGo to next message
Stefan Pflueger is currently offline  Stefan Pflueger
Messages: 99
Registered: February 2012
continuous participant
From: *superkabel.de
Hello Lu,

I uploaded a fix for this compilation error. Could you give it a test and let me know if it works? Thank you.

Cheers,
Stefan
Re: Fast Simulation- neutral particle [message #16202 is a reply to message #16196] Thu, 03 April 2014 14:06 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Hi Stefan,

it still has the same error during installing.

Best,
Lu
Re: Fast Simulation- neutral particle [message #16208 is a reply to message #16202] Thu, 03 April 2014 15:01 Go to previous messageGo to next message
Stefan Pflueger is currently offline  Stefan Pflueger
Messages: 99
Registered: February 2012
continuous participant
From: *superkabel.de
Hi,

@Lu:
ah sry, I just saw that you are using the mar14 release... That would explain the same error. I commited this fix into the trunk. I actually do not know if I'm allowed to commit to the release code.

@all:
How am I supposed to commit this change into the release branch? I tried to switch this file to the mar14 version, apply the fix and then commit this file. However the permission for the commit was denied. Is this the wrong way to fix it or am I simply missing the permissions.

Best regards,

Stefan
Re: Fast Simulation- neutral particle [message #16209 is a reply to message #16208] Thu, 03 April 2014 15:10 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.159.116*
I will do this, just give me the time to reach a PC.
Re: Fast Simulation- neutral particle [message #16210 is a reply to message #16208] Thu, 03 April 2014 15:45 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 93.48.231*
Now mar14 is updated, please try now.
Re: Fast Simulation- neutral particle [message #16213 is a reply to message #16210] Thu, 03 April 2014 16:11 Go to previous message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
The installation works now, and the example in scrut also works. Razz

Thanks to all helps.

Lu
Previous Topic: Fast Sim Forum
Next Topic: Vertex and kinematic fitter for full simulation.
Goto Forum:
  


Current Time: Mon Nov 18 10:14:30 CET 2019

Total time taken to generate the page: 0.02551 seconds