Skip to content

Commit

Permalink
Merge pull request #280 from VERITAS-Observatory/v491-dev8
Browse files Browse the repository at this point in the history 
v491 dev8 (minor plottting updates)
  • Loading branch information
@GernotMaier
GernotMaier authored Nov 19, 2024
2 parents 4916238 + d50f984 commit 2ac563c
Show file tree
Hide file tree
Showing 6 changed files with 56 additions and 13 deletions.
5 changes: 3 additions & 2 deletions inc/VPlotAnasumHistograms.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -101,7 +101,8 @@ class VPlotAnasumHistograms : public VAnalysisUtilities, public VPlotUtilities,
void plot_mscPlots( int irebin = 2, double xmin = -2., double xmax = 4., string mscwfile = "" );
void plot_qualityHistograms( double iSourceStrength = 1., bool bUpper = true, int iMethod = 0 );
TCanvas* plot_theta2( double t2min = 0., double t2max = 0.3, int irbin = 5,
double setYMax = -9999., double setYMin = -9999. );
double setYMax = -9999., double setYMin = -9999.,
bool containment_lines = false );
void plot_theta2Correction();
void plot_UncorrelatedSkyPlots();
void plot_CorrelatedSkyPlots();
Expand Down Expand Up @@ -144,7 +145,7 @@ class VPlotAnasumHistograms : public VAnalysisUtilities, public VPlotUtilities,
}
bool setRunNumber( int iRun ); // select run for plotting

ClassDef( VPlotAnasumHistograms, 16 );
ClassDef( VPlotAnasumHistograms, 17 );
};

#endif
1 change: 1 addition & 0 deletions inc/VPlotInstrumentResponseFunction.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -19,6 +19,7 @@
#include "VHistogramUtilities.h"
#include "VInstrumentResponseFunctionReader.h"
#include "VPlotUtilities.h"
#include "VUtilities.h"

using namespace std;

Expand Down
15 changes: 15 additions & 0 deletions inc/VUtilities.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -48,5 +48,20 @@ namespace VUtilities

return (*p == 0 ) ;
}

