// --------------------------------------------------------------------------
//
// Macro for standard transport simulation using UrQMD input and GEANT3
// CBM setup with MVD only
//
// V. Friese   06/02/2007
//
// --------------------------------------------------------------------------
    void transportTelescope(Int_t nEvents = 1)
//       void transportTelescope(Int_t nEvents = 3000)
//    void transportTelescope(Int_t nEvents = 51)
{
  // ========================================================================
  //          Adjust this part according to your requirements

  // Input file
  TString inDir   = gSystem->Getenv("VMCWORKDIR");

  TString workDir = inDir + "/macro/mvd/Data/";
  TString inFile  = inDir + "/input/urqmd.ftn14";
  
  // Number of events
  //Int_t   nEvents = 5000;

  // Output file name
  TString outDir  = "data";
  TString outFile = outDir + "telescope.box.root";

//   // Parameter file name
  TString parFile = outDir + "telescope.params.root";

  // Cave geometry
  TString caveGeom = "cave.geo";

  // Target geometry
  TString targetGeom = "target_au_250mu.geo";

  // Beam pipe geometry
  TString pipeGeom = "pipe_standard.geo";

  

  // MVD geometry
  
//     TString mvdGeom = "mvd/Beamtimesetupobgnum.root";
  TString mvdGeom = "mvd/MVD_012_Station_MA_Z_Digi_MISTRAL_150um.root";
  
  //TString mvdGeom = "beam/mvd_sensor.geo";
//  TString mvdGeom = "mvd/mvd_v08a.geo";

  // In general, the following parts need not be touched
  // ========================================================================




  // ----    Debug option   -------------------------------------------------
  gDebug = 0;
  // ------------------------------------------------------------------------


  // -----   Timer   --------------------------------------------------------
  TStopwatch timer;
  timer.Start();
  // ------------------------------------------------------------------------


  // ----  Load libraries   -------------------------------------------------
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libCbmBase");
  gSystem->Load("libCbmData");
  gSystem->Load("libCbmGenerators");
  gSystem->Load("libField");
  gSystem->Load("libGen");
  gSystem->Load("libPassive");
  //gSystem->Load("libEcal");
  //gSystem->Load("libKF");
  gSystem->Load("libMvd");
  gSystem->Load("libSts");
  //gSystem->Load("libLittrack");
  //gSystem->Load("libRich");
  //gSystem->Load("libTrd");
  //gSystem->Load("libTof");
  
  FairLogger *logger = FairLogger::GetLogger();
  logger->SetLogToScreen(kFALSE);
  
  // -----   Create simulation run   ----------------------------------------
  FairRunSim* fRun = new FairRunSim();
  fRun->SetName("TGeant3");              // Transport engine
  fRun->SetOutputFile(outFile);          // Output file
  FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
  // ------------------------------------------------------------------------
  

   

  // -----   Create media   -------------------------------------------------
  fRun->SetMaterials("media.geo");       // Materials
  // ------------------------------------------------------------------------


  // -----   Create geometry   ----------------------------------------------
  FairModule* cave= new CbmCave("CAVE");
  cave->SetGeometryFileName(caveGeom);
  fRun->AddModule(cave);

  FairModule* pipe= new CbmPipe("PIPE");
  pipe->SetGeometryFileName(pipeGeom);
  //fRun->AddModule(pipe);
  
//  FairModule* target= new CbmTarget("Target");
//  target->SetGeometryFileName(targetGeom);
//  fRun->AddModule(target);

//  FairModule* magnet= new CbmMagnet("MAGNET");
//  magnet->SetGeometryFileName(magnetGeom);
//  fRun->AddModule(magnet);
  
  FairDetector* mvd= new CbmMvd("MVD", kTRUE);
  mvd->SetGeometryFileName(mvdGeom); 
  
  fRun->AddModule(mvd);
  // ------------------------------------------------------------------------


  // ------------------------------------------------------------------------

  // -----   Create magnetic field   ----------------------------------------
  // Field Map Definition
 // --------------------
  // Constant Field
  CbmFieldConst *fMagField1=new CbmFieldConst();
  //fMagField1->SetField(0, 14 ,0 ); // values are in kG
  // MinX=-75, MinY=-40,MinZ=-12 ,MaxX=75, MaxY=40 ,MaxZ=124 );  // values are in cm
  //fMagField1->SetFieldRegion(-74, -39 ,-22 , 74, 39 , 160 );
  
  fMagField1->SetField(0, 0 ,0 ); // values are in kG
  fMagField1->SetFieldRegion(10, 10 ,0 , 11, 11 , 1 );
  
  fRun->SetField(fMagField1);
  
  // Store the partile trajectories

  fRun->SetStoreTraj(kTRUE);
  // ------------------------------------------------------------------------



  // -----   Create PrimaryGenerator   --------------------------------------
  FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
//  FairUrqmdGenerator*  urqmdGen = new FairUrqmdGenerator(inFile);
//   primGen->AddGenerator(urqmdGen);
  
  
  FairBoxGenerator* boxGenerator=new FairBoxGenerator(211,100);
// FairBoxGenerator* boxGenerator=new FairBoxGenerator(211,300);
  
   boxGenerator->SetXYZ(0,0,0);
  boxGenerator->SetPRange(120,121);
  boxGenerator->SetPhiRange(0,360);
  boxGenerator->SetThetaRange(0,360);
  //boxGenerator->SetBoxXYZ(-2,-2,2,2,-20)
   primGen->AddGenerator(boxGenerator);
//   primGen->SetTarget(0,1.); // set Target to position, with thickness

  primGen->SetEventTime(0); // starting time of the first event
  primGen->SetEventTimeInterval(15000,15000); // timewindow in which a new event appears
  
  fRun->SetGenerator(primGen);
  // ------------------------------------------------------------------------

    // -----   Create Tasks  --------------------------------------------------
  
  CbmMvdSensorBufferTask* buffer1 = new CbmMvdSensorBufferTask("MvdSensorBuffer1");
  buffer1-> SetWriteOutput(true);// Set true to enable writing Bufferoutput to file 
 
  
  fRun->AddTask(buffer1);


  // -----   Run initialisation   -------------------------------------------
  fRun->Init();
  
  // ------------------------------------------------------------------------

  
  
  // -----   Runtime database   ---------------------------------------------
  CbmFieldPar* fieldPar = (CbmFieldPar*) rtdb->getContainer("CbmFieldPar");
  fieldPar->SetParameters(fMagField1);
  fieldPar->setChanged();
  fieldPar->setInputVersion(fRun->GetRunId(),1);
  Bool_t kParameterMerged = kTRUE;
  FairParRootFileIo* parOut = new FairParRootFileIo(kParameterMerged);
  parOut->open(parFile.Data());
  rtdb->setOutput(parOut);
  rtdb->saveOutput();
  rtdb->print();
  // ------------------------------------------------------------------------
  

 
  // -----   Start run   ----------------------------------------------------
  fRun->Run(nEvents);
  // ------------------------------------------------------------------------
  
  
  // -----   Finish   -------------------------------------------------------
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  cout << endl << endl;
  cout << "Macro finished succesfully." << endl;
  cout << "Thank you for using transportTelescope.C." << endl;
  cout << "Output file is "    << outFile << endl;
  cout << "Parameter file is " << parFile << endl;
  cout << "Real time " << rtime << " s, CPU time " << ctime 
       << "s" << endl << endl;
  // ------------------------------------------------------------------------

  cout << " Test passed" << endl;
  cout << " All ok " << endl;
 // new TBrowser();
}