Herwig/Analysis/GammaGammaAnalysis.h

0001 // -*- C++ -*-
0002 //
0003 // GammaGammaAnalysis.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_GammaGammaAnalysis_H
0010 #define HERWIG_GammaGammaAnalysis_H
0011 //
0012 // This is the declaration of the GammaGammaAnalysis 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  * GammaGammaAnalysis is for the analysis of events with a pair of hard
0025  * photons produced.  These are selected as the two highest pt photons
0026  * in the final state of the event.  A topdrawer file with histograms
0027  * is written to the working directory.
0028  *
0029  * @see \ref GammaGammaAnalysisInterfaces "The interfaces"
0030  * defined for GammaGammaAnalysis.
0031  */
0032 class GammaGammaAnalysis: public AnalysisHandler {
0033
0034 public:
0035
0036   /**
0037    * The default constructor.
0038    */
0039   GammaGammaAnalysis();
0040
0041   /** @name Virtual functions required by the AnalysisHandler class. */
0042   //@{
0043   /**
0044    * Analyze a given Event. Note that a fully generated event
0045    * may be presented several times, if it has been manipulated in
0046    * between. The default version of this function will call transform
0047    * to make a lorentz transformation of the whole event, then extract
0048    * all final state particles and call analyze(tPVector) of this
0049    * analysis object and those of all associated analysis objects. The
0050    * default version will not, however, do anything on events which
0051    * have not been fully generated, or have been manipulated in any
0052    * way.
0053    * @param event pointer to the Event to be analyzed.
0054    * @param ieve the event number.
0055    * @param loop the number of times this event has been presented.
0056    * If negative the event is now fully generated.
0057    * @param state a number different from zero if the event has been
0058    * manipulated in some way since it was last presented.
0059    */
0060   virtual void analyze(tEventPtr event, long ieve, int loop, int state);
0061   //@}
0062
0063 public:
0064
0065   /**
0066    * The standard Init function used to initialize the interfaces.
0067    * Called exactly once for each class by the class description system
0068    * before the main function starts or
0069    * when this class is dynamically loaded.
0070    */
0071   static void Init();
0072
0073 protected:
0074
0075   /** @name Clone Methods. */
0076   //@{
0077   /**
0078    * Make a simple clone of this object.
0079    * @return a pointer to the new object.
0080    */
0081   virtual IBPtr clone() const {return new_ptr(*this);}
0082
0083   /** Make a clone of this object, possibly modifying the cloned object
0084    * to make it sane.
0085    * @return a pointer to the new object.
0086    */
0087   virtual IBPtr fullclone() const {return new_ptr(*this);}
0088   //@}
0089
0090 protected:
0091
0092   /** @name Standard Interfaced functions. */
0093   //@{
0094   /**
0095    * Finalize this object. Called in the run phase just after a
0096    * run has ended. Used eg. to write out statistics.
0097    */
0098   virtual void dofinish();
0099   //@}
0100
0101 private:
0102
0103   /**
0104    * The assignment operator is private and must never be called.
0105    * In fact, it should not even be implemented.
0106    */
0107   GammaGammaAnalysis & operator=(const GammaGammaAnalysis &) = delete;
0108
0109 private:
0110
0111   /**
0112    *   \f$p_T\f$ of the harder photon
0113    */
0114   Histogram _ptharder;
0115
0116   /**
0117    *   \f$p_T\f$ of the softer photon
0118    */
0119   Histogram _ptsofter;
0120
0121   /**
0122    *   \f$p_T\f$ of the photon pair
0123    */
0124   Histogram _ptpair;
0125
0126   /**
0127    *   Energy of the harder photon
0128    */
0129   Histogram _Eharder;
0130
0131   /**
0132    *   Energy of the softer photon
0133    */
0134   Histogram _Esofter;
0135
0136   /**
0137    *   Energy of the photon pair
0138    */
0139   Histogram _Epair;
0140
0141   /**
0142    *  Rapidity of the harder photon
0143    */
0144   Histogram _rapharder;
0145
0146   /**
0147    *  Rapidity of the softer photon
0148    */
0149   Histogram _rapsofter;
0150
0151   /**
0152    *  Rapidity of the photon pair
0153    */
0154   Histogram _rappair;
0155
0156   /**
0157    *  Azimuth of the harder photon
0158    */
0159   Histogram _phiharder;
0160
0161   /**
0162    *  Azimuth of the softer photon
0163    */
0164   Histogram _phisofter;
0165
0166   /**
0167    *  Azimuth of the photon pair
0168    */
0169   Histogram _deltaphi;
0170
0171   /**
0172    *  Invariant mass of the pair
0173    */
0174   Histogram _mpair;
0175 };
0176
0177 }
0178
0179 #endif /* HERWIG_GammaGammaAnalysis_H */