Tpc hit set reco

calibrations/calorimeter/calo_emc_pi0_tbt/CaloCalibEmc_Pi0.cc

@@ -5,15 +5,13 @@
 #include <calobase/RawClusterUtility.h>
 #include <calobase/RawTower.h>
 #include <calobase/RawTowerContainer.h>
+#include <calobase/RawTowerDefs.h>
 #include <calobase/RawTowerGeomContainer.h>
-#include <calobase/RawTowerGeom.h>
+#include <calobase/TowerInfoContainer.h>
+#include <calobase/TowerInfo.h>
 
-#include <calobase/TowerInfoContainerv1.h>
-#include <calobase/TowerInfov1.h>
-
-
-#include <g4vertex/GlobalVertex.h>
-#include <g4vertex/GlobalVertexMap.h>
+#include <globalvertex/GlobalVertex.h>
+#include <globalvertex/GlobalVertexMap.h>
 
 #include <fun4all/Fun4AllReturnCodes.h>
 #include <fun4all/SubsysReco.h>

calibrations/calorimeter/calo_emc_pi0_tbt/Makefile.am

@@ -13,7 +13,7 @@ AM_LDFLAGS = \
   -L$(OFFLINE_MAIN)/lib64 \
   -lHepMC \
   -lCLHEP \
-  -lg4vertex_io \
+  -lglobalvertex_io \
   -lcalo_io \
   -lSubsysReco
 

simulation/g4simulation/tpcresponse/Makefile.am → calibrations/tpc/TPCPedestalCalibration/Makefile.am

@@ -8,33 +8,28 @@ AM_CPPFLAGS = \
 AM_LDFLAGS = \
   -L$(libdir) \
   -L$(OFFLINE_MAIN)/lib \
-  -L$(OFFLINE_MAIN)/lib64 \
-  `root-config --libs`
+  -L$(OFFLINE_MAIN)/lib64
 
 pkginclude_HEADERS = \
-  TpcRS.h
+  TPCPedestalCalibration.h
 
 lib_LTLIBRARIES = \
-  libtpcresponse.la
+  libTPCPedestalCalibration.la
 
-libtpcresponse_la_SOURCES = \
-  TpcRS.cc
+libTPCPedestalCalibration_la_SOURCES = \
+  TPCPedestalCalibration.cc
 
-libtpcresponse_la_LIBADD = \
-  -lphg4hit \
+libTPCPedestalCalibration_la_LIBADD = \
   -lphool \
-  -lSubsysReco \
-  -ltpcrs \
-  -lGeom \
-  -lMathMore
+  -lSubsysReco
 
 BUILT_SOURCES = testexternals.cc
 
 noinst_PROGRAMS = \
   testexternals
 
 testexternals_SOURCES = testexternals.cc
-testexternals_LDADD   = libtpcresponse.la
+testexternals_LDADD   = libTPCPedestalCalibration.la
 
 testexternals.cc:
 	echo "//*** this is a generated file. Do not commit, do not edit" > $@

calibrations/tpc/TPCPedestalCalibration/TPCPedestalCalibration.cc

