From 7f90f66e19c73b1c2fb0323060b0d295c6c9bb3d Mon Sep 17 00:00:00 2001
From: Vardan Tadevosyan <tadevosn@jlab.org>
Date: Thu, 29 Jun 2017 16:05:05 -0400
Subject: [PATCH] Add SHMS calorimeter's total energy deposition divided by
momentum of the best track.
---
src/THcShower.cxx | 3 +++
src/THcShower.h | 7 ++++---
2 files changed, 7 insertions(+), 3 deletions(-)
diff --git a/src/THcShower.cxx b/src/THcShower.cxx
index 1d454d0..afc22d5 100644
--- a/src/THcShower.cxx
+++ b/src/THcShower.cxx
@@ -518,6 +518,7 @@ Int_t THcShower::DefineVariables( EMode mode )
{ "etracknorm", "Total energy divided by track momentum", "fEtrackNorm" },
{ "eprtrack", "Track Preshower energy", "fEPRtrack" },
{ "eprtracknorm", "Preshower energy divided by track momentum", "fEPRtrackNorm" },
+ { "etottracknorm", "Total energy divided by track momentum", "fETotTrackNorm" },
{ "ntracks", "Number of shower tracks", "fNtracks" },
{ 0 }
};
@@ -601,6 +602,7 @@ void THcShower::Clear(Option_t* opt)
fEtrackNorm = 0.;
fEPRtrack = 0.;
fEPRtrackNorm = 0.;
+ fETotTrackNorm = 0.;
fSizeClustArray = 0;
fNblockHighEnergy = 0.;
@@ -1142,6 +1144,7 @@ Int_t THcShower::FineProcess( TClonesArray& tracks )
fEtrackNorm=fEtrack/theTrack->GetP();
fEPRtrack=GetShEnergy(theTrack,1);
fEPRtrackNorm=fEPRtrack/theTrack->GetP();
+ fETotTrackNorm=fEtot/theTrack->GetP();
Xtr = -100.;
Ytr = -100.;
fNclustTrack = MatchCluster(theTrack, Xtr, Ytr);
diff --git a/src/THcShower.h b/src/THcShower.h
index 4b7c978..b727503 100644
--- a/src/THcShower.h
+++ b/src/THcShower.h
@@ -163,9 +163,10 @@ protected:
Double_t fEtot; // Total energy
Double_t fEtotNorm; // Total energy divided by spec central momentum
Double_t fEtrack; // Cluster energy associated to the best track
- Double_t fEtrackNorm; // Normalized to momentum cluster energy associated to the best track
- Double_t fEPRtrack; // Cluster's Preshower energy of the best track
- Double_t fEPRtrackNorm; // Normalized to momentum cluster's Preshower energy of the best track
+ Double_t fEtrackNorm; // Cluster energy divided by momentum for the best track
+ Double_t fEPRtrack; // Preshower part of cluster energy of the best track
+ Double_t fEPRtrackNorm; // Preshower part of cluster energy divided by momentum for the best track
+ Double_t fETotTrackNorm; // Total energy divided by momentum of the best track
THcShowerClusterList* fClusterList; // List of hit clusters
--
GitLab