diff --git a/Modules/HMPID/include/HMPID/HmpidTaskMatches.h b/Modules/HMPID/include/HMPID/HmpidTaskMatches.h index f50dbbcdda..22f643b7d3 100644 --- a/Modules/HMPID/include/HMPID/HmpidTaskMatches.h +++ b/Modules/HMPID/include/HMPID/HmpidTaskMatches.h @@ -65,6 +65,7 @@ class HmpidTaskMatches : public TaskInterface TH1F* mMatchInfoResidualsYTrackMIP[7]; TH1F* mMatchInfoChargeClusterMIP[7]; TH1F* mMatchInfoChargeClusterPhotons[7]; + TH2F* mMatchInfoClusterMIPMap[7]; TH2F* mMatchInfoThetaCherenkovVsMom[7]; std::vector mPublishedObjects; diff --git a/Modules/HMPID/src/HmpidTaskMatches.cxx b/Modules/HMPID/src/HmpidTaskMatches.cxx index 8f58dd7b4b..2c016536eb 100644 --- a/Modules/HMPID/src/HmpidTaskMatches.cxx +++ b/Modules/HMPID/src/HmpidTaskMatches.cxx @@ -50,6 +50,7 @@ HmpidTaskMatches::~HmpidTaskMatches() delete mMatchInfoResidualsYTrackMIP[iCh]; delete mMatchInfoChargeClusterMIP[iCh]; delete mMatchInfoChargeClusterPhotons[iCh]; + delete mMatchInfoClusterMIPMap[iCh]; delete mMatchInfoThetaCherenkovVsMom[iCh]; } } @@ -72,6 +73,9 @@ void HmpidTaskMatches::initialize(o2::framework::InitContext& /*ctx*/) getObjectsManager()->startPublishing(mMatchInfoResidualsYTrackMIP[iCh]); getObjectsManager()->startPublishing(mMatchInfoChargeClusterMIP[iCh]); getObjectsManager()->startPublishing(mMatchInfoChargeClusterPhotons[iCh]); + getObjectsManager()->startPublishing(mMatchInfoClusterMIPMap[iCh]); + getObjectsManager()->setDefaultDrawOptions(mMatchInfoClusterMIPMap[iCh], "colz"); + getObjectsManager()->setDisplayHint(mMatchInfoClusterMIPMap[iCh], "colz"); getObjectsManager()->startPublishing(mMatchInfoThetaCherenkovVsMom[iCh]); } @@ -123,8 +127,10 @@ void HmpidTaskMatches::monitorData(o2::framework::ProcessingContext& ctx) mMatchInfoResidualsXTrackMIP[ch]->Fill(xTrk - xMip); mMatchInfoResidualsYTrackMIP[ch]->Fill(yTrk - yMip); - if (matchHMP.getMipClusSize() < 20 && matchHMP.getMipClusSize() > 1) + if (matchHMP.getMipClusSize() < 20 && matchHMP.getMipClusSize() > 1) { mMatchInfoChargeClusterMIP[ch]->Fill(q); + } + mMatchInfoClusterMIPMap[ch]->Fill(xMip, yMip); mMatchInfoThetaCherenkovVsMom[ch]->Fill(TMath::Abs(matchHMP.getHmpMom()), matchHMP.getHMPsignal()); } // match info loop } @@ -148,6 +154,8 @@ void HmpidTaskMatches::BookHistograms() mMatchInfoResidualsYTrackMIP[iCh] = new TH1F(Form("Y Residuals chamber %i", iCh), Form("Y Residuals chamber %i; Y residuals (cm); Counts/1 cm", iCh), 100, -50., 50.); mMatchInfoChargeClusterMIP[iCh] = new TH1F(Form("MIP Cluster Charge chamber %i", iCh), Form("MIP Cluster Charge chamber %i; ADC charge; Counts/1 ADC", iCh), 2000, 200., 2200.); mMatchInfoChargeClusterPhotons[iCh] = new TH1F(Form("Photo-electron Cluster Charge chamber %i", iCh), Form("Photo-electron cluster charge chmaber %i; ADC charge; Counts/1 ADC", iCh), 180, 20., 200.); + mMatchInfoClusterMIPMap[iCh] = new TH2F(Form("MIP Cluster map chamber %i", iCh), Form("MIP Cluster map chamber %i; X (cm); Y (cm)", iCh), 133, 0, 133, 125, 0, 125); + mMatchInfoClusterMIPMap[iCh]->SetStats(0); mMatchInfoThetaCherenkovVsMom[iCh] = new TH2F(Form("ThetaCherenkov Vs Mom chamber %i", iCh), Form("ThetaCherenkov Vs Mom chamber %i; #it{p} (GeV/#it{c}); #theta_{Ckov} (rad);", iCh), 100, 0., 10., 1000, 0., 1.); } // @@ -164,6 +172,7 @@ void HmpidTaskMatches::reset() mMatchInfoResidualsYTrackMIP[iCh]->Reset(); mMatchInfoChargeClusterMIP[iCh]->Reset(); mMatchInfoChargeClusterPhotons[iCh]->Reset(); + mMatchInfoClusterMIPMap[iCh]->Reset(); mMatchInfoThetaCherenkovVsMom[iCh]->Reset(); } }