PndGemEmpty.cxx

Go to the documentation of this file.

/* $Id: */

// -------------------------------------------------------------------------
// -----                    PndGemEmpty source file         -----
// -----                  Created 15/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------

#include "PndGemEmpty.h"

#include "PndGemCluster.h"
#include "PndGemDigi.h"
#include "PndGemDigiPar.h"
#include "PndGemSensor.h"
#include "PndGemStation.h"

#include "PndDetectorList.h"

#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"

#include "TClonesArray.h"
#include "TMath.h"
#include "TRandom2.h"

#include <iomanip>

using std::cout;
using std::cerr;
using std::endl;
using std::flush;
using std::fixed;
using std::right;
using std::left;
using std::setw;
using std::setprecision;
using std::set;
using std::map;

// -----   Default constructor   ------------------------------------------
PndGemEmpty::PndGemEmpty()
  : FairTask("GEMEmpty", 0),
    fFunctor     (NULL),
    fInBranch    (NULL),
    fOutBranch   (NULL),
    fInBranchName("GEMDigi")
{
}
// -------------------------------------------------------------------------



// -----   Standard constructor   ------------------------------------------
PndGemEmpty::PndGemEmpty(Int_t iVerbose) 
  : FairTask("GEMEmpty", iVerbose),
    fFunctor     (NULL),
    fInBranch    (NULL),
    fOutBranch   (NULL),
    fInBranchName("GEMDigi") 
{
}
// -------------------------------------------------------------------------



// -----   Constructor with name   -----------------------------------------
PndGemEmpty::PndGemEmpty(const char* name, Int_t iVerbose) 
  : FairTask(name, iVerbose),
    fFunctor     (NULL),
    fInBranch    (NULL),
    fOutBranch   (NULL),
    fInBranchName("GEMDigi") {
}
// -------------------------------------------------------------------------



// -----   Destructor   ----------------------------------------------------
PndGemEmpty::~PndGemEmpty() {
  if ( 0!=fOutBranch ) {
    fOutBranch->Delete();
    delete fOutBranch;
  }
  if ( 0!=fFunctor ) delete fFunctor;
}
// -------------------------------------------------------------------------



// -----   Public method Exec   --------------------------------------------
void PndGemEmpty::Exec(Option_t*) {
  //  cout << endl << "======== PndGemEmpty::Exec(Event = " << fTNofEvents << " ) ====================" << endl;

  if ( fInBranchName.Contains("Sorted") ) {
    fInBranch->Clear();
    fInBranch = FairRootManager::Instance()->GetData(fInBranchName, fFunctor, 20.); //FairRootManager::Instance()->GetEventTime() +
  }
  else
    fInBranch = (TClonesArray*)FairRootManager::Instance()->GetObject(fInBranchName);

  fOutBranch = FairRootManager::Instance()->GetTClonesArray("GEMCluster");
  if ( ! fOutBranch ) Fatal("Exec", "No GEM Cluster Array");
  fOutBranch->Delete();

  std::vector<Int_t> clusterRefs;
  //PndGemCluster* cluster = NULL;
  for ( Int_t iout = 0 ; iout < fInBranch->GetEntries()/14 ; iout++ ) {
    new ((*fOutBranch)[iout]) PndGemCluster(585180824,
                                                      iout,
                                                      0,
                                                      iout*2.,
                                                      100,
                                                      0,
                                                      clusterRefs);
    
  }
  
  std::cout << "-I- PndGemEmpty::Exec. There are " << fInBranch->GetEntries() << " objects of " << fInBranchName.Data() << std::endl;
  
}
// -------------------------------------------------------------------------
    

// -----   Private method SetParContainers   -------------------------------
void PndGemEmpty::SetParContainers() {

  // Get run and runtime database
  FairRunAna* run = FairRunAna::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");

  FairRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("SetParContainers", "No runtime database");

}
// -------------------------------------------------------------------------




// -----   Private method Init   -------------------------------------------
InitStatus PndGemEmpty::Init() {

  // Get input array
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) Fatal("Init", "No FairRootManager");
  fInBranch = (TClonesArray*) ioman->GetObject(fInBranchName.Data());
  fOutBranch = ioman->Register("GEMCluster", "PndGemCluster", "PndGEM", kTRUE);

  fFunctor = new TimeGap();

  cout << "-I- " << fName.Data() << "::Init(). Initialization succesfull." << endl;  
  
  return kSUCCESS;
}
// -------------------------------------------------------------------------




// -----   Private method ReInit   -----------------------------------------
InitStatus PndGemEmpty::ReInit() {

  // Create sectorwise digi sets
  //  MakeSets();

  return kSUCCESS;
}
// -------------------------------------------------------------------------



// -----   Public method Finish   ------------------------------------------
void PndGemEmpty::Finish() {
  if ( fOutBranch ) fOutBranch->Clear();

  cout << "-------------------- " << fName.Data() << " : Summary -----------------------" << endl;
  cout << "-- DONE -- " << endl;
  cout << "---------------------------------------------------------------------" << endl; 

}
// -------------------------------------------------------------------------


ClassImp(PndGemEmpty)