bookkeep pt hard info when generating on-the-fly pythia events

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / AliEmcalJet.cxx

diff --git a/PWGJE/EMCALJetTasks/AliEmcalJet.cxx b/PWGJE/EMCALJetTasks/AliEmcalJet.cxx
index 450cd4b97617358ab87988b20b2cd4654fe2fa32..8ae3f3ca0e961da819c773bbb98cbbb157840141 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliEmcalJet.cxx
+++ b/PWGJE/EMCALJetTasks/AliEmcalJet.cxx
@@ -24,6 +24,7 @@ AliEmcalJet::AliEmcalJet() :
   fAreaPhi(0),       
   fAreaEmc(-1), 
   fAxisInEmcal(0), 
+  fFlavourTagging(0),
   fMaxCPt(0), 
   fMaxNPt(0), 
   fMCPt(0),
@@ -35,6 +36,8 @@ AliEmcalJet::AliEmcalJet() :
   fTrackIDs(),
   fMatched(2),
   fMatchingType(0),
+  fTaggedJet(0x0),
+  fTagStatus(-1),
   fPtSub(0),
   fPtVectSub(0),
   fTriggers(0)
@@ -60,6 +63,7 @@ AliEmcalJet::AliEmcalJet(Double_t px, Double_t py, Double_t pz) :
   fAreaPhi(0),       
   fAreaEmc(-1), 
   fAxisInEmcal(0),
+  fFlavourTagging(0),
   fMaxCPt(0), 
   fMaxNPt(0), 
   fMCPt(0),
@@ -71,6 +75,8 @@ AliEmcalJet::AliEmcalJet(Double_t px, Double_t py, Double_t pz) :
   fTrackIDs(),
   fMatched(2),
   fMatchingType(0),
+  fTaggedJet(0x0),
+  fTagStatus(-1),
   fPtSub(0),
   fPtVectSub(0),
   fTriggers(0)
@@ -102,6 +108,7 @@ AliEmcalJet::AliEmcalJet(Double_t pt, Double_t eta, Double_t phi, Double_t m) :
   fAreaPhi(0),       
   fAreaEmc(-1), 
   fAxisInEmcal(0),
+  fFlavourTagging(0),
   fMaxCPt(0), 
   fMaxNPt(0),
   fMCPt(0),
@@ -113,6 +120,8 @@ AliEmcalJet::AliEmcalJet(Double_t pt, Double_t eta, Double_t phi, Double_t m) :
   fTrackIDs(),
   fMatched(2),
   fMatchingType(0),
+  fTaggedJet(0x0),
+  fTagStatus(-1),
   fPtSub(0),
   fPtVectSub(0),
   fTriggers(0)
@@ -141,6 +150,7 @@ AliEmcalJet::AliEmcalJet(const AliEmcalJet &jet) :
   fAreaPhi(jet.fAreaPhi),       
   fAreaEmc(jet.fAreaEmc), 
   fAxisInEmcal(jet.fAxisInEmcal),
+  fFlavourTagging(jet.fFlavourTagging),
   fMaxCPt(jet.fMaxCPt), 
   fMaxNPt(jet.fMaxNPt), 
   fMCPt(jet.fMCPt),
@@ -152,6 +162,8 @@ AliEmcalJet::AliEmcalJet(const AliEmcalJet &jet) :
   fTrackIDs(jet.fTrackIDs),
   fMatched(jet.fMatched),
   fMatchingType(jet.fMatchingType),
+  fTaggedJet(jet.fTaggedJet),
+  fTagStatus(jet.fTagStatus),
   fPtSub(jet.fPtSub),
   fPtVectSub(jet.fPtVectSub),
   fTriggers(jet.fTriggers)
@@ -181,6 +193,7 @@ AliEmcalJet &AliEmcalJet::operator=(const AliEmcalJet &jet)
     fAreaPhi            = jet.fAreaPhi; 
     fAreaEmc            = jet.fAreaEmc; 
     fAxisInEmcal        = jet.fAxisInEmcal; 
+    fFlavourTagging     = jet.fFlavourTagging;
     fMaxCPt             = jet.fMaxCPt; 
     fMaxNPt             = jet.fMaxNPt;
     fMCPt               = jet.fMCPt;
@@ -195,6 +208,8 @@ AliEmcalJet &AliEmcalJet::operator=(const AliEmcalJet &jet)
     fClosestJetsDist[0] = jet.fClosestJetsDist[0];  
     fClosestJetsDist[1] = jet.fClosestJetsDist[1]; 
     fMatched            = jet.fMatched;
+    fTaggedJet          = jet.fTaggedJet;
+    fTagStatus          = jet.fTagStatus;
     fPtSub              = jet.fPtSub;
     fPtVectSub          = jet.fPtVectSub;
     fTriggers           = jet.fTriggers;
@@ -260,7 +275,8 @@ AliVParticle* AliEmcalJet::GetLeadingTrack(TClonesArray *tracks) const
   for (Int_t i = 0; i < GetNumberOfTracks(); i++) {
     AliVParticle *track = TrackAt(i, tracks);
     if (!track) {
-      AliError(Form("Unable to find jet track %d in collection %s (pos in collection %d, max %d)",i,tracks->GetName(),TrackAt(i),tracks->GetEntriesFast()));
+      AliError(Form("Unable to find jet track %d in collection %s (pos in collection %d, max %d)",
+                   i,tracks->GetName(),TrackAt(i),tracks->GetEntriesFast()));
       continue;
     }
     if (!maxTrack || track->Pt() > maxTrack->Pt()) 
@@ -277,7 +293,8 @@ AliVCluster* AliEmcalJet::GetLeadingCluster(TClonesArray *clusters) const
   for (Int_t i = 0; i < GetNumberOfClusters(); i++) {
     AliVCluster *cluster = ClusterAt(i, clusters);
     if (!cluster) {
-      AliError(Form("Unable to find jet cluster %d in collection %s (pos in collection %d, max %d)",i,clusters->GetName(),ClusterAt(i),clusters->GetEntriesFast()));
+      AliError(Form("Unable to find jet cluster %d in collection %s (pos in collection %d, max %d)",
+                   i,clusters->GetName(),ClusterAt(i),clusters->GetEntriesFast()));
       continue;
     }
     if (!maxCluster || cluster->E() > maxCluster->E())