The description of changes:

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

diff --git a/MUON/AliMUONRecoParam.cxx b/MUON/AliMUONRecoParam.cxx
index 4e99779ed710658b35add239d134e2c793143e0c..c98d54d3e9ad5b24d79e70767353e343d6a157d4 100644 (file)
--- a/MUON/AliMUONRecoParam.cxx
+++ b/MUON/AliMUONRecoParam.cxx
@@ -72,11 +72,19 @@ AliMUONRecoParam::AliMUONRecoParam()
   fRemoveConnectedTracksInSt12(kFALSE),
   fMaxTriggerTracks(0),
   fMaxTrackCandidates(0),
-  fSelectTrackOnSlope(kFALSE)
+  fSelectTrackOnSlope(kFALSE),
+  fMissingPadFractionLimit(0),
+  fFractionOfBuspatchOutsideOccupancyLimit(0)
 {
   /// Constructor
   
   SetNameTitle("Dummy","Dummy");
+  for (Int_t iCh = 0; iCh < 10; iCh++) {
+    fUseChamber[iCh] = kTRUE;
+    fDefaultNonBendingReso[iCh] = 0.;
+    fDefaultBendingReso[iCh] = 0.;
+  }
+  for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
   SetDefaultLimits();
 }
 
@@ -106,6 +114,7 @@ AliMUONRecoParam::GetCalibrationMode() const
   /// NOGAIN : only do pedestal subtraction
   /// GAIN : do pedestal subtraction, and apply gain correction, but with a
   ///        single capacitance value for all channels
+  /// INJECTIONGAIN : as GAIN, but with gain values taken as EMELEC factory values
   /// GAINCONSTANTCAPA : as GAIN, but with a channel-dependent capacitance value
   
   return fCalibrationMode.Data();
@@ -144,6 +153,18 @@ AliMUONRecoParam *AliMUONRecoParam::GetCosmicParam()
   return param;
 }
 
