delphes/scripts/dconfig.C

0001
0002 //
0003 // root -l dconfig.C
0004 //
0005
0006 void dconfig( void )
0007 {
0008   auto delphes = new DelphesConfigRICH("pfRICH");
0009
0010   // Define particle mass hypotheses in ascending mass order; yes, there is no
0011   // reason to overcomplicate things;
0012   delphes->AddMassHypothesis("pi+");
0013   delphes->AddMassHypothesis("K+");
0014   delphes->AddMassHypothesis("proton");
0015
0016   // Define radiator refractive index;
0017   delphes->SetRefractiveIndex(1.019);
0018   // Imitate tracker resolution;
0019   delphes->SetAdditionalSmearing(0.5);
0020
0021   // Define eta ranges in ascending order;
0022   {
0023     auto eta = delphes->AddEtaRange(-3.5, -3.0);
0024
0025     // Arguments [0..1]: momentum ranges for this eta in [GeV] in ascending order without gaps,
0026     // Arguments [2..N]: sigma values in [mrad] matching the mass hypotheses above, in the same
0027     // ascending order used to call delphes->AddMassHypothesis();
0028     eta->AddMomentumRange( 8., 10., 1.5, 1.5, 1.0);
0029     eta->AddMomentumRange(10., 11., 1.0, 1.0, 1.5);
0030   }
0031   {
0032     auto eta = delphes->AddEtaRange(-3.0, -1.7);
0033
0034     // Well, average bin momentum is used to calculate Cherenkov theta, so the bins
0035     // should not be too wide in the areas where dtheta/dp is high for at least one
0036     // of the hypotheses;
0037     eta->AddMomentumRange( 8., 10., 1.1, 1.2, 1.3);
0038     eta->AddMomentumRange(10., 11., 1.0, 1.0, 1.5);
0039   }
0040
0041   {
0042     auto eta = delphes->AddEtaRange(-1.7, -1.5);
0043
0044     eta->AddMomentumRange( 8.,  9., 1.1, 1.2, 1.3);
0045     eta->AddMomentumRange( 9., 10., 1.1, 1.2, 1.3);
0046     eta->AddMomentumRange(10., 11., 1.0, 1.0, 1.5);
0047   }
0048
0049   delphes->WriteTcl();
0050   exit(0);
0051 } // dconfig()