//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Implementation of class TpcLaserTask
//      see TpcLaserTask.h for details
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Felix Boehmer        TUM            (original author)
//
//
//-----------------------------------------------------------


// This Class' Header ---------------------------------------
#include "TpcSimpleEvtGen.h"

// Collaborating Class Headers --------
#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"
#include "TClonesArray.h"
#include "TpcPrimaryCluster.h"
#include "TVector3.h"
#include "TError.h"
#include "TpcDigiPar.h"
#include <iostream>
#include <string>
#include "TRandom3.h"


TpcSimpleEvtGen::TpcSimpleEvtGen()
  : FairTask("TPC Simple Event Gen")
{
  fmode = "iron55";
  fpos = TVector3(0.,0.,0.);
}

TpcSimpleEvtGen::~TpcSimpleEvtGen()
{
  ;  
}

InitStatus
TpcSimpleEvtGen::Init()
{
  //Get ROOT Manager
  FairRootManager* ioman= FairRootManager::Instance();

  if(ioman==0)
    {
      Error("TpcSimpleEvtGen::Init","RootManager not instantiated!");
      return kERROR;
    }
  // Get input collection
  fprimArray=(TClonesArray*) ioman->GetObject("TpcPrimaryCluster");
  
  if(fprimArray==0)
    {
      Error("TpcSimpleEvtGen::Init","PrimaryElectron-array not found!");
      return kERROR;
    }

  //clear input array of any real "physical" hits
  fprimArray->Delete();

  std::cerr<<"\n\nTpcSimpleEvtGen: running in mode *** "<<fmode.c_str()
	   <<" ***\n"<<std::endl;

   return kSUCCESS;
}



void 
TpcSimpleEvtGen::SetParContainers() {

  std::cout<<"TpcSimpleEvtGen::SetParContainers"<<std::endl;
  std::cout.flush();

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

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

  // Get Tpc digitisation parameter container
  fpar = (TpcDigiPar*) db->getContainer("TpcDigiPar");
  if (! fpar ) Fatal("SetParContainers", "TpcDigiPar not found");
}



void
TpcSimpleEvtGen::Exec(Option_t* opt)
{  
  if(fmode.compare("iron55") == 0) {
    Int_t nPrim = fprimArray->GetEntriesFast();
    //if(nPrim!=0) 
    //  Fatal("TpcSimpleEvtGen::Exec","Input-Array not empty!");
      
    std::cerr<<"\n\nTpcSimpleEvtGen: running in mode *** "<<fmode.c_str()
	     <<" ***\n"<<std::endl;
    std::cerr<<"TpcSimpleEvtGen: Position: ("<<fpos.X()<<", "
	     <<fpos.Y()<<", "<<fpos.Z()<<")"<<std::endl;
    
    TRandom3 randm;
    randm.SetSeed(1);
    
    for(unsigned int n=0; n<20; ++n) {
      TVector3 postemp = fpos+TVector3(randm.Uniform(-0.01,0.01),
				       randm.Uniform(-0.01,0.01),
				       randm.Uniform(-0.01,0.01));
				      
      nPrim = fprimArray->GetEntriesFast();
      TpcPrimaryCluster* clus =new((*fprimArray)[nPrim]) TpcPrimaryCluster(n,1,postemp,1,1);
      clus->setIndex(nPrim);
      //new((*fprimArray)[nPrim]) TpcPrimaryCluster(1,293,postemp,1,1);
    }
  }
}   
 
    
ClassImp(TpcSimpleEvtGen)