Changes for #94112: Standard AOD size (Jens)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Jun 2012 12:38:10 +0000 (12:38 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Jun 2012 12:38:10 +0000 (12:38 +0000)
ANALYSIS/AliAnalysisTaskESDfilter.cxx
STEER/AOD/AliAODPid.cxx
STEER/AOD/AliAODPid.h
STEER/AOD/AliAODTrack.cxx
STEER/AOD/AliAODTrack.h
STEER/AODLinkDef.h
STEER/STEERBase/AliTPCdEdxInfo.h

index 2ff038d..4987378 100644 (file)
@@ -2284,7 +2284,7 @@ void AliAnalysisTaskESDfilter::SetDetectorRawSignals(AliAODPid *aodpid, AliESDtr
   }
   aodpid->SetTOFpidResolution(tofRes);
 
-  aodpid->SetHMPIDsignal(track->GetHMPIDsignal());
+//  aodpid->SetHMPIDsignal(track->GetHMPIDsignal());
 
 }
 
index 5840380..0959b8c 100644 (file)
@@ -39,12 +39,10 @@ AliAODPid::AliAODPid():
     fTRDntls(0),
     fTRDslices(0x0),
     fTOFesdsignal(0),
-    fHMPIDsignal(0),
     fTPCdEdxInfo(0)
 {
   // default constructor
     for(Int_t i=0; i<kSPECIES; i++) fIntTime[i]   = 0; 
-    for(Int_t i=0; i<5; i++) fHMPIDprobs[i] = 0.;
     for(Int_t i=0; i<5; i++) fTOFpidResolution[i] = 0.;
     for(Int_t i=0; i<6; i++) {
       fTRDmomentum[i]      = 0.;
@@ -74,13 +72,11 @@ AliAODPid::AliAODPid(const AliAODPid& pid) :
   fTRDntls(pid.fTRDntls),
   fTRDslices(0x0),
   fTOFesdsignal(pid.fTOFesdsignal),
-  fHMPIDsignal(pid.fHMPIDsignal),
   fTPCdEdxInfo(0x0)
 {
   // Copy constructor
   SetTRDsignal(fTRDnSlices, pid.fTRDslices);
     for(Int_t i=0; i<kSPECIES; i++) fIntTime[i]=pid.fIntTime[i];
-    for(Int_t i=0; i<5; i++) fHMPIDprobs[i] = pid.fHMPIDprobs[i];
 
     for(Int_t i=0; i<6; i++){ 
       fTRDmomentum[i]=pid.fTRDmomentum[i];
@@ -130,10 +126,6 @@ AliAODPid& AliAODPid::operator=(const AliAODPid& pid)
     for (Int_t i = 0; i < 5; i++) fTOFpidResolution[i]=pid.fTOFpidResolution[i];
     for (Int_t i = 0; i < 5; i++) fIntTime[i]=pid.fIntTime[i];
     
-    fHMPIDsignal=pid.fHMPIDsignal;
-    
-    for(Int_t i = 0; i < 5; i++) fHMPIDprobs[i] = pid.fHMPIDprobs[i];
-    
      SetTPCdEdxInfo(pid.fTPCdEdxInfo);
   }
 
@@ -164,22 +156,6 @@ void AliAODPid::GetTOFpidResolution(Double_t tofRes[5]) const
 }
 
 //______________________________________________________________________________
-void AliAODPid::SetHMPIDprobs(Double_t hmpPid[5]) 
-{
-  //
-  // Set the HMPID PID probablities that are read from ESD
-  //  
-  for(Int_t i = 0; i < 5; i++ ) fHMPIDprobs[i] =  hmpPid[i];
-}
-//______________________________________________________________________________
-void AliAODPid::GetHMPIDprobs(Double_t *p) const
-{
-  //
-  // Set the HMPID PID probablities that are read from ESD
-  //  
-  for(Int_t i = 0; i < AliPID::kSPECIES; i++ ) p[i] =  fHMPIDprobs[i];
-}
-//______________________________________________________________________________
 void AliAODPid::SetITSdEdxSamples(const Double_t s[4])
 {
   //
@@ -210,4 +186,3 @@ void AliAODPid::SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo)
   (*fTPCdEdxInfo)=(*dEdxInfo);
 }
 
-
index 7fdf7b0..2f54c6f 100644 (file)
@@ -28,9 +28,9 @@ class AliAODPid : public TObject {
   void      SetITSsignal(Double_t its)                         {fITSsignal=its;}
   void      SetITSdEdxSamples(const Double_t s[4]);
   void      SetTPCsignal(Double_t tpc)                         {fTPCsignal=tpc;}
-  void      SetTPCsignalN(UShort_t tpcN)                       {fTPCsignalN=tpcN;}
+  void      SetTPCsignalN(UShort_t tpcN)                       {fTPCsignalN=(UChar_t)((tpcN<160)?tpcN:160);}
   void      SetTPCmomentum(Double_t tpcMom)                    {fTPCmomentum=tpcMom;}
-  void      SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo); 
+  void      SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo);
   inline void  SetTRDsignal(Int_t nslices, const Double_t * const trdslices);  
   void      SetTRDmomentum(Int_t nplane, Float_t trdMom)       {fTRDmomentum[nplane]=trdMom;}
   inline void  SetTRDncls(UChar_t ncls, Int_t layer = -1);
