obsolete/ak/gas-eval.C

0001
0002
0003 void gas_eval(const char *dfname, const char *cfname = 0)
0004 {
0005   auto *reco = new ReconstructionFactory(dfname, cfname, "pfRICH");
0006
0007   // Mark all pads hit by "calibration" (above the QE curve) photons in these
0008   // radiators as "useless" for the purposes of "GasVolume" generated photons;
0009   reco->AddBlackoutRadiator("QuartzWindow");
0010   //reco->AddBlackoutRadiator("Acrylic");
0011   //reco->AddBlackoutRadiator("BelleIIAerogel1");
0012
0013   // Carelessly remove (0x1 << n)x(0x1 << n) square area "around" these hits;
0014   reco->SetBlackoutBlowupValue(3);
0015
0016   // Whatever we want it;
0017   reco->SetSinglePhotonTimingResolution(0.050);
0018   // Bypass imaging;
0019   reco->SetExperimentalMode();
0020
0021   {
0022     CherenkovEvent *event;
0023
0024     // Event loop;
0025     while((event = reco->GetNextEvent())) {
0026       printf("%3d\n", (int)reco->GetDigitizedHits().size());
0027
0028       // Loop over digitized hits;
0029       for(unsigned iq=0; iq<reco->GetDigitizedHits().size(); iq++) {
0030     auto &hit = reco->GetDigitizedHits()[iq];
0031
0032     // From this point on see DigitizedHit.h class methods;
0033       } //for iq
0034     } //while
0035   }
0036 } // gas_eval()