+//_____________________________________________________________________________
+AliMUONRecoParam *AliMUONRecoParam::GetCalibrationParam() 
+{
+  /// Return default (dummy) reconstruction parameters for calibration environment
+  
+  AliMUONRecoParam *param = new AliMUONRecoParam();
+  param->SetCalibrationParam();
+  
+  return param;
+}
+
+
 //_____________________________________________________________________________
 void AliMUONRecoParam::SetLowFluxParam() 
 {
@@ -151,8 +172,8 @@ void AliMUONRecoParam::SetLowFluxParam()
   
   SetNameTitle("Low Flux","Low Flux");
   SetEventSpecie(AliRecoParam::kLowMult);
-  fMinBendingMomentum = 1.;
-  fMaxBendingMomentum = 3000.;
+  fMinBendingMomentum = 0.8;
+  fMaxBendingMomentum = 1.e10;
   fMaxNonBendingSlope = 0.3;
   fMaxBendingSlope = 0.4;
   fSelectTrackOnSlope = kFALSE;
@@ -162,10 +183,10 @@ void AliMUONRecoParam::SetLowFluxParam()
   fMaxBendingDistanceToTrack = 1.;
   fSigmaCutForTracking = 6.;
   fSigmaCutForImprovement = 5.;
-  fSigmaCutForTrigger = 8.;
+  fSigmaCutForTrigger = 4.;
   fStripCutForTrigger = 1.;
   fMaxStripAreaForTrigger = 3.;
-  fMaxNormChi2MatchTrigger = 16.;
+  fMaxNormChi2MatchTrigger = fSigmaCutForTrigger * fSigmaCutForTrigger;
   fCombinedClusterTrackReco = kFALSE;
   fTrackAllTracks = kTRUE;
   fRecoverTracks = kTRUE;
@@ -173,7 +194,7 @@ void AliMUONRecoParam::SetLowFluxParam()
   fMakeMoreTrackCandidates = kFALSE;
   fComplementTracks = kTRUE;
   fImproveTracks = kTRUE;
-  fRemoveConnectedTracksInSt12 = kTRUE;
+  fRemoveConnectedTracksInSt12 = kFALSE;
   fUseSmoother = kTRUE;
   for (Int_t iCh = 0; iCh < 10; iCh++) {
     fUseChamber[iCh] = kTRUE;
@@ -183,8 +204,7 @@ void AliMUONRecoParam::SetLowFluxParam()
   for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
   fBypassSt45 = 0;
   fMaxTriggerTracks = 100;
-  fMaxTrackCandidates = 10000;
-  
+  fMaxTrackCandidates = 10000;  
 }
 
 //_____________________________________________________________________________
@@ -194,8 +214,8 @@ void AliMUONRecoParam::SetHighFluxParam()
   
   SetNameTitle("High Flux","High Flux");
   SetEventSpecie(AliRecoParam::kHighMult);
-  fMinBendingMomentum = 1.;
-  fMaxBendingMomentum = 3000.;
+  fMinBendingMomentum = 0.8;
+  fMaxBendingMomentum = 1.e10;
   fMaxNonBendingSlope = 0.3;
   fMaxBendingSlope = 0.4;
   fSelectTrackOnSlope = kFALSE;
@@ -205,10 +225,10 @@ void AliMUONRecoParam::SetHighFluxParam()
   fMaxBendingDistanceToTrack = 1.;
   fSigmaCutForTracking = 6.;
   fSigmaCutForImprovement = 5.;
-  fSigmaCutForTrigger = 8.;
+  fSigmaCutForTrigger = 4.;
   fStripCutForTrigger = 1.;
   fMaxStripAreaForTrigger = 3.;
-  fMaxNormChi2MatchTrigger = 16.;
+  fMaxNormChi2MatchTrigger = fSigmaCutForTrigger * fSigmaCutForTrigger;
   fCombinedClusterTrackReco = kFALSE;
   fTrackAllTracks = kTRUE;
   fRecoverTracks = kTRUE;
@@ -237,21 +257,21 @@ void AliMUONRecoParam::SetCosmicParam()
   
   SetNameTitle("Cosmic","Cosmic");
   SetEventSpecie(AliRecoParam::kCosmic);
-  fMinBendingMomentum = 1.;
-  fMaxBendingMomentum = 10000000.;
-  fMaxNonBendingSlope = 0.4;
-  fMaxBendingSlope = 0.5;
+  fMinBendingMomentum = 0.8;
+  fMaxBendingMomentum = 1.e10;
+  fMaxNonBendingSlope = 0.3;
+  fMaxBendingSlope = 0.4;
   fSelectTrackOnSlope = kTRUE;
-  fNonBendingVertexDispersion = 200.;
-  fBendingVertexDispersion = 200.;
+  fNonBendingVertexDispersion = 170.;
+  fBendingVertexDispersion = 170.;
   fMaxNonBendingDistanceToTrack = 1.;
   fMaxBendingDistanceToTrack = 1.;
   fSigmaCutForTracking = 7.;
-  fSigmaCutForImprovement = 7.;
-  fSigmaCutForTrigger = 8.;
+  fSigmaCutForImprovement = 6.;
+  fSigmaCutForTrigger = 4.;
   fStripCutForTrigger = 1.5;
   fMaxStripAreaForTrigger = 3.;
-  fMaxNormChi2MatchTrigger = 16.;
+  fMaxNormChi2MatchTrigger = fSigmaCutForTrigger * fSigmaCutForTrigger;
   fPercentOfFullClusterInESD = 100.;
   fCombinedClusterTrackReco = kFALSE;
   fTrackAllTracks = kTRUE;
@@ -268,13 +288,37 @@ void AliMUONRecoParam::SetCosmicParam()
     fDefaultNonBendingReso[iCh] = 0.4;
     fDefaultBendingReso[iCh] = 0.4;
   }
-  for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
+  fRequestStation[0] = kTRUE;
+  fRequestStation[1] = kTRUE;
+  fRequestStation[2] = kTRUE;
+  fRequestStation[3] = kTRUE;
+  fRequestStation[4] = kTRUE;
   fBypassSt45 = 0;
   fPadGoodnessMask = 0x400BE80; // Ped Mean is Zero | Ped Mean Too Low | Ped Mean Too High | Ped Sigma Too Low | Ped Sigma Too High | Ped is missing | HV is missing | manu occupancy too high
   fMaxTriggerTracks = 100;
   fMaxTrackCandidates = 10000;
