Set of patches for report #61429: PID: Separating response functions from ESD (Marco)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Jan 2010 08:14:31 +0000 (08:14 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Jan 2010 08:14:31 +0000 (08:14 +0000)
20 files changed:
ANALYSIS/AliBackgroundSelection.cxx
ANALYSIS/AliBackgroundSelection.h
ANALYSIS/AliESDv0Cuts.h
ANALYSIS/AliKineTrackCuts.cxx
ANALYSIS/AliKineTrackCuts.h
CORRFW/AliCFPairIsPrimaryCuts.cxx
CORRFW/AliCFPairIsPrimaryCuts.h
CORRFW/AliCFPairPidCut.cxx
CORRFW/AliCFPairPidCut.h
CORRFW/AliCFPairQualityCuts.cxx
CORRFW/AliCFPairQualityCuts.h
CORRFW/AliCFV0TopoCuts.cxx
CORRFW/AliCFV0TopoCuts.h
HMPID/AliHMPIDReconstructor.cxx
HMPID/AliHMPIDTracker.cxx
PWG0/dNdPt/AlidNdPtAcceptanceCuts.h
PWG0/dNdPt/AlidNdPtEventCuts.h
PWG3/muon/AliESDMuonTrackCuts.h
STEER/AliITSPIDResponse.cxx
STEER/AliITSPIDResponse.h

index fdd5bb3..503c326 100644 (file)
@@ -54,7 +54,7 @@ AliBackgroundSelection::~AliBackgroundSelection() {
 \r
 }\r
 \r