@@ -38,8 +38,6 @@ class AliAODPid : public TObject {
   void      SetTOFsignal(Double_t tof)                         {fTOFesdsignal=tof;}
   void      SetTOFpidResolution(Double_t tofPIDres[5]);
   void      SetIntegratedTimes(Double_t timeint[5]);
-  void      SetHMPIDsignal(Double_t hmpid)                     {fHMPIDsignal=hmpid;}
-  void      SetHMPIDprobs(Double_t hmpPid[5]);
 
   Double_t  GetITSsignal()       const {return  fITSsignal;}
   void      GetITSdEdxSamples(Double_t *s) const;
@@ -48,19 +46,17 @@ class AliAODPid : public TObject {
     else return 0.;
   }
   Double_t  GetTPCsignal()       const {return  fTPCsignal;}
-  UShort_t  GetTPCsignalN()      const {return  fTPCsignalN;}
+  UShort_t  GetTPCsignalN()      const {return  (UShort_t)fTPCsignalN;}
   AliTPCdEdxInfo * GetTPCdEdxInfo()const{return fTPCdEdxInfo;}
 
   Double_t  GetTPCmomentum()     const {return  fTPCmomentum;}
   Int_t     GetTRDnSlices()      const {return  fTRDnSlices;}
   Double_t* GetTRDsignal()       const {return  fTRDslices;}
-  const Float_t*  GetTRDmomentum() const {return  fTRDmomentum;}
+  const Double_t*  GetTRDmomentum() const {return  fTRDmomentum;}
   UChar_t   GetTRDncls(UChar_t layer) const { if(layer > 5) return 0; return fTRDncls[layer];}
   inline UChar_t GetTRDncls() const;
   UChar_t   GetTRDntrackletsPID() const {return fTRDntls;}
   Double_t  GetTOFsignal()       const {return  fTOFesdsignal;}
-  Double_t  GetHMPIDsignal()     const {return  fHMPIDsignal;}
-  void      GetHMPIDprobs(Double_t *p) const;
 
   void      GetIntegratedTimes(Double_t timeint[5])  const; 
   void      GetTOFpidResolution (Double_t tofRes[5]) const;
@@ -68,23 +64,24 @@ class AliAODPid : public TObject {
  private :
   Double32_t  fITSsignal;        //[0.,0.,10] detector raw signal
   Double32_t  fITSdEdxSamples[4];//[0.,0.,10] ITS dE/dx samples
+
   Double32_t  fTPCsignal;        //[0.,0.,10] detector raw signal
-  UShort_t    fTPCsignalN;       // number of points used for TPC dE/dx
-  Double_t    fTPCmomentum;      // momentum at the inner wall of TPC;
+  UChar_t     fTPCsignalN;       // number of points used for TPC dE/dx
+  Double32_t  fTPCmomentum;      //[0.,0.,20] momentum at the inner wall of TPC;
+
   Int_t       fTRDnSlices;       // N slices used for PID in the TRD
   UChar_t     fTRDntls;          // number of tracklets used for PID calculation
   UChar_t     fTRDncls[6];       // number of clusters used for dE/dx calculation
-  Double32_t* fTRDslices;        //[fTRDnSlices]
-  Float_t     fTRDmomentum[6];   // momentum at the TRD layers
-  Double32_t  fTOFesdsignal;     // TOF signal - t0 (T0 interaction time)
-  Double32_t  fTOFpidResolution[5]; // TOF pid resolution for each mass hypotesys 
-  Double32_t  fIntTime[5];       // track time hypothesis
-  Double32_t  fHMPIDsignal;      // detector raw signal
-  Double32_t  fHMPIDprobs[5];    // detector pid probabilities
+  Double32_t* fTRDslices;        //[fTRDnSlices][0.,0.,10]
+  Double32_t  fTRDmomentum[6];   //[0.,0.,10]  momentum at the TRD layers
+
+  Double32_t  fTOFesdsignal;     //[0.,0.,20] TOF signal - t0 (T0 interaction time)
+  Double32_t  fTOFpidResolution[5]; //[0.,0.,20] TOF pid resolution for each mass hypotesys 
+  Double32_t  fIntTime[5];       //[0.,0.,20] track time hypothesis
  
   AliTPCdEdxInfo * fTPCdEdxInfo; // object containing dE/dx information for different pad regions
 
-  ClassDef(AliAODPid, 10);
+  ClassDef(AliAODPid, 11);
 };
 
 //_____________________________________________________________
index 854fd14..0362a98 100644 (file)
@@ -715,7 +715,7 @@ Double_t AliAODTrack::GetTRDmomentum(Int_t plane, Double_t */*sp*/) const
   // in TRD layer "plane".
 
   if (!fDetPid) return -1;
-  const Float_t *trdMomentum=fDetPid->GetTRDmomentum();
+  const Double_t *trdMomentum=fDetPid->GetTRDmomentum();
 
   if (!trdMomentum) {
     return -1.;
index ec0fce5..06f6d74 100644 (file)
@@ -275,14 +275,14 @@ class AliAODTrack : public AliVTrack {
   virtual AliTPCdEdxInfo* GetTPCdEdxInfo() const {return fDetPid?fDetPid->GetTPCdEdxInfo():0;}
   Double_t  GetTPCmomentum()     const { return fDetPid?fDetPid->GetTPCmomentum():0.;  }
   Double_t  GetTOFsignal()       const { return fDetPid?fDetPid->GetTOFsignal():0.;    }
-  Double_t  GetHMPIDsignal()     const { return fDetPid?fDetPid->GetHMPIDsignal():0.;  }
+  Double_t  GetHMPIDsignal()     const { return 0.;  } // TODO: To be implemented properly with the new HMPID object
   
   void      GetIntegratedTimes(Double_t *times) const {if (fDetPid) fDetPid->GetIntegratedTimes(times); }
   Double_t  GetTRDslice(Int_t plane, Int_t slice) const;
   Double_t  GetTRDmomentum(Int_t plane, Double_t */*sp*/=0x0) const;
   UChar_t   GetTRDncls(Int_t layer = -1) const;
   UChar_t   GetTRDntrackletsPID() const;
-  void      GetHMPIDpid(Double_t *p) const { if (fDetPid) fDetPid->GetHMPIDprobs(p); }
+  void      GetHMPIDpid(Double_t */*p*/) const { return; } // TODO: To be implemented properly with the new HMPID object
 
   const AliAODEvent* GetAODEvent(){return fAODEvent;}
   void SetAODEvent(const AliAODEvent* ptr){fAODEvent = ptr;}
index 674024a..f645446 100644 (file)
 #pragma link C++ class AliAODEvent+;
 #pragma link C++ class AliAODHeader+;
 #pragma link C++ class AliAODTrack+;
+
+#pragma read                                              \
+    sourceClass="AliAODPid"                               \
+    targetClass="AliAODPid"                               \
+    source="UShort_t fTPCsignalN; Double_t fTPCmomentum; Float_t fTRDmomentum[6]"  \
+    version="[-10]"                                       \
+    target="fTPCsignalN, fTPCmomentum, fTRDnSlices, fTRDmomentum"                                          \
+    code="{fTPCsignalN=(UChar_t)onfile.fTPCsignalN; fTPCmomentum=onfile.fTPCmomentum; for (Int_t i=0;i<6;++i) fTRDmomentum[i]=onfile.fTRDmomentum[i];}" 
+
 #pragma link C++ class AliAODPid+;
 #pragma link C++ class AliAODVertex+;
 #pragma link C++ class AliAODCluster+;
index 5c3e123..5a1b4f1 100644 (file)
@@ -15,12 +15,12 @@ public:
   
 private:
 
-  Double32_t  fTPCsignalRegion[4]; // TPC dEdx signal in 4 different regions - 0 - IROC, 1- OROC medium, 2 - OROC long, 3- OROC all, (default truncation used)  
+  Double32_t  fTPCsignalRegion[4]; //[0.,0.,10] TPC dEdx signal in 4 different regions - 0 - IROC, 1- OROC medium, 2 - OROC long, 3- OROC all, (default truncation used)  
   Char_t      fTPCsignalNRegion[3]; // number of clusters above threshold used in the dEdx calculation
   Char_t      fTPCsignalNRowRegion[3]; // number of crosed rows used in the dEdx calculation - signal below threshold included
 
   
-  ClassDef(AliTPCdEdxInfo,1)
+  ClassDef(AliTPCdEdxInfo,2)
 };
 
 #endif