Fork me on GitHub

Changes between Version 42 and Version 43 of WorkBook/LibraryInterface


Ignore:
Timestamp:
Mar 9, 2013, 6:56:36 PM (12 years ago)
Author:
Pavel Demin
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • WorkBook/LibraryInterface

    v42 v43  
    11== Library Interface ==
    22
    3 Here is an explanation of how Delphes' library can be called from a program. 
     3Here is an explanation of how Delphes' library can be called from a program.
    44
    55== Installation ==
     
    8989    modularDelphes = new Delphes("Delphes");
    9090    modularDelphes->SetConfReader(confReader);
    91    
     91
    9292    factory = modularDelphes->GetFactory();
    9393
     
    102102    {
    103103      modularDelphes->Clear();
    104       ConvertInput(factory, allParticleOutputArray, stableParticleOutputArray, partonOutputArray); 
     104      ConvertInput(factory, allParticleOutputArray, stableParticleOutputArray, partonOutputArray);
    105105      modularDelphes->ProcessTask();
    106106      ConvertOutput(modularDelphes);
     
    125125  TDatabasePDG *pdg;
    126126  TParticlePDG *pdgParticle;
     127  Int_t pdgCode;
     128
    127129  Int_t pid, status;
    128130  Double_t px, py, pz, e;
    129   Double_t x, y, z, t; 
     131  Double_t x, y, z, t;
    130132
    131133  pdg = TDatabasePDG::Instance();
    132  
     134
    133135  while(ReadParticle())
    134136  {
     
    136138
    137139    candidate->PID = pid;
     140    pdgCode = TMath::Abs(candidate->PID);
    138141
    139142    candidate->Status = status;
     
    155158      stableParticleOutputArray->Add(candidate);
    156159    }
    157     else if(status == 2)
     160    else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
    158161    {
    159162      partonOutputArray->Add(candidate);