diff --git a/src/THcShower.cxx b/src/THcShower.cxx index 1d454d0dde76f35b4905ee425b007711c5c8e1ab..afc22d5e03eac736707e229916bb5be1fea6be93 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 4b7c9786b84ea98b9869f32dac99aabcaf1260b7..b727503d8236713b2c07aa85dc4a87294c32e494 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