]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
Added a default customization function for most aspects of the standard cuts
authorpulvir <pulvir@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 31 Aug 2010 14:55:06 +0000 (14:55 +0000)
committerpulvir <pulvir@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 31 Aug 2010 14:55:06 +0000 (14:55 +0000)
PWG2/RESONANCES/AliRsnCutESD2010.cxx
PWG2/RESONANCES/AliRsnCutESD2010.h
PWG2/RESONANCES/AliRsnVAnalysisTaskME.cxx

index a6989fa9b0d85ea14ff2d59d8ad6b934eedb8802..1f3267c780f9aff25763d908ad3beba7cbdf3ddf 100644 (file)
@@ -120,16 +120,64 @@ AliRsnCutESD2010::AliRsnCutESD2010
 //
 // Main constructor.
 //
-
-  Initialize();
 }
 
 //_________________________________________________________________________________________________
-void AliRsnCutESD2010::Initialize()
+void AliRsnCutESD2010::InitializeToDefaults(Bool_t isSim)
 {
 //
 // Main constructor.
 //
+
+  // ----> set TPC range for PID and calibration
+  SetTPCrange(5.0, 3.0);
+  SetTPCpLimit(0.35);
+  
+  // ----> set ITS range for PID
+  SetITSband(4.0);
+  
+  // ----> set TPC calibration
+  if (isSim) SetTPCpar(2.15898 / 50.0, 1.75295E1, 3.40030E-9, 1.96178, 3.91720);
+  else       SetTPCpar(1.41543 / 50.0, 2.63394E1, 5.0411E-11, 2.12543, 4.88663);
+  
+  // ----> set standard quality cuts for TPC global tracks
+  fESDtrackCutsTPC.SetRequireTPCStandAlone(kTRUE); // require to have the projection at inner TPC wall
+  fESDtrackCutsTPC.SetMinNClustersTPC(70);
+  fESDtrackCutsTPC.SetMaxChi2PerClusterTPC(4);
+  fESDtrackCutsTPC.SetAcceptKinkDaughters(kFALSE);
+  fESDtrackCutsTPC.SetRequireTPCRefit(kTRUE);
+  fESDtrackCutsTPC.SetRequireITSRefit(kTRUE);
+  fESDtrackCutsTPC.SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
+  fESDtrackCutsTPC.SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); // DCA pt dependent: 7*(0.0050+0.0060/pt0.9)
+  fESDtrackCutsTPC.SetMaxDCAToVertexZ(1e6); // disabled
+  fESDtrackCutsTPC.SetDCAToVertex2D(kFALSE); // each DCA is checked separately
+  fESDtrackCutsTPC.SetRequireSigmaToVertex(kFALSE);
+  
+  // ----> set standard quality cuts for ITS standalone tracks
+  fESDtrackCutsITS.SetRequireITSStandAlone(kTRUE, kTRUE);
+  fESDtrackCutsITS.SetRequireITSRefit(kTRUE);
+  fESDtrackCutsITS.SetMinNClustersITS(4);
+  fESDtrackCutsITS.SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
+  fESDtrackCutsITS.SetMaxChi2PerClusterITS(1.);
+  fESDtrackCutsITS.SetMaxDCAToVertexXYPtDep("0.0595+0.0182/pt^1.55"); // DCA pt dependent
+  fESDtrackCutsITS.SetMaxDCAToVertexZ(1e6); // disabled
+  fESDtrackCutsITS.SetDCAToVertex2D(kFALSE); // each DCA is checked separately
+  
+  // ----> set the TOF calibration depending on type of input (sim/data)
+  SetTOFcorrectTExp(kTRUE);
+  SetTOFuseT0(kTRUE);
+  SetTOFresolution(100.0);
+  if (isSim)
+  {
+    SetTOFcalibrateESD(kFALSE);
+    SetTOFtuneMC(kTRUE);
+  }
+  else
+  {
+    //if (isPass2) SetTOFcalibrateESD(kFALSE); // potrebbe anche essere kFALSE
+    SetTOFcalibrateESD(kTRUE);
+    SetTOFtuneMC(kFALSE);
+  }
 }
 
 //_________________________________________________________________________________________________
@@ -271,12 +319,14 @@ Bool_t AliRsnCutESD2010::IsSelected(TObject *obj1, TObject* /*obj2*/)
         else
         {
           okTOF = kTRUE;
+          if (!fCheckTPC) okTOF = kFALSE;
           //cout << "TOF    -- not checked due to ZERO reference time" << endl;
         }
       }
       else
       {
         okTOF = kTRUE;
+        if (!fCheckTPC) okTOF = kFALSE;
         //cout << "TOF    -- not checked because TOF pid absent" << endl;
       }
     }
index de0411c7c15eb9c19c4fd95fad6cfa648074c599..b8fe6b90087cac33cafa3c54b13ab4d0e42e9bc3 100644 (file)
@@ -27,7 +27,7 @@ class AliRsnCutESD2010 : public AliRsnCut
     AliRsnCutESD2010(const char *name);
     virtual ~AliRsnCutESD2010() {;};
 
-    void             Initialize();
+    void             InitializeToDefaults(Bool_t isSim = kFALSE);
     AliESDtrackCuts* GetCutsTPC() {return &fESDtrackCutsTPC;}
     AliESDtrackCuts* GetCutsITS() {return &fESDtrackCutsITS;}
     virtual Bool_t   IsSelected(TObject *obj1, TObject *obj2 = 0x0);
index a010f21a38cb04169264b6157e33c94a35ac37ef..0a0d8ca2741376882488e0c71c54102632e87ffd 100644 (file)
@@ -174,7 +174,7 @@ Bool_t AliRsnVAnalysisTaskME::CheckAndPrintEvents()
 // or number of events is less or equal 1
 //
 
-  AliInfo(Form("Current Entry %d", Entry()));
+  AliInfo(Form("Current Entry %d", (Int_t)Entry()));
   Int_t nEvents = fInputHandler->GetBufferSize();
   if (nEvents <= 1) return kFALSE;
   fESDEvent = dynamic_cast<AliESDEvent*>(GetEvent(0));