sPhenix code displayed by LXR master/​analysis/​Prototype3/​HCAL/​ShowerCalib/​Proto3ShowerCalib.h	Source navigation Diff markup Identifier search General search
	Version: master Revert   Change

File indexing completed on 2025-08-05 08:14:44

 #ifndef __Proto3ShowerCalib_H__
 #define __Proto3ShowerCalib_H__
 
 #include <fun4all/SubsysReco.h>
 #include <TNtuple.h>
 #include <TFile.h>
 #include <string>
 #include <stdint.h>
 #include <fstream>
 
 class PHCompositeNode;
 class PHG4HitContainer;
 class Fun4AllHistoManager;
 class TH1F;
 class TF1;
 class TTree;
 class SvtxEvalStack;
 class PHG4Particle;
 class RawTowerGeom;
 class RawTowerContainer;
 class SvtxTrack;
 
 /// \class Proto3ShowerCalib
 class Proto3ShowerCalib : public SubsysReco
 {
 
 public:
 
   Proto3ShowerCalib(const std::string &filename = "ana.root");
   virtual
   ~Proto3ShowerCalib();
 
   int
   Init(PHCompositeNode *topNode);
   int
   InitRun(PHCompositeNode *topNode);
   int
   process_event(PHCompositeNode *topNode);
   int
   End(PHCompositeNode *topNode);
 
   int
   LoadRecalibMap(const std::string & file);
 
   class Eval_Run : public TObject
   {
   public:
     Eval_Run()
     {
       reset();
     }
     virtual
     ~Eval_Run()
     {
     }
 
     void
     reset()
     {
       run = -31454;
       event = -31454;
       beam_mom = -0;
       hodo_h = -31454;
       hodo_v = -31454;
       C2_sum = -31454;
       C1 = -31454;
       EoP = -31454;
       valid_hodo_v = false;
       valid_hodo_h = false;
       trigger_veto_pass = false;
       good_temp = false;
       good_e = false;
       good_h = false;
     }
 
     int run;
     int event;
     float beam_mom;
     int hodo_h;
     int hodo_v;
     float C2_sum;
     float C1;
     float EoP;
     bool valid_hodo_v;
     bool valid_hodo_h;
     bool trigger_veto_pass;
     bool good_temp;
     bool good_e;
     bool good_h;
 
   ClassDef(Eval_Run,2)
   };
 
   class HCAL_shower : public TObject
   {
    public:
     HCAL_shower()
     {
      reset();
     }
 
     virtual ~HCAL_shower(){}
 
     void reset()
     {
      emcal_e = 0.;
      hcalin_e = 0.;
      hcalout_e = 0.;
      sum_e = 0.;
      for(int itwr=0; itwr<16; itwr++)
      {
       hcalin_twr_e[itwr] = 0.;
       hcalout_twr_e[itwr] = 0.;
      }
      for(int itwr=0; itwr<64; itwr++)
       emcal_twr_e[itwr] = 0.;
 
      emcal_e_recal = 0.;
      hcalin_e_recal = 0.;
      hcalout_e_recal = 0.;
      sum_e_recal = 0.;
      hcal_asym = -999;
     }
 
    float emcal_e;
    float hcalin_e;
    float hcalout_e;
    float sum_e;
    float emcal_twr_e[64];
    float hcalin_twr_e[16];
    float hcalout_twr_e[16];
 
    float emcal_e_recal;
    float hcalin_e_recal;
    float hcalout_e_recal;
    float hcal_asym;
    float sum_e_recal;
   };
 
   class Scientillators : public TObject
   {
    public:
     Scientillators(){ reset(); }
     virtual ~Scientillators(){}
     void reset()
     {
       hcalin_edep = 0.;
       hcalout_edep = 0.;
       hcalin_lightyield = 0.;
       hcalout_lightyield = 0.;
     }
     float hcalin_edep;
     float hcalout_edep;
 
     float hcalin_lightyield;
     float hcalout_lightyield;
   };
 
   class Time_Samples : public TObject
   {
    public:
     Time_Samples(){ reset(); }
     virtual ~Time_Samples(){}
 
     void reset()
     {
      for(int itower=0; itower<16; itower++)
      {
       for(int isamp=0; isamp<24; isamp++)
       {
        hcalin_time_samples[itower][isamp] = -9999;
        hcalout_time_samples[itower][isamp] = -9999;
       }
      }
 
      for(int itower=0; itower<64; itower++)
      {
       for(int isamp=0; isamp<24; isamp++)
       {
        emcal_time_samples[itower][isamp] = -9999;
       }
      }
     }
    float hcalin_time_samples[16][24];
    float hcalout_time_samples[16][24];
    float emcal_time_samples[64][24];
   };
 
   void set_sim(const bool b)
   { _is_sim = b; }
 
   void fill_tower_samples(const bool);
   void fill_time_samples(const bool);
   void fill_slat_values(const bool);
 private:
 
   // calorimeter size
   enum
   {
     n_size_emcal = 8,
     n_size_hcalin = 4,
     n_size_hcalout = 4
   };
 
   Fun4AllHistoManager *
   get_HistoManager();
 
   std::string _filename;
 
   unsigned long _ievent;
 
   fstream fdata;
 
   Eval_Run _eval_run;
 
   HCAL_shower _shower;
 
   Time_Samples _time_samples;
 
   Scientillators _slats;
 
   //Add smearing
   TF1 *smearing;
   bool _is_sim;
   bool _fill_tower_samples;
   bool _fill_time_samples;
   bool _fill_slats;
 
   std::map<std::pair<int, int>, double> _emcal_recalib_const;
   std::map<std::pair<int, int>, double> _hcalin_recalib_const;
   std::map<std::pair<int, int>, double> _hcalout_recalib_const;
 
 };
 
 #endif // __Proto3ShowerCalib_H__

This page was automatically generated by the 2.3.7 LXR engine. The LXR team