CbmMCMatchLoaderTask.cxx

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CbmMCMatchLoaderTask.cxx

// -------------------------------------------------------------------------
// -----                CbmMCMatchLoaderTask source file             -----
// -----                  Created 18/07/08  by T.Stockmanns        -----
// -------------------------------------------------------------------------
#include "CbmMCMatchLoaderTask.h"

#include "CbmMCMatch.h"
#include "CbmDetectorList.h"

// framework includes
#include "FairRootManager.h"

// Root includes
#include "TClonesArray.h"

// libc includes
#include <iostream>
using std::cout;
using std::endl;

// -----   Default constructor   -------------------------------------------
CbmMCMatchLoaderTask::CbmMCMatchLoaderTask() 
  : FairTask("Creates CbmMCMatch"), 
    fMCLink(NULL),
    fEventNr(0),
    fMCMatch(NULL)
{
}
// -------------------------------------------------------------------------


// -----   Destructor   ----------------------------------------------------
CbmMCMatchLoaderTask::~CbmMCMatchLoaderTask()
{
}

// -----   Public method Init   --------------------------------------------
InitStatus CbmMCMatchLoaderTask::Init()
{
  fMCMatch = new CbmMCMatch("CbmMCMatch", "CbmMCMatch");

//  fMCMatch->InitStage(kMCTrack, "", "MCTrack");
  fMCMatch->InitStage(kStsPoint, "", "StsPoint");
  fMCMatch->InitStage(kStsDigi, "", "StsDigi");
  fMCMatch->InitStage(kStsCluster, "", "StsCluster");
  fMCMatch->InitStage(kStsHit, "", "StsHit");


  FairRootManager* ioman = FairRootManager::Instance();
    if (!ioman) {
      cout << "-E- CbmMCMatchLoaderTask::Init: "
          << "RootManager not instantiated!" << endl;
      return kFATAL;
    }

  fMCLink = (TClonesArray*)ioman->GetObject("MCLink");
  ioman->Register("MCMatch", "MCMatch", fMCMatch, kFALSE);

  return kSUCCESS;
}


// -------------------------------------------------------------------------
void CbmMCMatchLoaderTask::SetParContainers()
{
  // Get Base Container
//  FairRun* ana = FairRun::Instance();
//  FairRuntimeDb* rtdb=ana->GetRuntimeDb();

}


// -----   Public method Exec   --------------------------------------------
void CbmMCMatchLoaderTask::Exec(Option_t* opt)
{
/*  for (int trackIndex = 0; trackIndex < fTrack->GetEntriesFast(); trackIndex++){
    CbmTrack* myTrack = (CbmTrack*)fTrack->At(trackIndex);
    fMCMatch->AddElement(kTrack, trackIndex, kTrackCand, myTrack->GetRefIndex());
  }
*/

  if (!fMCLink) Fatal("Exec", "No fMCLink");
//  fMCLinkDet->Delete();
//  fMCLinkHit->Delete();
  fMCMatch->ClearMCList();

  fMCMatch->LoadInMCLists(fMCLink);
  fMCMatch->CreateArtificialStage(kMCTrack, "", "");

  fMCMatch->Print();
  cout << endl;
}

void CbmMCMatchLoaderTask::Finish()
{
}


ClassImp(CbmMCMatchLoaderTask);