EMCal-calibration/EMCalCalib/MCPhoton.h

0001 // $Id: $
0002
0003 /*!
0004  * \file MCPhoton.h
0005  * \brief
0006  * \author Jin Huang <jhuang@bnl.gov>
0007  * \version $Revision:   $
0008  * \date $Date: $
0009  */
0010
0011 #ifndef MCPhoton_H_
0012 #define MCPhoton_H_
0013
0014 #include <phool/PHObject.h>
0015 #include <phool/phool.h>
0016
0017 /*!
0018  * \brief MCPhoton
0019  */
0020 class MCPhoton : public PHObject
0021 {
0022 public:
0023   MCPhoton();
0024   virtual
0025   ~MCPhoton();
0026
0027   virtual void
0028   Clear(Option_t *option = "");
0029
0030   virtual void
0031   Reset();
0032
0033   int trackID;
0034   int charge;
0035   float quality;
0036   float chisq;
0037   float ndf;
0038   float nhits;
0039   unsigned int layers;
0040
0041   float dca2d;
0042   float dca2dsigma;
0043   float px;
0044   float py;
0045   float pz;
0046   float pcax;
0047   float pcay;
0048   float pcaz;
0049
0050   float presdphi;
0051   float presdeta;
0052   float prese3x3;
0053   float prese;
0054
0055   float cemcdphi;
0056   float cemcdeta;
0057   float cemce3x3;
0058   float cemce;
0059
0060   float hcalindphi;
0061   float hcalindeta;
0062   float hcaline3x3;
0063   float hcaline;
0064
0065   float hcaloutdphi;
0066   float hcaloutdeta;
0067   float hcaloute3x3;
0068   float hcaloute;
0069
0070   int gtrackID;
0071   int gflavor;
0072   int ng4hits;
0073
0074   float gpx;
0075   float gpy;
0076   float gpz;
0077   float gvx;
0078   float gvy;
0079   float gvz;
0080   float gfpx;
0081   float gfpy;
0082   float gfpz;
0083   float gfx;
0084   float gfy;
0085   float gfz;
0086   int gembed;
0087
0088   int nfromtruth;
0089
0090   enum
0091   {
0092     Max_N_Tower = 11
0093   };
0094
0095   float cemc_ieta[Max_N_Tower][Max_N_Tower];
0096   float cemc_iphi[Max_N_Tower][Max_N_Tower];
0097   float cemc_energy[Max_N_Tower][Max_N_Tower];
0098   float cemc_radius[Max_N_Tower][Max_N_Tower];
0099
0100   float hcalin_ieta[Max_N_Tower][Max_N_Tower];
0101   float hcalin_iphi[Max_N_Tower][Max_N_Tower];
0102   float hcalin_energy[Max_N_Tower][Max_N_Tower];
0103   float hcalin_radius[Max_N_Tower][Max_N_Tower];
0104
0105   float cemc_sum_energy;
0106   int cemc_sum_n_tower;
0107   float hcalin_sum_energy;
0108   int hcalin_sum_n_tower;
0109
0110   //! e/p based on CEMC and tracking
0111   double
0112   get_ep() const;
0113
0114   //! log-likelihood electron - ep matching
0115   float ll_ep_e;
0116
0117   //! log-likelihood hadron - ep matching
0118   float ll_ep_h;
0119
0120   //! log-likelihood electron - emcal + inner HCal energy deposition
0121   float ll_edep_e;
0122
0123   //! log-likelihood hadron - emcal + inner HCal energy deposition
0124   float ll_edep_h;
0125
0126   //! log-likelihood electron - shower shape
0127   float ll_shape_e;
0128
0129   //! log-likelihood hadron - shower shape
0130   float ll_shape_h;
0131
0132 ClassDef(MCPhoton,9)
0133 };
0134
0135 #endif /* MCPhoton_H_ */