@@ -82,12 +82,12 @@ void generateTPCCorrectionNTuple(const char* path = "InputSCDensityHistograms.ro
8282 const o2 ::gpu ::TPCFastTransformGeo & geo = fastTransform -> getGeometry ();
8383
8484 TFile * f = new TFile ("tpcCorrection.root" , "RECREATE" );
85- TNtuple * nt = new TNtuple ("dist" , "dist" , "slice :row:su:sv:dx:du:dv" );
85+ TNtuple * nt = new TNtuple ("dist" , "dist" , "sector :row:su:sv:dx:du:dv" );
8686
87- int32_t nSlices = 1 ; // fastTransform->getNumberOfSlices ();
88- // for( int32_t slice =0; slice<nSlices; slice ++){
89- for (int32_t slice = 0 ; slice < 1 ; slice ++ ) {
90- const o2 ::gpu ::TPCFastTransformGeo ::SliceInfo & sliceInfo = geo .getSliceInfo ( slice );
87+ int32_t nSectors = 1 ; // fastTransform->getNumberOfSectors ();
88+ // for( int32_t sector =0; sector<nSectors; sector ++){
89+ for (int32_t sector = 0 ; sector < 1 ; sector ++ ) {
90+ const o2 ::gpu ::TPCFastTransformGeo ::SectorInfo & sectorInfo = geo .getSectorInfo ( sector );
9191
9292 for (int32_t row = 0 ; row < geo .getNumberOfRows (); row ++ ) {
9393
@@ -96,12 +96,12 @@ void generateTPCCorrectionNTuple(const char* path = "InputSCDensityHistograms.ro
9696 for (float su = 0. ; su <= 1. ; su += 0.01 ) {
9797 for (float sv = 0. ; sv <= 1. ; sv += 0.01 ) {
9898 float u , v , y = 0 , z = 0 ;
99- geo .convScaledUVtoUV (slice , row , su , sv , u , v );
100- geo .convUVtoLocal (slice , u , v , y , z );
99+ geo .convScaledUVtoUV (sector , row , su , sv , u , v );
100+ geo .convUVtoLocal (sector , u , v , y , z );
101101
102102 // local 2 global
103103 float gx , gy , gz ;
104- geo .convLocalToGlobal (slice , x , y , z , gx , gy , gz );
104+ geo .convLocalToGlobal (sector , x , y , z , gx , gy , gz );
105105
106106 o2 ::tpc ::GlobalPosition3D positionCorrected (gx , gy , gz );
107107 sc -> correctElectron (positionCorrected );
@@ -111,15 +111,15 @@ void generateTPCCorrectionNTuple(const char* path = "InputSCDensityHistograms.ro
111111
112112 // global to local
113113 float x1 , y1 , z1 ;
114- geo .convGlobalToLocal (slice , gx , gy , gz , x1 , y1 , z1 );
114+ geo .convGlobalToLocal (sector , gx , gy , gz , x1 , y1 , z1 );
115115 float u1 = 0 , v1 = 0 ;
116- geo .convLocalToUV (slice , y1 , z1 , u1 , v1 );
116+ geo .convLocalToUV (sector , y1 , z1 , u1 , v1 );
117117
118118 float dx = x1 - x ;
119119 float du = u1 - u ;
120120 float dv = v1 - v ;
121- std ::cout << slice << " " << row << " " << su << " " << sv << " " << dx << " " << du << " " << dv << std ::endl ;
122- nt -> Fill (slice , row , su , sv , dx , du , dv );
121+ std ::cout << sector << " " << row << " " << su << " " << sv << " " << dx << " " << du << " " << dv << std ::endl ;
122+ nt -> Fill (sector , row , su , sv , dx , du , dv );
123123 }
124124 }
125125 }
0 commit comments