1 // ----------------------------------------------------//
3 // This macro reads the PMD clusters which //
4 // are stored in the file "PMD.RecPoints.root" //
6 // ----------------------------------------------------//
10 #include "TStopwatch.h"
12 extern AliRun *gAlice;
14 Int_t AliPMDRecpointRead(Int_t nevent = 1)
18 delete AliRunLoader::Instance();
19 delete gAlice;//if everything was OK here it is already NULL
22 AliRunLoader *fRunLoader = AliRunLoader::Open("galice.root","Event","update");
25 cerr<<"Can't load RunLoader"<<endl;
28 AliLoader *pmdloader = fRunLoader->GetLoader("PMDLoader");
29 Int_t nevent = fRunLoader->GetNumberOfEvents();
30 cout << " * *********** nevent = " << nevent << endl;
34 cerr<<" ===> Can not find PMD or PMDLoader <===\n";
39 pmdloader->LoadRecPoints("READ");
42 AliCDBManager *man = AliCDBManager::Instance();
43 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
46 TClonesArray *fRecpoints;
47 AliPMDUtility *cc = new AliPMDUtility();
51 TH2F *h2 = new TH2F("h2"," ",100,-100.,100.,100,-100.,100.);
53 FILE *fpw = fopen("junk_rec.dat","w");
55 for (Int_t ievt = 0; ievt < nevent; ievt++)
57 fRunLoader->GetEvent(ievt);
58 TTree *treeR = pmdloader->TreeR();
61 cout << " Can not get TreeR" << endl;
65 AliPMDrecpoint1 *pmdrecpoint;
66 TBranch *branch1 = treeR->GetBranch("PMDRecpoint");
67 branch1->SetAddress(&fRecpoints);
68 /**********************************************************************
69 * det : Detector, 0: PRE & 1:CPV *
70 * smn : Serial Module Number from 0 to 23 for both detector *
71 * xpos : x-position of the cluster *
72 * ypos : y-position of the cluster *
73 * THESE xpos & ypos are not the true xpos and ypos *
74 * for some of the unit modules. They are rotated. *
75 * adc : ADC contained in the cluster *
76 * ncell : Number of cells contained in the cluster *
77 * rad : radius of the cluster (1d fit) *
78 * xpad : TRUE x-position of the cluster *
79 * ypad : TRUE y-position of the cluster *
80 **********************************************************************/
83 Float_t xpos,ypos, xpad, ypad;
84 Float_t adc, ncell, sigx, sigy;
86 Int_t nmodules = branch1->GetEntries();
87 cout << " nmodules = " << nmodules << endl;
88 for (Int_t imodule = 0; imodule < nmodules; imodule++)
90 branch1->GetEntry(imodule);
91 Int_t nentries = fRecpoints->GetLast();
92 for(Int_t ient = 0; ient < nentries+1; ient++)
94 pmdrecpoint = (AliPMDrecpoint1*)fRecpoints->UncheckedAt(ient);
95 det = (Int_t) pmdrecpoint->GetDetector();
96 smn = (Int_t) pmdrecpoint->GetSMNumber();
97 xpos = pmdrecpoint->GetClusX();
98 ypos = pmdrecpoint->GetClusY();
99 adc = pmdrecpoint->GetClusADC();
100 ncell = pmdrecpoint->GetClusCells();
101 sigx = pmdrecpoint->GetClusSigmaX();
102 sigy = pmdrecpoint->GetClusSigmaY();
105 // User has to plug in his analysis code here
108 fprintf(fpw,"%d %d %d %d\n",
111 // Plot the cluster centroid to see the PMD geometry
112 // using the PMD Utility class
116 // Draw only for PRE plane
117 cc->RectGeomCellPos(smn,xpos,ypos,xx,yy);
122 // End of the User code