GSI Forum - RDF feed
https://forum.gsi.de/index.php
Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12439&th=3306#msg_12439
I have run today the new job on GRID for the tpc for pi+pi- analysis, with august11 release, and unfortunately the results are very strange:
- theta distribution peaks at 90°
- momentum distribution is around 0.5 GeV
- residual momentum peaks at 0.9
- for the invariant mass: only 40 events pass the selection on 100K events generated
I don't know what is the reason, but the results are very bad and strange,
Elisa
]]>Elisa Fioravanti2011-08-02T16:41:52-00:00Re: Run 952 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12441&th=3306#msg_12441
The run number if 952
Elisa]]>Elisa Fioravanti2011-08-02T20:46:40-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12446&th=3306#msg_12446
I see the same problem with eta_c channel, i.e. very low efficiency.
And the reason I suppose that outdated tpc reconstruction macro was used, at least if it is macro on grid /panda/user/p/pbarprod/tdr11/macros/pipi/run_reco_tpc_evt.C.
And it should be updated to trunk/macro/run/tdrct/npipi/run_reco_tpc_evt.C
with additional modification
the line
//tpcSPR->SetMCPid(); // use ideal particle identification
should be uncommented
and the line
tpcSPR->SetPDG(211);
should be commented.
The same is for other channels, except that for eta_c
Dima
]]>Dima Melnychuk2011-08-03T10:29:29-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12448&th=3306#msg_12448
this problem should not be related to the way of PID selection in the pattern recognition. We just tested both methods locally for the eta_c channel and still see the same efficiencies (approx. 28%).
Are the macros up-to date, so we can also test pi-pi+ locally? It seems like there is a bug somewhere in the code on the grid ...
Cheers
Felix]]>Felix Boehmer2011-08-03T11:36:45-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12449&th=3306#msg_12449
It's just what I have written, the old version of macro was used for TPC reconstruction on grid.
Running locally I also have close results for both PID selection in the pattern recognition.
So it's a question if/when Stefano will be able to correct macro and rerun.
Dima]]>Dima Melnychuk2011-08-03T11:47:02-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12450&th=3306#msg_12450
you can find the analysis macro in
/macro/run/tdrct/npipi/
- ana_invariantmass_2pi_tpc.C that I use to look at the invariant mass of two pions
- plot_invariant.C that I used to produce plots.
Concerning the reconstruction macro, It's not up-to date,
Elisa
]]>Elisa Fioravanti2011-08-03T11:54:54-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12451&th=3306#msg_12451
Would it be possible to modify all the tpc reco macros in svn, so that I have only to add them on the grid and submit the jobs? I do not have a computer in front of me right now...]]>StefanoSpataro2011-08-03T13:02:55-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12452&th=3306#msg_12452
we just tested eta_c with a completely clean trunk revision 12933, RKTrackRep revision 370 (grid version according to Johan). That is also true for the macros!
For both fixed PID and MCPID the efficiencies were around 30% (100 simulated events)...
So I don't understand what is going wrong there. Is the grid version really 12933 including all subfolders?
Cheers
Felix]]>Felix Boehmer2011-08-03T13:16:02-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12453&th=3306#msg_12453
I have activated
tpcSPR->SetMCPid();
in all tpc reco macros (run_reco_tpc_evt.C) in /trunk/macro/run/tdrct
So, they can be run on the grid.
Dima]]>Dima Melnychuk2011-08-03T13:16:30-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12454&th=3306#msg_12454
On the grid the macros are running not directly from
/macro/run/tdrct
but are copied to separate folder
/panda/user/p/pbarprod/tdr11/macros/
and in this folder they were outdated.
Dima]]>Dima Melnychuk2011-08-03T13:21:42-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12455&th=3306#msg_12455
thanks for the clarification, Dima
Still, if only the line about the MCPID was different ...
Anyway, with the consistent set of macros everything should be fine. Can you guys please inform us as soon as you see new results? We will try to test pi-pi+ locally in the meantime.
Cheers
Felix]]>Felix Boehmer2011-08-03T13:26:59-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12456&th=3306#msg_12456
StefanoSpataro2011-08-03T13:42:17-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12457&th=3306#msg_12457
From the modification time of macros in
/panda/user/p/pbarprod/tdr11/macros/
it seems that you could apply changes concerning
tpcSPR->SetMCPid();
but the macros didn't include previous modifications done by Johanes.
Dima]]>Dima Melnychuk2011-08-03T13:50:15-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12458&th=3306#msg_12458
we just took a look at the 2pion example for the npipi macros. From 100 simulated events we see 76 reconstructed entries in the invariant mass spectrum peaking at roughly 3.1 GeV, which is reasonable. All has been done with the macros from the trunk revision 12933.
Cheers
Felix]]>Felix Boehmer2011-08-03T15:51:05-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12459&th=3306#msg_12459
]]>StefanoSpataro2011-08-03T17:30:01-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12460&th=3306#msg_12460
*keeping fingers crossed*
Have a nice evening!
Felix]]>Felix Boehmer2011-08-03T17:50:22-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12463&th=3306#msg_12463
I have test the pi+pi-pi+pi- analysis with TPC data available on the GRID, and the results seem to be good and resonable, for example for the invariant mass I have an effiency of 73.61% for TPC
and for the STT data on the GRID 72.61%
calculated with the formula :
Efficiency = (Number of reconstructed events)/ (Number of generated events entering within geo. detector acceptance)
that we agreed in the last EVO meeting,
best
Elisa]]>Elisa Fioravanti2011-08-04T08:23:34-00:00Re: Run 951 GRID -TPC
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12464&th=3306#msg_12464
that is good to hear! Thanks to you and everyone involved!