- changing pT intervals for inv mass histos

- changing phi from degrees to radians

diff --git a/HLT/global/physics/AliHLTCaloHistoInvMass.cxx b/HLT/global/physics/AliHLTCaloHistoInvMass.cxx
index e60209adbbc073634939944edb7156711eaf6873..44e326d42b113d21729d40ea0916feba7a34ab5c 100644 (file)
--- a/HLT/global/physics/AliHLTCaloHistoInvMass.cxx
+++ b/HLT/global/physics/AliHLTCaloHistoInvMass.cxx
@@ -46,31 +46,31 @@ AliHLTCaloHistoInvMass::AliHLTCaloHistoInvMass(TString det) :
   fHistTwoClusterInvMass4(NULL)
 {
   // See header file for documentation
-  fHistTwoClusterInvMass0 = new TH1F(Form("%s fHistTwoClusterInvMass0", det.Data()), Form("Invariant mass of two clusters in %s, 0.8 GeV < E < 1.2", det.Data()), 200, 0, 1);
+  fHistTwoClusterInvMass0 = new TH1F(Form("%s fHistTwoClusterInvMass0", det.Data()), Form("Invariant mass of two clusters in %s, 1.5 GeV < E < 2.0", det.Data()), 200, 0, 1);
   fHistTwoClusterInvMass0->GetXaxis()->SetTitle("m_{#gamma#gamma} GeV");
   fHistTwoClusterInvMass0->GetYaxis()->SetTitle("Counts");
   fHistTwoClusterInvMass0->SetMarkerStyle(21);
   fHistArray->AddLast(fHistTwoClusterInvMass0);
 
-  fHistTwoClusterInvMass1 = new TH1F(Form("%s fHistTwoClusterInvMass1", det.Data()), Form("Invariant mass of two clusters in %s, 1.2 GeV < E < 1.6", det.Data()), 200, 0, 1);
+  fHistTwoClusterInvMass1 = new TH1F(Form("%s fHistTwoClusterInvMass1", det.Data()), Form("Invariant mass of two clusters in %s, 2.0 GeV < E < 2.5", det.Data()), 200, 0, 1);
   fHistTwoClusterInvMass1->GetXaxis()->SetTitle("m_{#gamma#gamma} GeV");
   fHistTwoClusterInvMass1->GetYaxis()->SetTitle("Counts");
   fHistTwoClusterInvMass1->SetMarkerStyle(21);
   fHistArray->AddLast(fHistTwoClusterInvMass1);
 
-  fHistTwoClusterInvMass2 = new TH1F(Form("%s fHistTwoClusterInvMass2", det.Data()), Form("Invariant mass of two clusters in %s, 1.6 GeV < E < 2.0", det.Data()), 200, 0, 1);
+  fHistTwoClusterInvMass2 = new TH1F(Form("%s fHistTwoClusterInvMass2", det.Data()), Form("Invariant mass of two clusters in %s, 2.5 GeV < E < 3.0", det.Data()), 200, 0, 1);
   fHistTwoClusterInvMass2->GetXaxis()->SetTitle("m_{#gamma#gamma} GeV");
   fHistTwoClusterInvMass2->GetYaxis()->SetTitle("Counts");
   fHistTwoClusterInvMass2->SetMarkerStyle(21);
   fHistArray->AddLast(fHistTwoClusterInvMass2);
 
-  fHistTwoClusterInvMass3 = new TH1F(Form("%s fHistTwoClusterInvMass3", det.Data()), Form("Invariant mass of two clusters in %s, 2.0 GeV < E < 4.0", det.Data()), 200, 0, 1);
+  fHistTwoClusterInvMass3 = new TH1F(Form("%s fHistTwoClusterInvMass3", det.Data()), Form("Invariant mass of two clusters in %s, 3.0 GeV < E < 5.0", det.Data()), 200, 0, 1);
   fHistTwoClusterInvMass3->GetXaxis()->SetTitle("m_{#gamma#gamma} GeV");
   fHistTwoClusterInvMass3->GetYaxis()->SetTitle("Counts");
   fHistTwoClusterInvMass3->SetMarkerStyle(21);
   fHistArray->AddLast(fHistTwoClusterInvMass3);
 
-  fHistTwoClusterInvMass4 = new TH1F(Form("%s fHistTwoClusterInvMass4", det.Data()), Form("Invariant mass of two clusters in %s E > 4.0 GeV", det.Data()), 200, 0, 1);
+  fHistTwoClusterInvMass4 = new TH1F(Form("%s fHistTwoClusterInvMass4", det.Data()), Form("Invariant mass of two clusters in %s E > 5.0 GeV", det.Data()), 200, 0, 1);
   fHistTwoClusterInvMass4->GetXaxis()->SetTitle("m_{#gamma#gamma} GeV");
   fHistTwoClusterInvMass4->GetYaxis()->SetTitle("Counts");
   fHistTwoClusterInvMass4->SetMarkerStyle(21);
@@ -147,7 +147,7 @@ Int_t AliHLTCaloHistoInvMass::FillInvariantMassHistograms(Int_t nc, Float_t cPos
   for(Int_t ic = 0; ic<(nc-1); ic++) { 
      
     //BALLE hardcoded variable
-    if(cEnergy[ic] < 0.4)
+    if(cEnergy[ic] < 0.5)
       continue;
 
     //Get momentum vector
@@ -159,7 +159,7 @@ Int_t AliHLTCaloHistoInvMass::FillInvariantMassHistograms(Int_t nc, Float_t cPos
     for(Int_t jc = ic+1; jc<nc; jc++) { 
      
     //BALLE hardcoded variable
-      if(cEnergy[jc] < 0.4)
+      if(cEnergy[jc] < 0.5)
        continue;
 
       
@@ -176,35 +176,38 @@ Int_t AliHLTCaloHistoInvMass::FillInvariantMassHistograms(Int_t nc, Float_t cPos
       //Fill histograms
       
       Float_t sum = cEnergy[ic]+cEnergy[ic];
-      if(sum > 1.2)
+      if(sum > 1.5)
       {
-        if(sum > 1.6)
+        if(sum > 2.0)
         {
-           if(sum > 2.0)
+           if(sum > 2.5)
            {
-              if(sum > 4.0)
+              if(sum > 3.0)
               {
-                 fHistTwoClusterInvMass4->Fill(m);
+                 if(sum > 5.0)
+                 {
+                    fHistTwoClusterInvMass4->Fill(m);
+                 }
+                 else
+                 {
+                    fHistTwoClusterInvMass3->Fill(m);
+                 }     
               }
               else
               {
-                 fHistTwoClusterInvMass3->Fill(m);
+                 fHistTwoClusterInvMass2->Fill(m);
               }
            }
            else
            {
-              fHistTwoClusterInvMass2->Fill(m);
+              fHistTwoClusterInvMass1->Fill(m);
            }
         }
         else
         {
-           fHistTwoClusterInvMass1->Fill(m);
+           fHistTwoClusterInvMass0->Fill(m);
         }
       }
-      else
-      {
-        fHistTwoClusterInvMass0->Fill(m);
-      }
     }
   }