GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » Update on eta_c reconstruction with STT and TPC
Update on eta_c reconstruction with STT and TPC [message #12112] Wed, 22 June 2011 00:13 Go to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *dip0.t-ipconnect.de
Hi all,

I have analyzed 10000 eta_c events with latest version of both STT and TPC reconstruction code.

First coming to STT. The low reconstruction efficiency of 10% persists after using corrected version of MVD, so the drop of efficiency have another origin and I have a suspect here.

But my main remark here is the low efficiency is not related to analysis macro but to small amount of event with 4 or more reconstructed charged candidates.

index.php?t=getfile&id=6467&private=0

And plot of invariant mass of phi and eta_c
index.php?t=getfile&id=6466&private=0

Coming to TPC I want to mention that I finally managed to run complete reconstruction chain till the rho analysis.

My version of TPC reconstruction macro /macro/run/tdrct/eta_c/run_reco_tpc_evt.C contains PndGFTrackToPndTrackConvertorTask task which convert GFTrack to PndTrack and I want here once more to ask to commit latest version of this task to svn (I don't have write access to GenfitTools directory). The file is attached. And to be more specific, maybe Stefano will commit it?

Obtained efficiency with TPC is even lower than with TPC (below 5%), which is also related to small number of events with 4 or more reconstructed charged candidates. And here I can quote the number of lost eta_c related to geometrical acceptance. In 67% of events at least one of the kaons is outside the range 140>theta>20 degree.

index.php?t=getfile&id=6469&private=0

And invariant mass resolution:

index.php?t=getfile&id=6468&private=0
looks comparable with STT results.

Trying to identify the source of low efficiency I compared the number of reconstructed track and number of charged candidates used for rho analysis. And this number appears to be lower by factor 2-3 (I checked it for TPC). So PndPidCorrelator, which creates charged candidates from PndTracks reject somehow significant amount of tracks. And I see that Stefano modified this class 5 days ago and approximately after that I observed drop of eta_c reconstruction efficiency with STT. So could Stefano look what could go wrong there?

Best regards,

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12113 is a reply to message #12112] Wed, 22 June 2011 00:32 Go to previous messageGo to next message
Johan Messchendorp is currently offline  Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen
first-grade participant

From: *xs4all.nl
Hi Dima,

PndGFTrackToPndTrackConvertorTask .... now updated in trunk.

j.
Re: Update on eta_c reconstruction with STT and TPC [message #12114 is a reply to message #12112] Wed, 22 June 2011 08:55 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 94.163.82*
Hi Dima,
Could you please try with a small sample and the previous correlator version?
Just to be sure, in theory my changes should have no effects on the reconstruction.
Re: Update on eta_c reconstruction with STT and TPC [message #12119 is a reply to message #12112] Wed, 22 June 2011 12:37 Go to previous messageGo to next message
Felix Boehmer is currently offline  Felix Boehmer
Messages: 149
Registered: May 2007
Location: Munich
first-grade participant

From: *gsi.de
Hi Dima,

I will try to look into the macros in the weekend to find out what goes wrong. The efficiencies you report are much too low, I will have to see what's going on there.


Cheers

Felix
Re: Update on eta_c reconstruction with STT and TPC [message #12125 is a reply to message #12112] Thu, 23 June 2011 01:07 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *play-internet.pl
Hi Stefano,

I should confirm that indeed previous correlator version doesn't change reconstruction efficiency.

But my another statement still remains true, that PndPidCorrelator produces for TPC 2-3 times less charged candidates than number of input tracks.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12143 is a reply to message #12125] Mon, 27 June 2011 13:57 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Hi,
I have tried to run some events using your reco macro w/o MCTrackAssociator (which does not work with tpc), but I have the following crash:

GEANE propagation failed
in line: 305 in file: /home/spataro/may11/pandaroot/GenfitTools/trackrep/GeaneTrackRep/GeaneTrackRep.cxx
with fatal flag 0
GFException thrown with excString:
GeaneTrackRep: PROTECT AGAINST LOW MOMENTA
in line: 178 in file: /home/spataro/may11/pandaroot/GenfitTools/trackrep/GeaneTrackRep/GeaneTrackRep.cxx
with fatal flag 0
RKTrackRep::RKutta ==> momentum too low: 2.38612 MeV
GFException thrown with excString:
RKTrackRep::Extrap ==> Runge Kutta propagation failed
in line: 997 in file: /home/spataro/may11/pandaroot/GenfitTools/trackrep/RKTrackRep/RKTrackRep.cxx
with fatal flag 0
PndTpcRiemannTrackingTask::Exec; Event Number: 131
RKTrackRep::setPDG particle -22 not known to TDatabasePDG -> abort



I don't know which code is using RKTrackRep, instead of GeaneTrackRep, and why -22 is used as pid hypothesis! I think we should use a standard hyp for the kalman and not the MC one, i.e. pion or muon.

I don't know hot to proceed after this crash, help from experts is required.

Re: Update on eta_c reconstruction with STT and TPC [message #12145 is a reply to message #12143] Mon, 27 June 2011 15:04 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Just to be more precise, the crash comes from PndTpcRiemannTrackingTask (I have commented out all the other tasks).
Re: Update on eta_c reconstruction with STT and TPC [message #12146 is a reply to message #12145] Mon, 27 June 2011 15:12 Go to previous messageGo to next message
Felix Boehmer is currently offline  Felix Boehmer
Messages: 149
Registered: May 2007
Location: Munich
first-grade participant

From: *natpool.mwn.de
Hi,

I will have a look at it asap.


Cheers

Felix
Re: Update on eta_c reconstruction with STT and TPC [message #12147 is a reply to message #12146] Mon, 27 June 2011 15:16 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
This is the output of the tpcriemanntask if i set the DEBUG to 1:

pdg: 211 charge: 3.000
MC pid pdg: -11 charge: 3.000
charge corrected MC pid pdg: -11
 center of track TVector3 A 3D physics vector (x,y,z)=(12.483792,3.663404,0.000000) (rho,theta,phi)=(13.010211,90.000000,16.354474)
 Radius of track [cm]: 3.778
 Dip of track [deg]:   46.17
 seed values: 
  start position: TVector3 A 3D physics vector (x,y,z)=(14.202656,7.027931,25.054993) (rho,theta,phi)=(29.645572,32.311851,26.327688)
  momentum [GeV]: 0.03149
  p_perp [GeV]:   0.02272
  direction: TVector3 A 3D physics vector (x,y,z)=(0.642457,-0.328217,-0.692476) (rho,theta,phi)=(1.000000,133.826418,-27.061496)
  winding: 1
  invertedTrack: 0
  pdg id: -11
Tracklet 3   nhits = 48
pdg: 211 charge: 3.000
MC pid pdg: 22 charge: 0.000
charge corrected MC pid pdg: -22
 center of track TVector3 A 3D physics vector (x,y,z)=(1.549914,37.385460,-0.000000) (rho,theta,phi)=(37.417574,90.000000,87.626011)
 Radius of track [cm]: 37.18
 Dip of track [deg]:   24.07
 seed values: 
  start position: TVector3 A 3D physics vector (x,y,z)=(-15.273078,4.230596,35.533201) (rho,theta,phi)=(38.907239,24.037374,164.517430)
  momentum [GeV]: 0.5482
  p_perp [GeV]:   0.2236
  direction: TVector3 A 3D physics vector (x,y,z)=(-0.363711,0.184549,0.913047) (rho,theta,phi)=(1.000000,24.070050,153.096469)
  winding: 1
  invertedTrack: 1
  pdg id: -22
RKTrackRep::setPDG particle -22 not known to TDatabasePDG -> abort



maybe it could help somehow.
Re: Update on eta_c reconstruction with STT and TPC [message #12150 is a reply to message #12112] Mon, 27 June 2011 17:54 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *baltnet.ru
Hi Stefano,

I want just to comment that from my side running 10 times 1000 event I didn't have such crush.

But I have a question here concerning MCTrackAssociator.
If it doesn't work with TPC can I still use in analysis macro MC based PID?

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12152 is a reply to message #12150] Mon, 27 June 2011 18:05 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
The MCTrackAssociator is not using the enum, while tpc reco code is still using it.
I can try to write a dummy MCTrackAssociatorEnum which uses the old enum. But this will be only a temporary fix.

I could also try to change tpc reco code in order to use the correct code and not the enum, to be compatible with the rest of the code, but for this I need to ask the OK from tpc developers.

Re: Update on eta_c reconstruction with STT and TPC [message #12153 is a reply to message #12152] Mon, 27 June 2011 18:14 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Hi Dima,
please update the package PndMCMatch and try to use PndMCTrackEnumAssociator. It should work... please let me know your feedback.

Re: Update on eta_c reconstruction with STT and TPC [message #12156 is a reply to message #12125] Tue, 28 June 2011 13:24 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Hi,
I have tried with a different sample, 500 events under eta_c folder, and I have seen:

1966 charged candidates
2865 PndTrackPostFitComplete tracks
2410 PndTrackPostFitComplete with flag>0.
1995 PndTrackPostFitComplete with flag>0. & p>0.1 & p<15

The remaining 1995-1966 missing candidates should be tracks were the back extrapolation to beam axis has failed.

Everything fine to me for what concerns the pidcorrelator.
Re: Update on eta_c reconstruction with STT and TPC [message #12352 is a reply to message #12112] Tue, 26 July 2011 16:20 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Hi,

I just noticed that the number of tracks in the "n charged" histogram is much too low.

I simulated 10 events and get 57 GFTracks. But in the "n charged" histogram there are only 15 tracks.

Also the PndPidCorrelator produces a lot of error messages, where propagation seems to go wrong.

I think we have problems with the PndPidCorrelator.

Could someone please take a look into it?

cheers,

Johannes
Re: Update on eta_c reconstruction with STT and TPC [message #12353 is a reply to message #12352] Tue, 26 July 2011 16:52 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.158.143*
Have you read my old message about low momentum values in tpcriemanntracks, two weeks ago? The correlator does not correlate tracks with p<100 MeV and p>10 GeV, for both tpc and stt. All the missing candidates are coming from very low momentum tracks.
Re: Update on eta_c reconstruction with STT and TPC [message #12354 is a reply to message #12112] Tue, 26 July 2011 17:34 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl

Concerning error messages from PndPidCorrelator, it really could be the place where tracks are lost.
With ideal pid in PndTpcRiemannTrackingTask I see the message from PndPidCorrelator

-W- PndPidCorrelator::GetTrackInfo :: Failed backward propagation

in average twice per event.

Whereas in case of specific particle hypothesis (pion or kaon) I see the same message 1 in 5 events.

This message is coming from PndPidTrackInfo.cxx and here tracks are rejected independent on momentum cut, actually this part of code is executed after cut on momentum.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12355 is a reply to message #12112] Tue, 26 July 2011 17:43 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Ah sorry, I haven't read this.

But I just fixed a little bug int the MVD correlator Task, there the RKTrackRep was not correctly copied and just GeaneTrackRep was used.

And I just tested with RK as cardinal Rep and there efficiencies seem to go up.

I just commited a new version of the RiemannTrackingTask, where Geane is no longer set as cardinal rep.

Dear Analysers, could you please check the efficiencies with the new versions?

regards,

Johannes
Re: Update on eta_c reconstruction with STT and TPC [message #12359 is a reply to message #12355] Wed, 27 July 2011 11:16 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Hello,

I would like to test the efficiency in the multipion analysis with the new version.
I run the macro run_sim_tpc_evt.C
in the folder /macro/run/tdrct/npipi/

but I have this error:

====== EMC 2):: ConstructASCIIGeometry() ===
=============================================
PndEmcReader: EMC geometry file == /home/fioravanti/fairsoft/pandaroot/geometry/emc_module12.dat
Emc module = 1
********
Emc module = 2
********

====== EMC:: ConstructRootGeometry() m3a ===
=============================================
File name = /home/fioravanti/fairsoft/pandaroot/geometry/emc_module3new.root

====== EMC:: ConstructRootGeometry() m4a ===
=============================================
fgeoName3:: /home/fioravanti/fairsoft/pandaroot/geometry/emc_module4_StraightGeo24.4 .root
File name Bw1= /home/fioravanti/fairsoft/pandaroot/geometry/emc_module4_StraightGeo24.4 .root
PndMdt::ConstructGeometry : No good MDT Barrel definition

And then root exits.
Could someone update the right macro to use?

thanks
elisa
Re: Update on eta_c reconstruction with STT and TPC [message #12361 is a reply to message #12112] Wed, 27 July 2011 11:42 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
Hi Elisa,

I assume that this problem is related to recent modifications of MDT classes. So if you use the trunk version of mdt you probably should use mdt geometry instead of
  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetBarrel("fast");
  Muo->SetEndcap("fast");
  Muo->SetMuonFilter("fast");
  Muo->SetMdtMagnet(kTRUE);
  Muo->SetMdtMFIron(kTRUE);
  fRun->AddModule(Muo);


the definition as in macro/mdt/sim_muo_dub.C
  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetMdtMagnet(kTRUE);
//  Muo->SetMdtMFIron(kFALSE);
  Muo->SetMdtCoil(kTRUE);
  Muo->SetBarrel("muon_TS_barrel_strip_v1_noGeo.root");
  Muo->SetEndcap("muon_TS_endcap_strip_v1_noGeo.root");
  Muo->SetForward("muon_Forward_strip_v1_noGeo.root");
  Muo->SetMuonFilter("muon_MuonFilter_strip_v1_noGeo.root");
  fRun->AddModule(Muo);


To be honest I didn't try it by myself, I use previously generated simulation data and repeat only from digitization stage.

But here I have a question to Johannes, if the modified padplane file comes to simulation or digitization stage?

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12363 is a reply to message #12112] Wed, 27 July 2011 12:12 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
Hi Elisa,

I have just thought that proper solution for the moment would be not to modify simulation macro but to modify the revision of mdt directory to the state before the update, i.e. 12752.

I.e. within mdt directory /pandaroot/mdt
type:
svn update -r 12752
and recompile.

Dima


Re: Update on eta_c reconstruction with STT and TPC [message #12364 is a reply to message #12354] Wed, 27 July 2011 12:19 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *19-87-r.retail.telecomitalia.it
Hi,
These tracks are skipped when geane backward extrapolation fails. This means that track parameters are wrong and cannot be back propagated... Maybe the track is not going to the inner part but outwards, this shoul happenin particular for secondaries... Not our case then. I think the output in momentum of theriemann tracks should be checked carefully.
Re: Update on eta_c reconstruction with STT and TPC [message #12365 is a reply to message #12112] Wed, 27 July 2011 12:26 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
And one more correction concerning mdt directory

The proper revision to go back is 12733,

i.e.
pandaroot/mdt$ svn update -r 12733

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12366 is a reply to message #12363] Wed, 27 July 2011 12:26 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Hello Dima,

thanks for your reply.
I have tried with svn update -r 12752 and I re-compiled but I have this error:

[ 57%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdt.cxx.o
[ 57%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtMagnet.cxx.o
[ 58%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtCoil.cxx.o
[ 58%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtMuonFilter.cxx.o
[ 58%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtMFIron.cxx.o
[ 58%] Building CXX object mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtGeoConstructorFast.cxx.o
/home/fioravanti/fairsoft/pandaroot/mdt/MdtMC/PndMdtGeoConstructorFast.c xx:34: error: no ‘void PndMdt::ConstructGeometryFast()’ member function declared in class ‘PndMdt’
make[2]: *** [mdt/CMakeFiles/Mdt.dir/MdtMC/PndMdtGeoConstructorFast.cxx.o] Error 1
make[1]: *** [mdt/CMakeFiles/Mdt.dir/all] Error 2
make: *** [all] Error 2


Instead, if I change the sim macro in this way:

PndMdt *Muo = new PndMdt("MDT",kTRUE);
Muo->SetMdtMagnet(kTRUE);
// Muo->SetMdtMFIron(kFALSE);
Muo->SetMdtCoil(kTRUE);
Muo->SetBarrel("muon_TS_barrel_strip_v1_noGeo.root");
Muo->SetEndcap("muon_TS_endcap_strip_v1_noGeo.root");
Muo->SetForward("muon_Forward_strip_v1_noGeo.root");
Muo->SetMuonFilter("muon_MuonFilter_strip_v1_noGeo.root");
fRun->AddModule(Muo);


everything works,
now, what is the best solution to run the simulation?

thanks a lot
Elisa
Re: Update on eta_c reconstruction with STT and TPC [message #12367 is a reply to message #12365] Wed, 27 July 2011 12:29 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Hello Dima,

now the simulation seems to work!

Thanks a lot
Elisa
Re: Update on eta_c reconstruction with STT and TPC [message #12368 is a reply to message #12361] Wed, 27 July 2011 12:56 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Hi Dima,

Quote:

But here I have a question to Johannes, if the modified padplane file comes to simulation or digitization stage?


the padplane file affects only digi and reco.

cheers

Johannes
Re: Update on eta_c reconstruction with STT and TPC [message #12369 is a reply to message #12112] Wed, 27 July 2011 15:26 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Hi,

I just run the eta_c macros for 100 events, and we have an efficiency of ~30% if we use only tpc (i.e. without MVD and GEM correlators).
If I add the MVD correlator I only find ~ 15 eta_cs.

So we definetly have a problem with our MVD and GEM correlators and are investigating!

cheers,

Johannes
Re: Update on eta_c reconstruction with STT and TPC [message #12370 is a reply to message #12367] Wed, 27 July 2011 15:27 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Hello everybody,

I would like to say that I'm running, for pi+pi- analysis, the reconstruction macro for the TPC, available into svn (update of pandaroot done this morning),
(macro/run/tdrct/npipi/run_reco_tpc_evt.C)

and I have al lot of this message:

in line: 305 in file:
/home/fioravanti/fairsoft/pandaroot/GenfitTools/trackrep/GeaneTrackRep/G eaneTrackRep.cxx
with fatal flag 0
GFException thrown with excString:
GEANE propagation failed

best,
Elisa
Re: Update on eta_c reconstruction with STT and TPC [message #12373 is a reply to message #12112] Thu, 28 July 2011 16:30 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Hi,

I worked on the MVD and GEM correlators and checked in new versions. The cuts are now in cm, so I updated all macros accordingly!

Please update (macros, tpc, externals(genfit and RKTrackRep)). Efficiencies for eta_c should now be around 30%.

regards,

Johannes
Re: Update on eta_c reconstruction with STT and TPC [message #12382 is a reply to message #12112] Thu, 28 July 2011 23:25 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *play-internet.pl
Hi Johannes,

I can only confirm that with 2000 events statistics eta_c reconstruction efficiency with TPC is around 30% as you have written. And mass distributions looks also good.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12383 is a reply to message #12382] Fri, 29 July 2011 00:03 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *21-79-r.retail.telecomitalia.it
Does it mean we could update july11 on the grid and rerun tpc jobs?
Which packages should be updated? Only tpc or maybe also something else?
Are there big differences with stt efficiency? Unfortunately I do not have these numbers with me.
Re: Update on eta_c reconstruction with STT and TPC [message #12384 is a reply to message #12112] Fri, 29 July 2011 00:36 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *play-internet.pl
Hi Stefano,

I also would like to clarify the question of data production on GRID.

But comming to TPC vs STT efficieny.
TPC efficiency = 28.9% (using vertex fit)
STT efficiency with GRID data =26.4% with all tracks and 23.7% with requirement of at least one STT hit in track.

STT efficiency with latest trunk version (2000 events) = 15.6%
(with requirement of at least one STT hit).
But I have no idea if the work on STT code is finished?
So I suppose that STT signal data should also be reproduced but the question is if to reproduce them together with TPC or later?

For TPC data reproduction in addition to tpc directory, genfit and RKTrackRep in GenfitTools/trackrep should be updated to rev 370 and updated tpc reco macros should be used.

And coming to the question of reconstruction of fully mixed events. I do not know how realistic is to expect to produce them by 5-7 August (the date when the GRID will be off), but without knowing when the code will be ready, there is still the question if the available DPM data should be used or they also needs to be reproduced. At least for TPC there were recent modifications related to clustering, which is the part of digitization, and I suppose TPC expert should define how important they are for DPM data which will be used for mixing.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12385 is a reply to message #12384] Fri, 29 July 2011 09:17 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *19-79-r.retail.telecomitalia.it
About Dpm events, they should be redone anyway.
Re: Update on eta_c reconstruction with STT and TPC [message #12386 is a reply to message #12385] Fri, 29 July 2011 10:24 Go to previous messageGo to next message
Johan Messchendorp is currently offline  Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen
first-grade participant

From: *xs4all.nl
Hi,

It would be good, if we could start some production over the weekend.
I will try to connect to the internet as much as I can (will travel to Barcelona today),
just let me know when (and what!) to update (in) the july11 release and the grid software.

Greets,

Johan.
Re: Update on eta_c reconstruction with STT and TPC [message #12387 is a reply to message #12386] Fri, 29 July 2011 10:47 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.158.180*
Probably it would be good f also the secondary problem with tpc detector would be fixed, in order to use ideal pid for tpc Riemann (at least, I suppose that is the problem)
Re: Update on eta_c reconstruction with STT and TPC [message #12389 is a reply to message #12384] Fri, 29 July 2011 14:33 Go to previous messageGo to next message
Gianluigi Boca is currently offline  Gianluigi Boca
Messages: 177
Registered: March 2004
first-grade participant
From: *gsi.de
Hi Dima,
I was still working on the Cleanup code for Stt in the last days.
Maybe you run with some non-final version and that's the cause of this sudden 15.6% efficiency.

So last night I put the final version (...yes, right Laughing Laughing Laughing ) of the Stt PR+Cleanup code in svn. I hope the efficiency for Stt returns to the previous values.

Also I updated all the analysis Macros (run_reco_stt_evt.C and
run_reco_stt_mix.C) for all the physics channels (eta_c,
llbar, npipi, psi3770) by adding the following line necessary
to activate the Cleanup code :

SttMvdTracking->Cleanup();

Please notice that in the non-mixing Macros (run_reco_stt_evt.C)
this line is COMMENTED OUT for the time being, because I would like to see if possible, if there are problems with the various analysis and no Cleanup.

If everything is fine, then activating the Cleanup should be beneficial even for the non-mixed events (it rejectes spurious tracks) even though I would expect a bit of a lowering in the reconstruction efficiency.

Lia has also updated the run_reco_stt_evt.C and run_reco_stt_mix.C adding the task of the extrapolation to the Gem region.
She is also checking the last details today. By tonight let's say at 20:00 all the Macro should be ready and can be used on the Grid. She will send aroung a message on the Forum when dona.

Tschuess an alle

Gianluigi

Dima Melnychuk wrote on Fri, 29 July 2011 00:36

Hi Stefano,

I also would like to clarify the question of data production on GRID.

But comming to TPC vs STT efficieny.
TPC efficiency = 28.9% (using vertex fit)
STT efficiency with GRID data =26.4% with all tracks and 23.7% with requirement of at least one STT hit in track.

STT efficiency with latest trunk version (2000 events) = 15.6%
(with requirement of at least one STT hit).
But I have no idea if the work on STT code is finished?
So I suppose that STT signal data should also be reproduced but the question is if to reproduce them together with TPC or later?

For TPC data reproduction in addition to tpc directory, genfit and RKTrackRep in GenfitTools/trackrep should be updated to rev 370 and updated tpc reco macros should be used.

And coming to the question of reconstruction of fully mixed events. I do not know how realistic is to expect to produce them by 5-7 August (the date when the GRID will be off), but without knowing when the code will be ready, there is still the question if the available DPM data should be used or they also needs to be reproduced. At least for TPC there were recent modifications related to clustering, which is the part of digitization, and I suppose TPC expert should define how important they are for DPM data which will be used for mixing.

Dima


Laughing
Re: Update on eta_c reconstruction with STT and TPC [message #12390 is a reply to message #12112] Fri, 29 July 2011 15:03 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
Hi,

I have several additional comments concerning eta_c reconstruction efficiency.

1. First comment is related to worse reconstruction efficiency for STT.
I used the rev. 12894, which includes the latest commit from Gianluigi available in svn from the 27.07 23:15.
It didn't include latest update from Lia, which only affected PndSecondaryTrackFinder.cxx.

For data which I produced yesterday with the trunk the mass distribution for eta_c and phi without any cuts look like:

index.php?t=getfile&id=6558&private=0
whereas for the data on grid

index.php?t=getfile&id=6557&private=0

So the left tail for eta_c distribution increased significantly, but it's still possible that it's a statistical fluctuation.

Especially looking at the final plots for new data:
index.php?t=getfile&id=6559&private=0
and grid data:
index.php?t=getfile&id=6560&private=0

they are very similar. So I could expect that with higher statistics, efficiency still could increase.

2. Second comment is related to ideal pid for tpc Riemann. I checked eta_c reconstruction with both pion and kaon hypothesis in riemann and results are very close (here the mass cuts are slightly tighter then for efficiency reported in previous post), i.e
Efficiency | sigma(eta_c) | sigma(phi)
pion:
27.04% | 11.74 MeV | 3.61 MeV
kaon:
27.29% | 11.94 MeV | 3.55 MeV

So using working ideal pid wouldn't change the results much.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12394 is a reply to message #12389] Fri, 29 July 2011 16:30 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *pv.infn.it
Quote:

Lia has also updated the run_reco_stt_evt.C and run_reco_stt_mix.C adding the task of the extrapolation to the Gem region.

Actually the task for the extrapolation was already there, I just set the additional flag to extrapolate to the GEM with the MC pdg code instead of the default one.
Dima, did you use this updated macro to produce your last plot? Could it be that the mc pdg code in the gem extrapolation worsen results? I would be very surprised (from my tests there was almost no change in efficiency).

Cheers,
Lia.
Re: Update on eta_c reconstruction with STT and TPC [message #12395 is a reply to message #12112] Fri, 29 July 2011 16:43 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
From: *fuw.edu.pl
Hi Lia,

yes I used the updated macro to produce the last results, i.e. with

PndSttMvdGemTracking * SttMvdGemTracking = new PndSttMvdGemTracking(0);
SttMvdGemTracking->SetPdgFromMC();
fRun->AddTask(SttMvdGemTracking);


So I will try to comment the line

SttMvdGemTracking->SetPdgFromMC();

and see how it will affect the efficiency, however it will take some time to run reconstruction on 2000 events (5-6 hours) and I suppose it's better to wait with grid reproduction of the signal with stt till I have the results.

But on the other hand I suppose that STT code on the grid can be updated till the latest version (if other modifications are not expected), since this point will affect only the switch in the macro.

Dima
Re: Update on eta_c reconstruction with STT and TPC [message #12397 is a reply to message #12395] Fri, 29 July 2011 17:04 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *pv.infn.it
Quote:

So I will try to comment the line

SttMvdGemTracking->SetPdgFromMC();

and see how it will affect the efficiency

Yes, I would appreciate this, thank you for testing it on a real channel.

Quote:

however it will take some time to run reconstruction on 2000 events (5-6 hours) and I suppose it's better to wait with grid reproduction of the signal with stt till I have the results.

Yes, I agree, it is better to see what happens first.

Quote:

But on the other hand I suppose that STT code on the grid can be updated till the latest version (if other modifications are not expected), since this point will affect only the switch in the macro.

I agree also on this: from STT side I can confirm we can update the stt and sttmvdtracking directories and we could wait for your results before updating the macros.

(For the secondary track finder I still have to work on it, so we cannot wait for it).

Cheers,
Lia.
Re: Update on eta_c reconstruction with STT and TPC [message #12398 is a reply to message #12395] Fri, 29 July 2011 19:02 Go to previous messageGo to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *57-82-r.retail.telecomitalia.it
Dima Melnychuk wrote on Fri, 29 July 2011 16:43


So I will try to comment the line

SttMvdGemTracking->SetPdgFromMC();

and see how it will affect the efficiency, however it will take some time to run reconstruction on 2000 events (5-6 hours) and I suppose it's better to wait with grid reproduction of the signal with stt till I have the results.



Considering the time scale we will do only a last grid production. Therefore, it is better to check everything before starting. The same fo the ideal pid of tpc, I have not well understood if this is finally fixed or not.
Previous Topic: Problem with secondary tracks
Next Topic: TpcRiemannTrackingTask - ideal PID now working - please rerun reconstruction
Goto Forum:
  


Current Time: Sun Nov 24 03:42:52 CET 2024

Total time taken to generate the page: 0.00943 seconds