eicsmear/erhic/EventRapgap.h

0001 /**
0002  \file
0003  Declaration of class erhic::EventRapgap.
0004
0005  \author    Thomas Burton
0006  \date      2011-07-07
0007  \copyright 2011 Brookhaven National Lab
0008  */
0009
0010 #ifndef INCLUDE_EICSMEAR_ERHIC_EVENTRAPGAP_H_
0011 #define INCLUDE_EICSMEAR_ERHIC_EVENTRAPGAP_H_
0012
0013 #include <string>
0014
0015 #include <Rtypes.h>
0016
0017 #include "eicsmear/erhic/EventMC.h"
0018
0019 namespace erhic {
0020
0021 /**
0022  Describes an event from the generator RAPGAP.
0023  */
0024 class EventRapgap : public EventMC {
0025  public:
0026   virtual bool Parse(const std::string&);
0027
0028
0029   /**  Returns a pointer to the scattered lepton in the event record.
0030        This is the first (only?) particle that matches the following:
0031        1) pdg code equals that of incident lepton beam.
0032        2) status code is 1 i.e. it's a stable/final-state particle.
0033        3) the parent is track 1 or 2
0034   */
0035   const ParticleMC* ScatteredLepton() const;
0036
0037   /**
0038    Returns a pointer to the exchanged boson.
0039    It would probably be the third track, but we'll go with the first status=21 boson
0040    that has particle 1 or 2 as parent
0041    */
0042   virtual const ParticleMC* ExchangeBoson() const;
0043
0044   Double32_t Get_cs() const {return cs;}
0045   Double32_t Get_sigma_cs() const {return sigma_cs;}
0046   Double32_t Get_s() const {return s;}
0047   Double32_t Get_q2() const {return q2;}
0048   Double32_t Get_xgam() const {return xgam;}
0049   Double32_t Get_xpr() const {return xpr;}
0050   Double32_t Get_Pt_h() const {return Pt_h;}
0051   Double32_t Get_t() const {return t;}
0052   Double32_t Get_x_pom() const {return x_pom;}
0053   Double32_t Get_sHat2() const {return sHat2;}
0054   Double32_t Get_z() const {return z;}
0055   Double32_t Get_x1() const {return x1;}
0056   Double32_t Get_phi1() const {return phi1;}
0057   Double32_t Get_pt2_hat() const {return pt2_hat;}
0058   Double32_t Get_sHat() const {return sHat;}
0059
0060  protected:
0061
0062   Int_t idir;
0063   Int_t idisdif;
0064   Int_t genevent;
0065
0066   Double32_t cs;
0067   Double32_t sigma_cs;
0068   Double32_t s;
0069   Double32_t q2;
0070   Double32_t xgam;
0071   Double32_t xpr;
0072   Double32_t Pt_h;
0073   Double32_t t;
0074   Double32_t x_pom;
0075   Double32_t sHat2;
0076   Double32_t z;
0077   Double32_t x1;
0078   Double32_t phi1;
0079   Double32_t pt2_hat;
0080   Double32_t sHat;  // Partonic centre-of-mass energy
0081
0082   ClassDef(erhic::EventRapgap, 2)
0083
0084 };
0085
0086 }  // namespace erhic
0087
0088 #endif  // INCLUDE_EICSMEAR_ERHIC_EVENTRAPGAP_H_