Protection against division by zero.

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

diff --git a/MUON/AliMUONRecoParam.cxx b/MUON/AliMUONRecoParam.cxx
index c70d7a3a38328f1e7c6f535a7bc3be85faeceb7b..87d1a7024ab7f6de2c21b7d78f523a7fb05d3896 100644 (file)
--- a/MUON/AliMUONRecoParam.cxx
+++ b/MUON/AliMUONRecoParam.cxx
@@ -69,7 +69,10 @@ AliMUONRecoParam::AliMUONRecoParam()
   fBypassSt45(0),
   fPadGoodnessMask(0),
   fChargeSigmaCut(4.0),
-  fRemoveConnectedTracksInSt12(kFALSE)
+  fRemoveConnectedTracksInSt12(kFALSE),
+  fMaxTriggerTracks(0),
+  fMaxTrackCandidates(0),
+  fSelectTrackOnSlope(kFALSE)
 {
   /// Constructor
   
@@ -152,8 +155,9 @@ void AliMUONRecoParam::SetLowFluxParam()
   fMaxBendingMomentum = 3000.;
   fMaxNonBendingSlope = 0.3;
   fMaxBendingSlope = 0.4;
-  fNonBendingVertexDispersion = 10.;
-  fBendingVertexDispersion = 10.;
+  fSelectTrackOnSlope = kFALSE;
+  fNonBendingVertexDispersion = 70.;
+  fBendingVertexDispersion = 70.;
   fMaxNonBendingDistanceToTrack = 1.;
   fMaxBendingDistanceToTrack = 1.;
   fSigmaCutForTracking = 6.;
@@ -178,6 +182,8 @@ void AliMUONRecoParam::SetLowFluxParam()
   }
   for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
   fBypassSt45 = 0;
+  fMaxTriggerTracks = 100;
+  fMaxTrackCandidates = 10000;
   
 }
 
@@ -192,8 +198,9 @@ void AliMUONRecoParam::SetHighFluxParam()
   fMaxBendingMomentum = 3000.;
   fMaxNonBendingSlope = 0.3;
   fMaxBendingSlope = 0.4;
-  fNonBendingVertexDispersion = 10.;
-  fBendingVertexDispersion = 10.;
+  fSelectTrackOnSlope = kFALSE;
+  fNonBendingVertexDispersion = 70.;
+  fBendingVertexDispersion = 70.;
   fMaxNonBendingDistanceToTrack = 1.;
   fMaxBendingDistanceToTrack = 1.;
   fSigmaCutForTracking = 6.;
@@ -218,6 +225,8 @@ void AliMUONRecoParam::SetHighFluxParam()
   }
   for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
   fBypassSt45 = 0;
+  fMaxTriggerTracks = 100;
+  fMaxTrackCandidates = 10000;
   
 }
 
@@ -230,15 +239,18 @@ void AliMUONRecoParam::SetCosmicParam()
   SetEventSpecie(AliRecoParam::kCosmic);
   fMinBendingMomentum = 1.;
   fMaxBendingMomentum = 10000000.;
-  fMaxNonBendingSlope = 0.3;
-  fMaxBendingSlope = 0.4;
-  fNonBendingVertexDispersion = 10.;
-  fBendingVertexDispersion = 10.;
-  fMaxNonBendingDistanceToTrack = 10.;
-  fMaxBendingDistanceToTrack = 10.;
-  fSigmaCutForTracking = 20.;
-  fSigmaCutForImprovement = 20.;
+  fMaxNonBendingSlope = 0.4;
+  fMaxBendingSlope = 0.5;
+  fSelectTrackOnSlope = kTRUE;
+  fNonBendingVertexDispersion = 200.;
+  fBendingVertexDispersion = 200.;
+  fMaxNonBendingDistanceToTrack = 1.;
+  fMaxBendingDistanceToTrack = 1.;
+  fSigmaCutForTracking = 7.;
+  fSigmaCutForImprovement = 7.;
   fSigmaCutForTrigger = 8.;
