/// \class AliMUONVTrackReconstructor
/// \brief Virtual class for the MUON track reconstruction
///
-/////////////////////////////////////////////
-/// Virtual MUON track reconstructor in ALICE
-/////////////////////////////////////////////
+// Author: Philippe Pillot
#include <TObject.h>
-class TClonesArray;
-class AliMUONData;
+class AliMUONTrack;
+class AliMUONTrackParam;
+class AliMUONVCluster;
class AliMUONTriggerTrack;
+class AliMUONTrackHitPattern;
+class AliMUONVClusterServer;
+class AliMUONVClusterStore;
+class AliMUONVTrackStore;
+class AliMUONVTriggerTrackStore;
+class AliMUONVTriggerStore;
+class AliMUONGeometryTransformer;
+class AliMUONDigitMaker;
+class AliMUONTriggerCircuit;
+class TClonesArray;
+class AliMUONRecoParam;
class AliMUONVTrackReconstructor : public TObject {
public:
- AliMUONVTrackReconstructor(AliMUONData* data); // default Constructor
+ AliMUONVTrackReconstructor(const AliMUONRecoParam* recoParam, AliMUONVClusterServer* clusterServer); // default Constructor
virtual ~AliMUONVTrackReconstructor(); // Destructor
- // Parameters for track reconstruction: public methods
- // Get and Set, Set to defaults
- /// Return minimum value (GeV/c) of momentum in bending plane
- Double_t GetMinBendingMomentum(void) const {return fMinBendingMomentum;}
- /// Return chamber resolution (cm) in bending plane
- Double_t GetBendingResolution(void) const {return fBendingResolution;}
- /// Return chamber resolution (cm) in non-bending plane
- Double_t GetNonBendingResolution(void) const {return fNonBendingResolution;}
- /// Return chamber thickness in number of radiation lengths
- Double_t GetChamberThicknessInX0(void) const {return fChamberThicknessInX0;}
-
- // Reconstructed tracks
- /// Return number of reconstructed tracks
- Int_t GetNRecTracks() const {return fNRecTracks;} // Number
- /// Set number of reconstructed tracks
- void SetNRecTracks(Int_t NRecTracks) {fNRecTracks = NRecTracks;}
- /// Return array of reconstructed tracks
- TClonesArray* GetRecTracksPtr(void) const {return fRecTracksPtr;} // Array
-
- // Functions
- void EventReconstruct(void);
- void EventReconstructTrigger(void);
- virtual void EventDump(void) = 0; // dump reconstructed event
- void EventDumpTrigger(void); // dump reconstructed trigger event
+
+ void EventReconstruct(AliMUONVClusterStore& clusterStore,
+ AliMUONVTrackStore& trackStore);
- // needed in MakeSegmentsPerStation !!!
- Double_t GetImpactParamFromBendingMomentum(Double_t BendingMomentum) const;
- Double_t GetBendingMomentumFromImpactParam(Double_t ImpactParam) const;
+ void EventReconstructTrigger(const AliMUONTriggerCircuit& triggerCircuit,
+ const AliMUONVTriggerStore& triggerStore,
+ AliMUONVTriggerTrackStore& triggerTrackStore);
- /// Return MUON data
- AliMUONData* GetMUONData() {return fMUONData;}
-
- /// Set trigger circuit
- void SetTriggerCircuit(TClonesArray* circuit) {fTriggerCircuit = circuit;}
-
+ void ValidateTracksWithTrigger(AliMUONVTrackStore& trackStore,
+ const AliMUONVTriggerTrackStore& triggerTrackStore,
+ const AliMUONVTriggerStore& triggerStore,
+ const AliMUONTrackHitPattern& trackHitPattern);
- protected:
-
- // Defaults parameters for reconstruction
- static const Double_t fgkDefaultMinBendingMomentum; ///< default min. bending momentum for reconstruction
- static const Double_t fgkDefaultMaxBendingMomentum; ///< default max. bending momentum for reconstruction
- static const Double_t fgkDefaultBendingResolution; ///< default bending coordinate resolution for reconstruction
- static const Double_t fgkDefaultNonBendingResolution; ///< default non bending coordinate resolution for reconstruction
- static const Double_t fgkDefaultMaxSigma2Distance; ///< default square of max. distance for window size
- // Simple magnetic field:
- // Value taken from macro MUONtracking.C: 0.7 T, hence 7 kG
- // Length and Position from reco_muon.F, with opposite sign:
- // Length = ZMAGEND-ZCOIL
- // Position = (ZMAGEND+ZCOIL)/2
- // to be ajusted differently from real magnetic field ????
- static const Double_t fgkDefaultSimpleBValue; ///< default value of magnetic field (dipole)
- static const Double_t fgkDefaultSimpleBLength; ///< default length of magnetic field (dipole)
- static const Double_t fgkDefaultSimpleBPosition; ///< default position of magnetic field (dipole)
+ /// Return reco parameters
+ const AliMUONRecoParam* GetRecoParam() const { return fkRecoParam; }
- // Parameters for track reconstruction
- Double_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane
- Double_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane
- Double_t fBendingResolution; ///< chamber resolution (cm) in bending plane
- Double_t fNonBendingResolution; ///< chamber resolution (cm) in non bending plane
- Double_t fMaxSigma2Distance; ///< maximum square distance in units of the variance (maximum chi2)
- Double_t fChamberThicknessInX0; ///< chamber thickness in number of radiation lengths
- Double_t fSimpleBValue; ///< simple magnetic field: value (kG)
- Double_t fSimpleBLength; ///< simple magnetic field: length (cm)
- Double_t fSimpleBPosition; ///< simple magnetic field: Z central position (cm)
+ /// Re-fit the given track
+ virtual Bool_t RefitTrack(AliMUONTrack &track, Bool_t enableImprovement = kTRUE) = 0;
- Double_t* fSegmentMaxDistBending; ///< maximum distance (cm) for segments in bending plane
- Double_t* fSegmentMaxDistNonBending; ///< maximum distance (cm) for segments in non bending plane
- // Hits for reconstruction (should be in AliMUON ????)
- TClonesArray* fHitsForRecPtr; ///< pointer to the array of hits for reconstruction
- Int_t fNHitsForRec; ///< number of hits for reconstruction
- // Information per chamber (should be in AliMUONChamber ????)
- Int_t* fNHitsForRecPerChamber; ///< number of HitsForRec
- Int_t* fIndexOfFirstHitForRecPerChamber; ///< index (0...) of first HitForRec
-
- // Segments inside a station
- TClonesArray** fSegmentsPtr; ///< array of pointers to the segments for each station
- Int_t* fNSegments; ///< number of segments for each station
-
- // Reconstructed tracks
+ protected:
+
TClonesArray *fRecTracksPtr; ///< pointer to array of reconstructed tracks
Int_t fNRecTracks; ///< number of reconstructed tracks
- // data container
- AliMUONData* fMUONData; ///< Data container for MUON subsystem
+ AliMUONVClusterServer* fClusterServer; ///< reference to our cluster server
+ const AliMUONRecoParam* fkRecoParam; ///< reference to reco parameters
+
+ Double_t* fMaxMCSAngle2; ///< maximum angle dispersion due to MCS
+
// Functions
AliMUONVTrackReconstructor (const AliMUONVTrackReconstructor& rhs); ///< copy constructor
AliMUONVTrackReconstructor& operator=(const AliMUONVTrackReconstructor& rhs); ///< assignment operator
- void SortHitsForRecWithIncreasingChamber();
- void MakeSegmentsPerStation(Int_t Station);
+ /// Make track candidates from clusters in stations(1..) 4 and 5
+ virtual Bool_t MakeTrackCandidates(AliMUONVClusterStore& clusterStore) = 0;
+ /// Make extra track candidates from clusters in stations(1..) 4 and 5
+ virtual Bool_t MakeMoreTrackCandidates(AliMUONVClusterStore& clusterStore) = 0;
+ /// Follow tracks in stations(1..) 3, 2 and 1
+ virtual Bool_t FollowTracks(AliMUONVClusterStore& clusterStore) = 0;
+ /// Complement the reconstructed tracks
+ virtual Bool_t ComplementTracks(const AliMUONVClusterStore& clusterStore) = 0;
+ void ImproveTracks();
+ /// Improve the given reconstructed track
+ virtual void ImproveTrack(AliMUONTrack &track) = 0;
+ void Finalize();
+ /// Finalize the given track
+ virtual Bool_t FinalizeTrack(AliMUONTrack &track) = 0;
+
+ Bool_t IsAcceptable(AliMUONTrackParam &trackParam);
+
+ TClonesArray* MakeSegmentsBetweenChambers(const AliMUONVClusterStore& clusterStore, Int_t ch1, Int_t ch2);
- virtual void AddHitsForRecFromRawClusters() = 0;
- virtual void MakeSegments(void) = 0;
- virtual void MakeTracks(void) = 0;
- virtual void MakeTrackCandidates(void) = 0;
- virtual void FollowTracks(void) = 0;
- virtual void RemoveDoubleTracks(void) = 0;
+ void RemoveUsedSegments(TClonesArray& segments);
+ void RemoveIdenticalTracks();
+ void RemoveDoubleTracks();
+ void RemoveConnectedTracks(Bool_t inSt345 = kFALSE);
- private:
+ void AskForNewClustersInStation(const AliMUONTrackParam &trackParam,
+ AliMUONVClusterStore& clusterStore, Int_t station);
+ void AskForNewClustersInChamber(const AliMUONTrackParam &trackParam,
+ AliMUONVClusterStore& clusterStore, Int_t chamber);
- AliMUONTriggerTrack* fTriggerTrack; ///< Trigger track structure
+ Double_t TryOneCluster(const AliMUONTrackParam &trackParam, AliMUONVCluster* cluster,
+ AliMUONTrackParam &trackParamAtCluster, Bool_t updatePropagator = kFALSE);
+ Bool_t TryOneClusterFast(const AliMUONTrackParam &trackParam, const AliMUONVCluster* cluster);
+ Double_t TryTwoClustersFast(const AliMUONTrackParam &trackParamAtCluster1, AliMUONVCluster* cluster2,
+ AliMUONTrackParam &trackParamAtCluster2);
+
+ Bool_t FollowLinearTrackInChamber(AliMUONTrack &trackCandidate, const AliMUONVClusterStore& clusterStore, Int_t nextChamber);
+ Bool_t FollowLinearTrackInStation(AliMUONTrack &trackCandidate, const AliMUONVClusterStore& clusterStore, Int_t nextStation);
- TClonesArray* fTriggerCircuit; //!< trigger circuit array
+ void CorrectThetaRange(Float_t&);
+
+ private:
// Functions
- void SetReconstructionParametersToDefaults(void);
-
- void ResetTracks(void);
- void ResetSegments(void);
- void ResetHitsForRec(void);
+ void ResetTracks();
- void ValidateTracksWithTrigger(void);
- Bool_t MakeTriggerTracks(void);
-
-
ClassDef(AliMUONVTrackReconstructor, 0) // MUON track reconstructor in ALICE
- };
+};
#endif