delphes/scripts/delphes-dirc.C

0001
0002 //
0003 // Units are [mrad], [mm], [T];
0004 //
0005 // root -l delphes-dirc.C
0006 //
0007
0008 //#define _E_PI_SEPARATION_MODE_
0009
0010 void delphes_dirc( void )
0011 {
0012   //printf("%f\n", (1.0 - erf(1.75/(sqrt(2.)*1.0)))/2);
0013
0014   auto dirc = new DelphesConfigDIRC("barrelDIRC");
0015   dirc->UsePtMode();
0016
0017   // Define particle mass hypotheses in ascending mass order; yes, there is no
0018   // reason to overcomplicate things;
0019 #ifdef _E_PI_SEPARATION_MODE_
0020   dirc->AddMassHypothesis(-11);
0021 #endif
0022   dirc->AddMassHypothesis("pi+");
0023 #ifndef _E_PI_SEPARATION_MODE_
0024   dirc->AddMassHypothesis("K+");
0025   dirc->AddMassHypothesis("proton");
0026 #endif
0027
0028   // "d(theta) ~ a/p + b" parameterization as taken from Wenqing's 11/23/21 slides,
0029   // roughly averaged over eta;
0030   dirc->SetTrackerAngularResolution(0.9, 0.1);
0031
0032   // eta and momentum range and binning;
0033   dirc->SetEtaRange     (-1.64, 1.25,  10);
0034   //dirc->SetEtaRange     ( 1.0,  1.01,   1);
0035   // Do not mind to use Pt rather than 1/Pt bins; [GeV/c];
0036 #ifdef _E_PI_SEPARATION_MODE_
0037   dirc->SetMomentumRange( 0.44, 3.00,  10);
0038   //dirc->SetMomentumRange( 1.20, 1.21,   1);
0039 #else
0040   dirc->SetMomentumRange( 0.44,10.00,  10);
0041 #endif
0042
0043   // Installation radius in [mm]; constant magnetic field in [T];
0044   dirc->SetInstallationRadius  (1060.0);
0045   dirc->SetMagneticField       ( 3.000);
0046
0047   dirc->SetParameterizationMap("./ctr_map_p1_0.95.root");
0048
0049   dirc->DoSigmaCalculations();
0050   // This is again some generic stuff;
0051   dirc->WriteTcl(false);
0052   exit(0);
0053 } // delphes_dirc()