print out the resp centre in the CheckChipResponseFile

Author: ChunzhengLab

Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C

@@ -24,6 +24,7 @@
 
 #define ENABLE_UPGRADES
 #include "ITSMFTSimulation/AlpideSimResponse.h"
+#include "ITS3Simulation/ChipSimResponse.h"
 
 #include "ITS3Base/SegmentationMosaix.h"
 #include "fairlogger/Logger.h"
@@ -34,40 +35,50 @@ using SegmentationMosaix = o2::its3::SegmentationMosaix;
 double um2cm(double um) { return um * 1e-4; }
 double cm2um(double cm) { return cm * 1e+4; }
 
-o2::itsmft::AlpideSimResponse *mAlpSimResp0 = nullptr,
-                              *mAlpSimResp1 = nullptr,
-                              *mAptSimResp1 = nullptr;
+std::unique_ptr<o2::its3::ChipSimResponse> mAlpSimResp0, mAlpSimResp1, mAptSimResp1;
 
-o2::itsmft::AlpideSimResponse* loadResponse(const std::string& fileName, const std::string& respName)
+std::unique_ptr<o2::its3::ChipSimResponse> loadResponse(const std::string& fileName, const std::string& respName)
 {
   TFile* f = TFile::Open(fileName.data());
   if (!f) {
     std::cerr << fileName << " not found" << std::endl;
     return nullptr;
   }
-  auto resp = (o2::itsmft::AlpideSimResponse*)f->Get(respName.data());
-  if (!resp)
+  auto base = f->Get<o2::itsmft::AlpideSimResponse>(respName.c_str());
+  if (!base) {
     std::cerr << respName << " not found in " << fileName << std::endl;
-  return resp;
+    return nullptr;
+  }
+  return std::make_unique<o2::its3::ChipSimResponse>(base);
 }
 
 void LoadRespFunc()
 {
   std::string AptsFile = "$(O2_ROOT)/share/Detectors/Upgrades/ITS3/data/ITS3ChipResponseData/APTSResponseData.root";
   std::string AlpideFile = "$(O2_ROOT)/share/Detectors/ITSMFT/data/AlpideResponseData/AlpideResponseData.root";
 
+  std::cout<<"=====================\n";
+  LOGP(info, "ALPIDE Vbb=0V response");
   mAlpSimResp0 = loadResponse(AlpideFile, "response0"); // Vbb=0V
-  LOG(info) << "ALPIDE Vbb=0V response" << std::endl;
+  mAlpSimResp0->computeCentreFromData();
   mAlpSimResp0->print();
+  LOGP(info, "Response Centre {}", mAlpSimResp0->getRespCentreDep());
+  std::cout<<"=====================\n";
+  LOGP(info, "ALPIDE Vbb=-3V response");
   mAlpSimResp1 = loadResponse(AlpideFile, "response1"); // Vbb=-3V
-  LOG(info) << "ALPIDE Vbb=-3V response" << std::endl;
+  mAlpSimResp1->computeCentreFromData();
   mAlpSimResp1->print();
+  LOGP(info, "Response Centre {}", mAlpSimResp1->getRespCentreDep());
+  std::cout<<"=====================\n";
+  LOGP(info, "APTS response");
   mAptSimResp1 = loadResponse(AptsFile, "response1"); // APTS
-  LOG(info) << "APTS response" << std::endl;
+  mAptSimResp1->computeCentreFromData();
   mAptSimResp1->print();
+  LOGP(info, "Response Centre {}", mAptSimResp1->getRespCentreDep());
+  std::cout<<"=====================\n";
 }
 
-std::vector<float> getCollectionSeediciencies(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getCollectionSeediciencies(o2::its3::ChipSimResponse* resp,
                                               const std::vector<float>& depths)
 {
   std::vector<float> seed;
@@ -81,7 +92,7 @@ std::vector<float> getCollectionSeediciencies(o2::itsmft::AlpideSimResponse* res
   return seed;
 }
 
-std::vector<float> getShareValues(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getShareValues(o2::its3::ChipSimResponse* resp,
                                   const std::vector<float>& depths)
 {
   std::vector<float> share;
@@ -103,7 +114,7 @@ std::vector<float> getShareValues(o2::itsmft::AlpideSimResponse* resp,
   return share;
 }
 
-std::vector<float> getEffValues(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getEffValues(o2::its3::ChipSimResponse* resp,
                                 const std::vector<float>& depths)
 {
   std::vector<float> all;
@@ -135,7 +146,7 @@ void CheckChipResponseFile()
 
   int colors[] = {kOrange + 7, kRed + 1, kAzure + 4};
   struct RespInfo {
-    o2::itsmft::AlpideSimResponse* resp;
+    std::unique_ptr<o2::its3::ChipSimResponse>& resp;
     std::string title;
     int color;
   };
@@ -154,11 +165,11 @@ void CheckChipResponseFile()
   for (auto& r : responses) {
     if (!r.resp)
       continue;
-    auto seed = getCollectionSeediciencies(r.resp, vecDepth);
-    auto shr = getShareValues(r.resp, vecDepth);
-    auto all = getEffValues(r.resp, vecDepth);
+    auto seed = getCollectionSeediciencies(r.resp.get(), vecDepth);
+    auto shr = getShareValues(r.resp.get(), vecDepth);
+    auto all = getEffValues(r.resp.get(), vecDepth);
 
-    TGraph* grSeed = new TGraph(vecDepth.size(), vecDepth.data(), seed.data());
+    auto grSeed = new TGraph(vecDepth.size(), vecDepth.data(), seed.data());
     grSeed->SetTitle(Form("%s seed", r.title.c_str()));
     grSeed->SetLineColor(r.color);
     grSeed->SetLineWidth(2);
@@ -168,7 +179,7 @@ void CheckChipResponseFile()
     grSeed->Draw("SAME LP");
     leg->AddEntry(grSeed, Form("%s seed", r.title.c_str()), "lp");
 
-    TGraph* grShare = new TGraph(vecDepth.size(), vecDepth.data(), shr.data());
+    auto grShare = new TGraph(vecDepth.size(), vecDepth.data(), shr.data());
     grShare->SetLineColor(r.color);
     grShare->SetLineWidth(2);
     grShare->SetMarkerColor(r.color);
@@ -177,7 +188,7 @@ void CheckChipResponseFile()
     grShare->Draw("SAME LP");
     leg->AddEntry(grShare, Form("%s share", r.title.c_str()), "p");
 
-    TGraph* grEff = new TGraph(vecDepth.size(), vecDepth.data(), all.data());
+    auto grEff = new TGraph(vecDepth.size(), vecDepth.data(), all.data());
     grEff->SetLineColor(r.color);
     grEff->SetLineWidth(2);
     grEff->SetMarkerColor(r.color);