Preparation for four momentum background subtraction, needs area vector

diff --git a/STEER/AliAODJet.cxx b/STEER/AliAODJet.cxx
index 83fc5f57d1c0ee8bbd411977839718d7954b1a74..55d13a0e3fa559d79da1c8fa6dd49a620b67c35f 100644 (file)
--- a/STEER/AliAODJet.cxx
+++ b/STEER/AliAODJet.cxx
@@ -32,6 +32,7 @@ AliAODJet::AliAODJet() :
     fNeutralFraction(0),
     fTrigger(0),
     fMomentum(0),
+    fVectorAreaCharged(0),
     fRefTracks(new TRefArray())
 {
   // constructor
@@ -47,6 +48,7 @@ AliAODJet::AliAODJet(Double_t px, Double_t py, Double_t pz, Double_t e):
     fNeutralFraction(0),
     fTrigger(0),
     fMomentum(0),
+    fVectorAreaCharged(0),
     fRefTracks(new TRefArray())
 {
   // constructor
@@ -63,6 +65,7 @@ AliAODJet::AliAODJet(TLorentzVector & p):
     fNeutralFraction(0),
     fTrigger(0),
     fMomentum(0),
+    fVectorAreaCharged(0),
     fRefTracks(new TRefArray())
 {
   // constructor
@@ -80,6 +83,7 @@ AliAODJet::~AliAODJet()
 {
   // destructor
     delete fMomentum;
+    delete fVectorAreaCharged;
     delete fRefTracks;
 }
 
@@ -89,6 +93,7 @@ AliAODJet::AliAODJet(const AliAODJet& jet) :
     fNeutralFraction(jet.fNeutralFraction),
     fTrigger(jet.fTrigger),
     fMomentum(0),
+    fVectorAreaCharged(0),
     fRefTracks(0)
 {
   // Copy constructor
@@ -99,6 +104,7 @@ AliAODJet::AliAODJet(const AliAODJet& jet) :
     fEffectiveAreaError[0] = jet.fEffectiveAreaError[0];
     fEffectiveAreaError[1] = jet.fEffectiveAreaError[1];
     fMomentum  = new TLorentzVector(*jet.fMomentum);
+    fVectorAreaCharged = new TLorentzVector(*jet.fVectorAreaCharged);
     fRefTracks = new TRefArray(*jet.fRefTracks);
 }
 
@@ -116,8 +122,8 @@ AliAODJet& AliAODJet::operator=(const AliAODJet& jet)
     fEffectiveAreaError[1] = jet.fEffectiveAreaError[1];
     fNeutralFraction = jet.fNeutralFraction;
     fTrigger = jet.fTrigger;
-    delete fMomentum;
-    fMomentum  = new TLorentzVector(*jet.fMomentum);
+    *fMomentum  = *jet.fMomentum;
+    *fVectorAreaCharged = *jet.fVectorAreaCharged;
     delete fRefTracks;
     fRefTracks = new TRefArray(*jet.fRefTracks);    
   }

diff --git a/STEER/AliAODJet.h b/STEER/AliAODJet.h
index 696d411d4f9d0c44fd85c0b796c6d2b30b1477ec..3534d5dafe511890468a8ca52ee7a5e8b1a2ab35 100644 (file)
--- a/STEER/AliAODJet.h
+++ b/STEER/AliAODJet.h
@@ -69,6 +69,14 @@ class AliAODJet : public AliVParticle {
 
     virtual Double_t   EffectiveAreaCharged()   const { return  fEffectiveArea[0];}
     virtual Double_t   EffectiveAreaNeutral()   const { return  fEffectiveArea[1];}
+    virtual void SetVectorAreaCharged(TLorentzVector *effVACh){
+      if(!fVectorAreaCharged)fVectorAreaCharged= new TLorentzVector(*effVACh);
+      else *fVectorAreaCharged = *effVACh;
+    }
+    virtual TLorentzVector*  VectorAreaCharged()   const {return fVectorAreaCharged;}
+
+
+
     virtual Double_t   ErrorEffectiveAreaCharged()   const { return  fEffectiveAreaError[0];}
     virtual Double_t   ErrorEffectiveAreaNeutral()   const { return  fEffectiveAreaError[1];}
     virtual Double_t   DeltaR(const AliVParticle* part);
@@ -99,9 +107,10 @@ class AliAODJet : public AliVParticle {
     Double32_t      fNeutralFraction;        //[0,1,12] Neutral fraction between 0 and 1 12 bit precision;
     UChar_t         fTrigger;                // Bit mask to flag jets triggered by a certain detector  
     TLorentzVector* fMomentum;               // Jet 4-momentum vector
+    TLorentzVector* fVectorAreaCharged;      // jet area four momentum 
     TRefArray*      fRefTracks;              // array of references to the tracks belonging to the jet
 
-    ClassDef(AliAODJet,6);
+    ClassDef(AliAODJet,7);
 
 };