example macros for EVE analysis using new ListAnalyser (Ben)

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

diff --git a/MUON/AliMUONRecoParam.cxx b/MUON/AliMUONRecoParam.cxx
index f18654f20cd6542a371063283d4561c85b449063..5780d4fae0a11a948fd05854075835d8fcf56382 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();
 }
 
@@ -145,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() 
 {
@@ -152,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;
@@ -184,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;  
 }
 
 //_____________________________________________________________________________
@@ -195,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;
@@ -238,17 +257,17 @@ 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.;
+  fSigmaCutForImprovement = 6.;
   fSigmaCutForTrigger = 8.;
   fStripCutForTrigger = 1.5;
   fMaxStripAreaForTrigger = 3.;
@@ -269,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
+
 }
 
 //_____________________________________________________________________________
@@ -444,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]);
@@ -473,7 +520,7 @@ AliMUONRecoParam::SetDefaultLimits()
        fHVSt345Limits[0]=1500;
        fHVSt345Limits[1]=2000;
 
-       fPedMeanLimits[0] = 50;
+       fPedMeanLimits[0] = 20;
        fPedMeanLimits[1] = 1024;
        
        fPedSigmaLimits[0] = 0.6;
@@ -493,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 % 
+
 }