fMomentumDist(0),
fUnfoldedMomentumDist(0),
fZenithDist(0),
- fPDist(0)
+ fPDist(0),
+ fNParticles(0)
{
//
// Default ctor.
fMomentumDist(0),
fUnfoldedMomentumDist(0),
fZenithDist(0),
- fPDist(0)
+ fPDist(0),
+ fNParticles(0)
{
//
// Standard ctor.
fOrigin[2] = 0.;
}
-//_____________________________________________________________________________
-AliGenACORDE::AliGenACORDE(const AliGenACORDE& gen)
- : AliGenerator(gen)
-{
- //
- // Copy constructor
- //
- gen.Copy(*this);
-}
-
-//_____________________________________________________________________________
-AliGenACORDE& AliGenACORDE::operator=(const AliGenACORDE& gen)
-{
- //
- // Asingment operator
- //
- gen.Copy(*this);
- return *this;
-}
-
//_____________________________________________________________________________
AliGenACORDE::~AliGenACORDE()
{
// Call the respective method inside the loop for the number
// of tracks per trigger.
- for (Int_t i = 0; i < fNpart; i++ ) {
+ for (Int_t i = 0; i < fNParticles; i++ ) {
if ( fCRMode == kMuonBundle ) {
this->GenerateOneMuonBundle();
// Initialize some internal methods.
//
+
+ printf("**************************************************************\n");
+ printf("<<< *** Starting the AliGenACORDE cosmic generator ******** >>>\n");
+ printf("<<< *** No. of muons generated at the surface of P2: %d, * >>>\n",fNParticles);
+ printf("**************************************************************\n");
+
// Determine some specific data members.
fPRange = TMath::Abs(fPMax-fPMin);
char buffer6[1024];
char buffer7[1024];
- sprintf(buffer1, h1Coef, y, y, y, y, y, y);
- sprintf(buffer2, h2Coef, y, y, y, y, y, y);
- sprintf(buffer3, h3Coef, y, y, y, y, y, y);
- sprintf(buffer4, s2Coef, y, y, y, y, y, y);
+ //sprintf(buffer1, h1Coef, y, y, y, y, y, y);
+ //sprintf(buffer2, h2Coef, y, y, y, y, y, y);
+ //sprintf(buffer3, h3Coef, y, y, y, y, y, y);
+ //sprintf(buffer4, s2Coef, y, y, y, y, y, y);
- sprintf(buffer5, h, buffer1, buffer2, buffer3, buffer4);
+ //sprintf(buffer5, h, buffer1, buffer2, buffer3, buffer4);
- sprintf(buffer6, flux, buffer5);
+ //sprintf(buffer6, flux, buffer5);
fMomentumDist = new TF1("fMomentumDist", buffer6, fPMin, fPMax);
- sprintf(buffer7, normalizedFlux, buffer5);
+ //sprintf(buffer7, normalizedFlux, buffer5);
fUnfoldedMomentumDist = new TF1("fUnfoldedMomentumDist", buffer7, fPMin, fPMax);
for (Int_t i = 0; i < 4; i++ ) {
fMomentumDist->SetParName(i, paramNames[i]);
Float_t weight = 0;
for ( Int_t i = 0; i < pEnd; i++ ) {
// Fill the distribution
- sprintf(name, "zenith%d", i+1);
- sprintf(title, "Zenith distribution, p=%f", fPMin+(Float_t)i);
+ //sprintf(name, "zenith%d", i+1);
+ //sprintf(title, "Zenith distribution, p=%f", fPMin+(Float_t)i);
zenith = new(mom[i]) TH1F(name, title, TMath::Abs(TMath::Nint(fZenithMax-fZenithMin)), TMath::Cos(fZenithMax*TMath::Pi()/180), TMath::Cos(fZenithMin*TMath::Pi()/180));
// Make a loop for the angle and fill the histogram for the weight