| Fast Simulation- neutral particle [message #16126] |
Mon, 31 March 2014 16:24  |
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 #16178 is a reply to message #16153] |
Wed, 02 April 2014 15:40  |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
The working environment looks fine, and ana_jpsi.C can work when I comment the line filling the MC list.
