Skip to content

HMP: Adding more plots and removing spaces from TaskMatches plot names #2142

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
knopers8 merged 2 commits into from
Feb 21, 2024
Merged

HMP: Adding more plots and removing spaces from TaskMatches plot names #2142

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
8cac6bd
Adding more QC plots
gvolpe79 Feb 19, 2024
8c17023
Minors
gvolpe79 Feb 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions Modules/HMPID/include/HMPID/HmpidTaskClusters.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,6 +55,8 @@ class HmpidTaskClusters : public TaskInterface
TH1F* hClusMultEv;
TH1F* hHMPIDchargeClus[7];
TH1F* hHMPIDchargeMipClus[7];
TH1F* hHMPIDclusX[7];
TH1F* hHMPIDclusY[7];
TH2F* hHMPIDpositionClus[7];

std::vector<TObject*> mPublishedObjects;
Expand Down
4 changes: 4 additions & 0 deletions Modules/HMPID/include/HMPID/HmpidTaskMatches.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -61,6 +61,10 @@ class HmpidTaskMatches : public TaskInterface
TH1F* mMatchInfoResidualsYTrackMIP[7];
TH1F* mMatchInfoChargeClusterMIP[7];
TH1F* mMatchInfoChargeClusterPhotons[7];
TH1F* mMatchInfoXMip[7];
TH1F* mMatchInfoYMip[7];
TH1F* mMatchInfoXTrack[7];
TH1F* mMatchInfoYTrack[7];
TH2F* mMatchInfoClusterMIPMap[7];
TH2F* mMatchInfoThetaCherenkovVsMom[7];

Expand Down
21 changes: 13 additions & 8 deletions Modules/HMPID/src/HmpidTaskClusters.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,6 +44,8 @@ HmpidTaskClusters::~HmpidTaskClusters()
for (Int_t iCh = 0; iCh < 7; iCh++) {
delete hHMPIDchargeClus[iCh];
delete hHMPIDchargeMipClus[iCh];
delete hHMPIDclusX[7];
delete hHMPIDclusY[7];
delete hHMPIDpositionClus[iCh];
}
}
Expand All @@ -64,6 +66,8 @@ void HmpidTaskClusters::initialize(o2::framework::InitContext& /*ctx*/)
getObjectsManager()->startPublishing(hHMPIDchargeClus[iCh]);
getObjectsManager()->startPublishing(hHMPIDchargeMipClus[iCh]);
getObjectsManager()->startPublishing(hHMPIDpositionClus[iCh]);
getObjectsManager()->startPublishing(hHMPIDclusX[iCh]);
getObjectsManager()->startPublishing(hHMPIDclusY[iCh]);
getObjectsManager()->setDefaultDrawOptions(hHMPIDpositionClus[iCh], "colz");
getObjectsManager()->setDisplayHint(hHMPIDpositionClus[iCh], "colz");
}
Expand All @@ -75,13 +79,11 @@ void HmpidTaskClusters::initialize(o2::framework::InitContext& /*ctx*/)

void HmpidTaskClusters::startOfActivity(const Activity& /*activity*/)
{
ILOG(Debug, Devel) << "startOfActivity" << ENDM;
reset();
}

void HmpidTaskClusters::startOfCycle()
{
ILOG(Debug, Devel) << "startOfCycle" << ENDM;
}

void HmpidTaskClusters::monitorData(o2::framework::ProcessingContext& ctx)
Expand Down Expand Up @@ -113,6 +115,8 @@ void HmpidTaskClusters::monitorData(o2::framework::ProcessingContext& ctx)
if (clusters[j].size() >= 3 && clusters[j].size() <= 7) {
hHMPIDchargeMipClus[chamber]->Fill(clusters[j].q());
}
hHMPIDclusX[chamber]->Fill(clusters[j].x());
hHMPIDclusY[chamber]->Fill(clusters[j].y());
hHMPIDpositionClus[chamber]->Fill(clusters[j].x(), clusters[j].y());
}
} // cluster loop
Expand All @@ -124,18 +128,15 @@ void HmpidTaskClusters::monitorData(o2::framework::ProcessingContext& ctx)

void HmpidTaskClusters::endOfCycle()
{

ILOG(Debug, Devel) << "endOfCycle" << ENDM;
}

