QA/run3/pidideal_complete.C File Reference

Go to the source code of this file.

Functions
int	pidideal_complete (Int_t nEvents=0)

Function Documentation

int pidideal_complete

(

Int_t

nEvents = 0

)

Definition at line 4 of file QA/run3/pidideal_complete.C.

References PndMasterRunAna::AddFriend(), assDisc, assDrc, assEMC, assMC, assMdt, assMvd, assStt, corr, PndMasterRunAna::Finish(), fRun, Geane, PndEmcMapper::Init(), nEvents, output, PndPidCorrelator::SetBarrelTrackBranch(), PndPidCorrelator::SetForwardTrackBranch(), PndPidCorrelator::SetIdeal(), PndPidCorrelator::SetIdealHyp(), PndMasterRunAna::SetInput(), PndMasterRunAna::SetOutput(), PndMasterRunAna::SetParamAsciiFile(), PndMasterRunAna::Setup(), and TString.

{
  //-----User Settings:------------------------------------------------------
  TString  parAsciiFile   = "all.par";
  TString  input          = "psi2s_Jpsi2pi_Jpsi_mumu.dec"; 
  TString  output         = "pidideal";
  TString  friend1        = "sim";
  TString  friend2        = "digi";
  TString  friend3        = "recoideal";
  TString  friend4        = "";
  
  // -----   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();
  
  // -----   Add tasks   ----------------------------------------------------
  FairGeane *Geane = new FairGeane();
  fRun->AddTask(Geane);

  PndPidCorrelator* corr = new PndPidCorrelator();
  //corr->SetVerbose();
  corr->SetBarrelTrackBranch("SttMvdGemGenTrack");
//  corr->SetInputIDBranch("SttMvdGemGenTrackID");
  corr->SetForwardTrackBranch("FtsIdealGenTrack");
//  corr->SetInputIDBranch2("FtsIdealGenTrackID");
  corr->SetIdeal(kTRUE);
  corr->SetIdealHyp(kTRUE);
  //corr->SetDebugMode(kTRUE);
  //corr->SetFast(kTRUE);
  //corr->SetBackPropagate(kFALSE);
  fRun->AddTask(corr);

  PndMcCloner *clone = new PndMcCloner();
  fRun->AddTask(clone);
 
  PndPidIdealAssociatorTask *assMC= new PndPidIdealAssociatorTask();
  fRun->AddTask(assMC);

  PndPidMvdAssociatorTask *assMvd= new PndPidMvdAssociatorTask();
  fRun->AddTask(assMvd);

  PndPidMdtHCAssociatorTask *assMdt= new PndPidMdtHCAssociatorTask();
  fRun->AddTask(assMdt);

  PndPidDrcAssociatorTask *assDrc= new PndPidDrcAssociatorTask();
  fRun->AddTask(assDrc);

  PndPidDiscAssociatorTask *assDisc= new PndPidDiscAssociatorTask();
  fRun->AddTask(assDisc);

  PndPidSttAssociatorTask *assStt= new PndPidSttAssociatorTask();
  fRun->AddTask(assStt);

  PndPidEmcBayesAssociatorTask *assEMC= new PndPidEmcBayesAssociatorTask();
  fRun->AddTask(assEMC);
  
  // -----   Intialise and run   --------------------------------------------
  PndEmcMapper::Init(1);
  fRun->Init();
  fRun->Run(0, nEvents);
  fRun->Finish();

  return 0;
}