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; 12/05/2022: still use the same, for the time being;
0030   dirc->SetTrackerAngularResolution(0.9, 0.1);
0031
0032   // eta and momentum range and binning;
0033   dirc->SetEtaRange     (-1.55, 1.79,  10);
0034   // Do not mind to use Pt rather than 1/Pt bins; [GeV/c];
0035 #ifdef _E_PI_SEPARATION_MODE_
0036   dirc->SetMomentumRange( 0.44, 3.00,  10);
0037   //dirc->SetMomentumRange( 1.20, 1.21,   1);
0038 #else
0039   dirc->SetMomentumRange( 0.44,10.00,  10);
0040 #endif
0041
0042   // Installation radius in [mm]; constant magnetic field in [T];
0043   dirc->SetInstallationRadius  (729.6);
0044   dirc->SetMagneticField       ( 1.700);
0045
0046   dirc->SetParameterizationMap("./ctr_map_p1_0.95.root");
0047
0048   dirc->DoSigmaCalculations();
0049   // This is again some generic stuff;
0050   dirc->WriteTcl(false);
0051   {
0052     auto fout = new TFile("barrelDIRC.root", "RECREATE");
0053     dirc->Write();
0054     fout->Close();
0055   }
0056   exit(0);
0057 } // delphes_dirc()