@@ -174,13 +174,23 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
174174 signal = new MCSignal (name, " Prompt jpsi (not from beauty)" , {prong}, {-1 });
175175 return signal;
176176 }
177+ if (!nameStr.compare (" eeFromNonpromptPsi2S" )) {
178+ MCProng prong (2 , {11 , 100443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false }, false , {503 }, {false });
179+ signal = new MCSignal (name, " ee pairs from non-prompt psi2s decays" , {prong, prong}, {1 , 1 }); // signal at pair level
180+ return signal;
181+ }
182+ if (!nameStr.compare (" eeFromPromptPsi2S" )) {
183+ MCProng prong (2 , {11 , 100443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false }, false , {503 }, {true });
184+ signal = new MCSignal (name, " ee pairs from prompt psi2s decays" , {prong, prong}, {1 , 1 }); // signal at pair level
185+ return signal;
186+ }
177187 if (!nameStr.compare (" eFromNonpromptPsi2S" )) {
178- MCProng prong (3 , {11 , 100443 , 503 }, {true , true , true }, {false , false , false }, {0 , 0 , 0 }, {0 , 0 , 0 }, {false , false , false });
188+ MCProng prong (2 , {11 , 100443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false }, false , { 503 }, { false });
179189 signal = new MCSignal (name, " Electrons from beauty psi2s decays" , {prong}, {-1 });
180190 return signal;
181191 }
182192 if (!nameStr.compare (" eFromPromptPsi2S" )) {
183- MCProng prong (3 , {11 , 100443 , 503 }, {true , true , true }, {false , false , true }, {0 , 0 , 0 }, {0 , 0 , 0 }, {false , false , false });
193+ MCProng prong (2 , {11 , 100443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false } , false , { 503 }, { true });
184194 signal = new MCSignal (name, " Electrons from prompt psi2s decays" , {prong}, {-1 });
185195 return signal;
186196 }
@@ -190,7 +200,7 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
190200 return signal;
191201 }
192202 if (!nameStr.compare (" nonPromptPsi2S" )) {
193- MCProng prong (2 , {100443 , 503 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false });
203+ MCProng prong (2 , {100443 , 503 }, {true , true }, {false , true }, {0 , 0 }, {0 , 0 }, {false , false });
194204 signal = new MCSignal (name, " Non-prompt psi2s" , {prong}, {-1 });
195205 return signal;
196206 }
@@ -487,8 +497,7 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
487497 return signal;
488498 }
489499 if (!nameStr.compare (" eFromPromptJpsi" )) {
490- MCProng prong (2 , {11 , 443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false });
491- prong.SetSourceBit (0 , MCProng::kPhysicalPrimary );
500+ MCProng prong (2 , {11 , 443 }, {true , true }, {false , false }, {0 , 0 }, {0 , 0 }, {false , false }, false , {503 }, {true });
492501 signal = new MCSignal (name, " Electrons from jpsi decays" , {prong}, {-1 });
493502 return signal;
494503 }
0 commit comments