////////////////////////////////////
#include "AliMUONVTrackReconstructor.h"
-#include "AliMUONData.h"
+#include "AliMUONRecData.h"
#include "AliMUONConstants.h"
#include "AliMUONHitForRec.h"
#include "AliMUONObjectPair.h"
#include "AliMUONTrack.h"
#include "AliMUONTrackParam.h"
#include "AliMUONTrackExtrap.h"
+#include "AliMUONTrackHitPattern.h"
#include "AliLog.h"
#include "AliTracker.h"
const Double_t AliMUONVTrackReconstructor::fgkDefaultMaxNormChi2MatchTrigger = 16.0;
//__________________________________________________________________________
-AliMUONVTrackReconstructor::AliMUONVTrackReconstructor(AliMUONData* data)
+AliMUONVTrackReconstructor::AliMUONVTrackReconstructor(AliMUONRecData* data)
: TObject(),
fMinBendingMomentum(fgkDefaultMinBendingMomentum),
fMaxBendingMomentum(fgkDefaultMaxBendingMomentum),
const AliMagF* kField = AliTracker::GetFieldMap();
if (!kField) AliFatal("No field available");
AliMUONTrackExtrap::SetField(kField);
+
+ fTrackHitPattern = new AliMUONTrackHitPattern(fMUONData);
}
//__________________________________________________________________________
delete [] fIndexOfFirstHitForRecPerChamber;
delete fTriggerTrack;
delete fHitsForRecPtr;
+ delete fTrackHitPattern;
}
//__________________________________________________________________________
MakeTracks();
if (fMUONData->IsTriggerTrackBranchesInTree())
ValidateTracksWithTrigger();
-
+
+ fTrackHitPattern->GetHitPattern(fRecTracksPtr);
+
// Add tracks to MUON data container
for(Int_t i=0; i<fNRecTracks; i++) {
AliMUONTrack * track = (AliMUONTrack*) fRecTracksPtr->At(i);
AliMUONTrack *track;
AliMUONTrackParam trackParam;
AliMUONTriggerTrack *triggerTrack;
+ AliMUONLocalTrigger *locTrg;
fMUONData->SetTreeAddress("RL");
fMUONData->GetRecTriggerTracks();
TClonesArray *recTriggerTracks = fMUONData->RecTriggerTracks();
-
- Bool_t matchTrigger;
+
+ fMUONData->SetTreeAddress("TC");
+ fMUONData->GetTrigger();
+ TClonesArray *localTrigger = fMUONData->LocalTrigger();
+
+ Int_t matchTrigger;
Int_t loTrgNum;
Double_t distTriggerTrack[3];
Double_t xTrack, yTrack, ySlopeTrack, chi2MatchTrigger, minChi2MatchTrigger, chi2;
track = (AliMUONTrack*) fRecTracksPtr->First();
while (track) {
- matchTrigger = kFALSE;
+ matchTrigger = -1;
chi2MatchTrigger = 0.;
loTrgNum = -1;
+ Int_t doubleMatch=-1; // Check if track matches 2 trigger tracks
+ Double_t doubleChi2 = -1.;
trackParam = *((AliMUONTrackParam*) (track->GetTrackParamAtHit()->Last()));
AliMUONTrackExtrap::ExtrapToZ(&trackParam, AliMUONConstants::DefaultChamberZ(10)); // extrap to 1st trigger chamber
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)
- if (chi2 < minChi2MatchTrigger && chi2 < fMaxNormChi2MatchTrigger) {
- minChi2MatchTrigger = chi2;
- matchTrigger = kTRUE;
- chi2MatchTrigger = chi2;
- loTrgNum=triggerTrack->GetLoTrgNum();
+ if (chi2 < fMaxNormChi2MatchTrigger) {
+ Bool_t isDoubleTrack = (TMath::Abs(chi2 - minChi2MatchTrigger)<1.);
+ if (chi2 < minChi2MatchTrigger && chi2 < fMaxNormChi2MatchTrigger) {
+ if(isDoubleTrack){
+ doubleMatch = loTrgNum;
+ doubleChi2 = chi2MatchTrigger;
+ }
+ minChi2MatchTrigger = chi2;
+ chi2MatchTrigger = chi2;
+ loTrgNum=triggerTrack->GetLoTrgNum();
+ locTrg = (AliMUONLocalTrigger*)localTrigger->UncheckedAt(loTrgNum);
+ matchTrigger=0;
+ if(locTrg->LoLpt()>0)matchTrigger=1;
+ if(locTrg->LoHpt()>0)matchTrigger=2;
+ }
+ else if(isDoubleTrack) {
+ doubleMatch = triggerTrack->GetLoTrgNum();
+ doubleChi2 = chi2;
+ }
}
triggerTrack = (AliMUONTriggerTrack*) recTriggerTracks->After(triggerTrack);
}
+ if(doubleMatch>=0){ // If two trigger tracks match, select the one passing more trigger cuts
+ AliDebug(1, Form("Two candidates found: %i and %i",loTrgNum,doubleMatch));
+ AliMUONLocalTrigger *locTrg1 = (AliMUONLocalTrigger*)localTrigger->UncheckedAt(doubleMatch);
+ if((locTrg1->LoLpt()>0 && matchTrigger<1) || (locTrg1->LoHpt() && matchTrigger<2)){
+ if(locTrg1->LoHpt()>0)matchTrigger=2;
+ else matchTrigger=1;
+ loTrgNum = doubleMatch;
+ chi2MatchTrigger=doubleChi2;
+ }
+ }
track->SetMatchTrigger(matchTrigger);
track->SetLoTrgNum(loTrgNum);
//__________________________________________________________________________
Bool_t AliMUONVTrackReconstructor::MakeTriggerTracks(void)
{
- // To make the trigger tracks from Local Trigger
+ /// To make the trigger tracks from Local Trigger
AliDebug(1, "Enter MakeTriggerTracks");
TTree* treeR;
Float_t x11 = 0.;
for (Int_t i=0; i<nlocals; i++) { // loop on Local Trigger
- locTrg = (AliMUONLocalTrigger*)localTrigger->UncheckedAt(i);
+ locTrg = (AliMUONLocalTrigger*)localTrigger->UncheckedAt(i);
+
+ Bool_t xTrig=kFALSE;
+ Bool_t yTrig=kFALSE;
+
+ if ( locTrg->LoSdev()==1 && locTrg->LoDev()==0 &&
+ locTrg->LoStripX()==0) xTrig=kFALSE; // no trigger in X
+ else xTrig=kTRUE; // trigger in X
+ if (locTrg->LoTrigY()==1 &&
+ locTrg->LoStripY()==15 ) yTrig = kFALSE; // no trigger in Y
+ else yTrig = kTRUE; // trigger in Y
- AliDebug(1, "AliMUONTrackReconstructor::MakeTriggerTrack using NEW trigger \n");
- AliMUONTriggerCircuit* circuit =
- (AliMUONTriggerCircuit*)fTriggerCircuit->At(locTrg->LoCircuit()-1); // -1 !!!
-
- y11 = circuit->GetY11Pos(locTrg->LoStripX());
- stripX21 = locTrg->LoStripX()+locTrg->LoDev()+1;
- y21 = circuit->GetY21Pos(stripX21);
- x11 = circuit->GetX11Pos(locTrg->LoStripY());
-
- AliDebug(1, Form(" MakeTriggerTrack %d %d %d %d %d %f %f %f \n",i,locTrg->LoCircuit(),
- locTrg->LoStripX(),locTrg->LoStripX()+locTrg->LoDev()+1,locTrg->LoStripY(),y11, y21, x11));
-
- Float_t thetax = TMath::ATan2( x11 , z11 );
- Float_t thetay = TMath::ATan2( (y21-y11) , (z21-z11) );
-
- fTriggerTrack->SetX11(x11);
- fTriggerTrack->SetY11(y11);
- fTriggerTrack->SetThetax(thetax);
- fTriggerTrack->SetThetay(thetay);
- fTriggerTrack->SetGTPattern(gloTrigPat);
- fTriggerTrack->SetLoTrgNum(i);
-
- fMUONData->AddRecTriggerTrack(*fTriggerTrack);
+ if (xTrig && yTrig) { // make Trigger Track if trigger in X and Y
+
+ AliDebug(1, "AliMUONTrackReconstructor::MakeTriggerTrack using NEW trigger \n");
+ AliMUONTriggerCircuit* circuit =
+ (AliMUONTriggerCircuit*)fTriggerCircuit->At(locTrg->LoCircuit()-1); // -1 !!!
+
+ y11 = circuit->GetY11Pos(locTrg->LoStripX());
+// need first to convert deviation to [0-30]
+// (see AliMUONLocalTriggerBoard::LocalTrigger)
+ Int_t deviation = locTrg->LoDev();
+ Int_t sign = 0;
+ if ( !locTrg->LoSdev() && deviation ) sign=-1;
+ if ( !locTrg->LoSdev() && !deviation ) sign= 0;
+ if ( locTrg->LoSdev() == 1 ) sign=+1;
+ deviation *= sign;
+ deviation += 15;
+ stripX21 = locTrg->LoStripX()+deviation+1;
+ y21 = circuit->GetY21Pos(stripX21);
+ x11 = circuit->GetX11Pos(locTrg->LoStripY());
+
+ AliDebug(1, Form(" MakeTriggerTrack %d %d %d %d %d %f %f %f \n",i,locTrg->LoCircuit(),
+ locTrg->LoStripX(),locTrg->LoStripX()+locTrg->LoDev()+1,locTrg->LoStripY(),y11, y21, x11));
+
+ Float_t thetax = TMath::ATan2( x11 , z11 );
+ Float_t thetay = TMath::ATan2( (y21-y11) , (z21-z11) );
+
+ fTriggerTrack->SetX11(x11);
+ fTriggerTrack->SetY11(y11);
+ fTriggerTrack->SetThetax(thetax);
+ fTriggerTrack->SetThetay(thetay);
+ fTriggerTrack->SetGTPattern(gloTrigPat);
+ fTriggerTrack->SetLoTrgNum(i);
+
+ fMUONData->AddRecTriggerTrack(*fTriggerTrack);
+ } // board is fired
} // end of loop on Local Trigger
return kTRUE;