void HmpidTaskClusters::endOfActivity(const Activity& /*activity*/)
{
ILOG(Debug, Devel) << "endOfActivity" << ENDM;
}

void HmpidTaskClusters::BookHistograms()
{
hClusMultEv = new TH1F("ClusMultEve", "HMPID Cluster multiplicity per event", 500, 0, 500);
hClusMultEv = new TH1F("ClusMultEve", "HMPID Cluster multiplicity per event", 1000, 0, 1000);

ThClusMult = new TProfile("ClusMult", "HMPID Cluster multiplicity per chamber;Chamber Id;# of clusters/event", 7, 0., 7.);
ThClusMult->Sumw2();
Expand All @@ -147,8 +148,10 @@ void HmpidTaskClusters::BookHistograms()
ThClusMult->SetStats(0);

for (Int_t iCh = 0; iCh < 7; iCh++) {
hHMPIDchargeClus[iCh] = new TH1F(Form("CluQ%i", iCh), Form("Cluster Charge in HMPID Chamber %i ; Entries/1 ADC; Q (ADC)", iCh), 2000, 0, 2000);
hHMPIDchargeMipClus[iCh] = new TH1F(Form("MipCluQ%i", iCh), Form("MIP Cluster Charge in HMPID Chamber %i;Entries/1 ADC;Q (ADC)", iCh), 2000, 200, 2200);
hHMPIDchargeClus[iCh] = new TH1F(Form("CluQ%i", iCh), Form("Cluster Charge in HMPID Chamber %i; Q (ADC); Entries/1 ADC", iCh), 2000, 0, 2000);
hHMPIDchargeMipClus[iCh] = new TH1F(Form("MipCluQ%i", iCh), Form("MIP Cluster Charge in HMPID Chamber %i; Q (ADC); Entries/1 ADC", iCh), 2000, 200, 2200);
hHMPIDclusX[iCh] = new TH1F(Form("ClusX%i", iCh), Form("X Cluster in HMPID Chamber %i; X (cm); Entries/1 cm", iCh), 133, 0, 133);
hHMPIDclusY[iCh] = new TH1F(Form("ClusY%i", iCh), Form("Y Cluster in HMPID Chamber %i; Y (cm); Entries/1 cm", iCh), 125, 0, 125);
hHMPIDpositionClus[iCh] = new TH2F(Form("ClusterPoistion%i", iCh), Form("Cluster Position in HMPID Chamber %i; X (cm); Y (cm)", iCh), 133, 0, 133, 125, 0, 125); // cmxcm
hHMPIDpositionClus[iCh]->SetStats(0);
}
Expand All @@ -165,6 +168,8 @@ void HmpidTaskClusters::reset()
for (Int_t iCh = 0; iCh < 7; iCh++) {
hHMPIDchargeClus[iCh]->Reset();
hHMPIDchargeMipClus[iCh]->Reset();
hHMPIDclusX[iCh]->Reset();
hHMPIDclusY[iCh]->Reset();
hHMPIDpositionClus[iCh]->Reset();
}
}
Expand Down
32 changes: 26 additions & 6 deletions Modules/HMPID/src/HmpidTaskMatches.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,10 @@ HmpidTaskMatches::~HmpidTaskMatches()
delete mMatchInfoResidualsYTrackMIP[iCh];
delete mMatchInfoChargeClusterMIP[iCh];
delete mMatchInfoChargeClusterPhotons[iCh];
delete mMatchInfoXMip[iCh];
delete mMatchInfoYMip[iCh];
delete mMatchInfoXTrack[iCh];
delete mMatchInfoYTrack[iCh];
delete mMatchInfoClusterMIPMap[iCh];
delete mMatchInfoThetaCherenkovVsMom[iCh];
}
Expand All @@ -66,6 +70,10 @@ void HmpidTaskMatches::initialize(o2::framework::InitContext& /*ctx*/)
getObjectsManager()->startPublishing(mMatchInfoClusterMIPMap[iCh]);
getObjectsManager()->setDefaultDrawOptions(mMatchInfoClusterMIPMap[iCh], "colz");
getObjectsManager()->setDisplayHint(mMatchInfoClusterMIPMap[iCh], "colz");
getObjectsManager()->startPublishing(mMatchInfoXMip[iCh]);
getObjectsManager()->startPublishing(mMatchInfoYMip[iCh]);
getObjectsManager()->startPublishing(mMatchInfoXTrack[iCh]);
getObjectsManager()->startPublishing(mMatchInfoXTrack[iCh]);
getObjectsManager()->startPublishing(mMatchInfoThetaCherenkovVsMom[iCh]);
}
}
Expand Down Expand Up @@ -112,6 +120,10 @@ void HmpidTaskMatches::monitorData(o2::framework::ProcessingContext& ctx)
if (matchHMP.getMipClusSize() < 20 && matchHMP.getMipClusSize() > 1) {
mMatchInfoChargeClusterMIP[ch]->Fill(q);
}
mMatchInfoXMip[ch]->Fill(xMip);
mMatchInfoYMip[ch]->Fill(yMip);
mMatchInfoXTrack[ch]->Fill(xTrk);
mMatchInfoYTrack[ch]->Fill(yTrk);
mMatchInfoClusterMIPMap[ch]->Fill(xMip, yMip);
mMatchInfoThetaCherenkovVsMom[ch]->Fill(TMath::Abs(matchHMP.getHmpMom()), matchHMP.getHMPsignal());
} // match info loop
Expand All @@ -129,13 +141,17 @@ void HmpidTaskMatches::BookHistograms()
{
for (int iCh = 0; iCh < 7; iCh++) {

mMatchInfoResidualsXTrackMIP[iCh] = new TH1F(Form("X Residuals chamber %i", iCh), Form("X Residuals chamber %i; X residuals (cm); Counts/1 cm", iCh), 100, -50., 50.);
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);
mMatchInfoResidualsXTrackMIP[iCh] = new TH1F(Form("XResCham%i", iCh), Form("X Residuals chamber %i; X residuals (cm); Entries/1 cm", iCh), 100, -50., 50.);
mMatchInfoResidualsYTrackMIP[iCh] = new TH1F(Form("YRescham%i", iCh), Form("Y Residuals chamber %i; Y residuals (cm); Entries/1 cm", iCh), 100, -50., 50.);
mMatchInfoChargeClusterMIP[iCh] = new TH1F(Form("MipClusQch%i", iCh), Form("MIP Cluster Charge chamber %i; ADC charge; Entries/1 ADC", iCh), 2000, 200., 2200.);
mMatchInfoChargeClusterPhotons[iCh] = new TH1F(Form("PhelQch%i", iCh), Form("Photo-electron cluster charge chmaber %i; ADC charge; Counts/1 ADC", iCh), 180, 20., 200.);
mMatchInfoClusterMIPMap[iCh] = new TH2F(Form("MipClusMap%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.);
mMatchInfoXMip[iCh] = new TH1F(Form("MipX%i", iCh), Form("X MIP in HMPID Chamber %i; X (cm); Entries/1 cm", iCh), 133, 0, 133);
mMatchInfoYMip[iCh] = new TH1F(Form("MipY%i", iCh), Form("Y MIP in HMPID Chamber %i; Y (cm); Entries/1 cm", iCh), 125, 0, 125);
mMatchInfoXTrack[iCh] = new TH1F(Form("TrackX%i", iCh), Form("X Track in HMPID Chamber %i; X (cm); Entries/1 cm", iCh), 133, 0, 133);
mMatchInfoYTrack[iCh] = new TH1F(Form("TrackY%i", iCh), Form("Y Track in HMPID Chamber %i; Y (cm); Entries/1 cm", iCh), 125, 0, 125);
mMatchInfoThetaCherenkovVsMom[iCh] = new TH2F(Form("ThetavsMom%i", iCh), Form("ThetaCherenkov Vs Mom chamber %i; #it{p} (GeV/#it{c}); #theta_{Ckov} (rad);", iCh), 100, 0., 10., 1000, 0., 1.);
}
//
}
Expand All @@ -152,6 +168,10 @@ void HmpidTaskMatches::reset()
mMatchInfoChargeClusterMIP[iCh]->Reset();
mMatchInfoChargeClusterPhotons[iCh]->Reset();
mMatchInfoClusterMIPMap[iCh]->Reset();
mMatchInfoXMip[iCh]->Reset();
mMatchInfoYMip[iCh]->Reset();
mMatchInfoXTrack[iCh]->Reset();
mMatchInfoYTrack[iCh]->Reset();
mMatchInfoThetaCherenkovVsMom[iCh]->Reset();
}
}
Expand Down