fMomentum(0.),
fBeta(0.),
fPmax (0.),
- fATarget (0.),
- fZTarget (0.),
fCharge(0),
fProtonDirection(1.),
fTemperatureG(0.),
fMomentum(0.),
fBeta(0.),
fPmax (10.),
- fATarget (208.),
- fZTarget (82.),
fCharge(1),
fProtonDirection(1.),
fTemperatureG(0.04),
// Initialization
//
Float_t kMass = TDatabasePDG::Instance()->GetParticle(kProton)->Mass();
- fMomentum = fCMS/2. * fZTarget / fATarget;
+ fMomentum = fCMS/2. * Float_t(fZTarget) / Float_t(fATarget);
fBeta = fMomentum / TMath::Sqrt(kMass * kMass + fMomentum * fMomentum);
if (fDebug) {
fDebugHist1 = new TH2F("DebugHist1", "nu vs N_slow", 100, 0., 100., 20, 0., 20.);
virtual void Generate();
virtual void SetPmax(Float_t pmax = 10.) {fPmax = pmax;}
virtual void SetNominalCmsEnergy(Float_t energy = 14000.) {fCMS = energy;}
- virtual void SetTarget(Float_t a=208, Float_t z=82) {fATarget = a; fZTarget = z;}
+ virtual void SetTarget(Int_t a = 208, Int_t z = 82) {fATarget = a; fZTarget = z;}
virtual void SetTarget(TString s, Int_t a, Int_t z) {AliGenerator::SetTarget(s, a, z);}
virtual void SetProtonDirection(Float_t dir = 1.);
virtual void SetCharge(Int_t c = 1) {fCharge = c;}
Float_t fMomentum; // Target nucleus momentum
Float_t fBeta; // Target nucleus beta
Float_t fPmax; // Maximum slow nucleon momentum
- Float_t fATarget; // Target nucleus mass number
- Float_t fZTarget; // Target nucleus charge number
Int_t fCharge; // Slow nucleon charge
Float_t fProtonDirection; // Direction of the proton
Float_t fTemperatureG; // Source Temperature for gray nucleons
void AliOmegaDalitz::
Rot(Double_t pin[3], Double_t pout[3], Double_t costheta, Double_t sintheta,
- Double_t cosphi, Double_t sinphi)
+ Double_t cosphi, Double_t sinphi) const
{
// Perform rotation
pout[0] = pin[0]*costheta*cosphi-pin[1]*sinphi+pin[2]*sintheta*cosphi;
virtual void Decay(TClonesArray* array);
virtual const TLorentzVector* Products() const {return fProducts;}
virtual void Copy(TObject&) const;
+
+ protected:
+ TH1F* fEPMass; // Histogram for electron pair mass
+ TH1F* fMPMass; // Histogram for muon pair mass
+ TLorentzVector fProducts[3]; // Decay products
+
private:
virtual void Rot(Double_t pin[3], Double_t pout[3],
Double_t costheta, Double_t sintheta,
- Double_t cosphi, Double_t sinphi);
+ Double_t cosphi, Double_t sinphi) const;
AliOmegaDalitz(const AliOmegaDalitz &dalitz);
AliOmegaDalitz & operator=(const AliOmegaDalitz & rhs);
- protected:
- TH1F* fEPMass; // Histogram for electron pair mass
- TH1F* fMPMass; // Histogram for muon pair mass
- TLorentzVector fProducts[3]; // Decay products
ClassDef(AliOmegaDalitz, 1) // AliDecayer implementation for omega Dalitz
};