ResoStruct getResponseResolutionMC( const std::string& outputdir, TTree* tree, float LYSF[], const std::string& name ) {
std::string fullVarName = getVarName(LYSF);
TH1D* h1 = new TH1D( name.c_str(), "", 500, 0., 500. );
tree->Project( name.c_str(), fullVarName.c_str() );
TF1* f1 = new TF1( Form("f1_%s", name.c_str()), "gaus", 100., 500.);
FitTools::doSingleFit( h1, f1, outputdir, name );
ResoStruct rs = getRespAndReso( f1, 0. );
return rs;
}
ResoStruct getResponseResolutionMC( const std::string& outputdir, TTree* tree, float LYSF[], const std::string& name ) {
std::string fullVarName = getVarName(LYSF);
//fullVarName += " + Ebgo";
//fullVarName = "Ebgo";
//fullVarName = "Eact";
TH1D* h1 = new TH1D( name.c_str(), "", 500, 0., 500. );
//tree->Project( name.c_str(), fullVarName.c_str(), "Escint>0.2");
tree->Project( name.c_str(), fullVarName.c_str(), "Escint>0.2 && Ehodo >0.1" );
TF1* f1 = new TF1( Form("f1_%s", name.c_str()), "gaus", 50., 500.);
FitTools::doSingleFit( h1, f1, outputdir, name,5.,1.5 );
ResoStruct rs = getRespAndReso( f1, 0. );
return rs;
}
int main( int argc, char* argv[]) {
std::string tag = "V00";
if( argc>1 ) {
std::string tag_str(argv[1]);
tag = tag_str;
}
DrawTools::setStyle();
TFile* file_mc = TFile::Open("EEShash_491MeV_10000ev_smear.root");
TFile* file_data = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_259_beam.root", tag.c_str()));
TTree* tree_data = (TTree*)file_data->Get("posTree");
TTree* tree_mc = (TTree*)file_mc->Get("EEShash");
std::string outputdir = "ResolutionStudiesPlots_"+tag;
std::string mkdir_command = "mkdir -p " + outputdir;
system( mkdir_command.c_str() );
TF1* f1_data = FitTools::fitSingleElectronPeak( outputdir, "data", tree_data );
ResoStruct rs_data = getRespAndReso(f1_data, 0.01);
float LYSF_ideal[] = {1., 1., 1., 1., 1., 1., 1., 1., 1., 1.};
float LYSF_real [] = {0.85, 0.94, 0.95, 0.98, 1.00, 1.02, 1.05, 1.05, 1.05, 0.74};
float LYSF_hole [] = {0.85, 0.94, 0.95, 0.98, 1.00, 0.00, 1.05, 1.05, 1.05, 0.74};
ResoStruct rs_ideal = getResponseResolutionMC( outputdir, tree_mc, LYSF_ideal, "ideal" );
ResoStruct rs_real = getResponseResolutionMC( outputdir, tree_mc, LYSF_real, "real" );
ResoStruct rs_hole = getResponseResolutionMC( outputdir, tree_mc, LYSF_hole, "hole" );
std::cout << std::endl;
std::cout << "* DATA: " << std::endl;
std::cout << "reso: " << rs_data.reso << " +/- " << rs_data.reso_error << std::endl;
std::cout << "Sres: " << rs_data.Sres << " +/- " << rs_data.Sres_error << std::endl;
std::cout << "MC reso ideal: "<< rs_ideal.reso << " +/- " << rs_ideal.reso_error << std::endl;
std::cout << "MC reso real: " << rs_real .reso << " +/- " << rs_real .reso_error << std::endl;
std::cout << "MC reso hole: " << rs_hole .reso << " +/- " << rs_hole .reso_error << std::endl;
std::cout << "MC Sres ideal: " << rs_ideal.Sres << " +/- " << rs_ideal.Sres_error << std::endl;
std::cout << "MC Sres real: " << rs_real .Sres << " +/- " << rs_real .Sres_error << std::endl;
std::cout << "MC Sres hole: " << rs_hole .Sres << " +/- " << rs_hole .Sres_error << std::endl;
std::cout << std::endl;
std::cout << "-> Adding photostatistics to MC (25% QE): " << std::endl;
ResoStruct rs_ideal_ps = addPhotoStatistics( rs_ideal, 0.25 );
ResoStruct rs_real_ps = addPhotoStatistics( rs_real, 0.25 );
ResoStruct rs_hole_ps = addPhotoStatistics( rs_hole, 0.25 );
std::cout << "MC reso ideal: "<< rs_ideal_ps.reso << " +/- " << rs_ideal_ps.reso_error << std::endl;
std::cout << "MC reso real: " << rs_real_ps .reso << " +/- " << rs_real_ps .reso_error << std::endl;
std::cout << "MC reso hole: " << rs_hole_ps .reso << " +/- " << rs_hole_ps .reso_error << std::endl;
std::cout << "MC Sres ideal: " << rs_ideal_ps.Sres << " +/- " << rs_ideal_ps.Sres_error << std::endl;
std::cout << "MC Sres real: " << rs_real_ps .Sres << " +/- " << rs_real_ps .Sres_error << std::endl;
std::cout << "MC Sres hole: " << rs_hole_ps .Sres << " +/- " << rs_hole_ps .Sres_error << std::endl;
std::cout << "-> Adding photostatistics to MC (15% QE): " << std::endl;
rs_ideal_ps = addPhotoStatistics( rs_ideal, 0.15 );
rs_real_ps = addPhotoStatistics( rs_real, 0.15 );
rs_hole_ps = addPhotoStatistics( rs_hole, 0.15 );
std::cout << "MC reso ideal: "<< rs_ideal_ps.reso << " +/- " << rs_ideal_ps.reso_error << std::endl;
std::cout << "MC reso real: " << rs_real_ps .reso << " +/- " << rs_real_ps .reso_error << std::endl;
std::cout << "MC reso hole: " << rs_hole_ps .reso << " +/- " << rs_hole_ps .reso_error << std::endl;
std::cout << "MC Sres ideal: " << rs_ideal_ps.Sres << " +/- " << rs_ideal_ps.Sres_error << std::endl;
std::cout << "MC Sres real: " << rs_real_ps .Sres << " +/- " << rs_real_ps .Sres_error << std::endl;
std::cout << "MC Sres hole: " << rs_hole_ps .Sres << " +/- " << rs_hole_ps .Sres_error << std::endl;
// FIRST: DIAGONAL13 SCAN
TFile* file_mc_3x3y = TFile::Open("EEShash_491MeV_10000ev_smear_3x3y.root");
TFile* file_mc_6x6y = TFile::Open("EEShash_491MeV_10000ev_smear_6x6y.root");
TFile* file_mc_9x9y = TFile::Open("EEShash_491MeV_10000ev_smear_9x9y.root");
TFile* file_mc_11p3x11p3y = TFile::Open("EEShash_491MeV_10000ev_smear_11p3x11p3y.root");
TTree* tree_mc_3x3y = (TTree*)file_mc_3x3y->Get("EEShash");
TTree* tree_mc_6x6y = (TTree*)file_mc_6x6y->Get("EEShash");
TTree* tree_mc_9x9y = (TTree*)file_mc_9x9y->Get("EEShash");
TTree* tree_mc_11p3x11p3y = (TTree*)file_mc_11p3x11p3y->Get("EEShash");
//TFile* file_data_3x3y = TFile::Open("AnalysisTrees_V00/Reco_BTF_141_20140430-183508_beam.root");
//TFile* file_data_6x6y = TFile::Open("AnalysisTrees_V00/Reco_BTF_143_20140430-191455_beam.root");
//TFile* file_data_9x9y = TFile::Open("AnalysisTrees_V00/Reco_BTF_167_20140430-210839_beam.root");
//TFile* file_data_11p3x11p3y = TFile::Open("AnalysisTrees_V00/Reco_BTF_219_20140501-092151_beam.root");
//TTree* tree_data_3x3y = (TTree*)file_data_3x3y->Get("recoTree");
//TTree* tree_data_6x6y = (TTree*)file_data_6x6y->Get("recoTree");
//TTree* tree_data_9x9y = (TTree*)file_data_9x9y->Get("recoTree");
//TTree* tree_data_11p3x11p3y = (TTree*)file_data_11p3x11p3y->Get("recoTree");
TFile* file_data_3x3y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_141_beam.root", tag.c_str()));
TFile* file_data_6x6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_143_beam.root", tag.c_str()));
TFile* file_data_9x9y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_167_beam.root", tag.c_str()));
TFile* file_data_11p3x11p3y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_219_beam.root", tag.c_str()));
TTree* tree_data_3x3y = (TTree*)file_data_3x3y->Get("posTree");
TTree* tree_data_6x6y = (TTree*)file_data_6x6y->Get("posTree");
TTree* tree_data_9x9y = (TTree*)file_data_9x9y->Get("posTree");
TTree* tree_data_11p3x11p3y = (TTree*)file_data_11p3x11p3y->Get("posTree");
std::vector<LateralScanStruct> lss_diag;
lss_diag.push_back( LateralScanStruct(0., tree_data, tree_mc) );
lss_diag.push_back( LateralScanStruct(3.*sqrt(2.), tree_data_3x3y, tree_mc_3x3y) );
lss_diag.push_back( LateralScanStruct(6.*sqrt(2.), tree_data_6x6y, tree_mc_6x6y) );
lss_diag.push_back( LateralScanStruct(9.*sqrt(2.), tree_data_9x9y, tree_mc_9x9y) );
lss_diag.push_back( LateralScanStruct(11.3*sqrt(2.), tree_data_11p3x11p3y, tree_mc_11p3x11p3y) );
std::string fullVarName_mc = getVarName(LYSF_hole);
drawLateralScan( outputdir, "diag13", lss_diag, "Diagonal", fullVarName_mc );
// SECOND: DIAGONAL02 SCAN
//TFile* file_data_d02_3x3y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_94_beam.root", tag.c_str()));
//TFile* file_data_d02_6x6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_96_beam.root", tag.c_str()));
//TFile* file_data_d02_9x9y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_98_beam.root", tag.c_str()));
//TFile* file_data_d02_m9xm9y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_171_beam.root", tag.c_str()));
//TFile* file_data_d02_m6xm6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_173_beam.root", tag.c_str()));
//TFile* file_data_d02_m3xm3y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_176_beam.root", tag.c_str()));
TFile* file_data_d02_9x9y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_171_beam.root", tag.c_str()));
TFile* file_data_d02_6x6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_173_beam.root", tag.c_str()));
TFile* file_data_d02_3x3y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_176_beam.root", tag.c_str()));
TTree* tree_data_d02_3x3y = (TTree*)file_data_d02_3x3y->Get("posTree");
TTree* tree_data_d02_6x6y = (TTree*)file_data_d02_6x6y->Get("posTree");
TTree* tree_data_d02_9x9y = (TTree*)file_data_d02_9x9y->Get("posTree");
//TTree* tree_data_d02_m3xm3y = (TTree*)file_data_d02_m3xm3y->Get("posTree");
//TTree* tree_data_d02_m6xm6y = (TTree*)file_data_d02_m6xm6y->Get("posTree");
//TTree* tree_data_d02_m9xm9y = (TTree*)file_data_d02_m9xm9y->Get("posTree");
std::vector<LateralScanStruct> lss_diag02;
lss_diag02.push_back( LateralScanStruct(0., tree_data, tree_mc) );
lss_diag02.push_back( LateralScanStruct(3.*sqrt(2.), tree_data_d02_3x3y, tree_mc_3x3y) );
lss_diag02.push_back( LateralScanStruct(6.*sqrt(2.), tree_data_d02_6x6y, tree_mc_6x6y) );
lss_diag02.push_back( LateralScanStruct(9.*sqrt(2.), tree_data_d02_9x9y, tree_mc_9x9y) );
//lss_diag02.push_back( LateralScanStruct(-3.*sqrt(2.), tree_data_d02_m3xm3y, tree_mc_3x3y) );
//lss_diag02.push_back( LateralScanStruct(-6.*sqrt(2.), tree_data_d02_m6xm6y, tree_mc_6x6y) );
//lss_diag02.push_back( LateralScanStruct(-9.*sqrt(2.), tree_data_d02_m9xm9y, tree_mc_9x9y) );
drawLateralScan( outputdir, "diag02", lss_diag02, "Diagonal", fullVarName_mc );
// THIRD: HORIZONTAL SCAN
TFile* file_mc_2x0y = TFile::Open("EEShash_491MeV_10000ev_smear_2x0y.root");
TFile* file_mc_4x0y = TFile::Open("EEShash_491MeV_10000ev_smear_4x0y.root");
TFile* file_mc_6x0y = TFile::Open("EEShash_491MeV_10000ev_smear_6x0y.root");
TFile* file_mc_8x0y = TFile::Open("EEShash_491MeV_10000ev_smear_8x0y.root");
TFile* file_mc_10x0y = TFile::Open("EEShash_491MeV_10000ev_smear_10x0y.root");
TFile* file_mc_12x0y = TFile::Open("EEShash_491MeV_10000ev_smear_12x0y.root");
TTree* tree_mc_2x0y = (TTree*)file_mc_2x0y ->Get("EEShash");
TTree* tree_mc_4x0y = (TTree*)file_mc_4x0y ->Get("EEShash");
TTree* tree_mc_6x0y = (TTree*)file_mc_6x0y ->Get("EEShash");
TTree* tree_mc_8x0y = (TTree*)file_mc_8x0y ->Get("EEShash");
TTree* tree_mc_10x0y = (TTree*)file_mc_10x0y->Get("EEShash");
TTree* tree_mc_12x0y = (TTree*)file_mc_12x0y->Get("EEShash");
TFile* file_data_12x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_144_20140430-200741_beam.root", tag.c_str()));
TFile* file_data_10x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_145_20140430-201024_beam.root", tag.c_str()));
TFile* file_data_8x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_146_20140430-201243_beam.root", tag.c_str()));
TFile* file_data_6x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_147_20140430-201509_beam.root", tag.c_str()));
TFile* file_data_4x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_148_20140430-201719_beam.root", tag.c_str()));
TFile* file_data_2x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_149_20140430-201913_beam.root", tag.c_str()));
TFile* file_data_m2x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_151_20140430-202639_beam.root", tag.c_str()));
TFile* file_data_m4x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_152_20140430-202852_beam.root", tag.c_str()));
TFile* file_data_m6x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_153_20140430-203054_beam.root", tag.c_str()));
TFile* file_data_m8x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_154_20140430-203255_beam.root", tag.c_str()));
TFile* file_data_m10x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_155_20140430-203551_beam.root", tag.c_str()));
TFile* file_data_m12x0y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_156_20140430-203818_beam.root", tag.c_str()));
TTree* tree_data_2x0y = (TTree*)file_data_2x0y ->Get("posTree");
TTree* tree_data_4x0y = (TTree*)file_data_4x0y ->Get("posTree");
TTree* tree_data_6x0y = (TTree*)file_data_6x0y ->Get("posTree");
TTree* tree_data_8x0y = (TTree*)file_data_8x0y ->Get("posTree");
TTree* tree_data_10x0y = (TTree*)file_data_10x0y ->Get("posTree");
TTree* tree_data_12x0y = (TTree*)file_data_12x0y ->Get("posTree");
TTree* tree_data_m2x0y = (TTree*)file_data_m2x0y ->Get("posTree");
TTree* tree_data_m4x0y = (TTree*)file_data_m4x0y ->Get("posTree");
TTree* tree_data_m6x0y = (TTree*)file_data_m6x0y ->Get("posTree");
TTree* tree_data_m8x0y = (TTree*)file_data_m8x0y ->Get("posTree");
TTree* tree_data_m10x0y = (TTree*)file_data_m10x0y->Get("posTree");
TTree* tree_data_m12x0y = (TTree*)file_data_m12x0y->Get("posTree");
std::vector<LateralScanStruct> lss_horiz;
lss_horiz.push_back( LateralScanStruct(0., tree_data, tree_mc) );
lss_horiz.push_back( LateralScanStruct(2., tree_data_2x0y, tree_mc_2x0y) );
lss_horiz.push_back( LateralScanStruct(4., tree_data_4x0y, tree_mc_4x0y) );
lss_horiz.push_back( LateralScanStruct(6., tree_data_6x0y, tree_mc_6x0y) );
lss_horiz.push_back( LateralScanStruct(8., tree_data_8x0y, tree_mc_8x0y) );
lss_horiz.push_back( LateralScanStruct(10., tree_data_10x0y, tree_mc_10x0y) );
lss_horiz.push_back( LateralScanStruct(12., tree_data_12x0y, tree_mc_12x0y) );
lss_horiz.push_back( LateralScanStruct(-2., tree_data_m2x0y, tree_mc_2x0y) );
lss_horiz.push_back( LateralScanStruct(-4., tree_data_m4x0y, tree_mc_4x0y) );
lss_horiz.push_back( LateralScanStruct(-6., tree_data_m6x0y, tree_mc_6x0y) );
lss_horiz.push_back( LateralScanStruct(-8., tree_data_m8x0y, tree_mc_8x0y) );
lss_horiz.push_back( LateralScanStruct(-10., tree_data_m10x0y, tree_mc_10x0y) );
lss_horiz.push_back( LateralScanStruct(-12., tree_data_m12x0y, tree_mc_12x0y) );
drawLateralScan( outputdir, "horiz", lss_horiz, "Horizontal", fullVarName_mc );
// FOURTH: VERTICAL SCAN
TFile* file_data_0x2y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_160_20140430-204719_beam.root", tag.c_str()));
TFile* file_data_0x4y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_159_20140430-204525_beam.root", tag.c_str()));
TFile* file_data_0x6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_158_20140430-204306_beam.root", tag.c_str()));
TFile* file_data_0x8y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_157_20140430-204053_beam.root", tag.c_str()));
TFile* file_data_0xm2y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_162_20140430-205129_beam.root", tag.c_str()));
TFile* file_data_0xm4y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_163_20140430-205337_beam.root", tag.c_str()));
TFile* file_data_0xm6y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_164_20140430-205542_beam.root", tag.c_str()));
TFile* file_data_0xm8y = TFile::Open(Form("PosAnTrees_%s/PosAn_BTF_165_20140430-205756_beam.root", tag.c_str()));
TTree* tree_data_0x2y = (TTree*)file_data_0x2y ->Get("posTree");
TTree* tree_data_0x4y = (TTree*)file_data_0x4y ->Get("posTree");
TTree* tree_data_0x6y = (TTree*)file_data_0x6y ->Get("posTree");
TTree* tree_data_0x8y = (TTree*)file_data_0x8y ->Get("posTree");
TTree* tree_data_0xm2y = (TTree*)file_data_0xm2y ->Get("posTree");
TTree* tree_data_0xm4y = (TTree*)file_data_0xm4y ->Get("posTree");
TTree* tree_data_0xm6y = (TTree*)file_data_0xm6y ->Get("posTree");
TTree* tree_data_0xm8y = (TTree*)file_data_0xm8y ->Get("posTree");
std::vector<LateralScanStruct> lss_vert;
lss_vert.push_back( LateralScanStruct(0., tree_data, tree_mc) );
lss_vert.push_back( LateralScanStruct(2., tree_data_0x2y, tree_mc_2x0y) );
lss_vert.push_back( LateralScanStruct(4., tree_data_0x4y, tree_mc_4x0y) );
lss_vert.push_back( LateralScanStruct(6., tree_data_0x6y, tree_mc_6x0y) );
lss_vert.push_back( LateralScanStruct(8., tree_data_0x8y, tree_mc_8x0y) );
lss_vert.push_back( LateralScanStruct(-2., tree_data_0xm2y, tree_mc_2x0y) );
lss_vert.push_back( LateralScanStruct(-4., tree_data_0xm4y, tree_mc_4x0y) );
lss_vert.push_back( LateralScanStruct(-6., tree_data_0xm6y, tree_mc_6x0y) );
lss_vert.push_back( LateralScanStruct(-8., tree_data_0xm8y, tree_mc_8x0y) );
drawLateralScan( outputdir, "vert", lss_vert, "Vertical", fullVarName_mc );
return 0;
}
