Added set method.

diff --git a/JETAN/AliJetEventParticles.cxx b/JETAN/AliJetEventParticles.cxx
index 67e9325faf89a34ab89a64f48a2490137058a810..e82f12029ab08f0487ca1e029ad4cf58d2536cbe 100644 (file)
--- a/JETAN/AliJetEventParticles.cxx
+++ b/JETAN/AliJetEventParticles.cxx
@@ -17,6 +17,8 @@
 ClassImp(AliJetEventParticles)
 
 AliJetEventParticles::AliJetEventParticles(Int_t size) :
+  TObject(), 
+  fHeader(),
   fNParticles(0),
   fParticles(new TClonesArray("AliJetParticle",size)),
   fVertexX(0.),
@@ -46,6 +48,7 @@ AliJetEventParticles::AliJetEventParticles(Int_t size) :
 
 AliJetEventParticles::AliJetEventParticles(const AliJetEventParticles& source) :
   TObject(source), 
+  fHeader(),
   fNParticles(source.fNParticles),
   fParticles(new TClonesArray("AliJetParticle",source.fNParticles)),
   fVertexX(source.GetVertexX()),
@@ -84,6 +87,48 @@ AliJetEventParticles::AliJetEventParticles(const AliJetEventParticles& source) :
   source.Hard(1,fHard[0][1],fHard[1][1],fHard[2][1],fHard[3][1],fHard[4][1]);
 }
 
+void AliJetEventParticles::Set(const AliJetEventParticles& source)
+{
+  source.Copy(*this); 
+  fHeader=source.GetHeader();
+  fNParticles=source.fNParticles;
+  if(fParticles) delete fParticles;
+  fParticles=new TClonesArray("AliJetParticle",fNParticles);
+  fVertexX=source.GetVertexX();
+  fVertexY=source.GetVertexY();
+  fVertexZ=source.GetVertexZ();
+  fTrials=source.Trials();
+  fNJets=source.NTriggerJets();
+  fNUQJets=source.NUQTriggerJets();
+  fXJet=source.GetXJet();
+  fYJet=source.GetXJet();
+  fImpact=source.GetImpact();
+  fNHardScatters=source.GetNhard();
+  fNwNwColl=source.GetNpart();
+
+  for(Int_t i =0; i<fNParticles; i++)
+    {
+      const AliJetParticle *kjp=(const AliJetParticle *)source.fParticles->At(i);
+      new((*fParticles)[i]) AliJetParticle(*(kjp));
+    }
+  for (Int_t i = 0; i < 4; i++) fZquench[i] = 0.;
+  for (Int_t i = 0; i < 10; i++) 
+    for (Int_t j = 0; j < 4; j++) {
+      fJets[j][i]=0;    // Trigger jets
+      fUQJets[j][i]=0;  // Unquenched trigger jets
+    }
+  source.GetZQuench(fZquench);
+  for (Int_t i = 0; i < NTriggerJets(); i++){
+    source.TriggerJet(i,fJets[0][i],fJets[1][i],fJets[2][i],fJets[3][i]);
+  }
+  for (Int_t i = 0; i < NUQTriggerJets(); i++){
+    source.UQJet(i,fUQJets[0][i],fUQJets[1][i],fUQJets[2][i],fUQJets[3][i]);
+  }
+
+  source.Hard(0,fHard[0][0],fHard[1][0],fHard[2][0],fHard[3][0],fHard[4][0]);
+  source.Hard(1,fHard[0][1],fHard[1][1],fHard[2][1],fHard[3][1],fHard[4][1]);
+}
+
 AliJetEventParticles::~AliJetEventParticles()
 {
   //destructor   
@@ -93,6 +138,7 @@ AliJetEventParticles::~AliJetEventParticles()
 
 void  AliJetEventParticles::Reset(Int_t size)
 {
+  fHeader="";
   //deletes all particles from the event
   if(fParticles) fParticles->Clear();
   if(size>=0) fParticles->Expand(size);
@@ -148,7 +194,6 @@ void AliJetEventParticles::AddSignal(const AliJetEventParticles& source)
   fXJet=source.GetXJet();
   fYJet=source.GetXJet();
 
-
   for (Int_t i = 0; i < 4; i++) fZquench[i] = 0.;
   for (Int_t i = 0; i < 10; i++) 
     for (Int_t j = 0; j < 4; j++) {

diff --git a/JETAN/AliJetEventParticles.h b/JETAN/AliJetEventParticles.h
index 69c6692a8243a522940570f506803642f68b6e50..ab0d67b0b0210da80f4039fb4463565563d13ce7 100644 (file)
--- a/JETAN/AliJetEventParticles.h
+++ b/JETAN/AliJetEventParticles.h
@@ -30,6 +30,7 @@ class AliJetEventParticles: public TObject
   void SetHeader(TString& s){fHeader=s;}
   void Reset(Int_t size=-1); //deletes all entries
   void Clear(Option_t *option="") {TObject::Clear(option);Reset();}
+  void Set(const AliJetEventParticles& source);
   
   //adds particle to the event
   void AddParticle(AliJetParticle* p);  
@@ -62,7 +63,7 @@ class AliJetEventParticles: public TObject
   Double_t GetXJet() const {return fXJet;}
   Double_t GetYJet() const {return fYJet;}    
   void     GetZQuench(Double_t z[4]) const;
-  TString  getHeader() const {return fHeader;}
+  TString  GetHeader() const {return fHeader;}
 
   Float_t  GetImpact() const {return fImpact;}
   Int_t    GetNhard()  const {return fNHardScatters;}