standard_tracking.C File Reference

Go to the source code of this file.

Functions
int	standard_tracking (Int_t nEvents=0)

Function Documentation

int standard_tracking

(

Int_t

nEvents = 0

)

Definition at line 4 of file standard_tracking.C.

References PndMasterRunAna::AddFriend(), PndMasterRunAna::Finish(), fRun, Geane, nEvents, PndTrkTracking2::NoMvdAloneTracking(), output, recoKalman, PndRecoKalmanTask2::SetBusyCut(), PndRecoKalmanTask::SetBusyCut(), PndMasterRunAna::SetInput(), PndTrkTracking2::SetInputBranchName(), PndMissingPzCleanerTask::SetInputTrackBranch(), PndMasterRunAna::SetOutput(), PndMasterRunAna::SetParamAsciiFile(), PndPersistencyTask::SetPersistency(), PndMissingPzCleanerTask::SetRemoveTrack(), PndRecoKalmanTask::SetTrackInBranchName(), PndRecoKalmanTask2::SetTrackInBranchName(), PndRecoKalmanTask2::SetTrackOutBranchName(), PndRecoKalmanTask::SetTrackOutBranchName(), PndRecoKalmanTask::SetTrackRep(), PndMasterRunAna::Setup(), SttMvdGemTracking, tracking, and TString.

{
  //-----User Settings:------------------------------------------------------
  TString  parAsciiFile   = "all.par";
  TString  prefix         = "evtcomplete";
  TString  input          = "psi2s_Jpsi2pi_Jpsi_mumu.dec";
  TString  output         = "standard";
  TString  friend1        = "sim";
  TString  friend2        = "digi";
  TString  friend3        = "";
  TString  friend4        = "";
  TString  fOptions       = "gf2"; // "gf2" for genfit 2

  // -----   Initial Settings   --------------------------------------------
  PndMasterRunAna *fRun= new PndMasterRunAna();
  fRun->SetInput(input);
  fRun->SetOutput(output);
  fRun->AddFriend(friend1);
  fRun->AddFriend(friend2);
  fRun->AddFriend(friend3);
  fRun->AddFriend(friend4);
  fRun->SetParamAsciiFile(parAsciiFile);
  fRun->Setup(prefix);

  FairGeane *Geane = new FairGeane();
    fRun->AddTask(Geane);


        PndTrkTracking2* tracking = NULL;
        fRun->AddTask(tracking = new PndTrkTracking2(0,false,false,true)); // 1
        tracking->SetInputBranchName("STTHit","MVDHitsPixel","MVDHitsStrip");
        tracking->NoMvdAloneTracking();
        tracking->SetPersistency(kFALSE);

    PndSttMvdGemTracking *SttMvdGemTracking = NULL;
    fRun->AddTask(SttMvdGemTracking = new PndSttMvdGemTracking(0));
    SttMvdGemTracking->SetPersistency(kTRUE);

    if (fOptions.Contains("gf1")){
                PndRecoKalmanTask* recoKalman = NULL;
                fRun->AddTask(recoKalman = new PndRecoKalmanTask());
                recoKalman->SetTrackInBranchName("SttMvdGemTrack");
                recoKalman->SetTrackOutBranchName("SttMvdGemGenTrack");

                recoKalman->SetBusyCut(50); // CHECK to be tuned
                //recoKalman->SetIdealHyp(kTRUE);
                //recoKalman->SetNumIterations(3);
                recoKalman->SetTrackRep(0); // 0 Geane (default), 1 RK
                //recoKalman->SetPropagateToIP(kFALSE);
        } else if (fOptions.Contains("gf2")){
                PndRecoKalmanTask2* recoKalman = NULL;
                fRun->AddTask(recoKalman = new PndRecoKalmanTask2());
                recoKalman->SetTrackInBranchName("SttMvdGemTrack");
                recoKalman->SetTrackOutBranchName("SttMvdGemGenTrack");

                recoKalman->SetBusyCut(50); // CHECK to be tuned
                //recoKalman->SetIdealHyp(kTRUE);
                //recoKalman->SetNumIterations(3);
                //recoKalman->SetTrackRep(0); // 0 Geane (default), 1 RK
                //recoKalman->SetPropagateToIP(kFALSE);
        }

        if (fOptions.Contains("filtered")){
                  PndMissingPzCleanerTask* cleaner = NULL;
                  fRun->AddTask(cleaner = new PndMissingPzCleanerTask()); //4
                  if ( (!fOptions.Contains("day1")) || (fOptions.Contains("gem")) )
                  {
                          cleaner->SetInputTrackBranch("SttMvdGemGenTrack");
                  } else {
                          cleaner->SetInputTrackBranch("SttMvdGenTrack");
                  }
                  cleaner->SetRemoveTrack(kTRUE);
        }


  // -----   Intialise and run   --------------------------------------------
  fRun->Init();
  fRun->Run(0, nEvents);
  fRun->Finish();
  return 0;
}