PndMasterRecoTask Class Reference

The default reconstruction tasks. More...

#include <PndMasterRecoTask.h>

Inheritance diagram for PndMasterRecoTask:

Classes
struct	recoTaskList
	struct of the task list This struct avoids to retrieve tasks using the integer value (misleading), but uses an easier scheme. The enum is "k" + the class name, i.e. "kPndTrkTracking2" for the class PndTracking2. The function PrintTaskList() can be used to check the list of the tasks and the corresponding number. More...

Public Member Functions
	PndMasterRecoTask (TString fOptions="")
	Default constructor with options This string can be: "" -> default settings full setup "day1" -> Setup for day1 experiment: no GEM, FTS1234, NO DISC, NO RICH "gem" (added to "day1") -> Setup for day1 experiment with 3 GEM planes "fts1256" (added to "day1") -> Setup for day1 experiment with FTS1256 insted of FTS1234 "mulitkalman" -> Run kalman tasks with all 5 hypothesis for a TrackCand Example: "day1+gem+fts1256" means day1 setup + GEM planes + fst1256. More...
virtual	~PndMasterRecoTask ()
	Destructor. More...
void	PrintTaskList ()
	It prints the list of tasks. More...
virtual void	SetPersistency (Bool_t pers=kTRUE)
	Set the persistency of all the tasks This function has to be implemented for each master task, since different tasks use different functions to change the persistency of the output TClonesArrays. More...
void	SetVerbose (Int_t iVerbose=1)
	Set the Verbosity to all the tasks. More...
void	SetVerbose (Int_t nTask, Int_t iVerbose)
	Set the Verbosity to a single task. More...
FairTask *	GetTask (Int_t nTask)
	Return the pointer to a single task. More...
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
virtual void	Exec (Option_t *opt)
Bool_t	GetPersistency ()

Protected Attributes
std::vector< PndPersistencyTask * >	fBranchTasks
std::vector< FairTask * >	fStandardTasks
std::map< PndPersistencyTask *, bool >	fFixedPersistency

Private Attributes
struct PndMasterRecoTask::recoTaskList	reco
TString	fOptions
	Options parsed to the reconstruction. More...

Detailed Description

The default reconstruction tasks.

Master Reconstruction Task Class

This class includes all the reconstruction tasks which need to be used in the default reconstruction macros.

Remarks

If you find some obsolete task which needs to be changed, contact the computing coordinator.

Author

Stefano Spataro spata.nosp@m.ro@t.nosp@m.o.inf.nosp@m.n.it, Torino University

Version

1.0

Date

Feb 1, 2016

Definition at line 20 of file PndMasterRecoTask.h.

Constructor & Destructor Documentation

PndMasterRecoTask::PndMasterRecoTask

(

TString

options = ""

)

Default constructor with options This string can be: "" -> default settings full setup "day1" -> Setup for day1 experiment: no GEM, FTS1234, NO DISC, NO RICH "gem" (added to "day1") -> Setup for day1 experiment with 3 GEM planes "fts1256" (added to "day1") -> Setup for day1 experiment with FTS1256 insted of FTS1234 "mulitkalman" -> Run kalman tasks with all 5 hypothesis for a TrackCand Example: "day1+gem+fts1256" means day1 setup + GEM planes + fst1256.

Default Constructor.

Reconstruction task list

