use Photonis 1 QE

DrcPidFast.cxx

@@ -9,7 +9,7 @@ DrcPidFast::DrcPidFast(){
   fMass[4]=0.938272;
 
   // read Cherenkov track resolution map
-  TFile* file = TFile::Open("ctr_map.root");
+  TFile* file = TFile::Open("ctr_map_p1_0.95.root");
   fTrrMap = new TH2F();
   file->GetObject("htrr", fTrrMap);  
 }

DrcPidFast.h

@@ -10,7 +10,7 @@
 #include "TRandom.h"
 #include <iostream>
 
-// prob - normalized to 1 probability for e,mu,pi,k,p
+// probability - normalized to 1 probability for e,mu,pi,k,p
 // sigma - deviation of the determined Cherenkov angle from expected in terms of Cherenkov track resolution
 struct DrcPidInfo {
     double probability[5];

README.md

@@ -1,4 +1,6 @@
 Fast reconstruction for EIC Barrel DIRC
 
+DrcPidFast::GetInfo returns normalized probability for e, mu, pi, k, p and deviation of the determined Cherenkov angle from the expected in terms of CTR.
+
 Usage example:
-root loadlib.C run.C
+root loadlib.C example.C

example.C

@@ -12,7 +12,7 @@ void example(int pdg=211,double mom=6){
 
   DrcPidInfo info;
   for(int i=0; i<1000; i++){
-    info = pid.GetInfo(pdg,momentum,0);
+    info = pid.GetInfo(pdg,momentum,0.5);
     hPi->Fill(info.sigma[2]);
     hK->Fill(info.sigma[3]);
   }

plot_map.C

@@ -16,8 +16,8 @@ void plot_map(int pdg=211){
     if(m>3){
       delete hK;
       delete hPi;
-      hK = new TH1F("hK","hK",200,-10,10);
-      hPi = new TH1F("hPi","hPi",200,-10,10);
+      hK = new TH1F("hK","hK",200,-15,15);
+      hPi = new TH1F("hPi","hPi",200,-15,15);
     }
     if(m<1){
       delete hK;
@@ -71,11 +71,13 @@ void plot_map(int pdg=211){
 
   gStyle->SetOptStat(0);
   prt_canvasAdd("map_sep",800,500);
-  hsep->GetYaxis()->SetRangeUser(0.6,10);
-  hsep->SetMaximum(10);
   gPad->SetLogz();
   hsep->Draw("colz");
+  hsep->GetYaxis()->SetRangeUser(0.6,10);
+  hsep->SetMaximum(10);
+  hsep->SetMinimum(1);
 
+
   prt_savepath ="data/plot_map";
-  prt_canvasSave();  
+  prt_canvasSave(0);  
 }