+  fStripCutForTrigger = 1.5;
+  fMaxStripAreaForTrigger = 3.;
   fMaxNormChi2MatchTrigger = 16.;
   fPercentOfFullClusterInESD = 100.;
   fCombinedClusterTrackReco = kFALSE;
@@ -253,11 +265,14 @@ void AliMUONRecoParam::SetCosmicParam()
   fSaveFullClusterInESD = kTRUE;
   for (Int_t iCh = 0; iCh < 10; iCh++) {
     fUseChamber[iCh] = kTRUE;
-    fDefaultNonBendingReso[iCh] = 0.144;
-    fDefaultBendingReso[iCh] = 0.01;
+    fDefaultNonBendingReso[iCh] = 0.4;
+    fDefaultBendingReso[iCh] = 0.4;
   }
   for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
   fBypassSt45 = 0;
+  fPadGoodnessMask = 0xC00BE80;
+  fMaxTriggerTracks = 100;
+  fMaxTrackCandidates = 10000;
   
 }
 
@@ -314,16 +329,14 @@ 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("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
-  
-  cout<<Form("Maximum non bending slope = %5.2f",fMaxNonBendingSlope)<<endl;
+  cout<<"Selection of track candidates:"<<endl;
+  if (fSelectTrackOnSlope) cout<<Form("\t- Non-bending slope < %5.2f",fMaxNonBendingSlope)<<endl;
+  else cout<<"\t- Impact parameter < 3 * vertex dispersion in the non-bending direction"<<endl;
+  cout<<Form("\t- if B!=0: Bending momentum > %5.2f",fMinBendingMomentum)<<endl;
+  if (fSelectTrackOnSlope) cout<<Form("\t  if B==0: Bending slope < %5.2f",fMaxBendingSlope)<<endl;
+  else cout<<"\t  if B==0: Impact parameter < 3 * vertex dispersion in the bending direction"<<endl;
   
-  cout<<Form("Maximum bending slope (used only if B=0) = %5.2f",fMaxBendingSlope)<<endl;
-  
-  if (strstr(fTrackingMode,"ORIGINAL"))
-    cout<<Form("Vertex dispertion = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
-  else if (strstr(option,"FULL"))
-    cout<<Form("Vertex dispertion (used for original tracking only) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
+  cout<<Form("Vertex dispersion (used to estimate initial bending momentum resolution) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
   
   cout<<Form("Maximum distance to track = (%5.2f,%5.2f)",fMaxNonBendingDistanceToTrack,fMaxBendingDistanceToTrack)<<endl;
   
@@ -433,6 +446,8 @@ void AliMUONRecoParam::Print(Option_t *option) const
   cout << "chamber bending resolution = |";
   for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultBendingReso[iCh]);
   cout << endl;
+  cout<<Form("maximum number of trigger tracks above which the tracking is cancelled = %d",fMaxTriggerTracks)<<endl;
+  cout<<Form("maximum number of track candidates above which the tracking is abandonned = %d",fMaxTrackCandidates)<<endl;
   
   cout<<"\t-----------------------------------------------------"<<endl<<endl;
   
@@ -453,7 +468,7 @@ AliMUONRecoParam::SetDefaultLimits()
        fPedMeanLimits[0] = 50;
        fPedMeanLimits[1] = 1024;
        
-       fPedSigmaLimits[0] = 0.1;
+       fPedSigmaLimits[0] = 0.6;
        fPedSigmaLimits[1] = 100;
 
        fGainA1Limits[0] = 0.1;
@@ -465,7 +480,8 @@ AliMUONRecoParam::SetDefaultLimits()
        fGainThresLimits[0] = 0;
        fGainThresLimits[1] = 4095;
        
-       fPadGoodnessMask = 0x8080;
+       fPadGoodnessMask = 0x8008080; // Ped is missing | HV is missing | killed
+
   
   fChargeSigmaCut = 4.0;
 }