Skip to content

Commit

Permalink
Merge pull request #23870 from EwaGlimos/master
Browse files Browse the repository at this point in the history 
PWGGA/ConvBase: add eta prime cuts
  • Loading branch information
@alibuild
alibuild authored Feb 17, 2025
2 parents 1f8287a + b7f5f2d commit 69ff443
Show file tree
Hide file tree
Showing 4 changed files with 101 additions and 6 deletions.
2 changes: 1 addition & 1 deletion PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiNeutralMeson.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -7791,7 +7791,7 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiNeutralMeson::CalculateBackground(Int_
* Mode 2 => pi+ and pi0 from same event
* Mode 3 => pi- and pi0 from same event
* Mode 4 => no pions from same event (default)
* Mode 5 => Ligesign mixing
* Mode 5 => Likesign mixing
* Mode 6 => Sideband
* Mode 7 => Rotation around pi0
*/
Expand Down
49 changes: 48 additions & 1 deletion PWGGA/GammaConv/macros/AddTask_GammaConvNeutralMesonPiPlPiMiNeutralMeson_CaloMode_pp.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -698,14 +698,61 @@ void AddTask_GammaConvNeutralMesonPiPlPiMiNeutralMeson_CaloMode_pp(
} else if( trainConfig == 264 ) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe30220000","32c51070a","0103603l00000000","0453503000000000"); // EG1

// same as default, but with mixed events
// same as default, but with mixed events, default
} else if(trainConfig == 275) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","i153503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 276){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","i153503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 277) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","i153503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with rotation
} else if(trainConfig == 278) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","ir53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 279){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","ir53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 280) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","ir53503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with mixed events, likesign
} else if(trainConfig == 281) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","ia53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 282){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","ia53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 283) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","ia53503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with mixed events, pi+pi- same event
} else if(trainConfig == 284) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","io53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 285){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","io53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 286) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","io53503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with mixed events, sideband right side of eta
} else if(trainConfig == 287) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","ie53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 288){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","ie53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 289) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","ie53503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with mixed events, sideband left side of eta
} else if(trainConfig == 290) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","if53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 291){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","if53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 292) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","if53503000000000"); // EG1, pT cut on one eta daugher

// same as default, but with mixed events, sideband both side of eta
} else if(trainConfig == 293) { // possible new defaults, EDC 13 TeV, 2 sigma for mass cut
cuts.AddCutHeavyMesonCalo("00010113","411790109fe309v0000","32c510796","h000613100000000","ig53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 294){
cuts.AddCutHeavyMesonCalo("0008e113","411790109fe309v0000","32c510796","h000693100000000","ig53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 295) {
cuts.AddCutHeavyMesonCalo("0008d113","411790109fe309v0000","32c510796","h0006a3100000000","ig53503000000000"); // EG1, pT cut on one eta daugher

// Eta Prime, 13 TeV, Systematics

Expand Down
49 changes: 48 additions & 1 deletion PWGGA/GammaConv/macros/AddTask_GammaConvNeutralMesonPiPlPiMiNeutralMeson_MixedMode_pp.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -732,15 +732,62 @@ AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","24466190wa01cc00000","32c510700","0103603l00000000","0453503000000000"); // INT7
cuts.AddCutHeavyMesonPCMCalo("00062113","0dm00009f9730000dge0404000","24466190wa01cc00000","32c510700","0103603l00000000","0453503000000000"); // PHI7

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed events
// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event (all particles from diff event, default)
} else if(trainConfig == 275) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","j153503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 276) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","j153503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 277) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","j153503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, rotation
} else if(trainConfig == 278) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","jr53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 279) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","jr53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 280) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","jr53503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event likesign
} else if(trainConfig == 281) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","ja53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 282) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","ja53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 283) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","ja53503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event, pi+pi- same event
} else if(trainConfig == 284) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","jo53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 285) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","jo53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 286) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","jo53503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event, sideband right side of eta
} else if(trainConfig == 287) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","je53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 288) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","je53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 289) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","je53503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event, sideband left side of eta
} else if(trainConfig == 290) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","jf53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 291) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","jf53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 292) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","jf53503000000000"); // EG1, pT cut on one eta daugher

// possible new default pp 13 TeV, 2 sigmas mass cut, mixed event, sideband both sides of eta
} else if(trainConfig == 293) {
cuts.AddCutHeavyMesonPCMCalo("00010113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h000643600000000","jf53503000000000"); // INT7, pT cut on eta
} else if(trainConfig == 294) {
cuts.AddCutHeavyMesonPCMCalo("0008e113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006c3600000000","jf53503000000000"); // EG2, pT cut on one eta daugher
} else if(trainConfig == 295) {
cuts.AddCutHeavyMesonPCMCalo("0008d113","0dm00009f9730000dge0404000","411790109fe309v0000","32c510796","h0006d3600000000","jf53503000000000"); // EG1, pT cut on one eta daugher

// Eta Prime, 13 TeV, Systematics

// based on cutstrings:
Expand Down
7 changes: 4 additions & 3 deletions PWGGA/GammaConvBase/AliConversionMesonCuts.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3813,6 +3813,7 @@ Bool_t AliConversionMesonCuts::SetBackgroundScheme(Int_t BackgroundScheme){
fDoBGProbability = kFALSE;
fBackgroundUseLikeSign = kTRUE;
fBackgroundUseSideband = kFALSE;
fBackgroundMode = 5;
break;
case 11: // b same event pi0 sideband candidates (right side of pi0 peak)
fBackgroundMode = 6;
Expand Down Expand Up @@ -3847,7 +3848,7 @@ Bool_t AliConversionMesonCuts::SetBackgroundScheme(Int_t BackgroundScheme){
fSidebandMixingRightLow = 0.180;
fSidebandMixingRightHigh = 0.220;
break;
case 14: //e same event with pi0 sideband candidates (right side of pi0 peak)
case 14: //e same event with eta sideband candidates (right side of eta peak)
fBackgroundMode = 6;
fUseRotationMethodInBG = kFALSE;
fUseTrackMultiplicityForBG = kFALSE;
Expand All @@ -3857,7 +3858,7 @@ Bool_t AliConversionMesonCuts::SetBackgroundScheme(Int_t BackgroundScheme){
fSidebandMixingLow = 0.600;
fSidebandMixingHigh = 0.650;
break;
case 15: //f same event with pi0 sideband candidates (left side of pi0 peak)
case 15: //f same event with eta sideband candidates (left side of eta peak)
fBackgroundMode = 6;
fUseRotationMethodInBG = kFALSE;
fUseTrackMultiplicityForBG = kFALSE;
Expand All @@ -3867,7 +3868,7 @@ Bool_t AliConversionMesonCuts::SetBackgroundScheme(Int_t BackgroundScheme){
fSidebandMixingLow = 0.42;
fSidebandMixingHigh = 0.47;
break;
case 16: //g same event with pi0 sideband candidates (both sides of pi0 peak)
case 16: //g same event with eta sideband candidates (both sides of eta peak)
fBackgroundMode = 6;
fUseRotationMethodInBG = kFALSE;
fUseTrackMultiplicityForBG = kFALSE;
Expand Down

0 comments on commit 69ff443

Please sign in to comment.