add a check properly the bad/exotic cluster in the corresponding array, add few histo...

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

diff --git a/MUON/AliMUONResponseTriggerV1.cxx b/MUON/AliMUONResponseTriggerV1.cxx
index 28c3e1e0c68daddc19349ce3f52e04b8fbd0e28e..ef4c81e3046942864d3b789b521a9e4fae9e3c17 100644 (file)
--- a/MUON/AliMUONResponseTriggerV1.cxx
+++ b/MUON/AliMUONResponseTriggerV1.cxx
@@ -71,24 +71,15 @@ namespace
                                
 
 //------------------------------------------------------------------   
-AliMUONResponseTriggerV1::AliMUONResponseTriggerV1() : AliMUONResponseTrigger(), fGenerCluster(0), fWorkCondition(2)
-{  
-  fHVvalues.Set(72);
-  SetHV();
-  
-  fBValues.Set(144);
-  SetBValues();
+AliMUONResponseTriggerV1::AliMUONResponseTriggerV1() : AliMUONResponseTrigger(), fGenerCluster(0), fHVvalues(), fBValues(), fWorkCondition(2)
+{
 }
 
 //------------------------------------------------------------------   
-AliMUONResponseTriggerV1::AliMUONResponseTriggerV1(Int_t mode) : AliMUONResponseTrigger(), fGenerCluster(0)
+AliMUONResponseTriggerV1::AliMUONResponseTriggerV1(Int_t mode) : AliMUONResponseTrigger(), fGenerCluster(0), fHVvalues(), fBValues(), fWorkCondition(mode)
 {
-  fWorkCondition = mode; // 1=streamer - 2=avalanche
-  
-  fHVvalues.Set(72);
+  //mode: 1=streamer - 2=avalanche
   SetHV();
-  
-  fBValues.Set(144);
   SetBValues();
 }
 
@@ -230,6 +221,10 @@ void AliMUONResponseTriggerV1::DisIntegrate(const AliMUONHit& hit, TList& digits
 //------------------------------------------------------------------
 void AliMUONResponseTriggerV1::SetHV()
 {
+  //
+  /// Set HV values from OCDB
+  //
+  fHVvalues.Set(72);
   TString side;
   Int_t newRPC=0,newPlane=0;
   
@@ -296,9 +291,12 @@ void AliMUONResponseTriggerV1::SetHV()
         voltage += (dcsValue->GetFloat())*(time2-time1);
       }
       
-      fHVvalues.AddAt(voltage/(timeend-timebegin)/1000,18*iPlane+iRPC); //voltage in kV, not in V
+      Double_t deltaTime = timeend - timebegin;
+      Double_t meanVoltage = ( deltaTime == 0. ) ? 0. : voltage/deltaTime/1000.;
+      fHVvalues.AddAt(meanVoltage,18*iPlane+iRPC); //voltage in kV, not in V
       
       voltage=0;
+      AliDebug(1,Form("HV value for MTR_%s_MT%d_RPC%d_HV.vEff = %g (kV)",side.Data(),newPlane,newRPC,meanVoltage));
     }
   }
 }
@@ -306,6 +304,12 @@ void AliMUONResponseTriggerV1::SetHV()
 //------------------------------------------------------------------  
 void AliMUONResponseTriggerV1::SetBValues()
 {
+  //
+  /// Set B values for cluster size function
+  //
+  
+  fBValues.Set(144);
+  
   Float_t bValues[2][4][18] =                                                                                                                        {{{1.97,2.47,2.47,2.47,2.97,2.97,2.47,2.47,1.97,2.22,1.97,2.47,1.97,2.97,2.97,2.47,2.47,1.97},  //MT11BP
     {2.22,2.22,1.97,2.47,2.97,2.97,1.97,2.47,1.97,1.97,1.97,2.47,1.97,2.97,2.97,1.97,1.97,1.97},  //MT12BP
     {2.22,2.22,2.47,2.47,2.97,2.97,2.47,2.47,2.22,1.97,1.97,2.47,1.97,2.97,2.97,1.97,1.97,1.97},  //MT21BP