From 43fe149177589c3fb82a28e2ced30407ef85061b Mon Sep 17 00:00:00 2001 From: Aahana Brahma Date: Tue, 8 Oct 2024 10:54:27 -0400 Subject: [PATCH 1/3] Adding the photon kinematics for the EMCal clusters QA type plots --- offline/QA/Jet/photonjetskinematics.cc | 184 +++++++++++++++++++++++++ offline/QA/Jet/photonjetskinematics.h | 72 ++++++++++ 2 files changed, 256 insertions(+) create mode 100644 offline/QA/Jet/photonjetskinematics.cc create mode 100644 offline/QA/Jet/photonjetskinematics.h diff --git a/offline/QA/Jet/photonjetskinematics.cc b/offline/QA/Jet/photonjetskinematics.cc new file mode 100644 index 0000000000..ccadf447a8 --- /dev/null +++ b/offline/QA/Jet/photonjetskinematics.cc @@ -0,0 +1,184 @@ +//____________________________________________________________________________.. + +#include "photonjetskinematics.h" +#include + +#include + +#include + +#include +#include +#include "TH1.h" +#include "TMath.h" +//#include + + +// Tower includes +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +// Cluster includes +#include +#include + + +//____________________________________________________________________________.. +photonjetskinematics::photonjetskinematics(const std::string &outputfilename): +SubsysReco("photonjetskinematics") +, outfilename(outputfilename) +{ + std::cout << "photonjetskinematics::photonjetskinematics(const std::string &name, const std::string&outputfilename) Calling ctor" << std::endl; +} +//____________________________________________________________________________.. +photonjetskinematics::~photonjetskinematics() +{ + std::cout << "photonjetskinematics::~photonjetskinematics() Calling dtor" << std::endl; +} + +//____________________________________________________________________________.. +int photonjetskinematics::Init(PHCompositeNode *topNode) +{ + std::cout << "photonjetskinematics::Init(PHCompositeNode *topNode) Initializing" << std::endl; + //outfile = new TFile(outfilename.c_str(), "RECREATE"); + PHTFileServer::get().open(outfilename, "RECREATE"); + + //initializing histograms + + h_emcal_cluster_chi2 = new TH1D("h_emcal_cluster_chi2", "h_emcal_cluster_chi2", 30, 0, 150); //1000,0,5e6 + h_emcal_cluster_chi2->GetXaxis()->SetTitle("#chi^{2}"); + + + h_emcal_cluster_energy = new TH1D("h_emcal_cluster_energy", "h_emcal_cluster_energy", 100, 0, 10); + h_emcal_cluster_energy->GetXaxis()->SetTitle("E_{T} [GeV]"); + + + h_emcal_cluster_eta_phi = new TH2F("h_emcal_cluster_eta_phi", "h_emcal_cluster_eta_phi", 140, -1.2, 1.2, 64, -M_PI, M_PI); + //eta used to be 24 increase to 100 + //TH2D changed to TH2F + h_emcal_cluster_eta_phi->GetXaxis()->SetTitle("#eta"); + h_emcal_cluster_eta_phi->GetYaxis()->SetTitle("#Phi"); + h_emcal_cluster_eta_phi->SetOption("COLZ"); + + + h_emcal_cluster_eta = new TH1D("h_emcal_cluster_eta", "Eta Distribution", 140, -1.2, 1.2); + h_emcal_cluster_eta->GetXaxis()->SetTitle("#eta"); + + h_emcal_cluster_phi = new TH1D("h_emcal_cluster_phi", "Phi Distribution", 64, -M_PI, M_PI); + h_emcal_cluster_phi->GetXaxis()->SetTitle("#phi"); + + return Fun4AllReturnCodes::EVENT_OK; + +} + +//____________________________________________________________________________.. +int photonjetskinematics::InitRun(PHCompositeNode *topNode) +{ + std::cout << "photonjetskinematics::InitRun(PHCompositeNode *topNode) Initializing for Run XXX" << std::endl; + return Fun4AllReturnCodes::EVENT_OK; +} + +//____________________________________________________________________________.. +int photonjetskinematics::process_event(PHCompositeNode *topNode) +{ + + RawClusterContainer* clusterContainer = findNode::getClass(topNode, "CLUSTERINFO_CEMC"); + if (!clusterContainer) + { + std::cout << PHWHERE << "funkyCaloStuff::process_event - Fatal Error - CLUSTER_CEMC node is missing. " << std::endl; + return 0; + } + + RawClusterContainer::ConstIterator clusterIter; + RawClusterContainer::ConstRange clusterEnd = clusterContainer->getClusters(); + + for (clusterIter = clusterEnd.first; clusterIter != clusterEnd.second; clusterIter++) + { + RawCluster* recoCluster = clusterIter->second; + CLHEP::Hep3Vector vertex(0, 0, 0); + CLHEP::Hep3Vector E_vec_cluster = RawClusterUtility::GetECoreVec(*recoCluster, vertex); + + float clusE = E_vec_cluster.mag(); + float clus_eta = E_vec_cluster.pseudoRapidity(); + float clus_phi = E_vec_cluster.phi(); + float clus_chisq = recoCluster->get_chi2(); + + + //Filling Histograms, only taking into account E vec cluster, not reco cluster + + h_emcal_cluster_chi2->Fill(clus_chisq); + h_emcal_cluster_energy->Fill(clusE); + h_emcal_cluster_eta_phi->Fill(clus_eta, clus_phi); + h_emcal_cluster_eta->Fill(clus_eta); // 1D eta plot + h_emcal_cluster_phi->Fill(clus_phi); // 1D phi plot + + // std::cout << "E vec phi = " << clus_phi << std::endl; + + // std::cout << "E vec energy = " << clusE << std::endl; + + // std::cout << "E vec chi2 = " << clus_chisq << std::endl; + + // std::cout << "E vec eta = " << clus_eta << std::endl; + + } + + // std::cout << "photonjetskinematics::process_event(PHCompositeNode *topNode) Processing Event" << std::endl; + return Fun4AllReturnCodes::EVENT_OK; +} + +//____________________________________________________________________________.. +int photonjetskinematics::ResetEvent(PHCompositeNode *topNode) +{ + // std::cout << "photonjetskinematics::ResetEvent(PHCompositeNode *topNode) Resetting internal structures, prepare for next event" << std::endl; + return Fun4AllReturnCodes::EVENT_OK; +} + +//____________________________________________________________________________.. +int photonjetskinematics::EndRun(const int runnumber) +{ + std::cout << "photonjetskinematics::EndRun(const int runnumber) Ending Run for Run " << runnumber << std::endl; + return Fun4AllReturnCodes::EVENT_OK; +} + +//____________________________________________________________________________.. +int photonjetskinematics::End(PHCompositeNode *topNode) +{ + std::cout << "photonjetskinematics::End - Output to " << outfilename << std::endl; + PHTFileServer::get().cd(outfilename); + + //outfile->cd("photonjetskinematics"); + + //Outputting the histograms + h_emcal_cluster_eta_phi->Write("h_emcal_cluster_eta_phi"); + h_emcal_cluster_energy->Write("h_emcal_cluster_energy"); + h_emcal_cluster_chi2->Write("h_emcal_cluster_chi2"); + h_emcal_cluster_eta->Write("h_emcal_cluster_eta"); // Save 1D eta plot + h_emcal_cluster_phi->Write("h_emcal_cluster_phi"); // Save 1D phi plot + + std::cout << "photonjetskinematics::End(PHCompositeNode *topNode) This is the End..." << std::endl; + // outfile->Write(); + // outfile->Close(); + return Fun4AllReturnCodes::EVENT_OK; + +} + +//____________________________________________________________________________.. +int photonjetskinematics::Reset(PHCompositeNode *topNode) +{ + std::cout << "photonjetskinematics::Reset(PHCompositeNode *topNode) being Reset" << std::endl; + return Fun4AllReturnCodes::EVENT_OK; +} + +//____________________________________________________________________________.. +void photonjetskinematics::Print(const std::string &what) const +{ + std::cout << "photonjetskinematics::Print(const std::string &what) const Printing info for " << what << std::endl; +} diff --git a/offline/QA/Jet/photonjetskinematics.h b/offline/QA/Jet/photonjetskinematics.h new file mode 100644 index 0000000000..08b2320547 --- /dev/null +++ b/offline/QA/Jet/photonjetskinematics.h @@ -0,0 +1,72 @@ +// Tell emacs that this is a C++ source +// -*- C++ -*-. +#ifndef PHOTONJETSKINEMATICS_H +#define PHOTONJETSKINEMATICS_H + +#include + +#include +#include +#include +#include + +class PHCompositeNode; +class TH1; +class TH2; +class TH3; + +class photonjetskinematics : public SubsysReco +{ + public: + + + photonjetskinematics(const std::string &outputfilename = "output"); + ~photonjetskinematics() override; + + /** Called during initialization. + Typically this is where you can book histograms, and e.g. + register them to Fun4AllServer (so they can be output to file + using Fun4AllServer::dumpHistos() method). + */ + int Init(PHCompositeNode *topNode) override; + + /** Called for first event when run number is known. + Typically this is where you may want to fetch data from + database, because you know the run number. A place + to book histograms which have to know the run number. + */ + int InitRun(PHCompositeNode *topNode) override; + + /** Called for each event. + This is where you do the real work. + */ + int process_event(PHCompositeNode *topNode) override; + + /// Clean up internals after each event. + int ResetEvent(PHCompositeNode *topNode) override; + + /// Called at the end of each run. + int EndRun(const int runnumber) override; + + /// Called at the end of all processing. + int End(PHCompositeNode *topNode) override; + + /// Reset + int Reset(PHCompositeNode * /*topNode*/) override; + + void Print(const std::string &what = "ALL") const override; + + private: +std::string outfilename; +TFile *outfile; + + +//Output histograms + TH1 *h_emcal_cluster_chi2 = nullptr; + TH1 *h_emcal_cluster_energy = nullptr; + TH2 *h_emcal_cluster_eta_phi = nullptr; + TH1 *h_emcal_cluster_eta = nullptr; // <-- Declare eta histogram + TH1 *h_emcal_cluster_phi = nullptr; // <-- Declare phi histogram +}; + +#endif // PHOTONJETSKINEMATICS_H From 3646d30039ff4d912afcfb7412b59222afc27d9a Mon Sep 17 00:00:00 2001 From: Aahana Brahma Date: Mon, 14 Oct 2024 17:54:43 -0400 Subject: [PATCH 2/3] Added the QA histogram manager and some other minor changes that was suggested by Derek. --- offline/QA/Jet/Makefile.am | 6 +- offline/QA/Jet/photonjetskinematics.cc | 141 +++++++++++++++---------- offline/QA/Jet/photonjetskinematics.h | 31 +++++- 3 files changed, 118 insertions(+), 60 deletions(-) diff --git a/offline/QA/Jet/Makefile.am b/offline/QA/Jet/Makefile.am index 1fb56cca53..c4180a1055 100644 --- a/offline/QA/Jet/Makefile.am +++ b/offline/QA/Jet/Makefile.am @@ -31,7 +31,8 @@ pkginclude_HEADERS = \ TrksInJetQAJetManager.h \ RhosinEvent.h \ JetSeedCount.h \ - JetKinematicCheck.h + JetKinematicCheck.h \ + photonjetskinematics.h libjetqa_la_SOURCES = \ ConstituentsinJets.cc\ @@ -47,7 +48,8 @@ libjetqa_la_SOURCES = \ TrksInJetQAJetManager.cc \ RhosinEvent.cc \ JetSeedCount.cc \ - JetKinematicCheck.cc + JetKinematicCheck.cc \ + photonjetskinematics.cc libjetqa_la_LIBADD = \ -L$(libdir) \ diff --git a/offline/QA/Jet/photonjetskinematics.cc b/offline/QA/Jet/photonjetskinematics.cc index ccadf447a8..0fa2fba144 100644 --- a/offline/QA/Jet/photonjetskinematics.cc +++ b/offline/QA/Jet/photonjetskinematics.cc @@ -3,6 +3,10 @@ #include "photonjetskinematics.h" #include +// Fun4all includes + +#include + #include #include @@ -30,13 +34,25 @@ #include #include +// QA includes +#include +// c++ includes +#include +#include +#include +#include //____________________________________________________________________________.. -photonjetskinematics::photonjetskinematics(const std::string &outputfilename): -SubsysReco("photonjetskinematics") -, outfilename(outputfilename) + +photonjetskinematics::photonjetskinematics(const std::string &m_modulename, const std::string &m_inputnode) : + SubsysReco(m_modulename) + , modulename(m_modulename) + , inputnode(m_inputnode) + , histtag("AllTrig") + , trgToSelect(JetQADefs::GL1::MBDNSJet1) + , doTrgSelect(false) { - std::cout << "photonjetskinematics::photonjetskinematics(const std::string &name, const std::string&outputfilename) Calling ctor" << std::endl; + if (Verbosity() > 1) std::cout << "photonjetskinematics::photonjetskinematics(const std::string &name, const std::string&outputfilename) Calling ctor" << std::endl; } //____________________________________________________________________________.. photonjetskinematics::~photonjetskinematics() @@ -45,34 +61,57 @@ photonjetskinematics::~photonjetskinematics() } //____________________________________________________________________________.. -int photonjetskinematics::Init(PHCompositeNode *topNode) +int photonjetskinematics::Init(PHCompositeNode* /*topNode*/) { - std::cout << "photonjetskinematics::Init(PHCompositeNode *topNode) Initializing" << std::endl; - //outfile = new TFile(outfilename.c_str(), "RECREATE"); - PHTFileServer::get().open(outfilename, "RECREATE"); + if (Verbosity() > 1) std::cout << "photonjetskinematics::Init(PHCompositeNode *topNode) Initializing" << std::endl; + manager = QAHistManagerDef::getHistoManager(); + if (!manager) + { + std::cerr << PHWHERE << "PANIC: couldn't grab histogram manager!" << std::endl; + assert(manager); + } + + // make sure module name is lower case + std::string smallModuleName = modulename; + std::transform( + smallModuleName.begin(), + smallModuleName.end(), + smallModuleName.begin(), + ::tolower); + // construct histogram names + std::vector vecHistNames = { + "emcal_cluster_chi2", + "emcal_cluster_energy", + "emcal_cluster_eta_phi", + "emcal_cluster_eta", + "emcal_cluster_phi"}; + + for (auto& histName : vecHistNames) + { + histName.insert(0, "h_" + smallModuleName + "_"); + if (!histtag.empty()) + { + histName.append("_" + histtag); + } + } //initializing histograms - h_emcal_cluster_chi2 = new TH1D("h_emcal_cluster_chi2", "h_emcal_cluster_chi2", 30, 0, 150); //1000,0,5e6 + h_emcal_cluster_chi2 = new TH1D(vecHistNames[0].data(), "h_emcal_cluster_chi2", 30, 0, 150); h_emcal_cluster_chi2->GetXaxis()->SetTitle("#chi^{2}"); - - h_emcal_cluster_energy = new TH1D("h_emcal_cluster_energy", "h_emcal_cluster_energy", 100, 0, 10); + h_emcal_cluster_energy = new TH1D(vecHistNames[1].data(), "h_emcal_cluster_energy", 100, 0, 10); h_emcal_cluster_energy->GetXaxis()->SetTitle("E_{T} [GeV]"); - - h_emcal_cluster_eta_phi = new TH2F("h_emcal_cluster_eta_phi", "h_emcal_cluster_eta_phi", 140, -1.2, 1.2, 64, -M_PI, M_PI); - //eta used to be 24 increase to 100 - //TH2D changed to TH2F + h_emcal_cluster_eta_phi = new TH2F(vecHistNames[2].data(), "h_emcal_cluster_eta_phi", 140, -1.2, 1.2, 64, -M_PI, M_PI); h_emcal_cluster_eta_phi->GetXaxis()->SetTitle("#eta"); h_emcal_cluster_eta_phi->GetYaxis()->SetTitle("#Phi"); h_emcal_cluster_eta_phi->SetOption("COLZ"); - - h_emcal_cluster_eta = new TH1D("h_emcal_cluster_eta", "Eta Distribution", 140, -1.2, 1.2); + h_emcal_cluster_eta = new TH1D(vecHistNames[3].data(), "Eta Distribution", 140, -1.2, 1.2); h_emcal_cluster_eta->GetXaxis()->SetTitle("#eta"); - h_emcal_cluster_phi = new TH1D("h_emcal_cluster_phi", "Phi Distribution", 64, -M_PI, M_PI); + h_emcal_cluster_phi = new TH1D(vecHistNames[4].data(), "Phi Distribution", 64, -M_PI, M_PI); h_emcal_cluster_phi->GetXaxis()->SetTitle("#phi"); return Fun4AllReturnCodes::EVENT_OK; @@ -80,9 +119,10 @@ int photonjetskinematics::Init(PHCompositeNode *topNode) } //____________________________________________________________________________.. -int photonjetskinematics::InitRun(PHCompositeNode *topNode) +int photonjetskinematics::InitRun(PHCompositeNode* /*topNode*/) { - std::cout << "photonjetskinematics::InitRun(PHCompositeNode *topNode) Initializing for Run XXX" << std::endl; + + if (Verbosity() > 1) std::cout << "photonjetskinematics::InitRun(PHCompositeNode *topNode) Initializing for Run XXX" << std::endl; return Fun4AllReturnCodes::EVENT_OK; } @@ -90,13 +130,23 @@ int photonjetskinematics::InitRun(PHCompositeNode *topNode) int photonjetskinematics::process_event(PHCompositeNode *topNode) { - RawClusterContainer* clusterContainer = findNode::getClass(topNode, "CLUSTERINFO_CEMC"); + RawClusterContainer* clusterContainer = findNode::getClass(topNode, inputnode); if (!clusterContainer) { std::cout << PHWHERE << "funkyCaloStuff::process_event - Fatal Error - CLUSTER_CEMC node is missing. " << std::endl; return 0; } + // if needed, check if trigger fired + if (doTrgSelect) + { + bool hasTrigger = JetQADefs::DidTriggerFire(trgToSelect, topNode); + if (!hasTrigger) + { + return Fun4AllReturnCodes::EVENT_OK; + } + } + RawClusterContainer::ConstIterator clusterIter; RawClusterContainer::ConstRange clusterEnd = clusterContainer->getClusters(); @@ -111,7 +161,6 @@ int photonjetskinematics::process_event(PHCompositeNode *topNode) float clus_phi = E_vec_cluster.phi(); float clus_chisq = recoCluster->get_chi2(); - //Filling Histograms, only taking into account E vec cluster, not reco cluster h_emcal_cluster_chi2->Fill(clus_chisq); @@ -119,23 +168,14 @@ int photonjetskinematics::process_event(PHCompositeNode *topNode) h_emcal_cluster_eta_phi->Fill(clus_eta, clus_phi); h_emcal_cluster_eta->Fill(clus_eta); // 1D eta plot h_emcal_cluster_phi->Fill(clus_phi); // 1D phi plot - - // std::cout << "E vec phi = " << clus_phi << std::endl; - - // std::cout << "E vec energy = " << clusE << std::endl; - - // std::cout << "E vec chi2 = " << clus_chisq << std::endl; - - // std::cout << "E vec eta = " << clus_eta << std::endl; - - } + } // std::cout << "photonjetskinematics::process_event(PHCompositeNode *topNode) Processing Event" << std::endl; - return Fun4AllReturnCodes::EVENT_OK; + return Fun4AllReturnCodes::EVENT_OK; } //____________________________________________________________________________.. -int photonjetskinematics::ResetEvent(PHCompositeNode *topNode) +int photonjetskinematics::ResetEvent(PHCompositeNode* /*topNode*/) { // std::cout << "photonjetskinematics::ResetEvent(PHCompositeNode *topNode) Resetting internal structures, prepare for next event" << std::endl; return Fun4AllReturnCodes::EVENT_OK; @@ -144,41 +184,36 @@ int photonjetskinematics::ResetEvent(PHCompositeNode *topNode) //____________________________________________________________________________.. int photonjetskinematics::EndRun(const int runnumber) { - std::cout << "photonjetskinematics::EndRun(const int runnumber) Ending Run for Run " << runnumber << std::endl; + if (Verbosity() > 1) std::cout << "photonjetskinematics::EndRun(const int runnumber) Ending Run for Run " << runnumber << std::endl; return Fun4AllReturnCodes::EVENT_OK; } //____________________________________________________________________________.. -int photonjetskinematics::End(PHCompositeNode *topNode) +int photonjetskinematics::End(PHCompositeNode* /*topNode*/) { - std::cout << "photonjetskinematics::End - Output to " << outfilename << std::endl; - PHTFileServer::get().cd(outfilename); - - //outfile->cd("photonjetskinematics"); - + if (Verbosity() > 1) std::cout << "photonjetskinematics::End - Output to " << outfilename << std::endl; + //Outputting the histograms - h_emcal_cluster_eta_phi->Write("h_emcal_cluster_eta_phi"); - h_emcal_cluster_energy->Write("h_emcal_cluster_energy"); - h_emcal_cluster_chi2->Write("h_emcal_cluster_chi2"); - h_emcal_cluster_eta->Write("h_emcal_cluster_eta"); // Save 1D eta plot - h_emcal_cluster_phi->Write("h_emcal_cluster_phi"); // Save 1D phi plot - - std::cout << "photonjetskinematics::End(PHCompositeNode *topNode) This is the End..." << std::endl; - // outfile->Write(); - // outfile->Close(); + manager->registerHisto(h_emcal_cluster_eta_phi); + manager->registerHisto(h_emcal_cluster_energy); + manager->registerHisto(h_emcal_cluster_chi2); + manager->registerHisto(h_emcal_cluster_eta); + manager->registerHisto(h_emcal_cluster_phi); + + if (Verbosity() > 1) std::cout << "photonjetskinematics::End(PHCompositeNode *topNode) This is the End..." << std::endl; return Fun4AllReturnCodes::EVENT_OK; } //____________________________________________________________________________.. -int photonjetskinematics::Reset(PHCompositeNode *topNode) +int photonjetskinematics::Reset(PHCompositeNode* /*topNode*/) { - std::cout << "photonjetskinematics::Reset(PHCompositeNode *topNode) being Reset" << std::endl; + if (Verbosity() > 1) std::cout << "photonjetskinematics::Reset(PHCompositeNode *topNode) being Reset" << std::endl; return Fun4AllReturnCodes::EVENT_OK; } //____________________________________________________________________________.. void photonjetskinematics::Print(const std::string &what) const { - std::cout << "photonjetskinematics::Print(const std::string &what) const Printing info for " << what << std::endl; + if (Verbosity() > 1) std::cout << "photonjetskinematics::Print(const std::string &what) const Printing info for " << what << std::endl; } diff --git a/offline/QA/Jet/photonjetskinematics.h b/offline/QA/Jet/photonjetskinematics.h index 08b2320547..9306209f32 100644 --- a/offline/QA/Jet/photonjetskinematics.h +++ b/offline/QA/Jet/photonjetskinematics.h @@ -9,7 +9,9 @@ #include #include #include +#include "JetQADefs.h" +class Fun4AllHistoManager; class PHCompositeNode; class TH1; class TH2; @@ -20,7 +22,7 @@ class photonjetskinematics : public SubsysReco public: - photonjetskinematics(const std::string &outputfilename = "output"); + photonjetskinematics(const std::string &modulename = "photonjetskinematics", const std::string &inputnode = "CLUSTERINFO_CEMC"); ~photonjetskinematics() override; /** Called during initialization. @@ -56,12 +58,31 @@ class photonjetskinematics : public SubsysReco void Print(const std::string &what = "ALL") const override; - private: -std::string outfilename; -TFile *outfile; + /// specifies a trigger to select + void SetTrgToSelect(const uint32_t trig = JetQADefs::GL1::MBDNSPhoton1) + { + doTrgSelect = true; + trgToSelect = trig; + } + /// set histogram tag + void SetHistTag(const std::string& tag) + { + histtag = tag; + } -//Output histograms + private: + std::string outfilename; + TFile *outfile; + // hist manager + Fun4AllHistoManager* manager; + std::string modulename; + std::string inputnode; + std::string histtag; + uint32_t trgToSelect; + bool doTrgSelect; + + ///Output histograms TH1 *h_emcal_cluster_chi2 = nullptr; TH1 *h_emcal_cluster_energy = nullptr; TH2 *h_emcal_cluster_eta_phi = nullptr; From 4a11793dcfe19ef80d6a217922717b8a5ae5b338 Mon Sep 17 00:00:00 2001 From: Aahana Brahma Date: Thu, 14 Nov 2024 10:08:41 -0500 Subject: [PATCH 3/3] Modified the outfile statements that Jenkins failed and put in the build report. --- offline/QA/Jet/photonjetskinematics.cc | 3 ++- offline/QA/Jet/photonjetskinematics.h | 5 +++-- 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/offline/QA/Jet/photonjetskinematics.cc b/offline/QA/Jet/photonjetskinematics.cc index 0fa2fba144..70cb304dc9 100644 --- a/offline/QA/Jet/photonjetskinematics.cc +++ b/offline/QA/Jet/photonjetskinematics.cc @@ -191,7 +191,8 @@ int photonjetskinematics::EndRun(const int runnumber) //____________________________________________________________________________.. int photonjetskinematics::End(PHCompositeNode* /*topNode*/) { - if (Verbosity() > 1) std::cout << "photonjetskinematics::End - Output to " << outfilename << std::endl; + // Commenting out the following line because Jenkins keeps failing the build test :( + // if (Verbosity() > 1) std::cout << "photonjetskinematics::End - Output to " << outfilename << std::endl; //Outputting the histograms manager->registerHisto(h_emcal_cluster_eta_phi); diff --git a/offline/QA/Jet/photonjetskinematics.h b/offline/QA/Jet/photonjetskinematics.h index 9306209f32..3821a9c7d7 100644 --- a/offline/QA/Jet/photonjetskinematics.h +++ b/offline/QA/Jet/photonjetskinematics.h @@ -72,9 +72,10 @@ class photonjetskinematics : public SubsysReco } private: - std::string outfilename; - TFile *outfile; + // std::string outfilename; + // TFile *outfile; // hist manager + Fun4AllHistoManager* manager; std::string modulename; std::string inputnode;