AliAnalysisTaskMuonTrackingEff.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKMUONTRACKINGEFF_H
#define ALIANALYSISTASKMUONTRACKINGEFF_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */ 

//Author: Nicolas LE BRIS - SUBATECH Nantes

#include "AliAnalysisTaskSE.h"
#include "AliMuonTrackCuts.h"
#include "TString.h"

class TList;
class TObjArray;
class AliCounterCollection;
class AliMUONGeometryTransformer;
class AliMUONTrack;
class AliMUONTrackParam;
class AliMpArea;
class AliMpPad;
class AliMUON2DMap;


class AliAnalysisTaskMuonTrackingEff : public AliAnalysisTaskSE
{
  
 public:
  
  AliAnalysisTaskMuonTrackingEff();
  AliAnalysisTaskMuonTrackingEff(TString name);
  virtual ~AliAnalysisTaskMuonTrackingEff();

  void SetDefaultStorage(const char* ocdbPath) { fOCDBpath = ocdbPath; }
  
  void SetAlignStorage(const char* ocdbPath) { fAlignOCDBpath = ocdbPath; }
  
  void SetRecoParamStorage(const char* ocdbPath) { fRecoParamOCDBpath = ocdbPath; }
  
  void SelectCentrality(Double_t min, Double_t max) {fCentMin = min; fCentMax = max;}
  
  // set cuts to select tracks to be considered
  void SetMuonTrackCuts(AliMuonTrackCuts &trackCuts);
  
  // set default cuts to select tracks to be considered
  void SetDefaultMuonTrackCuts(Bool_t isMC);
  
  // get cuts to select tracks to be considered
  AliMuonTrackCuts* MuonTrackCuts() {return fMuonTrackCuts;}
  
  void SetMuonPtCut(Double_t cut) {fPtCut = cut;}
  
  void UseMCLabel(Bool_t flag = kTRUE) { fUseMCLabel = flag; }
  
  void EnableDisplay(Bool_t flag = kTRUE) { fEnableDisplay = flag; }
  
  // Implementation of interface methods
  virtual void UserCreateOutputObjects();
  virtual void UserExec(Option_t *);
  virtual void NotifyRun();
  virtual void Terminate(Option_t *);
  
  
 private:
  
  AliAnalysisTaskMuonTrackingEff(const AliAnalysisTaskMuonTrackingEff& rhs);
  AliAnalysisTaskMuonTrackingEff& operator = (const AliAnalysisTaskMuonTrackingEff& rhs);
  
  // Identify clusters/chambers that can be removed from the track
  Bool_t TagRemovableClusters(AliMUONTrack &track, Bool_t removableChambers[10]);
  
  // Find which detection elements should have been hit and record the missing clusters
  void FindAndRecordMissingClusters(AliMUONTrackParam &param, Int_t chamber, Double_t trackInfo[6]);
  
  // Find the intersection point between the track (assuming straight line) and the DE in the global frame
  void Intersect(AliMUONTrackParam &param, Int_t deId, Double_t p[3]);
  
  // Check whether (global) area overlaps with the given DE
  Bool_t OverlapDE(AliMpArea &area, Int_t deId);
  
  // Register the cluster in the given stores
  void RecordCluster(Int_t chamber, Int_t deId, AliMpPad pad[2], Double_t trackInfo[6],
         TString clusterKey, TList *chamberHistList, Bool_t recordChamber);
  
  Bool_t FindPads(Int_t deId, Double_t pos[3], AliMpPad pad[2]);
  
  
private:
  
  static const Int_t fgkNofDE[11];  
  static const Int_t fgkNofBusPath; 
  static const Int_t fgkNofManu;    
  
  Bool_t   fOCDBLoaded;             
  TString  fOCDBpath;               
  TString  fAlignOCDBpath;          
  TString  fRecoParamOCDBpath;      
  Double_t fCentMin;                
  Double_t fCentMax;                
  AliMuonTrackCuts* fMuonTrackCuts; 
  Double_t fPtCut;                  
  Bool_t   fUseMCLabel;             
  Bool_t   fEnableDisplay;          

  AliMUONGeometryTransformer *fTransformer; 
  
  TObjArray *fDEPlanes; 
  
  AliCounterCollection* fClusters; 
  AliCounterCollection* fEvents; 
  TList* fChamberTDHistList; 
  TList* fChamberTTHistList; 
  TList* fChamberSDHistList; 
  TList* fExtraHistList;     


  ClassDef(AliAnalysisTaskMuonTrackingEff, 5)
  
};


//________________________________________________________________________
inline void AliAnalysisTaskMuonTrackingEff::SetMuonTrackCuts(AliMuonTrackCuts &trackCuts)
{
  delete fMuonTrackCuts;
  fMuonTrackCuts = new AliMuonTrackCuts(trackCuts);
}

//________________________________________________________________________
inline void AliAnalysisTaskMuonTrackingEff::SetDefaultMuonTrackCuts(Bool_t isMC)
{
  delete fMuonTrackCuts;
  fMuonTrackCuts = new AliMuonTrackCuts("stdCuts", "stdCuts");
  fMuonTrackCuts->SetAllowDefaultParams();
  fMuonTrackCuts->SetFilterMask(AliMuonTrackCuts::kMuMatchLpt | AliMuonTrackCuts::kMuEta |
                                AliMuonTrackCuts::kMuThetaAbs | AliMuonTrackCuts::kMuPdca);
  fMuonTrackCuts->SetIsMC(isMC);
}

#endif