@@ -0,0 +1,197 @@
+#include "TPCPedestalCalibration.h"
+
+#include <fun4all/Fun4AllReturnCodes.h>
+#include <phool/PHCompositeNode.h>
+#include <phool/PHIODataNode.h>    // for PHIODataNode
+#include <phool/PHNodeIterator.h>  // for PHNodeIterator
+#include <phool/PHObject.h>        // for PHObject
+#include <phool/getClass.h>
+#include <phool/phool.h>
+
+#include <Event/Event.h>
+#include <Event/EventTypes.h>
+#include <Event/packet.h>
+
+#include <TMath.h>
+#include <TFile.h>
+#include <TTree.h>
+
+#include <memory>
+#include <cassert>
+#include <iostream>
+
+/*************************************************************/
+/*                TPC Pedestal Calibration                   */
+/*               Thomas Marshall,Aditya Dash                 */
+/*        rosstom@ucla.edu,aditya55@physics.ucla.edu         */
+/*************************************************************/
+
+TPCPedestalCalibration::TPCPedestalCalibration(const std::string &name)
+ :SubsysReco("TPCPedestalCalibration")
+ , m_fname(name)
+{
+  // reserve memory for max ADC samples
+  m_adcSamples.resize(1024, 0);
+
+  for(int fee_no=0;fee_no<26;fee_no++)
+  {
+    for(int channel_no=0;channel_no<256;channel_no++)
+    {
+      m_aveADCFeeChannel[fee_no][channel_no]=0.0;
+      m_stdADCFeeChannel[fee_no][channel_no]=0.0;
+      m_countsADCFeeChannel[fee_no][channel_no]=0.0;
+      m_aliveArrayFeeChannel[fee_no][channel_no]=1;
+    }
+  }
+
+}
+
+int TPCPedestalCalibration::InitRun(PHCompositeNode *)
+{
+  m_file = TFile::Open(m_fname.c_str(), "recreate");
+  assert(m_file->IsOpen());
+
+  m_pedestalTree = new TTree("pedestalTree", "Each entry is one TPC Channel");
+
+  m_pedestalTree->Branch("isAlive",&m_isAlive,"isAlive/I");
+  m_pedestalTree->Branch("pedMean",&m_pedMean,"pedMean/F");
+  m_pedestalTree->Branch("pedStd",&m_pedStd,"pedStd/F");
+  m_pedestalTree->Branch("sector",&m_sector,"sector/I");
+  m_pedestalTree->Branch("fee",&m_outFEE,"fee/I");
+  m_pedestalTree->Branch("channel",&m_chan,"channel/I");
+  m_pedestalTree->Branch("module",&m_module,"module/I");
+  m_pedestalTree->Branch("slot",&m_slot,"slot/I");
+
+  return Fun4AllReturnCodes::EVENT_OK;
+}
+
+//____________________________________________________________________________..
+int TPCPedestalCalibration::process_event(PHCompositeNode *topNode)
+{
+  Event *_event = findNode::getClass<Event>(topNode, "PRDF");
+  if (_event == nullptr)
+  {
+    std::cout << "TPCRawDataTree::Process_Event - Event not found" << std::endl;
+    return -1;
+  }
+  if (_event->getEvtType() >= 8)  /// special events
+  {
+    return Fun4AllReturnCodes::DISCARDEVENT;
+  }
+
+  for (int packet : m_packets)
+  {
+    if (Verbosity())
+    {
+      std::cout << __PRETTY_FUNCTION__ << " : decoding packet " << packet << std::endl;
+    }
+
+    m_packet = packet;
+
+    std::unique_ptr<Packet> p (_event->getPacket(m_packet));
+    if (!p)
+    {
+      if (Verbosity())
+      {
+        std::cout << __PRETTY_FUNCTION__ << " : missing packet " << packet << std::endl;
+      }
+
+      continue;
+    }
+
+    m_nWaveormInFrame = p->iValue(0, "NR_WF");
+
+    for (int wf = 0; wf < m_nWaveormInFrame; wf++)
+    {
+      m_nSamples = p->iValue(wf, "SAMPLES");
+      m_fee = p->iValue(wf, "FEE");
+      m_Channel = p->iValue(wf, "CHANNEL");
+
+      if(m_nSamples==0)
+      {
+        m_aliveArrayFeeChannel[m_fee][m_Channel]=0;
+        continue;
+      }
+
+      assert(m_nSamples < (int) m_adcSamples.size());  // no need for movements in memory allocation
+      for (int s = 0; s < m_nSamples; s++)
+      {
+        m_adcSamples[s] = p->iValue(wf, s);
+      }
+
+      bool dead = false;
+      for(int adc_sam_no=0;adc_sam_no<m_nSamples;adc_sam_no++)
+      {
+        if (m_adcSamples[adc_sam_no] == 0 || TMath::IsNaN(float(m_adcSamples[adc_sam_no])))
+        {
+          m_aliveArrayFeeChannel[m_fee][m_Channel]=0;
+        }
+      }
+      if (dead) {continue;}
+
+      for(int adc_sam_no=0;adc_sam_no<m_nSamples;adc_sam_no++){
+        m_aveADCFeeChannel[m_fee][m_Channel]+=m_adcSamples[adc_sam_no];
+        m_stdADCFeeChannel[m_fee][m_Channel]+=pow(m_adcSamples[adc_sam_no],2);
+        m_countsADCFeeChannel[m_fee][m_Channel]+=1.0;
+      }
+    }
+  }  //   for (int packet : m_packets)
+
+  return Fun4AllReturnCodes::EVENT_OK;
+}
+
+int TPCPedestalCalibration::EndRun(const int runnumber)
+{
+  std::cout << "TPCPedestalCalibration::EndRun(const int runnumber) Ending Run for Run " << runnumber << std::endl;
+
+  for(int fee_no=0;fee_no<26;fee_no++)
+  {
+    for(int channel_no=0;channel_no<256;channel_no++)
+    {
+      if(m_countsADCFeeChannel[fee_no][channel_no] != 0.0)
+      {
+        float temp1 = m_aveADCFeeChannel[fee_no][channel_no]/m_countsADCFeeChannel[fee_no][channel_no];
+        float temp2 = m_stdADCFeeChannel[fee_no][channel_no]/m_countsADCFeeChannel[fee_no][channel_no];
+        m_aveADCFeeChannel[fee_no][channel_no] = temp1;
+        m_stdADCFeeChannel[fee_no][channel_no] = temp2;
+      }
+      else
+      {
+        m_aveADCFeeChannel[fee_no][channel_no] = 0.0;
+        m_stdADCFeeChannel[fee_no][channel_no] = 0.0;
+        m_aliveArrayFeeChannel[fee_no][channel_no]=0;
+      }
+
+      if(m_aveADCFeeChannel[fee_no][channel_no] > 200 || m_aveADCFeeChannel[fee_no][channel_no] < 10)
+      {
+        m_aliveArrayFeeChannel[fee_no][channel_no]=0;
+      }
+
+      m_pedMean=m_aveADCFeeChannel[fee_no][channel_no];
+      m_pedStd=sqrt(m_stdADCFeeChannel[fee_no][channel_no] - pow(m_aveADCFeeChannel[fee_no][channel_no],2));
+      m_isAlive=m_aliveArrayFeeChannel[fee_no][channel_no];
+      m_chan=channel_no;
+      m_outFEE=fee_no;
+      m_module=mod_arr[fee_no];
+      m_slot=slot_arr[fee_no];
+      m_pedestalTree->Fill();
+    }
+  }
+
+  return Fun4AllReturnCodes::EVENT_OK;
+}
+
+//____________________________________________________________________________..
+int TPCPedestalCalibration::End(PHCompositeNode *topNode)
+{
+  std::cout << "TPCPedestalCalibration::End(PHCompositeNode *topNode) This is the End..." << std::endl;
+
+  m_file->Write();
+
+  std::cout << __PRETTY_FUNCTION__ << " : completed saving to " << m_file->GetName() << std::endl;
+  m_file->ls();
+
+  m_file->Close();  
+
+  return Fun4AllReturnCodes::EVENT_OK;
+}