-Bool_t AliBackgroundSelection::IsSelected(TObject* obj){\r
+Bool_t AliBackgroundSelection::IsSelected(TObject* obj, TObject* /* event */){\r
 \r
   // reset fSelected\r
   SetSelected(kFALSE);\r
index a92c91a..da49742 100644 (file)
@@ -23,8 +23,8 @@ public:
   AliBackgroundSelection(const char* name, const char* title);
   AliBackgroundSelection(const AliBackgroundSelection& obj);  
   virtual ~AliBackgroundSelection();
-  virtual Bool_t IsSelected(TObject* obj);
-  virtual Bool_t IsSelected(TList*  ) {AliFatal("Not implemented");return 0;}
+  virtual Bool_t IsSelected(TObject* obj, TObject * /*event*/ = 0);
+  //virtual Bool_t IsSelected(TList*  ) {AliFatal("Not implemented");return 0;}
   virtual void   Init();
   virtual TList * GetOutput() {return fOutputHist;}
 
index 25bc86f..4dbd9af 100644 (file)
@@ -24,7 +24,7 @@ public:
   AliESDv0Cuts(const Char_t* name = "AliESDv0Cuts", const Char_t* title = "");
   virtual ~AliESDv0Cuts();
 
-  Bool_t IsSelected(TObject* /*obj*/) {return kTRUE;}
+  Bool_t IsSelected(TObject* /*obj*/, TObject* /* event */ = 0) {return kTRUE;}
   Bool_t IsSelected(TList* listObj);
   Bool_t IsSelected(TObject* obj1, TObject* obj2, TObject* obj3, TObject* obj4)
   {return AcceptV0((AliESDv0*) obj1, (AliESDtrack*) obj2, (AliESDtrack*) obj3, (const AliESDVertex*) obj4);}
index 60298ab..2c73476 100644 (file)
@@ -65,7 +65,7 @@ AliKineTrackCuts::AliKineTrackCuts(const Char_t* name, const Char_t* title) :
 
 
 //____________________________________________________________________
-Bool_t  AliKineTrackCuts::IsSelected(TObject* obj)
+Bool_t  AliKineTrackCuts::IsSelected(TObject* obj, TObject * /*event*/ )
 {
 
   TParticle * part = (TParticle *)obj;
index b3b3f6a..c9dba13 100644 (file)
@@ -16,7 +16,7 @@ public:
                       AliKineTrackCuts(const Char_t* name = "AliKineTrackCuts", const Char_t* title = "");
   virtual            ~AliKineTrackCuts(){;}
   
-              Bool_t  IsSelected(TObject* obj);
+              Bool_t  IsSelected(TObject* obj, TObject * /*event */ = 0);
              Bool_t  IsSelected(TList* /*list*/) {return kTRUE;}
             
               void    SetFinalParticles( Bool_t val=kTRUE )          { fOnlyFinalParticles = val; }
index d73a735..045423e 100755 (executable)
@@ -71,7 +71,7 @@ AliCFPairIsPrimaryCuts& AliCFPairIsPrimaryCuts::operator=(const AliCFPairIsPrima
 }
 
 //__________________________________________________________________________________
-Bool_t AliCFPairIsPrimaryCuts::IsSelected(TObject* obj) {
+Bool_t AliCFPairIsPrimaryCuts::IsSelected(TObject* obj, TObject * /* obj2 */ ) {
   //
   // loops over decisions of single cuts and returns if the track is accepted
   //
index 4c7115f..667d3b8 100755 (executable)
@@ -36,8 +36,8 @@ class AliCFPairIsPrimaryCuts : public AliCFCutBase
   AliCFPairIsPrimaryCuts& operator=(const AliCFPairIsPrimaryCuts& c) ;
   virtual ~AliCFPairIsPrimaryCuts() {delete fCutNeg; delete fCutPos; }
 
-  Bool_t IsSelected(TObject* obj) ;
-  Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+  Bool_t IsSelected(TObject* obj, TObject * /* obj2 */ = 0) ;
+  //Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
   virtual AliCFTrackIsPrimaryCuts* GetNegCut() const {return fCutNeg;}
   virtual AliCFTrackIsPrimaryCuts* GetPosCut() const {return fCutPos;}
 
index d18bfd5..c0d8950 100755 (executable)
@@ -69,7 +69,7 @@ AliCFPairPidCut& AliCFPairPidCut::operator=(const AliCFPairPidCut& c)
 }
 
 //__________________________________________________________________________________
-Bool_t AliCFPairPidCut::IsSelected(TObject* obj) {
+Bool_t AliCFPairPidCut::IsSelected(TObject* obj, TObject* /* obj2 */) {
   //
   // loops over decisions of single cuts and returns if the track is accepted
   //
index 74a1abc..8af42cf 100755 (executable)
@@ -43,8 +43,8 @@ class AliCFPairPidCut : public AliCFCutBase
   virtual AliCFTrackCutPid* GetNegCut() const {return fCutNeg;}
   virtual AliCFTrackCutPid* GetPosCut() const {return fCutPos;}
 
-  Bool_t IsSelected(TObject *obj); //boolean for detectors
-  Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+  Bool_t IsSelected(TObject *obj, TObject * /* obj2 */); //boolean for detectors
+  // Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
  private:
   AliCFTrackCutPid* fCutNeg; // PID cut on negative daughter
   AliCFTrackCutPid* fCutPos; // PID cut on positive daughter
index dbc18bd..8a8c46b 100755 (executable)
@@ -70,7 +70,7 @@ AliCFPairQualityCuts& AliCFPairQualityCuts::operator=(const AliCFPairQualityCuts
 }
 
 //__________________________________________________________________________________
-Bool_t AliCFPairQualityCuts::IsSelected(TObject* obj) {
+Bool_t AliCFPairQualityCuts::IsSelected(TObject* obj, TObject* /* obj2 */) {
   //
   // loops over decisions of single cuts and returns if the track is accepted
   //
index d7fc4f6..34e8954 100755 (executable)
@@ -38,8 +38,8 @@ class AliCFPairQualityCuts : public AliCFCutBase
   AliCFPairQualityCuts& operator=(const AliCFPairQualityCuts& c) ;
   virtual ~AliCFPairQualityCuts() {delete fCutNeg; delete fCutPos; }
 
-  Bool_t IsSelected(TObject* obj) ; 
-  Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+  Bool_t IsSelected(TObject* obj, TObject* /* obj2 */ = 0) ; 
+  //Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
 
   virtual AliCFTrackQualityCuts* GetNegCut() const {return fCutNeg;}
   virtual AliCFTrackQualityCuts* GetPosCut() const {return fCutPos;}
index 0214f9e..5df04ab 100755 (executable)
@@ -95,7 +95,7 @@ AliCFV0TopoCuts& AliCFV0TopoCuts::operator=(const AliCFV0TopoCuts& c)
 }
 
 //______________________________________________________________
-Bool_t AliCFV0TopoCuts::IsSelected(TObject *obj) {
+Bool_t AliCFV0TopoCuts::IsSelected(TObject *obj, TObject * /* event */) {
   //
   // computes V0 topological variables to cut on and return true 
   // in case the V0 is accepted
index 72b944a..0fbd7dc 100755 (executable)
@@ -42,8 +42,8 @@ class AliCFV0TopoCuts : public AliCFCutBase
   AliCFV0TopoCuts (const AliCFV0TopoCuts& c) ;
   AliCFV0TopoCuts& operator=(const AliCFV0TopoCuts& c) ;
   virtual ~AliCFV0TopoCuts() { } ;
-  Bool_t IsSelected(TObject* v0) ;
-  Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+  Bool_t IsSelected(TObject* v0, TObject * /*event*/ = 0) ;
+  //Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
   void   SetRecEventInfo(const TObject* evt) {fEvent = (AliVEvent*)evt;}
   void   SetMaxDcaDaughters (Double32_t dca)  {fMaxDcaDaughters = dca;}
   void   SetMinDcaNeg       (Double32_t dca)  {fMinDcaNeg = dca;}
index 0895ab0..6d6c495 100644 (file)
@@ -20,8 +20,8 @@
 #include "AliHMPIDReconstructor.h" //class header
 #include "AliHMPID.h"              //Reconstruct() 
 #include "AliHMPIDCluster.h"       //Dig2Clu()
-//#include "AliHMPIDPid.h"           //FillEsd() 
-//#include "AliHMPIDParam.h"         //FillEsd() 
+#include "AliHMPIDPid.h"           //FillEsd() 
+#include "AliHMPIDParam.h"         //FillEsd() 
 #include <AliCDBEntry.h>           //ctor
 #include <AliCDBManager.h>         //ctor
 #include <AliESDEvent.h>           //FillEsd()
@@ -190,7 +190,6 @@ void AliHMPIDReconstructor::FillESD(TTree */*digitsTree*/, TTree */*clustersTree
 // Fill ESD with all the infos from HMPID
 // Probability vector from AliHMPIDPid
 //...
-  /*
   AliHMPIDPid pID;
   Double_t prob[AliPID::kSPECIES];
   
@@ -199,5 +198,4 @@ void AliHMPIDReconstructor::FillESD(TTree */*digitsTree*/, TTree */*clustersTree
     pID.FindPid(pTrk,AliPID::kSPECIES,prob);
     pTrk->SetHMPIDpid(prob);
   }//ESD tracks loop
-  */
 }//FillESD()
index f9ce9ee..793eeef 100644 (file)
@@ -282,13 +282,11 @@ Int_t AliHMPIDTracker::Recon(AliESDEvent *pEsd,TObjArray *pClus,TObjArray *pNmea
 
     if(pTrk->GetHMPIDsignal()<0) continue;
         
-    /* Probabilies are set in AliReconstruction/AliPIDesd
     AliHMPIDPid pID;
     Double_t prob[5];
     pID.FindPid(pTrk,5,prob);
     pTrk->SetHMPIDpid(prob);
 //      Printf(" Prob e- %6.2f mu %6.2f pi %6.2f k %6.2f p %6.2f",prob[0]*100,prob[1]*100,prob[2]*100,prob[3]*100,prob[4]*100);
-    */
     delete hmpTrk;hmpTrk=0x0;
   }//iTrk
 
@@ -388,14 +386,12 @@ Int_t AliHMPIDTracker::ReconHiddenTrk(AliESDEvent *pEsd,TObjArray *pClus,TObjArr
       }
     }
     //
-    /*
     if(!reconHTA.CkovHiddenTrk(pTrk,(TClonesArray *)pClus->At(ipCh),indMip,nmean)) {                 //search for track parameters and Cerenkov angle of this track
       AliHMPIDPid pID;
       Double_t prob[5];
       pID.FindPid(pTrk,5,prob);
       pTrk->SetHMPIDpid(prob);
     }
-    */
 //      Printf(" Prob e- %6.2f mu %6.2f pi %6.2f k %6.2f p %6.2f",prob[0]*100,prob[1]*100,prob[2]*100,prob[3]*100,prob[4]*100);
   }//iTrk
 
index 1f34f59..2212872 100644 (file)
@@ -42,8 +42,8 @@ public:
   void Init();\r
 \r
   // check MC tracks\r
-  virtual Bool_t IsSelected(TObject *) {return kTRUE;}\r
-  virtual Bool_t IsSelected(TList *) {return kTRUE;}\r
+  virtual Bool_t IsSelected(TObject *, TObject * /* event */ = 0) {return kTRUE;}\r
+  //virtual Bool_t IsSelected(TList *) {return kTRUE;}\r
 \r
   //\r
   Bool_t AcceptTrack(AliESDtrack *track);\r
index 7125259..b24bca9 100644 (file)
@@ -68,8 +68,8 @@ public:
   void Init();\r
 \r
   // check MC tracks\r
-  Bool_t IsSelected(TObject *) {return kTRUE;}\r
-  Bool_t IsSelected(TList *) {return kTRUE;}\r
+  Bool_t IsSelected(TObject *, TObject * /*event */ = 0) {return kTRUE;}\r
+  //Bool_t IsSelected(TList *) {return kTRUE;}\r
 \r
   // accept event\r
   Bool_t AcceptEvent(AliESDEvent *event=0, AliMCEvent *mcEvent=0, const AliESDVertex *vtx=0);\r
index a19196a..6796092 100644 (file)
@@ -28,9 +28,9 @@ class AliESDMuonTrackCuts : public AliAnalysisCuts
 public:
   AliESDMuonTrackCuts(const Char_t* name = "AliESDMuonTrackCuts", const Char_t* title = "");
   virtual ~AliESDMuonTrackCuts();
-  Bool_t IsSelected(TObject* obj)
+  Bool_t IsSelected(TObject* obj, TObject* /* obj2 */ = 0)
        {return AcceptTrack((AliESDMuonTrack*)obj);}
-  Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
+  // Bool_t IsSelected(TList* /*list*/) {return kTRUE;}
   Bool_t AcceptTrack(AliESDMuonTrack* esdMuTrack);
   TObjArray* GetAcceptedTracks(AliESD* esd);
   Int_t CountAcceptedTracks(AliESD* esd);
index ee0e406..9f05c0b 100644 (file)
@@ -26,7 +26,6 @@
 #include "AliITSPidParams.h"
 #include "AliExternalTrackParam.h"
 
-Float_t AliITSPIDResponse::fgMIP = 82;
 ClassImp(AliITSPIDResponse)
 
 AliITSPIDResponse::AliITSPIDResponse(): 
@@ -41,7 +40,7 @@ AliITSPIDResponse::AliITSPIDResponse():
 
 //_________________________________________________________________________
 AliITSPIDResponse::AliITSPIDResponse(Double_t *param): 
-  fRes(param[1]),
+  fRes(param[0]),
   fKp1(15.77),
   fKp2(4.95),
   fKp3(0.312),
@@ -51,7 +50,6 @@ AliITSPIDResponse::AliITSPIDResponse(Double_t *param):
   //
   //  The main constructor
   //
-  fgMIP = param[0];
 }
 
 
@@ -62,7 +60,7 @@ Double_t AliITSPIDResponse::Bethe(Double_t p,Double_t mass) const {
   //
   Double_t bb=
     AliExternalTrackParam::BetheBlochAleph(p/mass,fKp1,fKp2,fKp3,fKp4,fKp5);
-  return bb*fgMIP;
+  return bb;
 }
 
 Double_t AliITSPIDResponse::GetResolution(Double_t bethe) const {
index 900a114..b46918a 100644 (file)
@@ -32,7 +32,6 @@ public:
 private:
 
   // Data members for truncated mean method
-  static Float_t fgMIP;     // dE/dx for MIP (relative units)
   Float_t fRes;             // relative dEdx resolution
   Double_t fKp1;             // ALEPH BB param 1
   Double_t fKp2;             // ALEPH BB param 2