- Added protection to prevent crash when the option "BypassSt45()"

[u/mrichter/AliRoot.git] / MUON / AliMUONRecoParam.cxx

diff --git a/MUON/AliMUONRecoParam.cxx b/MUON/AliMUONRecoParam.cxx
index 616022b9a6baa056d4e60e18f197a0e35274cade..6bfe46ee60a7bcd888e5dfc9e7b161fed09815b4 100644 (file)
--- a/MUON/AliMUONRecoParam.cxx
+++ b/MUON/AliMUONRecoParam.cxx
@@ -38,6 +38,7 @@ AliMUONRecoParam::AliMUONRecoParam()
 : AliDetectorRecoParam(),
   fClusteringMode("MLEM"),
   fTrackingMode("KALMAN"),
+  fMostProbBendingMomentum(0.),
   fMinBendingMomentum(0.),
   fMaxBendingMomentum(0.),
   fMaxNonBendingSlope(0.),
@@ -110,11 +111,23 @@ AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam()
   return param;
 }
 
+//_____________________________________________________________________________
+AliMUONRecoParam *AliMUONRecoParam::GetCosmicParam() 
+{
+  /// Return default reconstruction parameters for high flux environment
+  
+  AliMUONRecoParam *param = new AliMUONRecoParam();
+  param->SetCosmicParam();
+  
+  return param;
+}
+
 //_____________________________________________________________________________
 void AliMUONRecoParam::SetLowFluxParam() 
 {
   /// Set reconstruction parameters for low flux environment
   
+  fMostProbBendingMomentum = 2.;
   fMinBendingMomentum = 1.;
   fMaxBendingMomentum = 3000.;
   fMaxNonBendingSlope = 0.3;
@@ -145,6 +158,7 @@ void AliMUONRecoParam::SetHighFluxParam()
 {
   /// Set reconstruction parameters for high flux environment
   
+  fMostProbBendingMomentum = 2.;
   fMinBendingMomentum = 1.;
   fMaxBendingMomentum = 3000.;
   fMaxNonBendingSlope = 0.3;
@@ -170,6 +184,39 @@ void AliMUONRecoParam::SetHighFluxParam()
   
 }
 
+//_____________________________________________________________________________
+void AliMUONRecoParam::SetCosmicParam() 
+{
+  /// Set reconstruction parameters for high flux environment
+  
+  fMostProbBendingMomentum = 2.;
+  fMinBendingMomentum = 1.;
+  fMaxBendingMomentum = 10000000.;
+  fMaxNonBendingSlope = 0.3;
+  fNonBendingVertexDispersion = 10.;
+  fBendingVertexDispersion = 10.;
+  fMaxNonBendingDistanceToTrack = 10.;
+  fMaxBendingDistanceToTrack = 10.;
+  fSigmaCutForTracking = 20.;
+  fSigmaCutForImprovement = 20.;
+  fSigmaCutForTrigger = 8.;
+  fMaxNormChi2MatchTrigger = 16.;
+  fPercentOfFullClusterInESD = 100.;
+  fCombinedClusterTrackReco = kFALSE;
+  fTrackAllTracks = kTRUE;
+  fRecoverTracks = kTRUE;
+  fMakeTrackCandidatesFast = kFALSE;
+  fMakeMoreTrackCandidates = kFALSE;
+  fComplementTracks = kTRUE;
+  fImproveTracks = kTRUE;
+  fUseSmoother = kTRUE;
+  fSaveFullClusterInESD = kTRUE;
+  for (Int_t iCh = 0; iCh < 10; iCh++) fUseChamber[iCh] = kTRUE;
+  for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
+  fBypassSt45 = kFALSE;
+  
+}
+
 //_____________________________________________________________________________
 void AliMUONRecoParam::Print(Option_t *option) const
 {
@@ -190,6 +237,8 @@ void AliMUONRecoParam::Print(Option_t *option) const
   if (fSaveFullClusterInESD) cout<<Form("Save all cluster info in ESD for %5.2f %% of events",fPercentOfFullClusterInESD)<<endl;
   else cout<<"Save partial cluster info in ESD"<<endl;
   
+  cout<<Form("Most probable bending momentum (used only if B=0) = %5.2f",fMostProbBendingMomentum)<<endl;
+  
   cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
   
   cout<<Form("Maximum non bending slope = %5.2f",fMaxNonBendingSlope)<<endl;
@@ -278,14 +327,6 @@ void AliMUONRecoParam::Print(Option_t *option) const
   } while (++st < 5);
   if (discardedSt) cout<<endl;
   
-  if (strstr(option,"FULL")) {
-    cout<<"Bypass stations 4 and 5 (use trigger tracks as primary track candidates): ";
-    if (fBypassSt45) cout<<"ON"<<endl;
-    else cout<<"OFF"<<endl;
-  } else if (fBypassSt45)
-    cout<<"Bypass stations 4 and 5 (use trigger tracks as primary track candidates)"<<endl;
-  
-  
   cout<<"\t-------------------------------------"<<endl<<endl;
   
 }