New production requests (Nicole)

[u/mrichter/AliRoot.git] / prod / LHC09a16 / CheckESD.C

diff --git a/prod/LHC09a16/CheckESD.C b/prod/LHC09a16/CheckESD.C
new file mode 100644 (file)
index 0000000..fed6bc7
--- /dev/null
+++ b/prod/LHC09a16/CheckESD.C
@@ -0,0 +1,49 @@
+Bool_t CheckESD(const char* esdFileName = "AliESDs.root")\r
+{\r
+\r
+  AliCDBManager *cdb = AliCDBManager::Instance();\r
+  cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");\r
+\r
+  // open the ESD file\r
+  TFile* esdFile = TFile::Open(esdFileName);\r
+  if (!esdFile || !esdFile->IsOpen()) {\r
+    Error("CheckESD", "opening ESD file %s failed", esdFileName);\r
+    return kFALSE;\r
+  }\r
+  AliESDEvent * esd = new AliESDEvent;\r
+  TTree* tree = (TTree*) esdFile->Get("esdTree");\r
+  if (!tree) {\r
+    Error("CheckESD", "no ESD tree found");\r
+    return kFALSE;\r
+  }\r
+  esd->ReadFromTree(tree);\r
+\r
+  // loop over events\r
+  for (Int_t iEvent = 0; iEvent < tree->GetEntries(); iEvent++) {\r
+\r
+    // get the event summary data\r
+    tree->GetEvent(iEvent);\r
+    if (!esd) {\r
+      Error("CheckESD", "no ESD object found for event %d", iEvent);\r
+      return kFALSE;\r
+    }\r
+\r
+    Int_t nTracks = esd->GetNumberOfMuonTracks();\r
+    for(Int_t iTracks = 0; iTracks < nTracks; iTracks++) {\r
+      AliESDMuonTrack* muonTrack = esd->GetMuonTrack(iTracks);\r
+      if (muonTrack->ContainTrackerData()) {\r
+      }\r
+    }\r
+\r
+  }\r
+\r
+  delete esd;\r
+  esdFile->Close();\r
+  delete esdFile;\r
+\r
+  // result of check\r
+  Info("CheckESD", "check of ESD was successfull");\r
+  return kTRUE;\r
+\r
+}\r
+\r