-// Generating Pion Physics
- AliGenParam * genpion = new AliGenParam(1, AliGenMUONlib::kPion, "Vogt", "Pion");
- genpion->SetPtRange(0,100.);
- genpion->SetYRange(-8.,8);
- genpion->SetPhiRange(0.,360.);
- genpion->SetForceDecay(kPiToMu);
- genpion->SetTrackingFlag(1);
- Float_t ratiopion; // Ratio with respect to the reaction cross-section for the muonic channel in the kinematics limit of the MUONCocktail
- Float_t sigmapion = 0.93e-2; // Valerie presentation Clermont-16-jan-2004 and Alice-int-2002-06
- Float_t brpion = 0.9999; // Branching Ratio for Pion
- genpion->Init(); // Generating pT and Y parametrsation for the 4pi
- ratiopion = sigmapion * brpion * (0.80*fNumberOfParticipants+0.2*fNumberOfCollisions) / sigmaReaction * genpion->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
- genpion->SetPtRange(ptMin, ptMax);
- genpion->SetYRange(yMin, yMax);
- genpion->SetPhiRange(phiMin, phiMax);
- genpion->Init(); // Generating pT and Y parametrsation in the desired kinematic range
- printf(">>> ratio pion %f\n",ratiopion);
- AddGenerator(genpion,"Pion", ratiopion);
- fTotalRate+=ratiopion;
-
-// Generating Kaon Physics
- AliGenParam * genkaon = new AliGenParam(1, AliGenMUONlib::kKaon, "Vogt", "Kaon");
- genkaon->SetPtRange(0,100.);
- genkaon->SetYRange(-8.,8);
- genkaon->SetPhiRange(0.,360.);
- genkaon->SetForceDecay(kKaToMu);
- genkaon->SetTrackingFlag(1);
- Float_t ratiokaon; // Ratio with respect to the reaction cross-section for the muonic channel in the kinematics limit of the MUONCocktail
- Float_t sigmakaon = 1.23e-4; // Valerie presentation Clermont-16-jan-2004 and Alice-int-2002-06
- Float_t brkaon = 0.6351 ; // Branching Ratio for Kaon
- genkaon->Init(); // Generating pT and Y parametrsation for the 4pi
- ratiokaon = sigmakaon * brkaon * (0.80*fNumberOfParticipants+0.2*fNumberOfCollisions)/ sigmaReaction * genkaon->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
- genkaon->SetPtRange(ptMin, ptMax);
- genkaon->SetYRange(yMin, yMax);
- genkaon->SetPhiRange(phiMin, phiMax);
- genkaon->Init(); // Generating pT and Y parametrsation in the desired kinematic range
- printf(">>> ratio kaon %f\n",ratiokaon);
- AddGenerator(genkaon,"Kaon", ratiokaon);
- fTotalRate+=ratiokaon;
+ // Only if hadronic muons are included in the cocktail
+ if(fHadronicMuons) {
+ // Generating Pion Physics
+ // The scaling with Npart and Ncoll has been obtained to reproduced tha values presented by Valeri lors de presentatation
+ // a Clermont Ferrand http://pcrochet.home.cern.ch/pcrochet/files/valerie.pdf
+ // b range(fm) Ncoll Npart N_mu pT>0.4 GeV/c
+ // 0 - 3 1982 381 3.62
+ // 3 - 6 1388 297 3.07
+ // 6 - 9 674 177 1.76
+ // 9 - 12 188 71 0.655
+ // 12 - 16 15 10 0.086
+ // We found the hadronic muons scales quite well with the number of participants
+ AliGenParam * genpion = new AliGenParam(1, AliGenMUONlib::kPion, "default", "Pion");
+ genpion->SetPtRange(0,100.);
+ genpion->SetYRange(-8.,8);
+ genpion->SetPhiRange(0.,360.);
+ genpion->SetForceDecay(kPiToMu);
+ genpion->SetTrackingFlag(1);
+ Double_t ratiopion; // Ratio with respect to the reaction cross-section for the muonic channel in the kinematics limit of the MUONCocktail
+ Double_t sigmapion = 1.80e-2; // Just for reproducing Valeries's data
+ Double_t brpion = 0.9999; // Branching Ratio for Pion
+ genpion->Init(); // Generating pT and Y parametrsation for the 4pi
+ ratiopion = sigmapion * brpion * (0.93*fNumberOfParticipants+0.07*fNumberOfCollisions) / sigmaReaction * genpion->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
+ AliInfo(Form("Pseudo-Pion production cross-section in pp with shadowing %5.3g barns",sigmapion));
+ AliInfo(Form("Pion production probability per collisions in acceptance via the muonic channel %5.3g",ratiopion));
+ genpion->SetPtRange(ptMin, ptMax);
+ genpion->SetYRange(yMin, yMax);
+ genpion->SetPhiRange(phiMin, phiMax);
+ genpion->Init(); // Generating pT and Y parametrsation in the desired kinematic range
+ AddGenerator(genpion,"Pion", ratiopion);
+ fTotalRate+=ratiopion;
+
+ // Generating Kaon Physics
+ AliGenParam * genkaon = new AliGenParam(1, AliGenMUONlib::kKaon, "default", "Kaon");
+ genkaon->SetPtRange(0,100.);
+ genkaon->SetYRange(-8.,8);
+ genkaon->SetPhiRange(0.,360.);
+ genkaon->SetForceDecay(kKaToMu);
+ genkaon->SetTrackingFlag(1);
+ Double_t ratiokaon; // Ratio with respect to the reaction cross-section for the muonic channel in the kinematics limit of the MUONCocktail
+ Double_t sigmakaon = 2.40e-4; // Valerie presentation Clermont-16-jan-2004 and Alice-int-2002-06
+ Double_t brkaon = 0.6351 ; // Branching Ratio for Kaon
+ genkaon->Init(); // Generating pT and Y parametrsation for the 4pi
+ ratiokaon = sigmakaon * brkaon * (0.93*fNumberOfParticipants+0.07*fNumberOfCollisions)/ sigmaReaction * genkaon->GetRelativeArea(ptMin,ptMax,yMin,yMax,phiMin,phiMax);
+ AliInfo(Form("Pseudo-kaon production cross-section in pp with shadowing %5.3g barns",sigmakaon));
+ AliInfo(Form("Kaon production probability per collisions in acceptance via the muonic channel %5.3g",ratiokaon));
+ genkaon->SetPtRange(ptMin, ptMax);
+ genkaon->SetYRange(yMin, yMax);
+ genkaon->SetPhiRange(phiMin, phiMax);
+ genkaon->Init(); // Generating pT and Y parametrsation in the desired kinematic range
+ AddGenerator(genkaon,"Kaon", ratiokaon);
+ fTotalRate+=ratiokaon;
+ }