PWGPP-71 Move TPC noise filtering information form the AliESDevent to the AliESDHeader
authormivanov <marian.ivanov@cern.ch>
Sun, 2 Nov 2014 21:23:12 +0000 (22:23 +0100)
committermivanov <marian.ivanov@cern.ch>
Sun, 2 Nov 2014 21:23:12 +0000 (22:23 +0100)
STEER/ESD/AliESDEvent.cxx
STEER/ESD/AliESDEvent.h
STEER/ESD/AliESDHeader.cxx
STEER/ESD/AliESDHeader.h
TPC/Rec/AliTPCtracker.cxx

index 0a36525..3d61dae 100644 (file)
@@ -177,7 +177,6 @@ AliESDEvent::AliESDEvent():
   fDAQDetectorPattern(0xFFFF),
   fDAQAttributes(0xFFFF)
 {
-  for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=0;
 }
 //______________________________________________________________________________
 AliESDEvent::AliESDEvent(const AliESDEvent& esd):
@@ -273,7 +272,6 @@ AliESDEvent::AliESDEvent(const AliESDEvent& esd):
   AddObject(fMuonPads);
   GetStdContent();
   ConnectTracks();
-  for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=esd. fTPCNoiseFilterCounter[itype];
 
 }
 
index 92547da..8e3327e 100644 (file)
@@ -531,8 +531,6 @@ public:
   void SetDAQAttributes(UInt_t attributes) {fDAQAttributes = attributes;}
   UInt_t GetDAQDetectorPattern() const {return fDAQDetectorPattern;}
   UInt_t GetDAQAttributes() const {return fDAQAttributes;}
-  Char_t GetTPCNoiseFilterCounter(UInt_t index) {return fTPCNoiseFilterCounter[index%3];};
-  void SetTPCNoiseFilterCounter(UInt_t index,UChar_t value) {fTPCNoiseFilterCounter[index%3]=value;};
 protected:
   AliESDEvent(const AliESDEvent&);
   static Bool_t ResetWithPlacementNew(TObject *pObject);
@@ -600,7 +598,6 @@ protected:
   ULong64_t fDetectorStatus; // set detector event status bit for good event selection
   UInt_t fDAQDetectorPattern; // Detector pattern from DAQ: bit 0 is SPD, bit 4 is TPC, etc. See event.h
   UInt_t fDAQAttributes; // Third word of attributes from DAQ: bit 7 corresponds to HLT decision 
-  UChar_t fTPCNoiseFilterCounter[3];  // filter counter [0]=sector, [1]-timebin/sector, [2]-padrowsector 
 
   ClassDef(AliESDEvent,23)  //ESDEvent class 
 };
index 403fddc..09f3136 100644 (file)
@@ -59,6 +59,7 @@ AliESDHeader::AliESDHeader() :
   SetName("AliESDHeader");
   for(Int_t i = 0; i<kNMaxIR ; i++) fIRArray[i] = 0;
   fTriggerInputsNames.SetOwner(kTRUE);
+  for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=0;
 }
 
 AliESDHeader::~AliESDHeader() 
@@ -109,6 +110,7 @@ AliESDHeader::AliESDHeader(const AliESDHeader &header) :
     AliTriggerIR *ir = (AliTriggerIR*)((header.fIRBufferArray).At(i));
     if (ir) fIRBufferArray.Add(new AliTriggerIR(*ir));
   }
+  for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=header.fTPCNoiseFilterCounter[itype];
 }
 
 AliESDHeader& AliESDHeader::operator=(const AliESDHeader &header)
@@ -154,6 +156,7 @@ AliESDHeader& AliESDHeader::operator=(const AliESDHeader &header)
       AliTriggerIR *ir = (AliTriggerIR*)((header.fIRBufferArray).At(i));
       if (ir) fIRBufferArray.Add(new AliTriggerIR(*ir));
     }
+    for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=header.fTPCNoiseFilterCounter[itype];
   }
   return *this;
 }
@@ -201,7 +204,7 @@ void AliESDHeader::Reset()
    delete fIRArray[i];
    fIRArray[i]=0;
   }
-
+  for (Int_t itype=0; itype<3; itype++) fTPCNoiseFilterCounter[itype]=0;
   fIRBufferArray.Delete();
 }
 //______________________________________________________________________________
index 75c6a4c..c362ebf 100644 (file)
@@ -90,6 +90,8 @@ public:
   void      Print(const Option_t *opt=0) const;
 
   enum {kNTriggerInputs = 60};   //24 L0, 24 L1 and 12 L2 inputs
+  Char_t GetTPCNoiseFilterCounter(UInt_t index) {return fTPCNoiseFilterCounter[index%3];};
+  void SetTPCNoiseFilterCounter(UInt_t index,UChar_t value) {fTPCNoiseFilterCounter[index%3]=value;};
 
 private:
   void   SetIRInteractionMap() const;
