hpstr
The Heavy Photon Search Toolkit for Reconstruction (hpstr) provides an interface to physics data from the HPS experiment saved in the LCIO format and converts it into an ROOT based format.
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
RecoHitAnaHistos.cxx
Go to the documentation of this file.
1#include "RecoHitAnaHistos.h"
2#include <math.h>
3
7
8void RecoHitAnaHistos::FillTrackerHits(std::vector<TrackerHit*> *trkrHits, float weight ) {
9 int nHits = trkrHits->size();
10 Fill1DHisto("numTrkrHits_h", (float)nHits, weight);
11 for (int i=0; i < nHits; i++)
12 {
13 TrackerHit *hit = trkrHits->at(i);
14 Fill1DHisto("trkrHitEdep_h", hit->getCharge()*1000000.0, weight); // Scaled to MeV
15 }
16}
17
18void RecoHitAnaHistos::FillTracks(std::vector<Track*> *tracks, float weight ) {
19 int nTracks = tracks->size();
20 Fill1DHisto("numTracks_h", (float)nTracks, weight);
21}
22
23void RecoHitAnaHistos::FillEcalHits(std::vector<CalHit*> *ecalHits, float weight ) {
24 int nHits = ecalHits->size();
25 Fill1DHisto("numEcalHits_h", (float)nHits, weight);
26 for (int i=0; i < nHits; i++)
27 {
28 CalHit *hit = ecalHits->at(i);
29 Fill1DHisto("ecalHitEnergy_h", hit->getEnergy()*1000.0, weight); // Scaled to MeV
30 }
31}
32
33void RecoHitAnaHistos::FillEcalClusters(std::vector<CalCluster*> *ecalClusters, float weight ) {
34 int nClusters = ecalClusters->size();
35 Fill1DHisto("numEcalClusters_h", (float)nClusters, weight);
36 for (int i=0; i < nClusters; i++)
37 {
38 CalCluster *cluster = ecalClusters->at(i);
39 Fill1DHisto("ecalClusterEnergy_h", cluster->getEnergy()*1000.0, weight); // Scaled to MeV
40 }
41}
double getEnergy() const
Definition CalCluster.h:73
double getEnergy() const
Definition CalHit.h:36
void Fill1DHisto(const std::string &histoName, float value, float weight=1.)
description
void FillEcalClusters(std::vector< CalCluster * > *ecalClusters, float weight=1.)
description
virtual void Define1DHistos()
description
void FillEcalHits(std::vector< CalHit * > *ecalHits, float weight=1.)
description
void FillTrackerHits(std::vector< TrackerHit * > *trkrHits, float weight=1.)
description
void FillTracks(std::vector< Track * > *tracks, float weight=1.)
description
double getCharge() const
Definition TrackerHit.h:93