3 #include <calobase/RawTower.h>
4 #include <calobase/RawTowerContainer.h>
5 #include <calobase/RawTowerDefs.h>
6 #include <calobase/RawTowerGeom.h>
7 #include <calobase/RawTowerGeomContainer.h>
8 #include <calobase/RawTowerv1.h>
44 std::cout <<
"RetowerCEMC::process_event: entering" << std::endl;
48 RawTowerContainer *towersEM3 = findNode::getClass<RawTowerContainer>(topNode,
"TOWER_CALIB_CEMC");
51 RawTowerGeomContainer *geomIH = findNode::getClass<RawTowerGeomContainer>(topNode,
"TOWERGEOM_HCALIN");
55 std::cout <<
"RetowerCEMC::process_event: " << towersEM3->
size() <<
" TOWER_CALIB_CEMC towers" << std::endl;
62 _NETA = geomIH->get_etabins();
63 _NPHI = geomIH->get_phibins();
84 int this_IHetabin = geomIH->get_etabin(tower_geom->
get_eta());
85 int this_IHphibin = geomIH->get_phibin(tower_geom->
get_phi());
91 RawTowerContainer *emcal_retower = findNode::getClass<RawTowerContainer>(topNode,
"TOWER_CALIB_CEMC_RETOWER");
93 if (
Verbosity() > 0) std::cout <<
"RetowerCEMC::process_event: filling TOWER_CALIB_CEMC_RETOWER node, with initial size = " << emcal_retower->
size() << std::endl;
107 if (
Verbosity() > 0) std::cout <<
"RetowerCEMC::process_event: finished filling TOWER_CALIB_CEMC_RETOWER node, with final size = " << emcal_retower->
size() << std::endl;
109 if (
Verbosity() > 0) std::cout <<
"RetowerCEMC::process_event: exiting" << std::endl;
122 std::cout <<
PHWHERE <<
"DST Node missing, doing nothing." << std::endl;
130 std::cout <<
PHWHERE <<
"EMCal Node note found, doing nothing." << std::endl;
133 RawTowerContainer *test_emcal_retower = findNode::getClass<RawTowerContainer>(topNode,
"TOWER_CALIB_CEMC_RETOWER");
134 if (!test_emcal_retower)
136 if (
Verbosity() > 0) std::cout <<
"RetowerCEMC::CreateNode : creating TOWER_CALIB_CEMC_RETOWER node " << std::endl;
140 emcalNode->
addNode(emcalTowerNode);
144 std::cout <<
"RetowerCEMC::CreateNode : TOWER_CALIB_CEMC_RETOWER already exists! " << std::endl;