- default cuts are tuned
authorsgorbuno <sgorbuno@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 6 Nov 2009 16:42:16 +0000 (16:42 +0000)
committersgorbuno <sgorbuno@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 6 Nov 2009 16:42:16 +0000 (16:42 +0000)
- cosmetic

HLT/global/physics/AliHLTV0HistoComponent.cxx

index 4518e6f..c594d7b 100644 (file)
@@ -134,7 +134,7 @@ int AliHLTV0HistoComponent::DoInit( int argc, const char** argv )
   fAP->SetYTitle("p_{t}[GeV]");
   fAP->SetXTitle("(p^{+}_{L}-p^{-}_{L})/(p^{+}_{L}+p^{-}_{L})");
   fAP->SetStats(0);
-  fAP->SetOption("COLZ");
+  fAP->SetOption("CONT4 Z");
 
   fGammaXY = new TH2F("hGammaXY","HLT:  #gamma conversions",100,-100.,100.,100,-100.,100.);
   fGammaXY->SetMarkerStyle(8);
@@ -166,10 +166,10 @@ int AliHLTV0HistoComponent::DoInit( int argc, const char** argv )
   fGammaCuts[3] = 3.0;
   fGammaCuts[4] = 0.0;
   fGammaCuts[5] = 300.0;
-  fGammaCuts[6] = 100.;
+  fGammaCuts[6] = 3.5;
   fGammaCuts[7] = 0.05;
 
-  fAPCuts[0] = 0;
+  fAPCuts[0] = 60;
   fAPCuts[1] = 2.5;
   fAPCuts[2] = 3.5;
   fAPCuts[3] = 3.0;
@@ -185,7 +185,7 @@ int AliHLTV0HistoComponent::DoInit( int argc, const char** argv )
   fKsCuts[4] = 1.5;
   fKsCuts[5] = 50.0;
   fKsCuts[6] = 4.0;
-  fKsCuts[7] = 0.03;
+  fKsCuts[7] = 0.015;
 
   fLambdaCuts[0] = 60;   // [0] == 60  --- N clusters on each daughter track
   fLambdaCuts[1] = 3.0;  // [1] == 3.0 --- (daughter-primVtx)/sigma >= cut
@@ -286,10 +286,13 @@ int AliHLTV0HistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/,
        v0.GetMass(mass,error); 
        if( fGamma ) fGamma->Fill( mass );
 
-       if( TMath::Abs(mass)<=fGammaCuts[6]*error && TMath::Abs(mass)<=fGammaCuts[7] ){   
+       if( TMath::Abs(mass)<=fGammaCuts[6]*error || TMath::Abs(mass)<=fGammaCuts[7] ){   
          AliKFParticle gamma = v0;
          gamma.SetMassConstraint(0);
-         if( fGammaXY ) fGammaXY->Fill(gamma.GetX(), gamma.GetY());
+         if( fGammaXY
+             &&  t1->GetTPCNcls()>=60
+             && t2->GetTPCNcls()>=60
+             ) fGammaXY->Fill(gamma.GetX(), gamma.GetY());
          isGamma = 1;
          fNGammas++;
          vGammas.push_back( gamma );
@@ -355,7 +358,7 @@ int AliHLTV0HistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/,
        double mass, error;
        Ks.GetMass( mass, error);
        if( fKShort ) fKShort->Fill( mass );    
-       if( TMath::Abs( mass - kKsMass )<=fKsCuts[6]*error && TMath::Abs( mass - kKsMass )<=fKsCuts[7] ){  
+       if( TMath::Abs( mass - kKsMass )<=fKsCuts[6]*error || TMath::Abs( mass - kKsMass )<=fKsCuts[7] ){  
          isKs = 1;
          fNKShorts++;
        }
@@ -374,7 +377,7 @@ int AliHLTV0HistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/,
         && length >= fLambdaCuts[3]*sigmaLength
         && length >= fLambdaCuts[4]
         && r <= fLambdaCuts[5]
-        //&& TMath::Abs( ap )>.3
+        && TMath::Abs( ap )>.4
         ){
        
        AliKFParticle kP, kpi;
@@ -391,7 +394,7 @@ int AliHLTV0HistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/,
        double mass, error;
        lambda.GetMass( mass, error);
        if( fLambda ) fLambda->Fill( mass );
-       if( TMath::Abs( mass - kLambdaMass )<=fLambdaCuts[6]*error && TMath::Abs( mass - kLambdaMass )<=fLambdaCuts[7] ){
+       if( TMath::Abs( mass - kLambdaMass )<=fLambdaCuts[6]*error || TMath::Abs( mass - kLambdaMass )<=fLambdaCuts[7] ){
          fNLambdas++;
        }
       }