GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » Fairroot » General Discussion » genfit example macro
Re: genfit example macro [message #14404 is a reply to message #14319] Sun, 27 January 2013 13:41 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
Dear Alexander,
after your mail I realized that your FairRoot forum entry (which I read) was related to Panda...you were mentioning MDCs and TPC that we do not have then I supposed you were referring to some other experiment and I did not reply, sorry Smile

If you have the pandaroot trunk, go to macro/run and there you will find the standard full chain macros. You can execute:

sim_complete.C
digi_complete.C
reco_complete.C


In the reco you can leave only the fts tasks and comment out the previous parts.This means you have to leave only:

  PndFtsTrackerIdeal* trackFts = new PndFtsTrackerIdeal();
  trackFts->SetRelativeMomentumSmearing(0.05);
  trackFts->SetVertexSmearing(0.05, 0.05, 0.05);
  trackFts->SetTrackingEfficiency(1.);
  trackFts->SetTrackOutput("FtsIdealTrack");
  fRun->AddTask(trackFts);

  PndRecoKalmanTask* recoKalmanFwd = new PndRecoKalmanTask();
  recoKalmanFwd->SetTrackInBranchName("FtsIdealTrack");
  //recoKalmanFwd->SetTrackInIDBranchName("FtsIdealTrackID");
  recoKalmanFwd->SetTrackOutBranchName("FtsIdealGenTrack");
  recoKalmanFwd->SetBusyCut(50); // CHECK to be tuned
  //recoKalmanFwd->SetIdealHyp(kTRUE);
  //recoKalmanFwd->SetNumIterations(3);
  fRun->AddTask(recoKalmanFwd);

  PndMCTrackAssociator* trackMC3 = new PndMCTrackAssociator();
  trackMC3->SetTrackInBranchName("FtsIdealGenTrack");
  trackMC3->SetTrackOutBranchName("FtsIdealGenTrackID");
  fRun->AddTask(trackMC3);


From the first you can set the momentum and position smearing, the second is the kalman, the third is a simply match toMC.
By default PndFtsTrackerIdeal run mvd+gem+fts, but if you check the PndFtsTrackerIdeal hedes (in fts folder):

  void SetFtsActivity(Bool_t act=kTRUE){fBranchActive[0]=act;}
  void SetGemActivity(Bool_t act=kTRUE){fBranchActive[1]=act;}
  void SetMvdActivity(Bool_t act=kTRUE){fBranchActive[2]=act;}


you can switch on or off the different parts.

TPC is not supported anymore and there are no TPC classes anymore in the trunk.
If you find some problems or bugs, I would suggest to report on the PandaRoot forum and not in FairRoot, under bugs thread or also tracking if you prefer.

Good luck!!!
 
Read Message
Read Message
Previous Topic: configure csh script
Next Topic: FairMCPoint
Goto Forum:
  


Current Time: Sat Oct 16 19:44:36 CEST 2021

Total time taken to generate the page: 0.02307 seconds