Hi,
I've committed changes to
pc_app/dst/{fsDataCut.cxx,BrFsDstAna.cxx,BrFsDstAna.h}.
The changes were necessary to support the efficiency calculation for FS
spectrometer when
producing ntuples from dsts.
BrFsDstAna reads efficiency maps from
/brahms/u/ufstasze/effic/Effic5361_5972.root and
for each track calculates effective efficiencies (ffs and fs) and writes
to the ntuple.
Djamel tried this new stuff yesterday and claims that the spectra look
nicer.
By default the efficiency calculation is ignored. If you want to include
the efficiency calculation
then you should specified option -M 1 in the fsDataCut command lime.
By the way, we have also in hand the BrDcEnhancementModule. I tried it
and I got
91% in T5 at 3 deg and 97% in T5 at 12 deg. for runs with the T5/II
off!!!.
If any one by any reason wants to redo the global tracking I suggest to
add
this new module (specially for high field runs at 12 deg where the
statistics for bfs tracks is very important.)
The setup for the BrDcEnhanceModule should come just after
BrBfsTrackingModule,
and the BrBfsTrackingModule should have specified he following options
bfs->SetKeepMatched(kTRUE);
bfs->SetExtendD3Match(kTRUE);
bfs->SetExtendD4Match(kFALSE);
bfs->SetCleanUp(BrBfsTrackingModule::kCUT2, kTRUE);
bfs->SetCleanUp(BrBfsTrackingModule::kCUT3, kFALSE);
The setup for BrDcEn.... is:
gSystem->Load("/brahms/u/ufstasze/brahms/brahms_app/rk_app/bfs/dcenh/BrDcEnhancementModule.so");
BrDcEnhancementModule* enhT5 = new BrDcEnhancementModule("T5", "T5");
mainModule->AddModule(enhT5);
enhT5->SetRunNumber(runNo);
enhT5->SetVerbosity(verboseOption->GetValue());
enhT5->GetCombineT4T5()->SetFiducialCutDx(1.);
enhT5->GetCombineT4T5()->SetFiducialCutDy(1.);
enhT5->GetCombineT4T5()->SetSigmaCut(3.5);
enhT5->GetCombineT4T5()->SetMatchDalyOffset(0.0004);
enhT5->GetCombineT4T5()->SetMatchDangOffset(0.0011);
enhT5->GetCombineT4T5()->SetMatchDyOffset(0.075);
enhT5->GetCombineT4T5()->SetSigmaDaly(0.0014);
enhT5->GetCombineT4T5()->SetSigmaDang(0.0014);
enhT5->GetCombineT4T5()->SetSigmaDy(0.2);
That's all.
Regards Pawel.
--
----------------------------------------------------------------------
| Pawel Staszel |
| Niels Bohr Institute Tb 8 email: staszel@nbi.dk |
| Blegdamsvej 17 phone: (+45) 35 32 53 51 |
| København Ø FAX: (+45) 35 32 50 16 |
| Danmark |
----------------------------------------------------------------------
This archive was generated by hypermail 2b30 : Wed Jul 10 2002 - 08:46:05 EDT