Boost() moved to AliGenMC
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 Apr 2003 08:14:02 +0000 (08:14 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 Apr 2003 08:14:02 +0000 (08:14 +0000)
THijing/AliGenHijing.cxx
THijing/AliGenHijing.h

index 27a567d..70b924d 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.1  2003/03/15 14:45:57  morsch
+Classes imported from EVGEN
+
 Revision 1.47  2003/01/14 10:50:18  alibrary
 Cleanup of STEER coding conventions
 
@@ -216,8 +219,6 @@ AliGenHijing::AliGenHijing(Int_t npart)
 
     SetEnergyCMS();
     SetImpactParameterRange();
-    SetTarget();
-    SetProjectile();
     SetBoostLHC();
     SetJetEtaRange();
     SetJetPhiRange();
@@ -376,7 +377,9 @@ void AliGenHijing::Generate()
       if (fTrigger != kNoTrigger) {
          if (!CheckTrigger()) continue;
       }
-      if (fLHC) Boost();
+      Double_t dy    = - 0.5 * TMath::Log(Double_t(fZProjectile) * Double_t(fATarget) / 
+                                         (Double_t(fZTarget)    * Double_t(fAProjectile)));
+      if (fLHC) Boost(dy);
       
       
       Int_t np = fParticles->GetEntriesFast();
@@ -649,37 +652,6 @@ Bool_t AliGenHijing::Stable(TParticle*  particle)
 }
 
 
-void AliGenHijing::Boost()
-{
-//
-// Boost cms into LHC lab frame
-//
-    Double_t dy    = - 0.5 * TMath::Log(Double_t(fZProjectile) * Double_t(fATarget) / 
-                                     (Double_t(fZTarget)    * Double_t(fAProjectile)));
-    Double_t beta  = TMath::TanH(dy);
-    Double_t gamma = 1./TMath::Sqrt(1.-beta*beta);
-    Double_t gb    = gamma * beta;
-
-    printf("\n Boosting particles to lab frame %f %f %f", dy, beta, gamma);
-    
-    Int_t i;
-    Int_t np = fParticles->GetEntriesFast();
-    for (i = 0; i < np; i++) 
-    {
-       TParticle* iparticle = (TParticle*) fParticles->At(i);
-
-       Double_t e   = iparticle->Energy();
-       Double_t px  = iparticle->Px();
-       Double_t py  = iparticle->Py();
-       Double_t pz  = iparticle->Pz();
-
-       Double_t eb  = gamma * e -      gb * pz;
-       Double_t pzb =   -gb * e +   gamma * pz;
-
-       iparticle->SetMomentum(px, py, pzb, eb);
-    }
-}
-
 
 void AliGenHijing::MakeHeader()
 {
index e2b5db7..6fee35a 100644 (file)
@@ -34,10 +34,7 @@ class AliGenHijing : public AliGenMC
     virtual void    SetEnergyCMS(Float_t energy=5500) {fEnergyCMS=energy;}
     virtual void    SetReferenceFrame(TString frame="CMS")
        {fFrame=frame;}
-    virtual void    SetProjectile(TString proj="A", Int_t a=208, Int_t z=82)
-       {fProjectile = proj; fAProjectile = a; fZProjectile = z;}    
-    virtual void    SetTarget(TString tar="A", Int_t a=208, Int_t z=82)
-       {fTarget = tar; fATarget = a; fZTarget = z;}    
+
     virtual void    SetImpactParameterRange(Float_t bmin = 0, Float_t bmax = 15.)
        {fMinImpactParam=bmin; fMaxImpactParam=bmax;}
     virtual void    KeepFullEvent();
@@ -85,7 +82,6 @@ class AliGenHijing : public AliGenMC
     virtual Bool_t ProvidesCollisionGeometry() const {return kTRUE;}
     virtual AliCollisionGeometry* CollisionGeometry() const {return fCollisionGeometry;}
     virtual void EvaluateCrossSections();
-    virtual void Boost();
     virtual TGraph* CrossSection()     {return fDsigmaDb;}
     virtual TGraph* BinaryCollisions() {return fDnDb;}
     virtual Bool_t  CheckTrigger();
@@ -97,12 +93,7 @@ class AliGenHijing : public AliGenMC
 
  protected:
     TString     fFrame;         // Reference frame 
-    TString     fProjectile;    // Projectile
-    TString     fTarget;        // Target
-    Int_t       fAProjectile;    // Projectile A
-    Int_t       fZProjectile;    // Projectile Z
-    Int_t       fATarget;        // Target A
-    Int_t       fZTarget;        // Target Z
+
     Float_t     fMinImpactParam; // minimum impact parameter
     Float_t     fMaxImpactParam; // maximum impact parameter   
     Int_t       fKeep;           // Flag to keep full event information
@@ -138,8 +129,7 @@ class AliGenHijing : public AliGenMC
     Int_t      fSpecn;          // Num. of spectator neutrons
     Int_t      fSpecp;          // Num. of spectator protons
     Int_t       fLHC;            // Assume LHC as lab frame
-    TClonesArray* fParticles;    // Particle List
-    
+
  private:
     // adjust the weight from kinematic cuts
     void   AdjustWeights();