From ca8268ff6a7f979f98e31a6670eea8d9473378d9 Mon Sep 17 00:00:00 2001
From: Vardan Tadevosyan <tadevosn@jlab.org>
Date: Sat, 8 Mar 2014 15:06:51 +0400
Subject: [PATCH] Replace P versus Edep/P histogram by #DeltaP versus Edep/P
---
hcal_calib/THcShowerCalib.h | 13 ++++++++++---
hcal_calib/hcal_calib.cpp | 3 ++-
2 files changed, 12 insertions(+), 4 deletions(-)
diff --git a/hcal_calib/THcShowerCalib.h b/hcal_calib/THcShowerCalib.h
index 29008a9..3a4e027 100644
--- a/hcal_calib/THcShowerCalib.h
+++ b/hcal_calib/THcShowerCalib.h
@@ -39,7 +39,8 @@ class THcShowerCalib {
TH1F* hEunc;
TH1F* hEuncSel;
TH1F* hEcal;
- TH2F* hPvsEcal;
+ // TH2F* hPvsEcal;
+ TH2F* hDPvsEcal;
private:
Int_t fRunNumber;
@@ -226,7 +227,9 @@ void THcShowerCalib::Init() {
hEunc = new TH1F("hEunc", "Edep/P uncalibrated", 500, 0., 5.);
hEcal = new TH1F("hEcal", "Edep/P calibrated", 150, 0., 1.5);
- hPvsEcal = new TH2F("hPvsEcal", "P versus Edep/P ",150,0.,1.5, 100,0.7,0.9);
+ // hPvsEcal = new TH2F("hPvsEcal", "P versus Edep/P ",150,0.,1.5, 100,0.7,0.9);
+ hDPvsEcal = new TH2F("hDPvsEcal", "#DeltaP versus Edep/P ",
+ 150,0.,1.5, 250,-12.5,12.5);
// Initialize qumulative quantities.
@@ -728,7 +731,11 @@ void THcShowerCalib::FillHEcal() {
Double_t Enorm = trk.Enorm();
hEcal->Fill(Enorm);
- hPvsEcal->Fill(Enorm,P/1000.,1.);
+
+ Double_t delta;
+ fTree->SetBranchAddress("H.cal.trdelta",&delta);
+ // hPvsEcal->Fill(Enorm,P/1000.,1.);
+ hDPvsEcal->Fill(Enorm,delta,1.);
output << Enorm*P/1000. << " " << P/1000. << endl;
diff --git a/hcal_calib/hcal_calib.cpp b/hcal_calib/hcal_calib.cpp
index 4798bbb..b2cdcf8 100644
--- a/hcal_calib/hcal_calib.cpp
+++ b/hcal_calib/hcal_calib.cpp
@@ -43,6 +43,7 @@ void hcal_calib(Int_t RunNumber) {
// Momentum versus the calibrated energy deposition.
Canvas->cd(4);
- theShowerCalib.hPvsEcal->Draw();
+ // theShowerCalib.hPvsEcal->Draw();
+ theShowerCalib.hDPvsEcal->Draw();
}
--
GitLab