@@ -120,6 +122,7 @@ private:
   TObjArray    fIRBufferArray;// Array with interaction records before and after triggered event
   mutable TBits   fIRInt2InteractionsMap;  // map of the Int2 events (normally 0TVX) near the event, that's Int2Id-EventId within -90 +90 BXs
   mutable TBits   fIRInt1InteractionsMap;  // map of the Int1 events (normally V0A&V0C) near the event, that's Int1Id-EventId within -90 +90 BXs
+  UChar_t fTPCNoiseFilterCounter[3];  // filter counter [0]=sector, [1]-timebin/sector, [2]-padrowsector 
 
 
   ClassDef(AliESDHeader,12)
index 311b72d..7314ded 100644 (file)
@@ -343,8 +343,11 @@ Int_t AliTPCtracker::AcceptCluster(AliTPCseed * seed, AliTPCclusterMI * cluster)
     Int_t nclSeed=seed->GetNumberOfClusters();
     
     if (AliTPCReconstructor::StreamLevel()&kStreamErrParam) { // flag:stream in debug mode cluster and track extrapolation at given row together with error nad shape estimate
+      Int_t eventNr = fEvent->GetEventNumberInFile();
+       
     (*fDebugStreamer)<<"ErrParam"<<
       "iter="<<fIteration<<
+      "eventNr="<<eventNr<<
       "Cl.="<<cluster<<
       "nclSeed="<<nclSeed<<
       "T.="<<&param<<
@@ -1634,10 +1637,10 @@ void    AliTPCtracker::FilterOutlierClusters(){
   Int_t filteredSector= vecSectorOut9.Sum();                  // light version of export variable
   Int_t filteredSectorTime= vecMedianSectorTimeOut9.Sum();
   Int_t filteredSectorPadRow= vecMedianSectorPadRowOut9.Sum();
-  if (fEvent){
-    fEvent->SetTPCNoiseFilterCounter(0,TMath::Min(filteredSector,255));
-    fEvent->SetTPCNoiseFilterCounter(1,TMath::Min(filteredSectorTime,255));
-    fEvent->SetTPCNoiseFilterCounter(2,TMath::Min(filteredSectorPadRow,255));
+  if (fEvent) if (fEvent->GetHeader()){
+    fEvent->GetHeader()->SetTPCNoiseFilterCounter(0,TMath::Min(filteredSector,255));
+    fEvent->GetHeader()->SetTPCNoiseFilterCounter(1,TMath::Min(filteredSectorTime,255));
+    fEvent->GetHeader()->SetTPCNoiseFilterCounter(2,TMath::Min(filteredSectorPadRow,255));
   }
  
   //
@@ -3118,6 +3121,7 @@ void AliTPCtracker::DumpClusters(Int_t iter, TObjArray *trackArray)
     }
   }
   //
+  Int_t eventNr = fEvent->GetEventNumberInFile();
 
   for (Int_t sec=0;sec<fkNIS;sec++){
     for (Int_t row=0;row<fInnerSec->GetNRows();row++){
@@ -3126,6 +3130,7 @@ void AliTPCtracker::DumpClusters(Int_t iter, TObjArray *trackArray)
        AliTPCclusterMI* cl = (AliTPCclusterMI*)cla->At(icl);
        Float_t gx[3];  cl->GetGlobalXYZ(gx);
        (*fDebugStreamer)<<"clDump"<< 
+         "eventNr="<<eventNr<<
          "iter="<<iter<<
          "cl.="<<cl<<      
          "gx0="<<gx[0]<<
@@ -3138,6 +3143,7 @@ void AliTPCtracker::DumpClusters(Int_t iter, TObjArray *trackArray)
        AliTPCclusterMI* cl = (AliTPCclusterMI*)cla->At(icl);
        Float_t gx[3];  cl->GetGlobalXYZ(gx);
        (*fDebugStreamer)<<"clDump"<< 
+         "eventNr="<<eventNr<<
          "iter="<<iter<<
          "cl.="<<cl<<
          "gx0="<<gx[0]<<
@@ -3156,6 +3162,7 @@ void AliTPCtracker::DumpClusters(Int_t iter, TObjArray *trackArray)
        AliTPCclusterMI* cl = (AliTPCclusterMI*) cla->At(icl);
        cl->GetGlobalXYZ(gx);
        (*fDebugStreamer)<<"clDump"<< 
+         "eventNr="<<eventNr<<
          "iter="<<iter<<
          "cl.="<<cl<<
          "gx0="<<gx[0]<<
@@ -3169,6 +3176,7 @@ void AliTPCtracker::DumpClusters(Int_t iter, TObjArray *trackArray)
        AliTPCclusterMI* cl = (AliTPCclusterMI*) cla->At(icl);
        cl->GetGlobalXYZ(gx);
        (*fDebugStreamer)<<"clDump"<< 
+         "eventNr="<<eventNr<<
          "iter="<<iter<<
          "cl.="<<cl<<
          "gx0="<<gx[0]<<