Skip to content

Commit

Permalink
Added histograms for MC reconstructed
Browse files Browse the repository at this point in the history
  • Loading branch information
@antpaz
antpaz committed Apr 21, 2023
1 parent 06f76b6 commit 76ae5e8
Show file tree
Hide file tree
Showing 2 changed files with 46 additions and 19 deletions.
63 changes: 44 additions & 19 deletions PWGMM/UE/UeKNO/AliAnalysisTaskFlatenicity.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -118,20 +118,24 @@ ClassImp(AliAnalysisTaskFlatenicity) // classimp: necessary for root
fMultSelection(0x0), hPtPrimIn(0), hPtPrimOut(0), hPtSecOut(0), hPtOut(0),
hFlatV0vsFlatTPC(0), hFlatenicityBefore(0), hFlatenicity(0),
hFlatenicityMC(0), hFlatResponse(0), hFlatVsPt(0), hFlatVsPtMC(0),
hActivityV0DataSectBefore(0), hActivityV0DataSect(0), hV0vsVtxz(0),
hActivityV0McSect(0), hFlatVsNchMC(0), hFlatVsV0M(0), hFlatMCVsV0M(0),
hEtamc(0), hEtamcAlice(0), hCounter(0), hCountProduV0m(0),
hCountAuthV0m(0), hCountProdu_FlatMC(0), hCountAuth_FlatMC(0),
hMultMCmVsV0M(0), hMultMCaVsV0M(0), hMultMCcVsV0M(0), hMultmVsV0M(0),
hMultmVsV0Malice(0), hMultaVsV0M(0), hMultcVsV0M(0), hV0MBadruns(0),
hChgProdu_All_pt(0), hChgAuth_All_pt(0), hChgProdu_pt_V0(0),
hChgAuth_pt_V0(0), hChgProdu_pt_Flat(0), hChgAuth_pt_Flat(0) {
hRecFlatVsPtMC(0), hActivityV0DataSectBefore(0), hActivityV0DataSect(0),
hV0vsVtxz(0), hActivityV0McSect(0), hFlatVsNchMC(0), hFlatVsV0M(0),
hFlatMCVsV0M(0), hEtamc(0), hEtamcAlice(0), hCounter(0),
hCountProduV0m(0), hCountAuthV0m(0), hCountProdu_FlatMC(0),
hCountAuth_FlatMC(0), hMultMCmVsV0M(0), hMultMCaVsV0M(0),
hMultMCcVsV0M(0), hMultmVsV0M(0), hMultmVsV0Malice(0), hMultaVsV0M(0),
hMultcVsV0M(0), hV0MBadruns(0), hChgProdu_All_pt(0), hChgAuth_All_pt(0),
hChgProdu_pt_V0(0), hChgAuth_pt_V0(0), hChgProdu_pt_Flat(0),
hChgAuth_pt_Flat(0) {
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hFlatVsPtV0M[i_c] = 0;
}
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hFlatVsPtV0MMC[i_c] = 0;
}
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hRecFlatVsPtV0MMC[i_c] = 0;
}
for (Int_t i_d = 0; i_d < nDet; ++i_d) {
hComponentsMult[i_d] = 0;
}
Expand Down Expand Up @@ -168,14 +172,15 @@ AliAnalysisTaskFlatenicity::AliAnalysisTaskFlatenicity(const char *name)
fMultSelection(0x0), hPtPrimIn(0), hPtPrimOut(0), hPtSecOut(0), hPtOut(0),
hFlatV0vsFlatTPC(0), hFlatenicityBefore(0), hFlatenicity(0),
hFlatenicityMC(0), hFlatResponse(0), hFlatVsPt(0), hFlatVsPtMC(0),
hActivityV0DataSectBefore(0), hActivityV0DataSect(0), hV0vsVtxz(0),
hActivityV0McSect(0), hFlatVsNchMC(0), hFlatVsV0M(0), hFlatMCVsV0M(0),
hEtamc(0), hEtamcAlice(0), hCounter(0), hCountProduV0m(0),
hCountAuthV0m(0), hCountProdu_FlatMC(0), hCountAuth_FlatMC(0),
hMultMCmVsV0M(0), hMultMCaVsV0M(0), hMultMCcVsV0M(0), hMultmVsV0M(0),
hMultmVsV0Malice(0), hMultaVsV0M(0), hMultcVsV0M(0), hV0MBadruns(0),
hChgProdu_All_pt(0), hChgAuth_All_pt(0), hChgProdu_pt_V0(0),
hChgAuth_pt_V0(0), hChgProdu_pt_Flat(0), hChgAuth_pt_Flat(0)
hRecFlatVsPtMC(0), hActivityV0DataSectBefore(0), hActivityV0DataSect(0),
hV0vsVtxz(0), hActivityV0McSect(0), hFlatVsNchMC(0), hFlatVsV0M(0),
hFlatMCVsV0M(0), hEtamc(0), hEtamcAlice(0), hCounter(0),
hCountProduV0m(0), hCountAuthV0m(0), hCountProdu_FlatMC(0),
hCountAuth_FlatMC(0), hMultMCmVsV0M(0), hMultMCaVsV0M(0),
hMultMCcVsV0M(0), hMultmVsV0M(0), hMultmVsV0Malice(0), hMultaVsV0M(0),
hMultcVsV0M(0), hV0MBadruns(0), hChgProdu_All_pt(0), hChgAuth_All_pt(0),
hChgProdu_pt_V0(0), hChgAuth_pt_V0(0), hChgProdu_pt_Flat(0),
hChgAuth_pt_Flat(0)

