timeordered/eventDisplay_timebased.C

Go to the documentation of this file.

int eventDisplay_timebased()
{
    //-----User Settings:-----------------------------------------------
  TString  SimEngine      ="TGeant3"; 
  TString  InputFile     ="dpm_simulation.root";
  //------------------------------------------------------------------


// Load basic libraries
//  gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");
//  rootlogon();
//  gSystem->Load("libEve");
//  gSystem->Load("libEventDisplay");
//  gSystem->Load("libPndEventDisplay");
  FairLogger* logger = FairLogger::GetLogger();
//  logger->SetLogToScreen(kTRUE);
  logger->SetLogVerbosityLevel("LOW");
  logger->SetLogScreenLevel("INFO");


  PndFileNameCreator creator(InputFile.Data());
  TString digiFile = creator.GetDigiFileName("timebased");
  TString recoFile = creator.GetRecoFileName("timebased");
  TString trackF = creator.GetTrackFindingFileName("");
  TString ParFile = creator.GetParFileName();
                                     
  // -----   Reconstruction run   -------------------------------------------
  FairRunAna *fRun= new FairRunAna();
  fRun->SetInputFile(InputFile.Data());
  fRun->AddFriend(recoFile.Data());
  fRun->AddFriend(digiFile.Data());
//  fRun->AddFriend(trackF.Data());
  fRun->SetOutputFile("tst.root");

  fRun->RunWithTimeStamps();

  FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
  FairParRootFileIo* parInput1 = new FairParRootFileIo();
  parInput1->open(ParFile.Data());
  rtdb->setFirstInput(parInput1);
  FairEventManager *fMan= new FairEventManager();

  FairGeane *Geane = new FairGeane();
  fRun->AddTask(Geane);
//
 //     fRun->Init();
 
 // SetPalette(1);

 //----------------------Traks and points -------------------------------------
  FairMCTracks *Track =  new FairMCTracks ("Monte-Carlo Tracks");
  FairMCPointDraw *MvdPoints =   new FairMCPointDraw ("MVDPoint",kBlue,  kFullSquare);
//  FairMCPointDraw *EMCPoints =   new FairMCPointDraw ("EmcHit",kOrange,  kFullSquare);
//  FairMCPointDraw *TofPoint =    new FairMCPointDraw ("TofPoint",kYellow,  kFullSquare);
//  FairMCPointDraw *TofSciFPoint= new FairMCPointDraw ("TofSciFPoint",kTeal, kFullSquare);
//  FairMCPointDraw *MdtPoint =    new FairMCPointDraw ("MdtPoint",kAzure, kFullSquare);
//  FairMCPointDraw *PndDrcPoint = new FairMCPointDraw ("PndDrcPoint",kViolet, kFullSquare);
//  FairMCPointDraw *PndDchPoint = new FairMCPointDraw ("PndDchPoint",kPink, kFullSquare);
//  FairMCPointDraw *PndTpcPoint = new FairMCPointDraw ("PndTpcPoint",kCyan,  kFullSquare);
  FairMCPointDraw *PndSTTPoint = new FairMCPointDraw ("STTPoint",kMagenta, kFullSquare);
//  FairMCPointDraw *PndGEMPoint = new FairMCPointDraw ("GEMPoint",kRed, kFullSquare);
//  FairMCPointDraw *PndDskPoint = new FairMCPointDraw ("DskCerenkov",kGreen, kFullSquare);
//  FairHitDraw *EMCRecoHit = new FairHitDraw("EmcRecoHit");
//
  PndTrackDraw* genfitTrack = new PndTrackDraw("BarrelTrack", kTRUE);
  genfitTrack->SetTimeWindowPlus(5);
  genfitTrack->SetTimeWindowMinus(5);
  fMan->AddTask(genfitTrack);

  PndTrackDraw* PndTrackRiemann = new PndTrackDraw("MVDTrack");
  PndTrackRiemann->SetTimeWindowPlus(10);
  PndTrackRiemann->SetTimeWindowMinus(10);

//  PndRiemannTrackDraw* RiemannTrack = new PndRiemannTrackDraw("MVDRiemannTrack");
//  PndMvdDigiPixelDraw* MvdDigiPixel = new PndMvdDigiPixelDraw("MVDPixelDigis");
                                                            
  FairHitDraw *MvdRecoHit =   new FairHitDraw ("MVDHitsPixel");
  MvdRecoHit->SetTimeWindowPlus(5);
  MvdRecoHit->SetTimeWindowMinus(5);
  FairHitDraw *MvdRecoStrip = new FairHitDraw ("MVDHitsStrip");
  MvdRecoStrip->SetTimeWindowPlus(5);
  MvdRecoStrip->SetTimeWindowMinus(5);
  FairHitDraw *STTHits = new FairHitDraw ("STTSortedHits");
  STTHits->SetTimeWindowPlus(200);
  STTHits->SetTimeWindowMinus(1);
  fMan->AddTask(Track);
  fMan->AddTask(MvdPoints);

  PndSttIsochroneDraw* STTIsochrone = new PndSttIsochroneDraw("STTHit");
  STTIsochrone->SetTimeWindowPlus(250);
  STTIsochrone->SetTimeWindowMinus(1);
  STTIsochrone->UseIsochroneTime();
  fMan->AddTask(STTIsochrone);

//  fMan->AddTask(MvdDigiPixel);
//  fMan->AddTask(EMCPoints);
//  fMan->AddTask(TofPoint);
//  fMan->AddTask( TofSciFPoint);
//  fMan->AddTask( MdtPoint);
//  fMan->AddTask( PndDrcPoint);
//  fMan->AddTask( PndDchPoint);
//  fMan->AddTask( PndTpcPoint);
  fMan->AddTask( PndSTTPoint);
//  fMan->AddTask( PndGEMPoint);
//  fMan->AddTask( PndDskPoint);
//
//  fMan->AddTask(EMCRecoHit);
  fMan->AddTask(MvdRecoHit);
  fMan->AddTask(MvdRecoStrip);
  fMan->AddTask(STTHits);
  
//  fMan->AddTask(RiemannCand);
//  fMan->AddTask(PndTrackRiemann);

  fMan->Init();

  return 0;
}