Herwig/Analysis/GammaJetAnalysis.h

0001 // -*- C++ -*-
0002 //
0003 // GammaJetAnalysis.h is a part of Herwig - A multi-purpose Monte Carlo event generator
0004 // Copyright (C) 2002-2019 The Herwig Collaboration
0005 //
0006 // Herwig is licenced under version 3 of the GPL, see COPYING for details.
0007 // Please respect the MCnet academic guidelines, see GUIDELINES for details.
0008 //
0009 #ifndef HERWIG_GammaJetAnalysis_H
0010 #define HERWIG_GammaJetAnalysis_H
0011 //
0012 // This is the declaration of the GammaJetAnalysis class.
0013 //
0014
0015 #include "ThePEG/Repository/CurrentGenerator.h"
0016 #include "ThePEG/Handlers/AnalysisHandler.h"
0017 #include "Herwig/Utilities/Histogram.h"
0018
0019 namespace Herwig {
0020
0021 using namespace ThePEG;
0022
0023 /** \ingroup Analysis
0024  * GammaJetAnalysis selects the photon with the hightest pt in the final
0025  * state and books a number of histograms from its four momentum.  The
0026  * results are witten in topdrawer format to the working directory.
0027  *
0028  * @see \ref GammaJetAnalysisInterfaces "The interfaces"
0029  * defined for GammaJetAnalysis.
0030  */
0031 class GammaJetAnalysis: public AnalysisHandler {
0032
0033 public:
0034
0035   /**
0036    * The default constructor.
0037    */
0038   GammaJetAnalysis();
0039
0040   /** @name Virtual functions required by the AnalysisHandler class. */
0041   //@{
0042   /**
0043    * Analyze a given Event. Note that a fully generated event
0044    * may be presented several times, if it has been manipulated in
0045    * between. The default version of this function will call transform
0046    * to make a lorentz transformation of the whole event, then extract
0047    * all final state particles and call analyze(tPVector) of this
0048    * analysis object and those of all associated analysis objects. The
0049    * default version will not, however, do anything on events which
0050    * have not been fully generated, or have been manipulated in any
0051    * way.
0052    * @param event pointer to the Event to be analyzed.
0053    * @param ieve the event number.
0054    * @param loop the number of times this event has been presented.
0055    * If negative the event is now fully generated.
0056    * @param state a number different from zero if the event has been
0057    * manipulated in some way since it was last presented.
0058    */
0059   virtual void analyze(tEventPtr event, long ieve, int loop, int state);
0060   //@}
0061
0062 public:
0063
0064   /**
0065    * The standard Init function used to initialize the interfaces.
0066    * Called exactly once for each class by the class description system
0067    * before the main function starts or
0068    * when this class is dynamically loaded.
0069    */
0070   static void Init();
0071
0072 protected:
0073
0074   /** @name Clone Methods. */
0075   //@{
0076   /**
0077    * Make a simple clone of this object.
0078    * @return a pointer to the new object.
0079    */
0080   virtual IBPtr clone() const {return new_ptr(*this);}
0081
0082   /** Make a clone of this object, possibly modifying the cloned object
0083    * to make it sane.
0084    * @return a pointer to the new object.
0085    */
0086   virtual IBPtr fullclone() const {return new_ptr(*this);}
0087   //@}
0088
0089 protected:
0090
0091   /** @name Standard Interfaced functions. */
0092   //@{
0093   /**
0094    * Finalize this object. Called in the run phase just after a
0095    * run has ended. Used eg. to write out statistics.
0096    */
0097   virtual void dofinish();
0098
0099 private:
0100
0101   /**
0102    * The assignment operator is private and must never be called.
0103    * In fact, it should not even be implemented.
0104    */
0105   GammaJetAnalysis & operator=(const GammaJetAnalysis &) = delete;
0106
0107 private:
0108   /**
0109    *   \f$p_T\f$ of the photon
0110    */
0111   Histogram _ptg;
0112   Histogram _ptgZoom;
0113
0114   /**
0115    *   Energy of the photon
0116    */
0117   Histogram _Eg;
0118
0119   /**
0120    *  Rapidity of the photon
0121    */
0122   Histogram _rapg;
0123
0124   /**
0125    *  Azimuth of the photon
0126    */
0127   Histogram _phig;
0128
0129 };
0130
0131 }
0132
0133 #endif /* HERWIG_GammaJetAnalysis_H */