Просмотр исходного кода

Merge branch 'master' of devel.mephi.ru:PEParfenov/ModelPicoDstFlowCheck

Peter лет назад: 3
Родитель
Сommit
59a3d4234a
2 измененных файлов с 14 добавлено и 2 удалено
  1. 7 1
      get_flow_model.C
  2. 7 1
      get_flow_pico.C

+ 7 - 1
get_flow_model.C

@@ -19,7 +19,7 @@
 
 // R__LOAD_LIBRARY(libPicoDst.so)
 
-const std::vector<float> vResTpc     = {0.295456, 0.451264, 0.471794, 0.432848, 0.3665, 0.297251, 0.260145, 0.269551};
+const std::vector<float> vResTpc     = {0.224964, 0.349906, 0.368124, 0.336364, 0.28126, 0.224188, 0.190761, 0.19178};
 
 void get_flow_model(TString inputFileName, TString outputFileName)
 {
@@ -122,6 +122,12 @@ void get_flow_model(TString inputFileName, TString outputFileName)
       if (pt > pt_max_cut) continue;
       if (abs(eta) > eta_cut) continue;
       if (abs(eta) < eta_gap) continue;
+      
+      // PID-related cut (or to cut out neutral particles)
+      auto particle = (TParticlePDG*) TDatabasePDG::Instance()->GetParticle(pdg[iTr]);
+      if (!particle) continue;
+      float charge = 1./3.*particle->Charge();
+      if (charge == 0) continue;
 
       // TPC Left EP
       if (eta < 0)

+ 7 - 1
get_flow_pico.C

@@ -22,7 +22,7 @@
 
 // R__LOAD_LIBRARY(libPicoDst.so)
 
-const std::vector<float> vResMcTpc   = {0.295164, 0.450308, 0.470348, 0.432475, 0.366723, 0.298105, 0.260313, 0.266459};
+const std::vector<float> vResMcTpc   = {0.226827, 0.347749, 0.369009, 0.336583, 0.281697, 0.224447, 0.18955, 0.186935};
 const std::vector<float> vResRecoTpc = {0.211611, 0.329277, 0.350632, 0.319347, 0.269741, 0.21128, 0.181021, 0.176676};
 
 void get_flow_pico(TString inputFileName, TString outputFileName)
@@ -105,6 +105,12 @@ void get_flow_pico(TString inputFileName, TString outputFileName)
       if (abs(eta) > eta_cut) continue;
       if (abs(eta) < eta_gap) continue;
 
+      // PID-related cut (or to cut out neutral particles)
+      auto particle = (TParticlePDG*) TDatabasePDG::Instance()->GetParticle(mcTrack->GetPdg());
+      if (!particle) continue;
+      float charge = 1./3.*particle->Charge();
+      if (charge == 0) continue;
+
       // Mc-specific track cuts
       if (mcTrack->GetMotherId() != motherId_cut) continue;