-  SetManuOccupancyLimits(-1,0.5); // reject manu above occ=50%
+  SetPedMeanLimits(20, 700);
+  SetManuOccupancyLimits(-1.,0.01); // reject manu above occ=1%
+
+  SetBuspatchOccupancyLimits(-1,0.01);  
+  SetMissingPadFractionLimit(0.1); // 10 %   
+  SetFractionOfBuspatchOutsideOccupancyLimit(0.05); // 5 %   
+}
+
+
+//_____________________________________________________________________________
+void AliMUONRecoParam::SetCalibrationParam() 
+{
+  /// Set (dummy) reconstruction parameters for calibration environment
   
+  SetNameTitle("Calibration","Calibration");
+  SetEventSpecie(AliRecoParam::kCalib);
+
+  fPedMeanLimits[0] = 5000;
+  fPedMeanLimits[1] = 0;
+
+  fPadGoodnessMask = 0x8C00; // Pedestal is missing | is too low | too high
+
 }
 
 //_____________________________________________________________________________
@@ -443,9 +487,13 @@ void AliMUONRecoParam::Print(Option_t *option) const
   
   cout << "Occupancy limits are :" << endl;
   
-  cout << Form("%7.2f <= Manu occupancy < %7.2f",ManuOccupancyLowLimit(),ManuOccupancyHighLimit()) << endl;
-  cout << Form("%7.2f <= Buspatch occupancy < %7.2f",BuspatchOccupancyLowLimit(),BuspatchOccupancyHighLimit()) << endl;
-  cout << Form("%7.2f <= DE occupancy < %7.2f",DEOccupancyLowLimit(),DEOccupancyHighLimit()) << endl;
+  cout << Form("%e <= Manu occupancy < %7.2f",ManuOccupancyLowLimit(),ManuOccupancyHighLimit()) << endl;
+  cout << Form("%e <= Buspatch occupancy < %7.2f",BuspatchOccupancyLowLimit(),BuspatchOccupancyHighLimit()) << endl;
+  cout << Form("%e <= DE occupancy < %7.2f",DEOccupancyLowLimit(),DEOccupancyHighLimit()) << endl;
+  
+  cout << "'QAChecker' limits" << endl;  
+  cout << Form("MissingPadFractionLimit = %5.2f %%",MissingPadFractionLimit()*100.0) << endl;  
+  cout << Form("FractionOfBuspatchOutsideOccupancyLimit = %5.2f %%",FractionOfBuspatchOutsideOccupancyLimit()*100.0) << endl;
   
   cout << "chamber non bending resolution = |";
   for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultNonBendingReso[iCh]);
@@ -472,7 +520,7 @@ AliMUONRecoParam::SetDefaultLimits()
        fHVSt345Limits[0]=1500;
        fHVSt345Limits[1]=2000;
 
-       fPedMeanLimits[0] = 50;
+       fPedMeanLimits[0] = 20;
        fPedMeanLimits[1] = 1024;
        
        fPedSigmaLimits[0] = 0.6;
@@ -492,12 +540,16 @@ AliMUONRecoParam::SetDefaultLimits()
   fManuOccupancyLimits[0] = -1.0; 
   fManuOccupancyLimits[1] = 1.0;
 
-  fBuspatchOccupancyLimits[0] = -1.0; 
+  fBuspatchOccupancyLimits[0] = 1E-6; 
   fBuspatchOccupancyLimits[1] = 1.0;
 
   fDEOccupancyLimits[0] = -1.0; 
   fDEOccupancyLimits[1] = 1.0;
 
   fChargeSigmaCut = 4.0;
+  
+  fMissingPadFractionLimit = 0.1; // 10 % 
+  fFractionOfBuspatchOutsideOccupancyLimit = 0.05; // 5 % 
+
 }