macro/old/analysis_xsecQ.C

0001 // SPDX-License-Identifier: LGPL-3.0-or-later
0002 // Copyright (C) 2023 Christopher Dilks
0003
0004 R__LOAD_LIBRARY(EpicAnalysis)
0005 #include "AnalysisDelphes.h"
0006
0007 // cross section in Q2 bins
0008 void analysis_xsecQ(
0009     TString configFile="datarec/arc/crossCheck*.root", /* delphes tree(s) */
0010     TString outfilePrefix="xsecQ" /* output filename prefix*/
0011 ) {
0012
0013   // setup analysis ========================================
0014   AnalysisDelphes *A = new AnalysisDelphes(
0015       configFile,
0016       outfilePrefix
0017       );
0018
0019
0020   // set binning scheme ====================================
0021
0022   // bin 1 (cf Harut slide 11)
0023   A->BinScheme("x")->BuildBin("CenterDelta", 0.3, 0.05 );
0024   A->BinScheme("z")->BuildBin("CenterDelta", 0.7, 0.05 );
0025   A->BinScheme("pt")->BuildBin("CenterDelta", 0.5, 0.05 );
0026
0027   // bin 2
0028   A->BinScheme("x")->BuildBin("CenterDelta", 0.3, 0.05 );
0029   A->BinScheme("z")->BuildBin("CenterDelta", 0.4, 0.05 );
0030   A->BinScheme("pt")->BuildBin("CenterDelta", 0.15, 0.05 );
0031
0032   // diagonalization of (pt,x,z) bins
0033   A->diagonalPtXZ = true;
0034
0035   // Q2 bins
0036   A->BinScheme("q2")->BuildBins(10,1,121,false);
0037
0038   // y minima
0039   A->BinScheme("y")->BuildBin("Min",0.03);
0040   A->BinScheme("y")->BuildBin("Min",0.05);
0041   A->BinScheme("y")->BuildBin("Min",0.10);
0042
0043
0044   // perform the analysis ==================================
0045   A->Execute();
0046 };