From c618cc1b50e02de004dce9f79912d56609eb9449 Mon Sep 17 00:00:00 2001
From: adrien-matta <a.matta@surrey.ac.uk>
Date: Wed, 30 Apr 2014 11:33:45 +0100
Subject: [PATCH] * Fixing Geometrical efficiency macro
---
NPAnalysis/macros/GeometricalEfficiency.C | 21 ++++++++++-----------
1 file changed, 10 insertions(+), 11 deletions(-)
diff --git a/NPAnalysis/macros/GeometricalEfficiency.C b/NPAnalysis/macros/GeometricalEfficiency.C
index abcf77661..06ed0c122 100644
--- a/NPAnalysis/macros/GeometricalEfficiency.C
+++ b/NPAnalysis/macros/GeometricalEfficiency.C
@@ -65,20 +65,18 @@ void GeometricalEfficiency(const char * fname = "myResult"){
tree->SetBranchStatus("InteractionCoordinates", true);
// Prepare histograms
- TH1F *hDetecTheta = new TH1F("hDetecTheta", "DetecTheta", 90, 0, 180);
- TH1F *hDetecThetaCM = new TH1F("hDetecThetaCM", "hDetecThetaCM", 90, 0, 180);
- TH1F *hEmittTheta = new TH1F("hEmittTheta", "EmittTheta", 90, 0, 180);
- TH1F *hEmittThetaCM = new TH1F("hEmittThetaCM", "hEmittThetaCM", 90, 0, 180);
+ TH1F *hDetecTheta = new TH1F("hDetecTheta", "DetecTheta", 180,0,180);
+ TH1F *hDetecThetaCM = new TH1F("hDetecThetaCM", "hDetecThetaCM", 180,0,180);
+ TH1F *hEmittTheta = new TH1F("hEmittTheta", "EmittTheta", 180,0,180);
+ TH1F *hEmittThetaCM = new TH1F("hEmittThetaCM", "hEmittThetaCM", 180,0,180);
// Read the TTree
- Int_t nentries = tree->GetEntries();
- // cout << "TTree contains " << nentries << " events" << endl;
- for (Int_t i = 0; i < nentries; i++) {
- //if (i%1000 == 0) cout << "Entry " << i << endl;
+ int nentries = tree->GetEntries();
+ for (int i = 0; i < nentries; i++) {
tree->GetEntry(i);
// Fill histos
hEmittTheta->Fill(initCond->GetThetaLab_WorldFrame(0));
- hEmittTheta->Fill(initCond->GetThetaCM(0));
+ hEmittThetaCM->Fill(initCond->GetThetaCM(0));
if (interCoord->GetDetectedMultiplicity() > 0){
hDetecTheta->Fill(interCoord->GetDetectedAngleTheta(0));
@@ -87,13 +85,14 @@ void GeometricalEfficiency(const char * fname = "myResult"){
}
TCanvas* c4 = new TCanvas("c4", "CM Frame");
+
TH1F* SolidACM = new TH1F(*hDetecThetaCM);
SolidACM->Sumw2();
- TF1* CCM = new TF1("CCM",Form("%i /(4*%f)",nentries,M_PI),0,180);
TF1* C = new TF1("C",Form("%i /(4*%f)",nentries,M_PI),0,180);
SolidACM->Divide(C,1);
+ SolidACM->Divide(hEmittThetaCM);
SolidACM->Draw();
- TF1* f = new TF1("f",Form("2 * %f * sin(x*%f/180.) *2*%f/180.",M_PI,M_PI,M_PI),0,180);
+ TF1* f = new TF1("f",Form("2 * %f * sin(x*%f/180.) *1*%f/180.",M_PI,M_PI,M_PI),0,180);
f->Draw("SAME");
f->Draw("SAME");
SolidACM->GetXaxis()->SetTitle("#theta_{CM} (deg)");
--
GitLab