fix in the method Gain
authorbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 17 Mar 2006 15:25:07 +0000 (15:25 +0000)
committerbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 17 Mar 2006 15:25:07 +0000 (15:25 +0000)
PMD/AliPMDDigitizer.cxx
PMD/AliPMDDigits2Recpoints.C
PMD/AliPMDRecpointRead.C
PMD/AliPMDtracker.cxx

index 12eaf51..5dd212d 100644 (file)
@@ -640,7 +640,7 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   TrackAssignment2Cell();
   ResetCell();
 
-  Float_t gain1,gain2;
+  Float_t gain1;
   Float_t adc;
   Float_t deltaE = 0.;
   Int_t detno = 0;
@@ -655,17 +655,16 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
                {
                  if (idet == 0)
                    {
-              gain1 =Gain(idet,ism,jrow,kcol);
-  //           cout<<"The gain Factor is"<<gain1<<endl;
+                     gain1 = Gain(idet,ism,jrow,kcol);
+
                      deltaE = fPRE[ism][jrow][kcol]*gain1;
                      trno   = fPRETrackNo[ism][jrow][kcol];
                      detno = 0;
                    }
                  else if (idet == 1)
                    {
-              gain2 =Gain(idet,ism,jrow,kcol);
-   //          cout<<"The Gain factor is"<<gain2<<endl;
-                     deltaE = fCPV[ism][jrow][kcol]*gain2;
+                     gain1 = Gain(idet,ism,jrow,kcol);
+                     deltaE = fCPV[ism][jrow][kcol]*gain1;
                      trno   = fCPVTrackNo[ism][jrow][kcol];
                      detno = 1;
                    }
@@ -683,7 +682,7 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   
   fPMDLoader->WriteDigits("OVERWRITE");
   ResetCellADC();
-
+  
 }
 //____________________________________________________________________________
 
@@ -1239,7 +1238,7 @@ Float_t AliPMDDigitizer::Gain(Int_t det, Int_t smn, Int_t row, Int_t col) const
 
   if(!fCalibData) {
     AliError("No calibration data loaded from CDB!!!");
-    return -1;
+    return 1;
   }
 
   Float_t GainFact;
index fe7d43a..24688df 100644 (file)
@@ -38,12 +38,12 @@ void AliPMDDigits2Recpoints(Int_t nevt=1)
 
   // Create the PMD Cluster Finder 
   AliPMDClusterFinder *clus = new AliPMDClusterFinder(fRunLoader);
-  clus->SetDebug(1);
-  if (itype == 0)
+
+  if (itype == 1)
     {
       clus->Load();
     }
-  else if (itype == 1)
+  else if (itype == 0)
     {
       clus->LoadClusters();
     }  
@@ -51,23 +51,23 @@ void AliPMDDigits2Recpoints(Int_t nevt=1)
 
   for (Int_t ievt = 0; ievt < nevt; ievt++)
     {
-      if (itype == 0)
+      if (itype == 1)
        {
          // from digits data
          clus->Digits2RecPoints(ievt);
        }
-      else if (itype == 1)
+      else if (itype == 0)
        {
          // from raw data
          AliRawReaderFile reader(ievt);
          clus->Digits2RecPoints(ievt, &reader);
        }
     }
-  if (itype == 0)
+  if (itype == 1)
     {
       clus->UnLoad();
     }
-  else if (itype == 1)
+  else if (itype == 0)
     {
       clus->UnLoadClusters();
     }
index 8740686..15992af 100644 (file)
@@ -54,8 +54,8 @@ Int_t AliPMDRecpointRead(Int_t nevent = 1)
        }
 
       AliPMDrecpoint1  *pmdrecpoint;
-      TBranch *branch = treeR->GetBranch("PMDRecpoint");
-      branch->SetAddress(&fRecpoints);
+      TBranch *branch1 = treeR->GetBranch("PMDRecpoint");
+      branch1->SetAddress(&fRecpoints);
       /**********************************************************************
        *    det   : Detector, 0: PRE & 1:CPV                                *
        *    smn   : Serial Module Number from which Super Module Number     *
@@ -75,13 +75,13 @@ Int_t AliPMDRecpointRead(Int_t nevent = 1)
       Int_t   ism, ium;
       Int_t   det,smn;
       Float_t xpos,ypos, xpad, ypad;
-      Float_t adc, ncell, rad;
+      Float_t adc, ncell, sigx, sigy;
       Float_t xx, yy;
-      Int_t   nmodules = treeR->GetEntries();
+      Int_t   nmodules = branch1->GetEntries();
       cout << " nmodules = " << nmodules << endl;
       for (Int_t imodule = 0; imodule < nmodules; imodule++)
        {
-         treeR->GetEntry(imodule); 
+         branch1->GetEntry(imodule); 
          Int_t nentries = fRecpoints->GetLast();
          for(Int_t ient = 0; ient < nentries+1; ient++)
            {
@@ -92,7 +92,8 @@ Int_t AliPMDRecpointRead(Int_t nevent = 1)
              ypos  = pmdrecpoint->GetClusY();
              adc   = pmdrecpoint->GetClusADC();
              ncell = pmdrecpoint->GetClusCells();
-             rad   = pmdrecpoint->GetClusRadius();
+             sigx  = pmdrecpoint->GetClusSigmaX();
+             sigy  = pmdrecpoint->GetClusSigmaY();
              //
              // Now change the xpos and ypos to its original values
              // for the unit modules which are earlier changed.
@@ -119,8 +120,8 @@ Int_t AliPMDRecpointRead(Int_t nevent = 1)
              // User has to plug in his analysis code here
              //
 
-             fprintf(fpw,"%d %d %d %d %f %f %f %f %f\n",
-                     det,smn,ism,ium,xpad,ypad,adc,ncell,rad);
+             fprintf(fpw,"%d %d %d %d %f %f %f %f %f %f\n",
+                     det,smn,ism,ium,xpad,ypad,adc,ncell,sigx,sigy);
              //
              // Plot the cluster centroid to see the PMD geometry
              // using the PMD Utility class
index 2539f45..4426afd 100644 (file)
@@ -146,14 +146,13 @@ void AliPMDtracker::Clusters2Tracks(AliESD *event)
   AliPMDDiscriminator *pmddiscriminator = new AliPMDEmpDiscriminator();
   pmddiscriminator->Discrimination(fPMDcontin,fPMDcontout);
 
-  const Float_t kzpos0 = 361.5;    // for PREshower plane BKN
-  const Float_t kzpos1 = 361.5;    // for CPV plane
+  const Float_t kzpos = 361.5;    // middle of the PMD
   Int_t   ism =0, ium=0;
   Int_t   det,smn;
   Float_t xpos,ypos;
   Float_t xpad = 0, ypad = 0;
   Float_t adc, ncell, rad;
-  Float_t xglobal, yglobal, zglobal;
+  Float_t xglobal, yglobal, zglobal = 0;
   Float_t pid;
 
 
@@ -217,11 +216,11 @@ void AliPMDtracker::Clusters2Tracks(AliESD *event)
 
       if (det == 0)
        {
-         zglobal = kzpos0 + 0.5; // to be found out
+         zglobal = kzpos + 1.6; // PREshower plane
        }
       else if (det == 1)
        {
-         zglobal = kzpos1 - 0.5; // to be found out BKN
+         zglobal = kzpos - 1.7; // CPV plane
        }