Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
g4hitshift.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file g4hitshift.cc
1 #include <g4main/PHG4Hit.h>
2 #include <g4main/PHG4HitContainer.h>
3 
4 #include "g4hitshift.h"
5 
6 #include <fun4all/Fun4AllReturnCodes.h>
7 
8 #include <fun4all/Fun4AllHistoManager.h>
9 
10 #include <phool/PHCompositeNode.h>
11 
12 // G4Hits includes
13 #include <g4main/PHG4Hit.h>
14 #include <g4main/PHG4HitContainer.h>
15 #include <g4main/PHG4TruthInfoContainer.h>
16 #include <g4main/PHG4Particle.h>
17 #include <g4main/PHG4HitDefs.h> // for hit_idbits
18 
19 // G4Cells includes
20 #include <g4detectors/PHG4Cell.h>
21 #include <g4detectors/PHG4CellContainer.h>
22 #include <g4detectors/PHG4CellDefs.h> // for genkey, keytype
23 
24 // Tower includes
25 #include <calobase/RawTower.h>
26 #include <calobase/RawTowerContainer.h>
27 #include <calobase/RawTowerGeom.h>
28 #include <calobase/RawTowerGeomContainer.h>
29 #include <calobase/RawTowerDefs.h>
30 // Cluster includes
31 #include <calobase/RawCluster.h>
32 #include <calobase/RawClusterContainer.h>
33 
34 #include <fun4all/Fun4AllHistoManager.h>
35 #include <fun4all/Fun4AllReturnCodes.h>
36 
37 
38 #include "g4detectors/PHG4CylinderGeomContainer.h"
39 #include "g4detectors/PHG4CylinderGeom_Spacalv1.h" // for PHG4CylinderGeom_Spaca...
40 #include "g4detectors/PHG4CylinderGeom_Spacalv3.h"
41 #include "g4detectors/PHG4CylinderCellGeomContainer.h"
42 #include "g4detectors/PHG4CylinderCellGeom_Spacalv1.h"
43 
44 
45 #include <phool/getClass.h>
46 #include <TProfile.h>
47 #include <TFile.h>
48 #include <TNtuple.h>
49 #include <TMath.h>
50 #include <cassert>
51 #include <sstream>
52 #include <string>
53 #include <TF1.h>
54 
55 
56 //____________________________________________________________________________..
57 g4hitshift::g4hitshift(const std::string &name):
58  SubsysReco(name)
59 {
60 }
61 
62 //____________________________________________________________________________..
63 g4hitshift::~g4hitshift()
64 {
65 }
66 
67 //____________________________________________________________________________..
68 int g4hitshift::Init(PHCompositeNode *topNode)
69 {
70  return Fun4AllReturnCodes::EVENT_OK;
71 }
72 
73 //____________________________________________________________________________..
74 int g4hitshift::InitRun(PHCompositeNode *topNode)
75 {
76  return Fun4AllReturnCodes::EVENT_OK;
77 }
78 
79 //____________________________________________________________________________..
80 int g4hitshift::process_event(PHCompositeNode *topNode)
81 {
82  PHG4HitContainer* hits = findNode::getClass<PHG4HitContainer>(topNode, "G4HIT_CEMC");
83  if (hits)
84  {
85  //-----------------------------------------------------------------------
86  //Loop over G4Hits
87  //-----------------------------------------------------------------------
88  PHG4HitContainer::ConstRange hit_range = hits->getHits();
89  for (PHG4HitContainer::ConstIterator hit_iter = hit_range.first; hit_iter != hit_range.second; hit_iter++)
90  {
91  int scint_id = hit_iter->second->get_scint_id();
92 
93  PHG4CylinderGeom_Spacalv3::scint_id_coder decoder(scint_id);
94 
95  int sectornumber = decoder.sector_ID;
96  int shiftedsector = sectornumber + 8;
97  if (shiftedsector > 31)
98  {
99  shiftedsector = shiftedsector - 31;
100  }
101  PHG4CylinderGeom_Spacalv3::scint_id_coder encoder(shiftedsector,decoder.tower_ID,decoder.fiber_ID);
102 
103  hit_iter->second->set_scint_id(encoder.scint_ID);
104  }
105  }
106  return Fun4AllReturnCodes::EVENT_OK;
107 }
108 
109 //____________________________________________________________________________..
110 int g4hitshift::ResetEvent(PHCompositeNode *topNode)
111 {
112  return Fun4AllReturnCodes::EVENT_OK;
113 }
114 
115 //____________________________________________________________________________..
116 int g4hitshift::EndRun(const int runnumber)
117 {
118  return Fun4AllReturnCodes::EVENT_OK;
119 }
120 
121 //____________________________________________________________________________..
122 int g4hitshift::End(PHCompositeNode *topNode)
123 {
124  return Fun4AllReturnCodes::EVENT_OK;
125 }
126 
127 //____________________________________________________________________________..
128 int g4hitshift::Reset(PHCompositeNode *topNode)
129 {
130  return Fun4AllReturnCodes::EVENT_OK;
131 }
132 
133 //____________________________________________________________________________..
134 void g4hitshift::Print(const std::string &what) const
135 {
136 }