obsolete/jd/pfrich.C

0001
0002 // ? export LD_LIBRARY_PATH=/home/ayk/eic/EicSandbox/build/lib:${LD_LIBRARY_PATH}
0003 //   export LD_LIBRARY_PATH=/home/ayk/eic/irt/build/lib:${LD_LIBRARY_PATH}
0004 //
0005 //   root -l './pfrich.C("rich.root")'
0006 //
0007
0008 void pfrich(const char *dfname, const char *cfname = 0)
0009 {
0010   auto fcfg  = new TFile(cfname ? cfname : dfname);
0011   auto geometry = dynamic_cast<CherenkovDetectorCollection*>(fcfg->Get("CherenkovDetectorCollection"));
0012   auto fdata = new TFile(dfname);
0013   TTree *t = dynamic_cast<TTree*>(fdata->Get("t"));
0014   auto event = new CherenkovEvent();
0015   t->SetBranchAddress("e", &event);
0016
0017   int nEvents = t->GetEntries();
0018
0019   auto hxy = new TH2D("hxy", "", 650, -650., 650., 650, -650.0, 650.);
0020
0021   for(unsigned ev=0; ev<nEvents; ev++) {
0022     t->GetEntry(ev);
0023
0024     //printf("@@@ %2lu\n", event->ChargedParticles().size());
0025     for(auto particle: event->ChargedParticles()) {
0026       //if (particle->GetPDG() != 211) continue;
0027
0028       for(auto rhistory: particle->GetRadiatorHistory()) {
0029     auto history  = particle->GetHistory (rhistory);
0030
0031     for(auto photon: history->Photons()) {
0032       //printf("Here! %d\n", photon->WasDetected());
0033       if (!photon->WasDetected() ) continue;
0034
0035       TVector3 phx = photon->GetDetectionPosition();
0036       hxy->Fill(phx.X(), phx.Y());
0037     } //for photon
0038       } //for rhistory
0039     } //for particle
0040   } //for ev
0041
0042   gStyle->SetOptStat(0);
0043   auto cv = new TCanvas("cv", "", 1000, 1000);
0044   hxy->GetXaxis()->SetTitle("Sensor plane X, [mm]");
0045   hxy->GetYaxis()->SetTitle("Sensor plane Y, [mm]");
0046   hxy->Draw("COL");
0047 } // pfrich()