calibrations/tpc/TPCPedestalCalibration/TPCPedestalCalibration.h

@@ -0,0 +1,74 @@
+// Tell emacs that this is a C++ source
+//  -*- C++ -*-.
+#ifndef TPCPedestalCalibration_H
+#define TPCPedestalCalibration_H
+
+#include <fun4all/SubsysReco.h>
+
+#include <string>
+#include <vector>
+
+class PHCompositeNode;
+class TFile;
+class TTree;
+
+class TPCPedestalCalibration : public SubsysReco
+{
+ public:
+  explicit TPCPedestalCalibration(const std::string &name = "TPCPedestalCalibration.root");
+
+  ~TPCPedestalCalibration() override {}
+
+  int InitRun(PHCompositeNode *topNode) override;
+
+  int process_event(PHCompositeNode *topNode) override;
+
+  int EndRun(const int runnumber) override;
+
+  int End(PHCompositeNode *topNode) override;
+
+  void AddPacket(int packet)
+  {
+    m_packets.push_back(packet);
+  }
+
+  void SetSector(int sectorNum)
+  {
+    m_sector = sectorNum;
+  }
+
+ protected:
+  //! which packet to decode
+  std::vector<int> m_packets{1001};
+
+ private:
+  std::string m_fname;
+  TFile * m_file = nullptr;
+  TTree * m_pedestalTree = nullptr;
+
+  int m_packet = 0;
+  int m_nWaveormInFrame = 0;
+  int m_nSamples = 0;
+  int m_fee = 0;
+  int m_Channel = 0;
+  std::vector<unsigned short> m_adcSamples; 
+
+  float m_aveADCFeeChannel[26][256];
+  float m_stdADCFeeChannel[26][256];
+  float m_countsADCFeeChannel[26][256];
+  int m_aliveArrayFeeChannel[26][256];
+
+  int m_isAlive = 1;
+  float m_pedMean = 0;
+  float m_pedStd = 0;
+  int m_sector = 0;
+  int m_outFEE = 99;
+  int m_chan = 999;
+  int m_module = 9;
+  int m_slot = 99;
+
+  int mod_arr[26]={2,2,1,1,1,3,3,3,3,3,3,2,2,1,2,2,1,1,2,2,3,3,3,3,3,3};
+  int slot_arr[26] = {5,6,1,3,2,12,10,11,9,8,7,1,2,4,8,7,6,5,4,3,1,3,2,4,6,5};
+};
+
+#endif // TPCPedestalCalibration_H

simulation/g4simulation/tpcresponse/configure.ac → calibrations/tpc/TPCPedestalCalibration/configure.ac

@@ -1,4 +1,4 @@
-AC_INIT(tpcresponse,[1.00])
+AC_INIT(tpcpedestalcalibration,[1.00])
 AC_CONFIG_SRCDIR([configure.ac])
 
 AM_INIT_AUTOMAKE
@@ -9,7 +9,7 @@ LT_INIT([disable-static])
 dnl   no point in suppressing warnings people should 
 dnl   at least see them, so here we go for g++: -Wall
 if test $ac_cv_prog_gxx = yes; then
-   CXXFLAGS="$CXXFLAGS -Wall -Werror -Wextra"
+   CXXFLAGS="$CXXFLAGS -Wall -Werror"
 fi
 
 AC_CONFIG_FILES([Makefile])