{
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
Expand All @@ -184,6 +189,9 @@ AliAnalysisTaskFlatenicity::AliAnalysisTaskFlatenicity(const char *name)
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hFlatVsPtV0MMC[i_c] = 0;
}
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hRecFlatVsPtV0MMC[i_c] = 0;
}
for (Int_t i_d = 0; i_d < nDet; ++i_d) {
hComponentsMult[i_d] = 0;
}
Expand Down Expand Up @@ -352,6 +360,11 @@ void AliAnalysisTaskFlatenicity::UserCreateOutputObjects() {
min_flat, max_flat, 100, -0.5, 99.5);
fOutputList->Add(hFlatVsNchMC);

hRecFlatVsPtMC = new TH2D("hRecFlatVsPtMC",
"MC rec; Flatenicity; #it{p}_{T} (GeV/#it{c})",
nbins_flat, min_flat, max_flat, nPtbins, Ptbins);
fOutputList->Add(hRecFlatVsPtMC);

for (Int_t i_d = 0; i_d < nDet; ++i_d) {
hComponentsMultmc[i_d] = new TH2D(Form("hTrueMult_%s", DetName[i_d]), "",
600, -0.5, 599.0, 200, -0.5, 199.5);
Expand All @@ -369,6 +382,7 @@ void AliAnalysisTaskFlatenicity::UserCreateOutputObjects() {
499.5, 500, -0.5, 499.5);
fOutputList->Add(hRmCombinedMult[i_c]);
}

for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hFlatVsPtV0MMC[i_c] =
new TH2D(Form("hFlatVsPtV0MMC_c%d", i_c),
Expand All @@ -378,6 +392,15 @@ void AliAnalysisTaskFlatenicity::UserCreateOutputObjects() {
nbins_flat, min_flat, max_flat, nPtbins, Ptbins);
fOutputList->Add(hFlatVsPtV0MMC[i_c]);
}
for (Int_t i_c = 0; i_c < nCent; ++i_c) {
hRecFlatVsPtV0MMC[i_c] =
new TH2D(Form("hRecFlatVsPtV0MMC_c%d", i_c),
Form("Measured %1.0f-%1.0f%%V0M; rec Flatenicity; "
"#it{p}_{T} (GeV/#it{c})",
centClass[i_c], centClass[i_c + 1]),
nbins_flat, min_flat, max_flat, nPtbins, Ptbins);
fOutputList->Add(hRecFlatVsPtV0MMC[i_c]);
}

