/* $Id$ */
-/// \class AliMUONVTrackReconstructor
////////////////////////////////////
///
+/// \class AliMUONVTrackReconstructor
/// Virtual MUON track reconstructor in ALICE (class renamed from AliMUONEventReconstructor)
///
/// This class contains as data:
/// * EventReconstruct to build the muon tracks
/// * EventReconstructTrigger to build the trigger tracks
///
+/// \author Philippe Pillot
+///
////////////////////////////////////
-#include <stdlib.h>
-#include <Riostream.h>
-#include <TClonesArray.h>
-#include <TMath.h>
-#include <TTree.h>
-
#include "AliMUONVTrackReconstructor.h"
#include "AliMUONData.h"
#include "AliMUONConstants.h"
#include "AliMUONTrack.h"
#include "AliMUONTrackParam.h"
#include "AliMUONTrackExtrap.h"
-#include "AliMagF.h"
+
#include "AliLog.h"
#include "AliTracker.h"
+#include <Riostream.h>
+#include <TClonesArray.h>
+#include <TMath.h>
+
+/// \cond CLASSIMP
ClassImp(AliMUONVTrackReconstructor) // Class implementation in ROOT context
+/// \endcond
//************* Defaults parameters for reconstruction
const Double_t AliMUONVTrackReconstructor::fgkDefaultMinBendingMomentum = 3.0;
fBendingVertexDispersion(fgkDefaultBendingVertexDispersion),
fNonBendingVertexDispersion(fgkDefaultNonBendingVertexDispersion),
fMaxNormChi2MatchTrigger(fgkDefaultMaxNormChi2MatchTrigger),
- fSegmentMaxDistBending(0x0),
- fSegmentMaxDistNonBending(0x0),
fHitsForRecPtr(0x0),
fNHitsForRec(0),
fNHitsForRecPerChamber(0x0),
fTriggerCircuit(0x0)
{
/// Constructor for class AliMUONVTrackReconstructor
- fSegmentMaxDistBending = new Double_t[AliMUONConstants::NTrackingSt()];
- fSegmentMaxDistNonBending = new Double_t[AliMUONConstants::NTrackingSt()];
fNHitsForRecPerChamber = new Int_t[AliMUONConstants::NTrackingCh()];
fIndexOfFirstHitForRecPerChamber = new Int_t[AliMUONConstants::NTrackingCh()];
- SetReconstructionParametersToDefaults();
-
// Memory allocation for the TClonesArray of hits for reconstruction
// Is 10000 the right size ????
fHitsForRecPtr = new TClonesArray("AliMUONHitForRec", 10000);
AliMUONVTrackReconstructor::~AliMUONVTrackReconstructor(void)
{
/// Destructor for class AliMUONVTrackReconstructor
- delete [] fSegmentMaxDistBending;
- delete [] fSegmentMaxDistNonBending;
delete [] fNHitsForRecPerChamber;
delete [] fIndexOfFirstHitForRecPerChamber;
delete fTriggerTrack;
delete fHitsForRecPtr;
}
- //__________________________________________________________________________
-void AliMUONVTrackReconstructor::SetReconstructionParametersToDefaults(void)
-{
- /// Set reconstruction parameters for making segments to default values
- // Would be much more convenient with a structure (or class) ????
-
- // ******** Parameters for making segments
- // should be parametrized ????
- // according to interval between chambers in a station ????
- // Maximum distance in non bending plane
- // 5 * 0.22 just to remember the way it was made in TRACKF_STAT
- // SIGCUT*DYMAX(IZ)
- for (Int_t st = 0; st < AliMUONConstants::NTrackingSt(); st++)
- fSegmentMaxDistNonBending[st] = 5. * 0.22;
- // Maximum distance in bending plane:
- // values from TRACKF_STAT, corresponding to (J psi 20cm),
- // scaled to the real distance between chambers in a station
- fSegmentMaxDistBending[0] = TMath::Abs( 1.5 *
- (AliMUONConstants::DefaultChamberZ(1) - AliMUONConstants::DefaultChamberZ(0)) / 20.0);
- fSegmentMaxDistBending[1] = TMath::Abs( 1.5 *
- (AliMUONConstants::DefaultChamberZ(3) - AliMUONConstants::DefaultChamberZ(2)) / 20.0);
- fSegmentMaxDistBending[2] = TMath::Abs( 3.0 *
- (AliMUONConstants::DefaultChamberZ(5) - AliMUONConstants::DefaultChamberZ(4)) / 20.0);
- fSegmentMaxDistBending[3] = TMath::Abs( 6.0 *
- (AliMUONConstants::DefaultChamberZ(7) - AliMUONConstants::DefaultChamberZ(6)) / 20.0);
- fSegmentMaxDistBending[4] = TMath::Abs( 6.0 *
- (AliMUONConstants::DefaultChamberZ(9) - AliMUONConstants::DefaultChamberZ(8)) / 20.0);
-
- return;
-}
-
//__________________________________________________________________________
void AliMUONVTrackReconstructor::EventReconstruct(void)
{
- // To reconstruct one event
+ /// To reconstruct one event
AliDebug(1,"Enter EventReconstruct");
+
ResetTracks(); //AZ
ResetHitsForRec(); //AZ
AddHitsForRecFromRawClusters();
AliMUONHitForRec *hit1Ptr, *hit2Ptr;
AliMUONObjectPair *segment;
- Double_t bendingSlope, distBend, distNonBend, extBendCoor, extNonBendCoor, extrapFact;
- Double_t impactParam = 0., bendingMomentum = 0.; // to avoid compilation warning
+ Double_t bendingSlope = 0, impactParam = 0., bendingMomentum = 0.; // to avoid compilation warning
// first and second chambers (0...) in the station
Int_t ch1 = 2 * station;
Int_t ch2 = ch1 + 1;
hit1++) {
// pointer to the HitForRec
hit1Ptr = (AliMUONHitForRec*) ((*fHitsForRecPtr)[hit1]);
- // extrapolation, on the straight line joining the HitForRec to the vertex (0,0,0),
- // to the Z of the HitForRec in the second chamber of the station
// Loop over HitsForRec's in the second chamber of the station
for (Int_t hit2 = fIndexOfFirstHitForRecPerChamber[ch2];
hit2 < fIndexOfFirstHitForRecPerChamber[ch2] + fNHitsForRecPerChamber[ch2];
hit2++) {
// pointer to the HitForRec
hit2Ptr = (AliMUONHitForRec*) ((*fHitsForRecPtr)[hit2]);
- // absolute values of distances, in bending and non bending planes,
- // between the HitForRec in the second chamber
- // and the previous extrapolation
- extrapFact = hit2Ptr->GetZ()/ hit1Ptr->GetZ();
- extBendCoor = extrapFact * hit1Ptr->GetBendingCoor();
- extNonBendCoor = extrapFact * hit1Ptr->GetNonBendingCoor();
- distBend = TMath::Abs(hit2Ptr->GetBendingCoor() - extBendCoor);
- distNonBend = TMath::Abs(hit2Ptr->GetNonBendingCoor() - extNonBendCoor);
- // bending slope
if ( hit1Ptr->GetZ() - hit2Ptr->GetZ() != 0. ) {
+ // bending slope
bendingSlope = (hit1Ptr->GetBendingCoor() - hit2Ptr->GetBendingCoor()) / (hit1Ptr->GetZ() - hit2Ptr->GetZ());
// impact parameter
impactParam = hit1Ptr->GetBendingCoor() - hit1Ptr->GetZ() * bendingSlope;
AliWarning("hit1Ptr->GetZ() = hit2Ptr->GetZ(): no segment created");
continue;
}
- // check for distances not too large,
- // and impact parameter not too big if stations downstream of the dipole.
- // Conditions "distBend" and "impactParam" correlated for these stations ????
- if ((distBend < fSegmentMaxDistBending[station]) && (distNonBend < fSegmentMaxDistNonBending[station]) &&
- (bendingMomentum < fMaxBendingMomentum) && (bendingMomentum > fMinBendingMomentum)) {
+ // check for bending momentum within tolerances
+ if ((bendingMomentum < fMaxBendingMomentum) && (bendingMomentum > fMinBendingMomentum)) {
// make new segment
segment = new ((*segments)[segments->GetLast()+1]) AliMUONObjectPair(hit1Ptr, hit2Ptr, kFALSE, kFALSE);
if (AliLog::GetGlobalDebugLevel() > 1) {
- cout << "segmentIndex(0...): " << segments->GetLast()
- << " distBend: " << distBend
- << " distNonBend: " << distNonBend
- << endl;
+ cout << "segmentIndex(0...): " << segments->GetLast() << endl;
segment->Dump();
cout << "HitForRec in first chamber" << endl;
hit1Ptr->Dump();
distTriggerTrack[0] = (triggerTrack->GetX11()-xTrack)/kDistSigma[0];
distTriggerTrack[1] = (triggerTrack->GetY11()-yTrack)/kDistSigma[1];
distTriggerTrack[2] = (TMath::Tan(triggerTrack->GetThetay())-ySlopeTrack)/kDistSigma[2];
- loTrgNum=triggerTrack->GetLoTrgNum();
chi2 = 0.;
for (Int_t iVar = 0; iVar < 3; iVar++) chi2 += distTriggerTrack[iVar]*distTriggerTrack[iVar];
chi2 /= 3.; // Normalized Chi2: 3 degrees of freedom (X,Y,slopeY)
minChi2MatchTrigger = chi2;
matchTrigger = kTRUE;
chi2MatchTrigger = chi2;
+ loTrgNum=triggerTrack->GetLoTrgNum();
}
triggerTrack = (AliMUONTriggerTrack*) recTriggerTracks->After(triggerTrack);
}
track->SetMatchTrigger(matchTrigger);
track->SetLoTrgNum(loTrgNum);
track->SetChi2MatchTrigger(chi2MatchTrigger);
-
track = (AliMUONTrack*) fRecTracksPtr->After(track);
}