From bcf83487279470a9a9ed15e8f75ac6b5a708abce Mon Sep 17 00:00:00 2001 From: cblume Date: Mon, 28 Feb 2000 19:06:16 +0000 Subject: [PATCH] Cleanup --- TRD/digitsAna.C | 115 --------------------------------------------- TRD/digitsCreate.C | 57 ---------------------- 2 files changed, 172 deletions(-) delete mode 100644 TRD/digitsAna.C delete mode 100644 TRD/digitsCreate.C diff --git a/TRD/digitsAna.C b/TRD/digitsAna.C deleted file mode 100644 index a31cb4aa1e7..00000000000 --- a/TRD/digitsAna.C +++ /dev/null @@ -1,115 +0,0 @@ -void digitsAna () { - -///////////////////////////////////////////////////////////////////////// -// -// Example macro for the analysis of the TRD digits and the use -// of the AliTRDmatrix class. -// -///////////////////////////////////////////////////////////////////////// - - // Dynamically link some shared libs - if (gClassTable->GetID("AliRun") < 0) { - gROOT->LoadMacro("loadlibs.C"); - loadlibs(); - } - - // Input file name - Char_t *alifile = "galice_v1.root"; - - // Event number - Int_t nEvent = 0; - - // Define the objects - AliTRDv1 *TRD; - TClonesArray *TRDDigits; - AliTRDdigit *OneTRDDigit; - - // Connect the AliRoot file containing Geometry, Kine, Hits, and Digits - TFile *gafl = (TFile*) gROOT->GetListOfFiles()->FindObject(alifile); - if (!gafl) { - cout << "Open the ALIROOT-file " << alifile << endl; - gafl = new TFile(alifile); - } - else { - cout << alifile << " is already open" << endl; - } - - // Get AliRun object from file or create it if not on file - if (!gAlice) { - gAlice = (AliRun*) gafl->Get("gAlice"); - if (gAlice) - cout << "AliRun object found on file" << endl; - else - gAlice = new AliRun("gAlice","Alice test program"); - } - - // Import the Trees for the event nEvent in the file - Int_t nparticles = gAlice->GetEvent(nEvent); - if (nparticles <= 0) break; - - // Get the pointer to the hit-tree - TTree *DigitsTree = gAlice->TreeD(); - - // Get the pointer to the detector classes - TRD = (AliTRDv1*) gAlice->GetDetector("TRD"); - // Get the pointer to the hit container - if (TRD) TRDDigits = TRD->Digits(); - - // Define the detector matrix for one chamber (Sector 6, Chamber 3, Plane 1) - const Int_t iSec = 6; - const Int_t iCha = 3; - const Int_t iPla = 1; - Int_t rowMax = TRD->GetRowMax(iPla,iCha); - Int_t colMax = TRD->GetColMax(iPla); - Int_t timeMax = TRD->GetTimeMax(); - cout << " rowMax = " << rowMax - << " colMax = " << colMax - << " timeMax = " << timeMax << endl; - AliTRDmatrix *TRDMatrix = new AliTRDmatrix(rowMax,colMax,timeMax,iSec,iCha,iPla); - - Int_t nEntries = DigitsTree->GetEntries(); - cout << "Number of entries in digits tree = " << nEntries << endl; - - // Loop through all entries in the tree - Int_t nbytes; - for (Int_t iEntry = 0; iEntry < nEntries; iEntry++) { - - cout << "iEntry = " << iEntry << endl; - - // Import the tree - gAlice->ResetDigits(); - nbytes += DigitsTree->GetEvent(iEntry); - - // Get the number of digits in the detector - Int_t nTRDDigits = TRDDigits->GetEntriesFast(); - cout << " nTRDDigits = " << nTRDDigits << endl; - - // Loop through all TRD digits - for (Int_t iTRDDigits = 0; iTRDDigits < nTRDDigits; iTRDDigits++) { - - // Get the information for this digit - OneTRDDigit = (AliTRDdigit*) TRDDigits->UncheckedAt(iTRDDigits); - Int_t signal = OneTRDDigit->fSignal; - Int_t sector = OneTRDDigit->fSector; - Int_t chamber = OneTRDDigit->fChamber; - Int_t plane = OneTRDDigit->fPlane; - Int_t row = OneTRDDigit->fRow; - Int_t col = OneTRDDigit->fCol; - Int_t time = OneTRDDigit->fTime; - - // Fill the detector matrix - if (signal > 1) { - TRDMatrix->SetSignal(row,col,time,signal); - } - - } - - } - - // Display the detector matrix - TRDMatrix->Draw(); - TRDMatrix->DrawRow(18); - TRDMatrix->DrawCol(58); - TRDMatrix->DrawTime(20); - -} diff --git a/TRD/digitsCreate.C b/TRD/digitsCreate.C deleted file mode 100644 index d9c283d8dc8..00000000000 --- a/TRD/digitsCreate.C +++ /dev/null @@ -1,57 +0,0 @@ -void digitsCreate () { - -///////////////////////////////////////////////////////////////////////// -// -// Creates the digits from the hit information. An additional hit-tree -// is added to the input file. -// -///////////////////////////////////////////////////////////////////////// - - // Dynamically link some shared libs - if (gClassTable->GetID("AliRun") < 0) { - gROOT->LoadMacro("loadlibs.C"); - loadlibs(); - } - - // Input (and output) file name - Char_t *alifile = "galice_v1.root"; - - // Event number - Int_t nEvent = 0; - - // Connect the AliRoot file containing Geometry, Kine, and Hits - TFile *gafl = (TFile*) gROOT->GetListOfFiles()->FindObject(alifile); - if (!gafl) { - cout << "Open the ALIROOT-file " << alifile << endl; - gafl = new TFile(alifile,"UPDATE"); - } - else { - cout << alifile << " is already open" << endl; - } - - // Get AliRun object from file or create it if not on file - if (!gAlice) { - gAlice = (AliRun*) gafl->Get("gAlice"); - if (gAlice) - cout << "AliRun object found on file" << endl; - else - gAlice = new AliRun("gAlice","Alice test program"); - } - - // Import the Trees for the event nEvent in the file - Int_t nparticles = gAlice->GetEvent(nEvent); - if (nparticles <= 0) break; - - // Get the pointer to the detector class - AliTRDv2 *TRD = (AliTRDv2*) gAlice->GetDetector("TRD"); - - // Create the digitd and fill the digits-tree - TRD->Hits2Digits(); - - // Write the new tree into the input file - cout << "Entries in hit tree = " << gAlice->TreeD()->GetEntries()) << endl; - Char_t treeName[7]; - sprintf(treeName,"TreeD%d",nEvent); - gAlice->TreeD()->Write(treeName); - -} -- 2.39.3