read_particles_from_TAUOLA.h

00001 #ifndef _READ_PARTICLES_FROM_TAUOLA_H_
00002 #define _READ_PARTICLES_FROM_TAUOLA_H_
00003 
00004 // HepMC headers
00005 #include "HepMC/GenEvent.h"
00006 #include "HepMC/GenParticle.h"
00007 #include "HepMC/IO_GenEvent.h"
00008 
00009 #include "TauSpinner/SimpleParticle.h"
00010 #include <vector>
00011 using namespace std;
00012 using namespace TauSpinner;
00013 
00014 /** Read HepMC::IO_GenEvent.
00015 
00016   Read HepMC event from data file
00017   and return particles needed for tau spin weight calculation.
00018   
00019   This routine is prepared for use with files generated by Pythia8.
00020   Fills:
00021   
00022   'X'              - Heavy particle (W+/-, H+/-, H, Z)
00023   'tau'            - first tau
00024   'tau2'           - second tau or nu_tau, if 'X' decays to one tau only
00025   'tau_daughters'  - daughters of 'tau'
00026   'tau2_daughters' - daughters of 'tau2' or empty list, if 'tau2' is nu_tau.
00027   
00028   Returns:
00029   0 - no more events to read               (finished processing the file)
00030   1 - no decay found in the event          (finished processing current event)
00031   2 - decay found and processed correctly.
00032       Event will continue to be processed
00033       with next function call. */
00034 int readParticlesFromTAUOLA_HepMC(HepMC::IO_GenEvent &input_file, SimpleParticle &X, SimpleParticle &tau, SimpleParticle &tau2, vector<SimpleParticle> &tau_daughters, vector<SimpleParticle> &tau2_daughters);
00035 
00036 /** Get daughters of HepMC::GenParticle
00037 
00038   Recursively searches for final-state daughters of 'x' */
00039 vector<SimpleParticle> *getDaughters(HepMC::GenParticle *x);
00040 
00041 #endif