Update to store the accumulated charge of two time windows for PID

diff --git a/TRD/AliTRDmcmSim.cxx b/TRD/AliTRDmcmSim.cxx
index 9f070f703fffb04213717455c1ee009c5eaf1606..1c8ca81a3ca525536311920a19b996205167723d 100644 (file)
--- a/TRD/AliTRDmcmSim.cxx
+++ b/TRD/AliTRDmcmSim.cxx
@@ -1595,6 +1595,8 @@ void AliTRDmcmSim::FitTracklet()
         }
         new ((*fTrackletArray)[fTrackletArray->GetEntriesFast()]) AliTRDtrackletMCM((UInt_t) fMCMT[cpu], fDetector*2 + fRobPos%2, fRobPos, fMcmPos);
         ((AliTRDtrackletMCM*) (*fTrackletArray)[fTrackletArray->GetEntriesFast()-1])->SetLabel(mcLabel);
+        ((AliTRDtrackletMCM*) (*fTrackletArray)[fTrackletArray->GetEntriesFast()-1])->SetQ0(q0);
+        ((AliTRDtrackletMCM*) (*fTrackletArray)[fTrackletArray->GetEntriesFast()-1])->SetQ1(q1);
       }
     }
   }

diff --git a/TRD/AliTRDtrackletMCM.cxx b/TRD/AliTRDtrackletMCM.cxx
index 907f19a25a5543077f4bb2dd4f9ede4412ad2e8c..e9597623a8a715f3ed2337a503f349eea3f112f3 100644 (file)
--- a/TRD/AliTRDtrackletMCM.cxx
+++ b/TRD/AliTRDtrackletMCM.cxx
@@ -35,6 +35,8 @@ AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord) :
   fTrackletWord(trackletWord), 
   fMCM(-1), 
   fROB(-1), 
+  fQ0(0),
+  fQ1(0),
   fLabel(-1)
 { 
     fGeo = new AliTRDgeometry();
@@ -47,6 +49,8 @@ AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord, Int_t hcid) :
   fTrackletWord(trackletWord), 
   fMCM(-1),
   fROB(-1),
+  fQ0(0),
+  fQ1(0),
   fLabel(-1)
 { 
     fGeo = new AliTRDgeometry();
@@ -59,6 +63,8 @@ AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord, Int_t hcid, Int_t rob,
   fTrackletWord(trackletWord), 
   fMCM(mcm),
   fROB(rob),
+  fQ0(0),
+  fQ1(0),
   fLabel(-1)
 { 
     fGeo = new AliTRDgeometry();
@@ -71,6 +77,8 @@ AliTRDtrackletMCM::AliTRDtrackletMCM(const AliTRDtrackletMCM &rhs) :
   fTrackletWord(rhs.fTrackletWord),
   fMCM(rhs.fMCM),
   fROB(rhs.fROB),
+  fQ0(rhs.fQ0),
+  fQ1(rhs.fQ1),
   fLabel(rhs.fLabel)
 {
     fGeo = new AliTRDgeometry();

diff --git a/TRD/AliTRDtrackletMCM.h b/TRD/AliTRDtrackletMCM.h
index aaf126fbff40f8f445afdc758548638652fc0efb..e9bb0d603bc93ece195bc7c37c07c92a2c0c8882 100644 (file)
--- a/TRD/AliTRDtrackletMCM.h
+++ b/TRD/AliTRDtrackletMCM.h
@@ -45,6 +45,9 @@ class AliTRDtrackletMCM : public AliTRDtrackletBase {
   Float_t GetY() const { return (GetYbin() * 160e-4); }
   Float_t GetZ() const { return fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId / 12) % 5)->GetRowPos( 4 * (fROB / 2) + fMCM / 4); }
 
+  Int_t GetQ0() const { return fQ0; }
+  Int_t GetQ1() const { return fQ1; }
+
   UInt_t GetTrackletWord() const { return fTrackletWord; }
   void SetTrackletWord(UInt_t trackletWord) { fTrackletWord = trackletWord; }
 
@@ -53,6 +56,8 @@ class AliTRDtrackletMCM : public AliTRDtrackletBase {
   void SetMCM(Int_t mcm) { fMCM = mcm; }
   void SetROB(Int_t rob) { fROB = rob; }
   void SetLabel(Int_t label) { fLabel = label; }
+  void SetQ0(Int_t charge) { fQ0 = charge; }
+  void SetQ1(Int_t charge) { fQ1 = charge; }
 
  protected:
   AliTRDgeometry *fGeo; //! TRD geometry
@@ -63,12 +68,15 @@ class AliTRDtrackletMCM : public AliTRDtrackletBase {
   Int_t fMCM; // MCM no. in which the tracklet was found
   Int_t fROB; // ROB no. on which the tracklet was found
 
+  Int_t fQ0; // accumulated charge in the first time window
+  Int_t fQ1; // accumulated charge in the second time window
+
   Int_t fLabel; // label for MC track
 
  private:
   AliTRDtrackletMCM& operator=(const AliTRDtrackletMCM &rhs);   // not implemented
 
-  ClassDef(AliTRDtrackletMCM, 1);
+  ClassDef(AliTRDtrackletMCM, 2);
 };
 
 #endif