#include "AliMUONTracker.h"
-#include "AliMUONReconstructor.h"
+#include "AliCodeTimer.h"
+#include "AliESDEvent.h"
+#include "AliESDMuonTrack.h"
+#include "AliESDVertex.h"
+#include "AliLog.h"
+#include "AliMUONClusterStoreV2.h"
+#include "AliMUONESDInterface.h"
+#include "AliMUONLegacyClusterServer.h"
#include "AliMUONRecoParam.h"
+#include "AliMUONReconstructor.h"
#include "AliMUONTrack.h"
#include "AliMUONTrackExtrap.h"
#include "AliMUONTrackHitPattern.h"
#include "AliMUONTrackParam.h"
-#include "AliMUONVCluster.h"
-#include "AliMUONVClusterServer.h"
-#include "AliMUONVDigitStore.h"
#include "AliMUONTrackReconstructor.h"
#include "AliMUONTrackReconstructorK.h"
#include "AliMUONTrackStoreV1.h"
#include "AliMUONTriggerTrackStoreV1.h"
-#include "AliMUONClusterStoreV2.h"
+#include "AliMUONTriggerTrack.h"
+#include "AliMUONLocalTrigger.h"
+#include "AliMUONVCluster.h"
+#include "AliMUONVClusterServer.h"
+#include "AliMUONVDigitStore.h"
#include "AliMUONVTriggerStore.h"
-#include "AliMUONESDInterface.h"
-
-#include "AliESDEvent.h"
-#include "AliESDVertex.h"
-#include "AliESDMuonTrack.h"
-#include "AliLog.h"
-#include "AliCodeTimer.h"
-
#include <Riostream.h>
-#include <TTree.h>
#include <TRandom.h>
+#include <TTree.h>
/// \cond CLASSIMP
ClassImp(AliMUONTracker)
//_____________________________________________________________________________
-AliMUONTracker::AliMUONTracker(AliMUONVClusterServer& clusterServer,
- const AliMUONVDigitStore& digitStore,
+AliMUONTracker::AliMUONTracker(AliMUONVClusterServer* clusterServer,
+ const AliMUONVDigitStore& digitStore,
const AliMUONDigitMaker* digitMaker,
const AliMUONGeometryTransformer* transformer,
const AliMUONTriggerCircuit* triggerCircuit)
fTrackReco(0x0),
fClusterStore(0x0),
fTriggerStore(0x0),
- fClusterServer(clusterServer), // not owner
- fDigitStore(digitStore) // not owner
+ fClusterServer(clusterServer),
+ fIsOwnerOfClusterServer(kFALSE),
+ fDigitStore(digitStore), // not owner
+ fInputClusterStore(0x0),
+ fTriggerTrackStore(0x0)
{
/// constructor
if (fTransformer && fDigitMaker)
fTrackHitPatternMaker = new AliMUONTrackHitPattern(*fTransformer,*fDigitMaker);
- TIter next(fDigitStore.CreateIterator());
- fClusterServer.UseDigits(next);
+ if (!fClusterServer)
+ {
+ AliInfo("No cluster server given. Will use AliMUONLegacyClusterServer");
+ fIsOwnerOfClusterServer = kTRUE;
+ }
+ else
+ {
+ TIter next(fDigitStore.CreateIterator());
+ fClusterServer->UseDigits(next);
+
+ SetupClusterServer(*fClusterServer);
+ }
}
//_____________________________________________________________________________
delete fTrackHitPatternMaker;
delete fClusterStore;
delete fTriggerStore;
+ if ( fIsOwnerOfClusterServer ) delete fClusterServer;
+ delete fInputClusterStore;
+ delete fTriggerTrackStore;
}
//_____________________________________________________________________________
return fClusterStore;
}
+//_____________________________________________________________________________
+AliMUONVTriggerTrackStore*
+AliMUONTracker::TriggerTrackStore() const
+{
+ /// Return (and create if necessary) the trigger track container
+ if (!fTriggerTrackStore)
+ {
+ fTriggerTrackStore = new AliMUONTriggerTrackStoreV1;
+ }
+ return fTriggerTrackStore;
+}
+
//_____________________________________________________________________________
Int_t AliMUONTracker::LoadClusters(TTree* clustersTree)
{
/// Load triggerStore from clustersTree
- ClusterStore()->Clear();
-
delete fTriggerStore;
+ delete fInputClusterStore;
+ fInputClusterStore=0x0;
if ( ! clustersTree ) {
AliFatal("No clustersTree");
return 2;
}
- ClusterStore()->Connect(*clustersTree,kFALSE);
+ if ( fIsOwnerOfClusterServer )
+ {
+ fInputClusterStore = AliMUONVClusterStore::Create(*clustersTree);
+ if ( fInputClusterStore )
+ {
+ AliInfo(Form("Created %s from cluster tree",fInputClusterStore->ClassName()));
+ fInputClusterStore->Clear();
+ fInputClusterStore->Connect(*clustersTree,kFALSE);
+ }
+ delete fClusterServer;
+ fClusterServer = new AliMUONLegacyClusterServer(*fTransformer,fInputClusterStore,
+ AliMUONReconstructor::GetRecoParam()->BypassSt4(),
+ AliMUONReconstructor::GetRecoParam()->BypassSt5());
+ SetupClusterServer(*fClusterServer);
+ }
+
fTriggerStore->Connect(*clustersTree,kFALSE);
clustersTree->GetEvent(0);
AliDebug(1,"");
AliCodeTimerAuto("")
- if (!fTrackReco) CreateTrackReconstructor();
+ if (!fTrackReco)
+ {
+ fTrackReco = CreateTrackReconstructor(AliMUONReconstructor::GetRecoParam()->GetTrackingMode(),fClusterServer);
+ }
// if the required tracking mode does not exist
if (!fTrackReco) return 1;
return 3;
}
+ // Make trigger tracks
+ if ( fTriggerCircuit )
+ {
+ TriggerTrackStore()->Clear();
+ fTrackReco->EventReconstructTrigger(*fTriggerCircuit,*fTriggerStore,*(TriggerTrackStore()));
+ }
+
+ if ( ( AliMUONReconstructor::GetRecoParam()->BypassSt4() ||
+ AliMUONReconstructor::GetRecoParam()->BypassSt5() ) &&
+ TriggerTrackStore()->GetSize() > 5 )
+ {
+ // Hard cut to reject shower events
+
+ AliCodeTimerAuto("MUON Shower events");
+
+ AliWarning(Form("Probably got a shower event (%d trigger tracks). Will not reconstruct tracks.",
+ TriggerTrackStore()->GetSize()));
+
+ return 0;
+ }
+
// Make tracker tracks
AliMUONVTrackStore* trackStore = new AliMUONTrackStoreV1;
fTrackReco->EventReconstruct(*(ClusterStore()),*trackStore);
- // Make trigger tracks
- AliMUONVTriggerTrackStore* triggerTrackStore(0x0);
- if ( fTriggerCircuit ) {
- triggerTrackStore = new AliMUONTriggerTrackStoreV1;
- fTrackReco->EventReconstructTrigger(*fTriggerCircuit,*fTriggerStore,*triggerTrackStore);
- }
-
// Match tracker/trigger tracks
- if ( triggerTrackStore && fTrackHitPatternMaker ) {
- fTrackReco->ValidateTracksWithTrigger(*trackStore,*triggerTrackStore,*fTriggerStore,*fTrackHitPatternMaker);
+ if ( fTrackHitPatternMaker )
+ {
+ fTrackReco->ValidateTracksWithTrigger(*trackStore,*(TriggerTrackStore()),*fTriggerStore,*fTrackHitPatternMaker);
}
// Fill ESD
// cleanup
delete trackStore;
- delete triggerTrackStore;
return 0;
}
AliDebug(1,"");
AliCodeTimerAuto("")
- AliMUONTrack* track;
- AliESDMuonTrack esdTrack;
- Double_t vertex[3] = {0., 0., 0.};
- TIter next(trackStore.CreateIterator());
-
// get ITS vertex
+ Double_t vertex[3] = {0., 0., 0.};
const AliESDVertex* esdVert = esd->GetVertex();
if (esdVert->GetNContributors()) {
esdVert->GetXYZ(vertex);
}
// fill ESD event including all info in ESD cluster if required and only for the given fraction of events
+ AliMUONTrack* track;
+ AliMUONLocalTrigger* locTrg;
+ AliESDMuonTrack esdTrack;
+ TIter next(trackStore.CreateIterator());
if (AliMUONReconstructor::GetRecoParam()->SaveFullClusterInESD() &&
gRandom->Uniform(100.) <= AliMUONReconstructor::GetRecoParam()->GetPercentOfFullClusterInESD()) {
while ( ( track = static_cast<AliMUONTrack*>(next()) ) ) {
- AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex, &fDigitStore);
+
+ if (track->GetMatchTrigger() > 0) {
+ locTrg = static_cast<AliMUONLocalTrigger*>(fTriggerStore->FindLocal(track->LoCircuit()));
+ AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex, &fDigitStore, locTrg);
+ } else AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex, &fDigitStore);
+
esd->AddMuonTrack(&esdTrack);
}
} else {
while ( ( track = static_cast<AliMUONTrack*>(next()) ) ) {
- AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex);
+
+ if (track->GetMatchTrigger() > 0) {
+ locTrg = static_cast<AliMUONLocalTrigger*>(fTriggerStore->FindLocal(track->LoCircuit()));
+ AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex, 0x0, locTrg);
+ } else AliMUONESDInterface::MUONToESD(*track, esdTrack, vertex);
+
esd->AddMuonTrack(&esdTrack);
}
}
+ // fill the local trigger decisions not matched with tracks (associate them to "ghost" tracks)
+ UInt_t ghostId = 0xFFFFFFFF - 1;
+ Bool_t matched = kFALSE;
+ AliMUONTriggerTrack *triggerTrack;
+ TIter itTriggerTrack(fTriggerTrackStore->CreateIterator());
+ while ( ( triggerTrack = static_cast<AliMUONTriggerTrack*>(itTriggerTrack()) ) ) {
+
+ locTrg = static_cast<AliMUONLocalTrigger*>(fTriggerStore->FindLocal(triggerTrack->GetLoTrgNum()));
+
+ // check if this local trigger has already been matched
+ TIter itTrack(trackStore.CreateIterator());
+ while ( ( track = static_cast<AliMUONTrack*>(itTrack()) ) ) {
+ matched = (track->LoCircuit() == locTrg->LoCircuit());
+ if (matched) break;
+ }
+ if (matched) continue;
+
+ AliMUONESDInterface::MUONToESD(*locTrg, esdTrack, ghostId);
+
+ esd->AddMuonTrack(&esdTrack);
+ ghostId -= 1;
+ }
+
}
//_____________________________________________________________________________
-void AliMUONTracker::CreateTrackReconstructor()
+AliMUONVTrackReconstructor* AliMUONTracker::CreateTrackReconstructor(const char* trackingMode, AliMUONVClusterServer* clusterServer)
{
/// Create track reconstructor, depending on tracking mode set in RecoParam
- TString opt(AliMUONReconstructor::GetRecoParam()->GetTrackingMode());
+ AliMUONVTrackReconstructor* trackReco(0x0);
+
+ TString opt(trackingMode);
opt.ToUpper();
if (strstr(opt,"ORIGINAL"))
{
- fTrackReco = new AliMUONTrackReconstructor(fClusterServer);
+ trackReco = new AliMUONTrackReconstructor(clusterServer);
}
else if (strstr(opt,"KALMAN"))
{
- fTrackReco = new AliMUONTrackReconstructorK(fClusterServer);
+ trackReco = new AliMUONTrackReconstructorK(clusterServer);
}
else
{
- AliError(Form("tracking mode \"%s\" does not exist",opt.Data()));
- return;
+ AliErrorClass(Form("tracking mode \"%s\" does not exist",opt.Data()));
+ return 0x0;
}
- AliInfo(Form("Will use %s for tracking",fTrackReco->ClassName()));
+ AliInfoClass(Form("Will use %s for tracking",trackReco->ClassName()));
+
+ return trackReco;
}
//_____________________________________________________________________________
{
/// Clear internal clusterStore
- ClusterStore()->Clear();
+ delete fInputClusterStore;
+ fInputClusterStore = 0x0;
+}
+
+
+//_____________________________________________________________________________
+void
+AliMUONTracker::SetupClusterServer(AliMUONVClusterServer& clusterServer)
+{
+ /// Setup the cluster server
+
+ if ( AliMUONReconstructor::GetRecoParam()->BypassSt4() ||
+ AliMUONReconstructor::GetRecoParam()->BypassSt5() )
+ {
+ Bool_t ok = clusterServer.UseTriggerTrackStore(TriggerTrackStore());
+
+ TString msg1;
+ TString msg2;
+
+ if ( AliMUONReconstructor::GetRecoParam()->BypassSt45() )
+ {
+ msg1 = "STATIONS 4 AND 5";
+ msg2 = "THOSE TWO STATIONS";
+ }
+ else if ( AliMUONReconstructor::GetRecoParam()->BypassSt4() )
+ {
+ msg1 = "STATION 4";
+ msg2 = "THAT STATION";
+ }
+ else if ( AliMUONReconstructor::GetRecoParam()->BypassSt5() )
+ {
+ msg1 = "STATION 5";
+ msg2 = "THAT STATION";
+ }
+
+ if ( ok )
+ {
+ AliWarning(Form("WILL USE TRIGGER TRACKS TO GENERATE CLUSTERS IN %s, "
+ "THUS BYPASSING REAL CLUSTERS IN %s!!!",msg1.Data(),msg2.Data()));
+ }
+ else
+ {
+ AliWarning("BYPASSING OF ST4 AND/OR 5 REQUESTED, BUT CLUSTERSERVER DOES NOT SEEM TO SUPPORT IT !!!");
+ }
+ }
}
+
+