Here all the reconstruction tasks are added to the task, with the standard settings. A check is done after each task if the tasklist enum is broken or not. At the end the event counter is added (each 100 events), and the verbosity is set to 0 to all the tasks (it can be changed afterwards with SetVerbosity() functions.

Definition at line 31 of file PndMasterRecoTask.cxx.

References Add(), PndIdealTrackFinder::AddBranchName(), PndTrackSmearTask::AddInputTrackBranch(), PndMasterTask::fBranchTasks, PndMasterTask::fFixedPersistency, fOptions, PndTrkTracking2::NoMvdAloneTracking(), recoKalman, PndRecoMultiKalmanTask::SetBusyCut(), PndRecoKalmanTask::SetBusyCut(), PndRecoMultiKalmanTask::SetFitHypotheses(), PndTrkTracking2::SetInputBranchName(), PndMissingPzCleanerTask::SetInputTrackBranch(), PndIdealTrackFinder::SetOutputBranchName(), PndIdealTrackFinder::SetRelativeMomentumSmearing(), PndMissingPzCleanerTask::SetRemoveTrack(), PndRecoKalmanTask::SetTrackInBranchName(), PndRecoKalmanTask2::SetTrackInBranchName(), PndRecoMultiKalmanTask2::SetTrackInBranchName(), PndRecoMultiKalmanTask::SetTrackInBranchName(), PndIdealTrackFinder::SetTrackingEfficiency(), PndRecoMultiKalmanTask2::SetTrackOutBranchName(), PndRecoKalmanTask2::SetTrackOutBranchName(), PndRecoKalmanTask::SetTrackOutBranchName(), PndRecoMultiKalmanTask::SetTrackOutBranchName(), PndRecoMultiKalmanTask::SetTrackRep(), PndRecoKalmanTask::SetTrackRep(), PndIdealTrackFinder::SetTrackSelector(), PndMasterTask::SetVerbose(), PndIdealTrackFinder::SetVertexSmearing(), SttMvdGemTracking, trackFts, tracking, and TString.

                                                    :
  PndMasterTask("Master Reconstruction Task"), fOptions(options)
{

  std::array<TString,5> hypoName= {"Electron","Muon","Pion","Kaon","Proton"};
  // -----  MVD + STT Pattern Recognition -----------------------------------
  //  use the constructor with input :
  //      printout flag (int) , plotting flag (bool), MC comparison flag (bool), SciTil.
  PndTrkTracking2* tracking = NULL;
  fBranchTasks.push_back(tracking = new PndTrkTracking2(0, false, false, true));
  tracking->SetInputBranchName("STTHit", "MVDHitsPixel", "MVDHitsStrip");
  // tracking->SetInputBranchName("STTHitMix","MVDHitsPixelMix","MVDHitsStripMix");
  //  don't do the Pattern Recognition second part, starting from the Mvd;
  tracking->NoMvdAloneTracking();
  // do Cleanup only when there is Mixing;
  // tracking->Cleanup();
  fFixedPersistency[tracking] = kFALSE;

  // ----- MVD + STT + GEM Pattern Recognition --------------
  if ((!fOptions.Contains("nogem")) && (!fOptions.Contains("gem0"))) {
    PndSttMvdGemTracking *SttMvdGemTracking = NULL;
    fBranchTasks.push_back(SttMvdGemTracking = new PndSttMvdGemTracking(0));
    fFixedPersistency[SttMvdGemTracking] = kFALSE;
  }

  // ----- Barrel Kalman Task     ----------------------------
  if (!fOptions.Contains("multikalman"))
  { //no multikalman
    PndRecoKalmanTask* recoKalman = NULL;
    fBranchTasks.push_back(recoKalman = new PndRecoKalmanTask());
    if ((!fOptions.Contains("nogem")) && (!fOptions.Contains("gem0"))) {
      recoKalman->SetTrackInBranchName("SttMvdGemTrack");
      //      recoKalman->SetTrackInIDBranchName("SttMvdGemTrackID");
      recoKalman->SetTrackOutBranchName("SttMvdGemGenTrack"); // PID appendix follows;
    } else {
      recoKalman->SetTrackInBranchName("SttMvdTrack");
      //      recoKalman->SetTrackInIDBranchName("SttMvdTrackID");
      recoKalman->SetTrackOutBranchName("SttMvdGenTrack"); // PID appendix follows
    }
    recoKalman->SetBusyCut(50); // CHECK to be tuned
    //recoKalman->SetIdealHyp(kTRUE);
    //recoKalman->SetNumIterations(3);
    recoKalman->SetTrackRep(0); // 0 Geane (default), 1 RK
    //recoKalman->SetPropagateToIP(kFALSE);
  }
  else //yes, multikalman
  {
    PndRecoMultiKalmanTask* recoKalman = NULL;
    fBranchTasks.push_back(recoKalman = new PndRecoMultiKalmanTask());
    recoKalman->SetFitHypotheses(fOptions);  //Walter added
    if ((!fOptions.Contains("nogem")) && (!fOptions.Contains("gem0"))) {
      //if ((!fOptions.Contains("day1")) || (fOptions.Contains("gem"))) {
      recoKalman->SetTrackInBranchName("SttMvdGemTrack");
      //      recoKalman->SetTrackInIDBranchName("SttMvdGemTrackID");
      recoKalman->SetTrackOutBranchName("SttMvdGemGenTrack"); // PID appendix follows
    } else {
      recoKalman->SetTrackInBranchName("SttMvdTrack");
      //      recoKalman->SetTrackInIDBranchName("SttMvdTrackID");
      recoKalman->SetTrackOutBranchName("SttMvdGenTrack"); // PID appendix follows
    }
    recoKalman->SetBusyCut(50); // CHECK to be tuned
    //recoKalman->SetIdealHyp(kTRUE);
    //recoKalman->SetNumIterations(3);
    recoKalman->SetTrackRep(0); // 0 Geane (default), 1 RK
    //recoKalman->SetPropagateToIP(kFALSE);
  }

  if (fOptions.Contains("filtered")) {
    PndMissingPzCleanerTask* cleaner = NULL;
    fBranchTasks.push_back(cleaner = new PndMissingPzCleanerTask());
    if ((!fOptions.Contains("nogem")) && (!fOptions.Contains("gem0"))) {
      cleaner->SetInputTrackBranch("SttMvdGemGenTrack"); // _filtered appendix follows
    } else {
      cleaner->SetInputTrackBranch("SttMvdGenTrack");// _filtered appendix follows
    }
    cleaner->SetRemoveTrack(kTRUE);
  }


  // -----  FTS Ideal Tracking    ----------------------------
  TString ftstrkname="FtsIdealTrack", ftstrkgenname="FtsIdealGenTrack";
  if(fOptions.Contains("ftsca")) {
    ftstrkname="FtsCaTrack";
    ftstrkgenname="FtsCaGenTrack";
  }
  PndIdealTrackFinder* trackFts = NULL;
  fBranchTasks.push_back(trackFts = new PndIdealTrackFinder());
  trackFts->SetTrackSelector("FtsTrackFunctor");
  trackFts->AddBranchName("FTSHit");
  trackFts->AddBranchName("MVDHitsPixel");
  trackFts->AddBranchName("MVDHitsStrip");
  trackFts->SetRelativeMomentumSmearing(0.05);
  trackFts->SetVertexSmearing(0.05, 0.05, 0.05);
  trackFts->SetTrackingEfficiency(1.);
  trackFts->SetOutputBranchName("FtsIdealTrack");
  fFixedPersistency[trackFts] = kFALSE;

  // ----- Forward Kalman Task     ---------------------------
  if (!fOptions.Contains("multikalman")) { // no multikalman
    if(fOptions.Contains("genfit2")) {
      PndRecoKalmanTask2* recoKalman = NULL;
      fBranchTasks.push_back(recoKalman = new PndRecoKalmanTask2());
      recoKalman->SetTrackInBranchName(ftstrkname);
      recoKalman->SetTrackOutBranchName(ftstrkgenname);
    } else {
      PndRecoKalmanTask* recoKalman = NULL;
      fBranchTasks.push_back(recoKalman = new PndRecoKalmanTask());

      recoKalman->SetTrackInBranchName(ftstrkname);
      recoKalman->SetTrackOutBranchName(ftstrkgenname);
    }
  } else { //yes, multikalman!
    if(fOptions.Contains("genfit2")) {
      PndRecoMultiKalmanTask2* recoKalman = NULL;
      fBranchTasks.push_back(recoKalman = new PndRecoMultiKalmanTask2());
      recoKalman->SetTrackInBranchName(ftstrkname);
      recoKalman->SetTrackOutBranchName(ftstrkgenname);
    } else {
      PndRecoMultiKalmanTask* recoKalman = NULL;
      fBranchTasks.push_back(recoKalman = new PndRecoMultiKalmanTask());
      recoKalman->SetTrackInBranchName(ftstrkname);
      recoKalman->SetTrackOutBranchName(ftstrkgenname);
    }
  }

  // --- Online Emulator smearing
  if(fOptions.Contains("fakeonline"))
  {
    PndTrackSmearTask* smearer=NULL;
    fBranchTasks.push_back(smearer=new PndTrackSmearTask());
    TString barrelbranchname="SttMvdGemGenTrack";
    TString ftsbranchname="FtsIdealGenTrack";
    if (fOptions.Contains("nogem")||fOptions.Contains("gem0")) {
      barrelbranchname="SttMvdGenTrack";
    }

    if(fOptions.Contains("filtered")) barrelbranchname+="_filtered";

    if(fOptions.Contains("multikalman")) { //yes, multikalman
      for(int k=0; k<5; ++k) {
        smearer->AddInputTrackBranch(Form("%s%s",barrelbranchname.Data(),hypoName[k].Data()));
        smearer->AddInputTrackBranch(Form("%s%s",ftsbranchname.Data(),hypoName[k].Data()));
      }
    } else { // no multikalman
      smearer->AddInputTrackBranch(barrelbranchname);
      smearer->AddInputTrackBranch(ftsbranchname);
    }
  }

  std::for_each(fBranchTasks.begin(), fBranchTasks.end(),
  [this](const FairTask* task) {
    Add((TTask*)task);
  });

  SetVerbose(0);
}

PndMasterRecoTask::~PndMasterRecoTask ( )

virtual

Destructor.

Definition at line 190 of file PndMasterRecoTask.cxx.

                                      {
}

Member Function Documentation

void PndBlackBoxTask::Exec ( Option_t *  opt )

virtualinherited

Reimplemented in PndMvdQATask, and PndBufferTestTask.

Definition at line 50 of file PndBlackBoxTask.cxx.

{
  return;
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

FairTask * PndMasterTask::GetTask ( Int_t  nTask )

inherited

Return the pointer to a single task.

This function allows to retrieve the pointer of a task inside the task list, when you want to use some particular setter or change someting in the task.

Parameters

nTask

Index of the task you want to modify. Instead of the integer you can use the enum "k" + task name (i.e. kMvdDigiTask, kSttHitProducerRealFast, etc...)

Returns

A pointer to the corresponding task

Retrieve the pointer to a Task in the list

Definition at line 75 of file PndMasterTask.cxx.

{
  TList* thistasks = this->GetListOfTasks();

  return ((FairTask*)thistasks->At(nTask));
}

InitStatus PndBlackBoxTask::Init ( )

virtualinherited

Reimplemented in PndMvdQATask, and PndBufferTestTask.

Definition at line 43 of file PndBlackBoxTask.cxx.

{
  return kSUCCESS;
}

void PndMasterTask::PrintTaskList ( )

inherited

It prints the list of tasks.

Print the list of the task which are included in the list

This function print the task number, its title and its name. This can be important when it is needed to modify a particular task and the task number is needed. In any case, the use of enum should be preferred.

Remarks

The task names are not well defined in the classes, some of them are missing. We should define a better scheme.

Definition at line 24 of file PndMasterTask.cxx.

References counter, and next.

{
  TList* thistasks = this->GetListOfTasks();

  Int_t counter = 0;
  TIter next(thistasks->MakeIterator());
  while (FairTask *task = (FairTask*)next())
    {
      if (counter < 10)
        std::cout << "Task  #" << counter << "\tTitle: " << task->GetTitle() << "\tName: " << task->GetName() << std::endl;
      else
        std::cout << "Task #" << counter << "\tTitle: " << task->GetTitle() << "\tName: " << task->GetName() << std::endl;
      counter++;
    }

  return;
}

InitStatus PndBlackBoxTask::ReInit ( )

virtualinherited

Definition at line 25 of file PndBlackBoxTask.cxx.

{
  return kSUCCESS;
}

void PndBlackBoxTask::SetParContainers ( )

virtualinherited

Reimplemented in PndMvdQATask.

Definition at line 21 of file PndBlackBoxTask.cxx.

{
}

void PndMasterTask::SetPersistency ( Bool_t  pers = kTRUE )

virtualinherited

Set the persistency of all the tasks This function has to be implemented for each master task, since different tasks use different functions to change the persistency of the output TClonesArrays.

Parameters

fPersistency

Persistency level: 0 no TCA, 1 all TCA

Set the Persistency of all the tasks in the same way

Reimplemented in PndMasterLocalRecoTask, PndMasterMultiPidTask, PndMasterPidTask, PndMasterDigiTask, PndMasterRecoIdealTask, and PndMasterSimTask.

Definition at line 66 of file PndMasterTask.cxx.

References PndMasterTask::fBranchTasks, and PndMasterTask::fFixedPersistency.

Referenced by PndMasterRunAna::AddDigiOnlyTasks(), and PndMasterRunAna::AddRecoTasks().

{
          std::for_each(fBranchTasks.begin(), fBranchTasks.end(), [pers](PndPersistencyTask* task){task->SetPersistency(pers);});
          std::for_each(fFixedPersistency.begin(), fFixedPersistency.end(), [](std::pair<PndPersistencyTask*, bool> task){task.first->SetPersistency(task.second);});

          return;
}

void PndMasterTask::SetVerbose ( Int_t  iVerbose = 1 )

inherited

Set the Verbosity to all the tasks.

Parameters

iVerbose

Verbosity level: 0 no messages, the higher more messages

Set the Verbosity of all the tasks at the same number

Definition at line 43 of file PndMasterTask.cxx.

References next.

Referenced by PndMasterDigiOnlyTask::PndMasterDigiOnlyTask(), PndMasterDigiTask::PndMasterDigiTask(), PndMasterLocalRecoTask::PndMasterLocalRecoTask(), PndMasterMultiPidTask::PndMasterMultiPidTask(), PndMasterPidTask::PndMasterPidTask(), PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask(), PndMasterSimTask::PndMasterSimTask(), and PndMasterTask::SetVerbose().

{
  TList* thistasks = this->GetListOfTasks();

  TIter next(thistasks->MakeIterator());
  while (FairTask *task = (FairTask*)next())
    {
      task->SetVerbose(iVerbose);
    }

  return;
}

void PndMasterTask::SetVerbose ( Int_t  nTask, Int_t  iVerbose  )

inherited

Set the Verbosity to a single task.

This function has to be used when you want to change the verbosity to a single task

Parameters

nTask	Index of the task you want to modify. Instead of the integer you can use the enum k + task name (i.e. kMvdDigiTask, hSttHitProducerRealFast
iVerbose	Verbosity. 0 -> no messages; the higher -> more messages

Set the Verbosity of a single task

Definition at line 57 of file PndMasterTask.cxx.

References PndMasterTask::SetVerbose().

{
  TList* thistasks = this->GetListOfTasks();
  ((FairTask*)thistasks->At(nTask))->SetVerbose(iVerbose);

  return;
}

Member Data Documentation

std::vector<PndPersistencyTask*> PndMasterTask::fBranchTasks

protectedinherited

Definition at line 80 of file PndMasterTask.h.

Referenced by PndMasterDigiOnlyTask::PndMasterDigiOnlyTask(), PndMasterRecoTask(), and PndMasterTask::SetPersistency().

std::map<PndPersistencyTask*, bool> PndMasterTask::fFixedPersistency

protectedinherited

Definition at line 82 of file PndMasterTask.h.

Referenced by PndMasterDigiOnlyTask::PndMasterDigiOnlyTask(), PndMasterRecoTask(), and PndMasterTask::SetPersistency().

TString PndMasterRecoTask::fOptions

private

Options parsed to the reconstruction.

Definition at line 58 of file PndMasterRecoTask.h.

Referenced by PndMasterRecoTask().

std::vector<FairTask*> PndMasterTask::fStandardTasks

protectedinherited

Definition at line 81 of file PndMasterTask.h.

Referenced by PndMasterPidTask::PndMasterPidTask().

struct PndMasterRecoTask::recoTaskList PndMasterRecoTask::reco

private

---

The documentation for this class was generated from the following files:

• PndMasterRecoTask.h
• PndMasterRecoTask.cxx