waveformexists added. testing in progress

Control/CalypsoExample/ClusterReadExample/python/ClusterReadExampleConfig.py

@@ -33,7 +33,7 @@ if __name__ == "__main__":
     ConfigFlags.IOVDb.GlobalTag = "OFLCOND-XXXX-XXX-XX"          # Needed to bypass autoconfig, only the "OFLCOND" matters at the moment
     ConfigFlags.GeoModel.FaserVersion     = "FASER-01"           # Default FASER geometry
     #ConfigFlags.Detector.SimulateFaserSCT = True
-    ConfigFlags.Input.Files = ["my.COSMICHITS.pool.root"]
+    ConfigFlags.Input.Files = ["run001332.ESD.pool.root"]#["my.COSMICHITS.pool.root"]
     ConfigFlags.lock()
 
 # Configure components

Control/CalypsoExample/ClusterReadExample/src/ClusterReadAlg.cxx

@@ -22,7 +22,16 @@
 #include <TTree.h>
 #include <TBranch.h>
 
+/*
+To add as new branches:
+POOLCollectiontree - IsSimulation (0 or 1)
+- event number
+- event time
+check issimulation before continuing
+- coords, energy deposit, etc
+- wtf is this doing that isnt done by clusterfit again...?
 
+*/
 
 
 //static const std::string moduleFailureReason{"ClusterReadAlg: Exceeds max clusters"};
@@ -37,6 +46,7 @@ ClusterReadAlg::ClusterReadAlg(const std::string& name, ISvcLocator* pSvcLocator
 // Initialize method:
 StatusCode ClusterReadAlg::initialize() {
   ATH_MSG_INFO("ClusterReadAlg::initialize()");
+  std::cout<<m_numberOfEvents;
 
   ATH_CHECK(m_faserTriggerDataKey.initialize());
   ATH_CHECK(m_clusterContainerKey.initialize());
@@ -56,6 +66,7 @@ StatusCode ClusterReadAlg::initialize() {
   m_tree->Branch("pull", &m_pull, "pull/F");
   m_tree->Branch("chi2", &m_refitChi2, "chi2/F");
   m_tree->Branch("mc_Event", &m_mcEventKey2, "chi2/F");
+  m_tree->Branch("waveformExists", &m_waveformExists, "waveform/B");
 
 
   ATH_CHECK(histSvc()->regHist("/HIST/chi2", m_chi2));
@@ -76,6 +87,7 @@ StatusCode ClusterReadAlg::initialize() {
 // Execute method:
 StatusCode ClusterReadAlg::execute(const EventContext& ctx) const
 {
+
   ++m_numberOfEvents;  
   setFilterPassed(false, ctx);
 
@@ -109,6 +121,8 @@ StatusCode ClusterReadAlg::execute(const EventContext& ctx) const
   SG::ReadHandle<FaserSCT_RDO_Container> h_sctRDO(m_faserRdoKey);
   //SG::ReadHandle<TrackerSimDataCollection> h_collectionMap(m_sctMap);
 
+  m_waveformExists = true;
+
   m_tree->Fill();
 
   // Done

Control/CalypsoExample/ClusterReadExample/src/ClusterReadAlg.h

@@ -107,6 +107,7 @@ class ClusterReadAlg : public AthReentrantAlgorithm, AthHistogramming
     mutable float m_pull;
     mutable float m_refitChi2;
     mutable float m_mcEventKey2;
+    mutable bool  m_waveformExists;
 
     mutable std::atomic<int> m_numberOfEvents{0};
     mutable std::atomic<int> m_numberOfTriggeredEvents{0};