hFlatMCVsV0M = new TH2D("hFlatMCVsV0M", "", nCent, centClass, nbins_flat,
min_flat, max_flat);
Expand Down Expand Up @@ -581,6 +604,7 @@ void AliAnalysisTaskFlatenicity::UserExec(Option_t *) {
vector<Float_t> ptRecon;
vector<Int_t> idRecon;
vector<Int_t> isprimRecon;

fnRecon = FillArray(ptRecon, idRecon, isprimRecon);

// Trigger selection
Expand All @@ -606,6 +630,7 @@ void AliAnalysisTaskFlatenicity::UserExec(Option_t *) {

vector<Float_t> ptDetMC;
vector<Int_t> idDetMC;

if (isGoodVtxPosMC) {
fnDetec = FillMCarray(ptDetMC, idDetMC);
GetMCchargedDetDists(fnDetec, ptDetMC, idDetMC);
Expand Down Expand Up @@ -820,10 +845,13 @@ void AliAnalysisTaskFlatenicity::MakeMCanalysis() {
if (esdtrack->Pt() < fPtMin)
continue;
hPtOut->Fill(esdtrack->Pt());

Int_t mcLabel = -1;
mcLabel = TMath::Abs(esdtrack->GetLabel());
if (fMC->IsPhysicalPrimary(mcLabel)) {
hPtPrimOut->Fill(esdtrack->Pt());
hRecFlatVsPtMC->Fill(fFlatMC, esdtrack->Pt());
hRecFlatVsPtV0MMC[fV0Mindex]->Fill(fFlatMC, esdtrack->Pt());
} else {
hPtSecOut->Fill(esdtrack->Pt());
}
Expand Down Expand Up @@ -1205,7 +1233,6 @@ Double_t AliAnalysisTaskFlatenicity::GetFlatenicityMC() {
}
sRho_tmp /= (1.0 * nCells * nCells);
Float_t sRho = TMath::Sqrt(sRho_tmp);

if (mRho > 0) {
if (fRemoveTrivialScaling) {
flatenicity = TMath::Sqrt(1.0 * nMult) * sRho / mRho;
Expand Down Expand Up @@ -1600,7 +1627,6 @@ void AliAnalysisTaskFlatenicity::GetMCchargedTrueDists(

if (multGen < 1)
return;
//! ap The INEL>0 condition

// hCountEvent->Fill(1.0);
hCountProduV0m->Fill(fv0mpercentile);
Expand All @@ -1619,7 +1645,6 @@ void AliAnalysisTaskFlatenicity::GetMCchargedDetDists(

if (multRec < 1)
return;
//! ap The INEL>0 condition
hCountAuthV0m->Fill(fv0mpercentile);
hCountAuth_FlatMC->Fill(fFlatAltMC);

Expand Down
2 changes: 2 additions & 0 deletions PWGMM/UE/UeKNO/AliAnalysisTaskFlatenicity.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -119,6 +119,7 @@ class AliAnalysisTaskFlatenicity : public AliAnalysisTaskSE {
TH2D *hFlatResponse;
TH2D *hFlatVsPt;
TH2D *hFlatVsPtMC;
TH2D *hRecFlatVsPtMC;
TProfile *hActivityV0DataSectBefore;
TProfile *hActivityV0DataSect;
TProfile *hV0vsVtxz;
Expand Down Expand Up @@ -150,6 +151,7 @@ class AliAnalysisTaskFlatenicity : public AliAnalysisTaskSE {
TH2D *hChgAuth_pt_Flat;
TH2D *hFlatVsPtV0M[9];
TH2D *hFlatVsPtV0MMC[9];
TH2D *hRecFlatVsPtV0MMC[9];
TH2D *hComponentsMult[4];
TH2D *hCombinedMult[3];
TH2D *hComponentsMultmc[4];
Expand Down

0 comments on commit 76ae5e8

Please sign in to comment.