// friendlier colors
inline int color_id(unsigned int index)
{
// Array containing the specified ROOT color ID
int colors[] = {12, 633, 9, 418, 801, 881, 900, 32, 423, 393, 798, 616, 409};
unsigned int size = sizeof(colors) / sizeof(colors[0]);

if (index >= size)
{
return index;
}
return colors[index];
}

}
#endif
38 changes: 30 additions & 8 deletions src/VPlotAnasumHistograms.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -645,7 +645,7 @@ void VPlotAnasumHistograms::plot_CorrelatedSkyPlots()
*
*
*/
TCanvas* VPlotAnasumHistograms::plot_theta2( double t2min, double t2max, int irbin, double setYMax, double setYMin )
TCanvas* VPlotAnasumHistograms::plot_theta2( double t2min, double t2max, int irbin, double setYMax, double setYMin, bool containment_lines )
{
// int iPlotPSF = 0;

Expand Down Expand Up @@ -690,33 +690,43 @@ TCanvas* VPlotAnasumHistograms::plot_theta2( double t2min, double t2max, int irb
htheta2_off->Draw( "hist e" );
htheta2_on->Draw( "hist e same" );

// get 68% containment radius (up to theta2 = 0.05deg2)
// get 68% containment radius (up to theta2 = 0.2deg2)
double nt2 = 0.;
for( int i = 1; i < htheta2_diff->GetXaxis()->FindBin( 0.05 ); i++ )
double max_integration_radius = 0.1;
for( int i = 1; i < htheta2_diff->GetXaxis()->FindBin( max_integration_radius ); i++ )
{
nt2 += htheta2_diff->GetBinContent( i );
}
double nt2_68 = 0.;
double t_68 = -1.;
double t_95 = 0.;
if( nt2 > 0. )
{
for( int i = 1; i < htheta2_diff->GetXaxis()->FindBin( 0.05 ); i++ )
for( int i = 1; i < htheta2_diff->GetXaxis()->FindBin( max_integration_radius ); i++ )
{
nt2_68 += htheta2_diff->GetBinContent( i );
if( nt2_68 / nt2 > 0.68 )
if( nt2_68 / nt2 > 0.68 && t_68 < 0. )
{
cout << "Theta2 containment radius (68%, binning dependent): " << htheta2_diff->GetXaxis()->GetBinLowEdge( i ) << " deg2" << endl;
t_68 = htheta2_diff->GetXaxis()->GetBinLowEdge( i );
}
if( nt2_68 / nt2 > 0.95 )
{
t_95 = htheta2_diff->GetXaxis()->GetBinLowEdge( i );
break;
}
}
}
cout << "Theta2 containment radius (68%, binning dependent): " << t_68 << " deg2" << endl;
cout << "Theta2 containment radius (95%, binning dependent): " << t_95 << " deg2" << endl;

c_t2_diff->cd();
c_t2_diff->SetLeftMargin( 0.14 );
htheta2_diff->SetFillColor( 418 );
htheta2_diff->SetXTitle( "#Theta^{2} [deg^{2}]" );
htheta2_diff->SetTitle( "" );
sprintf( hname, "Number of events / %.3f deg^{2}", htheta2_diff->GetXaxis()->GetBinWidth( 2 ) );
htheta2_diff->SetYTitle( hname );
htheta2_diff->GetYaxis()->SetTitleOffset( 1.6 );
htheta2_diff->GetYaxis()->SetTitleOffset( 3.0 );
setHistogramPlottingStyle( htheta2_diff, 418, 2, 1, 1, irbin, 1001 );
htheta2_diff->SetAxisRange( t2min, t2max );
if( setYMax > -999. )
Expand All @@ -730,6 +740,18 @@ TCanvas* VPlotAnasumHistograms::plot_theta2( double t2min, double t2max, int irb

htheta2_diff->Draw( "hist e" );

if( containment_lines )
{
TLine *iL_68 = new TLine( t_68, 0., t_68, htheta2_diff->GetMaximum()*0.5);
iL_68->SetLineStyle(2);
iL_68->SetLineColor(2);
iL_68->Draw();
TLine *iL_95 = new TLine( t_95, 0., t_95, htheta2_diff->GetMaximum()*0.5);
iL_95->SetLineStyle(2);
iL_95->SetLineColor(2);
iL_95->Draw();
}


if( fDebug )
{
Expand Down Expand Up @@ -1132,7 +1154,7 @@ TCanvas* VPlotAnasumHistograms::plot_radec( int sPlot, double rmax, double zmin,
{
sprintf( hname, "c_skysig_%d_%d", fRunNumber, sPlot );
sprintf( htitle, "sky map, run %d", fRunNumber );
c_skysig = new TCanvas( hname, htitle, 10, 10, 400, 400 );
c_skysig = new TCanvas( hname, htitle, 10, 10, 800, 800 );
c_skysig->Draw();
c_skysig->SetRightMargin( 0.14 );
c_skysig->SetLeftMargin( 0.11 );
Expand Down
8 changes: 6 additions & 2 deletions src/VPlotInstrumentResponseFunction.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -102,7 +102,7 @@ bool VPlotInstrumentResponseFunction::addInstrumentResponseData( string iFile, d
iTempIRFReader->setDebug( fDebug );
if( iColor < 0 )
{
iColor = fData.size() + 1;
iColor = VUtilities::color_id(fData.size() + 1);
}
if( iLineStyle < 0 )
{
Expand Down Expand Up @@ -1574,7 +1574,9 @@ TCanvas* VPlotInstrumentResponseFunction::plotPSF( vector< double > i_Energy_TeV
sprintf( hname, "Theta_ID_%d", iCumulative );
}
TCanvas* c = new TCanvas( hname, hname, 10, 10, 600, 600 );
c->Divide( TMath::Nint( sqrt( i_Energy_TeV_lin.size() ) ), TMath::Nint( sqrt( i_Energy_TeV_lin.size() ) ) );
cout << "PSF plotting: " << i_Energy_TeV_lin.size() << " energies (";
cout << TMath::Ceil( sqrt( i_Energy_TeV_lin.size() ) ) << ")" << endl;
c->Divide( TMath::Ceil( sqrt( i_Energy_TeV_lin.size() ) ), TMath::Ceil( sqrt( i_Energy_TeV_lin.size() ) ) );
for( unsigned int j = 0; j < i_Energy_TeV_lin.size(); j++ )
{
c->cd( j + 1 );
Expand Down Expand Up @@ -1632,10 +1634,12 @@ TCanvas* VPlotInstrumentResponseFunction::plotPSF( vector< double > i_Energy_TeV
TH1D* hCumu = get_Cumulative_Histogram( h, true, true );
if( iCumulative )
{
hCumu->SetLineColor( VUtilities::color_id(i+1) );
hCumu->Draw( "same" );
}
else
{
h->SetLineColor( VUtilities::color_id(i+1) );
// rebin
if( iPlotTheta2 )
{
Expand Down
2 changes: 1 addition & 1 deletion src/printCrabSensitivity.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@ void print_sensitivity( string anasum_file, double alpha = 1. / 6. )

t->GetEntry( t->GetEntries() - 1 );

cout << "Rates (on/Off): " << Rate << "\t" << RateOff << endl;
cout << "Rates (on/off): " << Rate << "\t" << RateOff << endl;

vector< double > fSourceStrength;
fSourceStrength.push_back( 1. );
Expand Down

0 comments on commit 2ac563c

Please sign in to comment.