Skip to content
Snippets Groups Projects
Commit c0bd3970 authored by Savannah Rose Shively's avatar Savannah Rose Shively
Browse files

cleanup. energyloss attempt

parent e2a77da1
No related branches found
No related tags found
No related merge requests found
Pipeline #2767621 failed
......@@ -44,7 +44,6 @@ 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_mcEventKey.initialize() );
ATH_CHECK( m_faserSiHitKey.initialize() );
//TH_CHECK( m_faserRdoKey.initialize());
......@@ -56,31 +55,22 @@ StatusCode ClusterReadAlg::initialize() {
// Get the SCT ID helper
ATH_CHECK(detStore()->retrieve(m_idHelper, "FaserSCT_ID"));
// Histograms
m_chi2 = new TH1F("chi2", "Fit #chi^2;#chi^2;Events/bin", 100, 0, 50);
// Histograms
// Tree
m_tree = new TTree("residTree","Cosmics residuals");
m_tree->Branch("id", &m_hash, "id/I");
//m_tree->Branch("residual", &m_residual, "residual/F");
//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("mc_Event", &m_mcEventKey, "mcEvent/F");
m_tree->Branch("waveformExists", &m_waveformExists, "waveform/B");
m_tree->Branch("station", &m_station, "station/I");
//hist test
m_hist = new TH1D("Stations", "Stations hit", 4, 0, 3);
ATH_CHECK(histSvc()->regHist("/HIST/myhist", m_hist));
m_hist2 = new TH1D("EnergyLoss", "Energy",100,0,0.1);
ATH_CHECK(histSvc()->regHist("/HIST/Ehist", m_hist2));
//ATH_CHECK(histSvc()->regHist("/HIST/chi2", m_chi2));
ATH_CHECK(histSvc()->regTree("/HIST/residTree", m_tree));
//rdo
//m_hist = new TH1D("GroupSize", "RDO Group Size", 8, 0, 8); //first string is root object name, second is histogram title
//m_hprof = new TProfile("IncAngleGroup", "Mean Group Size vs Incident Angle", 10, -.1, .1 ,0,5);
//m_incAnglHist = new TH1D("IncAngleHist", "Incident Angle Count", 10, -.1, .1);
//ATH_CHECK(histSvc()->regHist("/HIST/myhist", m_hist));
//ATH_CHECK(histSvc()->regHist("/HIST/myhistprof", m_hprof));
//ATH_CHECK(histSvc()->regHist("/HIST/myhistAngl", m_incAnglHist));
return StatusCode::SUCCESS;
......@@ -131,9 +121,12 @@ StatusCode ClusterReadAlg::execute(const EventContext& ctx) const
int row =hit.getRow();
int module = hit.getModule();
int sensor = hit.getSensor();
std::cout<<"St:"<<station<<"\tPl:"<<plane<<"\tMod:"<<module<<"\n";
m_station[station]++;
int energyloss= hit.energyLoss();
//std::cout<<"St:"<<station<<"\tPl:"<<plane<<"\tMod:"<<module<<"\n";
std::cout<<"\nENERGY LOSS"<<energyloss<<'\n';
++m_station[station];
m_hist->Fill(station);
if (energyloss!=0) m_hist2->Fill(energyloss);
}
//SG::ReadHandle<FaserSCT_RDO_Container> h_sctRDO(m_faserRdoKey);
......@@ -155,8 +148,8 @@ StatusCode ClusterReadAlg::finalize()
ATH_MSG_INFO( m_numberOfTriggeredEvents << " triggered events processed" );
ATH_MSG_INFO( m_numberOfClusterCollection<< " cluster collections processed" );
ATH_MSG_INFO( m_numberOfCluster<< " cluster processed" );
ATH_MSG_INFO( m_numberOfFits << " fits performed" );
ATH_MSG_INFO("station distribution" << m_station[0] << m_station[1]<< m_station[2]<< m_station[3]);
......
......@@ -105,22 +105,16 @@ class ClusterReadAlg : public AthReentrantAlgorithm, AthHistogramming
//SG::ReadHandleKey<TrackerSimDataCollection> m_sctMap {this, "TrackerSimDataCollection", "SCT_SDO_Map"};
mutable TH1* m_chi2;
mutable TTree* m_tree;
mutable int m_hash;
mutable float m_residual;
mutable float m_pull;
mutable float m_refitChi2;
mutable float m_mcEventKey2;
mutable bool m_waveformExists;
mutable int m_station[5]={};
mutable int m_station[5]={}; //replace with coordinate vector later
mutable float m_energyLoss;
mutable std::atomic<int> m_numberOfEvents{0};
mutable std::atomic<int> m_numberOfTriggeredEvents{0};
mutable std::atomic<int> m_numberOfClusterCollection{0};
mutable std::atomic<int> m_numberOfCluster{0};
mutable std::atomic<int> m_numberOfFits{0};
//isvalids
mutable bool m_triggerisValid = false;
......@@ -128,7 +122,7 @@ class ClusterReadAlg : public AthReentrantAlgorithm, AthHistogramming
//RDO
TH1* m_hist; // Example histogram
TH1* m_incAnglHist;
TH1* m_hist2;
TProfile* m_hprof;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment