Skip to content
Merged
Changes from all commits
Commits
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
Original file line number Diff line number Diff line change
Expand Up @@ -81,7 +81,6 @@ TH1F* fhPtNegA[kIdBfNoOfSpecies] = {nullptr};
TH2F* fhNPosNegA[kIdBfNoOfSpecies] = {nullptr};
TH1F* fhDeltaNA[kIdBfNoOfSpecies] = {nullptr};


TH1F* fhEtaB = nullptr;
TH1F* fhEtaA = nullptr;

Expand Down Expand Up @@ -617,7 +616,7 @@ struct IdentifiedBfFilterTracks {
fhEtaB = new TH1F("fHistEtaB", "#eta distribution for reconstructed before;#eta;counts", 40, -2.0, 2.0);
fhEtaA = new TH1F("fHistEtaA", "#eta distribution for reconstructed;#eta;counts", etabins, etalow, etaup);
fhPhiB = new TH1F("fHistPhiB", "#phi distribution for reconstructed before;#phi;counts", 360, 0.0, constants::math::TwoPI);
fhdEdxB = new TH2F("fHistdEdxB", "dE/dx vs P before; dE/dx (a.u.); P (GeV/c)",1000, 0.0, 1000.0, 100, 0.0, 15.0);
fhdEdxB = new TH2F("fHistdEdxB", "dE/dx vs P before; dE/dx (a.u.); P (GeV/c)", 1000, 0.0, 1000.0, 100, 0.0, 15.0);
fhPhiA = new TH1F("fHistPhiA", "#phi distribution for reconstructed;#phi;counts", 360, 0.0, constants::math::TwoPI);
fhDCAxyB = new TH1F("DCAxyB", "DCA_{xy} distribution for reconstructed before;DCA_{xy} (cm);counts", 1000, -4.0, 4.0);
fhDCAxyA = new TH1F("DCAxyA", "DCA_{xy} distribution for reconstructed;DCA_{xy} (cm);counts", 1000, -4., 4.0);
Expand Down Expand Up @@ -653,10 +652,9 @@ struct IdentifiedBfFilterTracks {
fhDeltaNA[sp] = new TH1F(TString::Format("fhDeltaNA_%s", speciesName[sp]).Data(),
TString::Format("N(%s^{#plus}) #minus N(%s^{#minus}) distribution for reconstructed;N(%s^{#plus}) #minus N(%s^{#minus})", speciesTitle[sp], speciesTitle[sp], speciesTitle[sp], speciesTitle[sp]).Data(),
79, -39.5, 39.5);
fhdEdxA[sp] = new TH2F(TString::Format("fhdEdxA_%s", speciesName[sp]).Data(),
TString::Format("dE/dx vs P reconstructed %s; dE/dx (a.u.); P (GeV/c)",speciesTitle[sp]).Data(),
1000, 0.0, 1000.0, ptbins, ptlow, ptup);

fhdEdxA[sp] = new TH2F(TString::Format("fhdEdxA_%s", speciesName[sp]).Data(),
TString::Format("dE/dx vs P reconstructed %s; dE/dx (a.u.); P (GeV/c)", speciesTitle[sp]).Data(),
1000, 0.0, 1000.0, ptbins, ptlow, ptup);
}

/* add the hstograms to the output list */
Expand Down Expand Up @@ -690,7 +688,6 @@ struct IdentifiedBfFilterTracks {
fOutputList->Add(fhNPosNegA[sp]);
fOutputList->Add(fhDeltaNA[sp]);
fOutputList->Add(fhdEdxA[sp]);

}
}

Expand Down Expand Up @@ -1224,7 +1221,7 @@ void IdentifiedBfFilterTracks::fillTrackHistosBeforeSelection(TrackObject const&
fhPtB->Fill(track.pt());
fhEtaB->Fill(track.eta());
fhPhiB->Fill(track.phi());
fhdEdxB->Fill(track.tpcSignal(),track.p());
fhdEdxB->Fill(track.tpcSignal(), track.p());
if (track.sign() > 0) {
fhPtPosB->Fill(track.pt());
} else {
Expand Down Expand Up @@ -1253,7 +1250,7 @@ void IdentifiedBfFilterTracks::fillTrackHistosAfterSelection(TrackObject const&
}
fhPA[sp]->Fill(track.p());
fhPtA[sp]->Fill(track.pt());
fhdEdxA[sp]->Fill(track.tpcSignal(),track.p());
fhdEdxA[sp]->Fill(track.tpcSignal(), track.p());
if (track.sign() > 0) {
fhPtPosA[sp]->Fill(track.pt());
} else {
Expand Down