Skip to content

support for NUA, part 2 (the custom case) #5935

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
victor-gonzalez merged 2 commits into from
May 2, 2024
Merged

support for NUA, part 2 (the custom case) #5935

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
a2b2303
support for NUA, part 2 (the custom case)
abilandz May 1, 2024
86425ce
fix for MegaLinter
abilandz May 1, 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 PWGCF/MultiparticleCorrelations/Core/MuPa-Configurables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -120,6 +120,8 @@ struct : ConfigurableGroup {
struct : ConfigurableGroup {
Configurable<vector<int>> cfApplyNUAPDF{"cfApplyNUAPDF", {0, 0, 0}, "Apply (1) or do not apply (0) NUA on variable, ordering is the same as in enum eNUAPDF (phi, pt, eta)"};
Configurable<vector<int>> cfUseDefaultNUAPDF{"cfUseDefaultNUAPDF", {1, 1, 1}, "Use (1) or do not use (0) default NUA profile, ordering is the same as in enum eNUAPDF (phi, pt, eta)"};
Configurable<vector<string>> cfCustomNUAPDFHistNames{"cfCustomNUAPDFHistNames", {"a", "bb", "ccc"}, "the names of histograms holding custom NUA in an external file."};
Configurable<string> cfFileWithCustomNUA{"cfFileWithCustomNUA", "/home/abilandz/DatasetsO2/customNUA.root", "path to external ROOT file which holds all histograms with custom NUA"}; // for AliEn file prepend "/alice/cern.ch/", for CCDB prepend "/alice-ccdb.cern.ch"
} cf_nua;

// *) Internal validation:
Expand Down
16 changes: 9 additions & 7 deletions PWGCF/MultiparticleCorrelations/Core/MuPa-DataMembers.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -168,13 +168,15 @@ struct NestedLoops {

// *) Toy NUA (can be applied both in real data analysis and in analysis 'on-the-fly'):
struct NUA {
TList* fNUAList = NULL; // list to hold all NUA objects
TProfile* fNUAFlagsPro = NULL; // profile to hold all flags for NUA objects
Bool_t fApplyNUAPDF[eNUAPDF_N] = {kFALSE}; // apply NUA to particular kine variable (see the corresponding enum eNUAPDF)
Bool_t fUseDefaultNUAPDF[eNUAPDF_N] = {kTRUE}; // by default, use simple hardcoded expressions for NUA acceptance profile
TF1* fDefaultNUAPDF[eNUAPDF_N] = {NULL}; // default distributions used as pdfs to simulate events on-the-fly
TH1D* fCustomNUAPDF[eNUAPDF_N] = {NULL}; // custom, user-supplied distributions used to simulate NUA
Double_t fMaxValuePDF[eNUAPDF_N] = {0.}; // see algorithm used in Accept(...). I implemented it as a data member, so that it is not calculated again and again at each particle call
TList* fNUAList = NULL; // list to hold all NUA objects
TProfile* fNUAFlagsPro = NULL; // profile to hold all flags for NUA objects
Bool_t fApplyNUAPDF[eNUAPDF_N] = {kFALSE}; // apply NUA to particular kine variable (see the corresponding enum eNUAPDF)
Bool_t fUseDefaultNUAPDF[eNUAPDF_N] = {kTRUE}; // by default, use simple hardcoded expressions for NUA acceptance profile
TF1* fDefaultNUAPDF[eNUAPDF_N] = {NULL}; // default distributions used as pdfs to simulate events on-the-fly
TH1D* fCustomNUAPDF[eNUAPDF_N] = {NULL}; // custom, user-supplied distributions used to simulate NUA
TString* fCustomNUAPDFHistNames[eNUAPDF_N] = {NULL}; // these are the names of histograms holding custom NUA in an external file. There is a configurable for this one.
TString fFileWithCustomNUA = ""; // path to external ROOT file which holds all histograms with custom NUA
Double_t fMaxValuePDF[eNUAPDF_N] = {0.}; // see algorithm used in Accept(...). I implemented it as a data member, so that it is not calculated again and again at each particle call
} nua;

// *) Internal validation:
Expand Down
Loading