/[PAMELA software]/calo/flight/CaloEnergy/src/CaloEnergy.cpp
ViewVC logotype

Diff of /calo/flight/CaloEnergy/src/CaloEnergy.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.9 by mocchiut, Thu Aug 13 15:56:57 2009 UTC revision 1.19 by mocchiut, Thu Dec 3 10:43:19 2009 UTC
# Line 35  CaloEnergy::CaloEnergy(PamLevel2 *l2p, B Line 35  CaloEnergy::CaloEnergy(PamLevel2 *l2p, B
35    //    //
36  }  }
37    
38    void CaloEnergy::Delete(){
39      Clear();
40      delete this;
41    }
42    
43    void CaloEnergy::UseLevel2(){
44      if ( clong ){
45        delete clong;
46        clong = new CaloLong(L2);
47        clong->SplitInto(0,22);
48        //    clong->HeavyTail(true);
49        clong->HeavyTail(false);
50      };
51      if ( cp ) delete cp;
52      cp = NULL;
53    }
54    
55  void CaloEnergy::UseCaloPreSampler(){  void CaloEnergy::UseCaloPreSampler(){
56    //    //
57    // use the presampler setting forcefitmode to 1000 means to force the DV routine to find the track inside the calorimeter using the "shower" approach developed for electrons    // use the presampler setting forcefitmode to 1000 means to force the DV routine to find the track inside the calorimeter using the "shower" approach developed for electrons
# Line 46  void CaloEnergy::UseCaloPreSampler(){ Line 63  void CaloEnergy::UseCaloPreSampler(){
63  //   cp->ForceCaloFit();  //   cp->ForceCaloFit();
64  //   cp->SetDebug(true);  //   cp->SetDebug(true);
65  //   cp->Process();  //   cp->Process();
66    if ( clong ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer());    if ( clong ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer());
67  }  }
68    
69    
# Line 55  void CaloEnergy::UseLongFit(){ Line 72  void CaloEnergy::UseLongFit(){
72    fLong = true;    fLong = true;
73    if ( !clong ){    if ( !clong ){
74      clong = new CaloLong(L2);      clong = new CaloLong(L2);
75      if ( cp ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer());      if ( cp ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer());
76      clong->SplitInto(0,22);      clong->SplitInto(0,22);
77        //    clong->HeavyTail(true);
78        clong->HeavyTail(false);
79    };    };
80    //    //
81  }  }
# Line 79  void CaloEnergy::Set(){ Line 98  void CaloEnergy::Set(){
98    //    //
99    indep = false;    indep = false;
100    //    //
101      fAllpl = true;
102    fLong = false;    fLong = false;
103    fPl = 1;    fPl = 1;
104    fRad = -1;    fRad = -1;
# Line 86  void CaloEnergy::Set(){ Line 106  void CaloEnergy::Set(){
106    clong = NULL;    clong = NULL;
107    x0max = -1.;    x0max = -1.;
108    //    //
109      multicol = false;
110      //
111    this->DefineGeometry();    this->DefineGeometry();
112    fXosel =true;    fXosel =true;
113    fXesel = true;    fXesel = true;
# Line 99  void CaloEnergy::Set(){ Line 121  void CaloEnergy::Set(){
121    fXemin = 1000;    fXemin = 1000;
122    fYomin = 1000;    fYomin = 1000;
123    fYemin = 1000;    fYemin = 1000;
124      //
125      this->UseCaloPreSampler(); // use it by default, to go back to "standard" mode use CaloEnergy::UseLevel2().
126      //
127  }  }
   
 void CaloEnergy::SetMinimumContainment(TString section, Int_t plane){  
   section.ToUpper();  
   if ( section.Contains("XO") ) fXomin = plane;  
   if ( section.Contains("XE") ) fXemin = plane;  
   if ( section.Contains("YO") ) fYomin = plane;  
   if ( section.Contains("YE") ) fYemin = plane;  
 }  
   
 void CaloEnergy::SetMinimumContainment(Int_t plane){  
   this->SetMinimumContainment("XEXOYEYO",plane);  
 }  
   
 Int_t CaloEnergy::GetMinimumContainment(TString section){  
   section.ToUpper();  
   if ( section.Contains("XO") ) return(fXomin);  
   if ( section.Contains("XE") ) return(fXemin);  
   if ( section.Contains("YE") ) return(fYemin);  
   if ( section.Contains("YO") ) return(fYomin);  
   printf(" ERROR: section not recognized \n");  
   return(-1000);  
 }  
   
 void CaloEnergy::SetConversionFactor(TString section, Float_t conv){  
   section.ToUpper();  
   if ( section.Contains("XO") ) fConv_rxo = conv;  
   if ( section.Contains("XE") ) fConv_rxe = conv;  
   if ( section.Contains("YO") ) fConv_ryo = conv;  
   if ( section.Contains("YE") ) fConv_rye = conv;  
 }  
   
 void CaloEnergy::SetConversionFactor(Float_t conv){  
   this->SetConversionFactor("XEXOYEYO",conv);  
 }  
   
 Float_t CaloEnergy::GetConversionFactor(TString section){  
   section.ToUpper();  
   if ( section.Contains("XO") ) return(fConv_rxo);  
   if ( section.Contains("XE") ) return(fConv_rxe);  
   if ( section.Contains("YO") ) return(fConv_ryo);  
   if ( section.Contains("YE") ) return(fConv_rye);  
   printf(" ERROR: section not recognized \n");  
   return(-1000.);  
 }  
   
 Int_t CaloEnergy::GetMaxplane(TString section){  
   section.ToUpper();  
   if ( section.Contains("XO") ) return fMax_planexo;  
   if ( section.Contains("XE") ) return fMax_planexe;  
   if ( section.Contains("YO") ) return fMax_planeyo;  
   if ( section.Contains("YE") ) return fMax_planeye;  
   return(-1);  
 }  
   
 Float_t CaloEnergy::GetEnergyAtMaxplane(TString section){  
   section.ToUpper();  
   if ( section.Contains("XO") ) return xomax_en;  
   if ( section.Contains("XE") ) return xemax_en;  
   if ( section.Contains("YO") ) return yomax_en;  
   if ( section.Contains("YE") ) return yemax_en;  
   return(-1);  
 }  
   
 Float_t CaloEnergy::GetMaxEnergy(TString section){  
   section.ToUpper();  
   if ( fLong ){  
     this->Process(section);      
     return fXOen_maxplane;  
   } else {  
     if ( section.Contains("XO") ) return fXOen_maxplane;  
     if ( section.Contains("XE") ) return fXEen_maxplane;  
     if ( section.Contains("YO") ) return fYOen_maxplane;  
     if ( section.Contains("YE") ) return fYEen_maxplane;  
   };  
   return(-1);  
 }  
   
 Float_t CaloEnergy::GetMaxEnergy(){    
   if ( fLong ){  
     if ( debug ) printf(" oh! call process! with asntr %s and sntr %s \n",asntr.Data(),sntr.Data());  
     this->Process(asntr);  
   };  
   return((fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane));  
 }  
   
128  void CaloEnergy::DefineGeometry(){  void CaloEnergy::DefineGeometry(){
129    //    //
130    // Use CaloStrip to determine once the position of border strips for each section    // Use CaloStrip to determine once the position of border strips for each section
131    //    //
132    //    //  fM = 2. + 0.096; // real position from cbar BUG the 0.096 is already taken into account in the border calculation made by Giovanna
133    fM = 2. + 0.096; // real position from cbar    fM = 2. ; // real position from cbar
134    fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc.    //  fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc. BUG! this way we count from the silicon border not from the silicon sensitive area
135      fM1 = 2. - 0.122 - 0.096 + 0.096; // due to calculation of xe1 etc.
136    if ( fM1 < 0. ) fM1 = 0.;    if ( fM1 < 0. ) fM1 = 0.;
137    //    //
138    CaloStrip *cs = new CaloStrip(fSimu);    CaloStrip *cs = new CaloStrip(fSimu);
# Line 199  void CaloEnergy::DefineGeometry(){ Line 140  void CaloEnergy::DefineGeometry(){
140    // view y plane 0 strip 0    // view y plane 0 strip 0
141    cs->Set(1,0,0);    cs->Set(1,0,0);
142    xe1= cs->GetY();    xe1= cs->GetY();
   //  z1= cs->GetZ();  
143    // view y plane 0 strip 31    // view y plane 0 strip 31
144    cs->Set(1,0,31);    cs->Set(1,0,31);
145    xe2= cs->GetY();    xe2= cs->GetY();
# Line 215  void CaloEnergy::DefineGeometry(){ Line 155  void CaloEnergy::DefineGeometry(){
155    // view y plane 0 strip 95    // view y plane 0 strip 95
156    cs->Set(1,0,95);    cs->Set(1,0,95);
157    xe6= cs->GetY();    xe6= cs->GetY();
   
158    // view x plane 0 strip 0    // view x plane 0 strip 0
159    cs->Set(0,0,0);    cs->Set(0,0,0);
160    yo1= cs->GetX();    yo1= cs->GetX();
   //  z2= cs->GetZ();  
161    // view x plane 0 strip 31    // view x plane 0 strip 31
162    cs->Set(0,0,31);    cs->Set(0,0,31);
163    yo2= cs->GetX();    yo2= cs->GetX();
# Line 235  void CaloEnergy::DefineGeometry(){ Line 173  void CaloEnergy::DefineGeometry(){
173    // view x plane 0 strip 95    // view x plane 0 strip 95
174    cs->Set(0,0,95);    cs->Set(0,0,95);
175    yo6= cs->GetX();    yo6= cs->GetX();
   
     
176    // view y plane 1 strip 0    // view y plane 1 strip 0
177    cs->Set(1,1,0);    cs->Set(1,1,0);
178    xo1= cs->GetY();    xo1= cs->GetY();
   //  z3= cs->GetZ();  
179    // view y plane 1 strip 31    // view y plane 1 strip 31
180    cs->Set(1,1,31);    cs->Set(1,1,31);
181    xo2= cs->GetY();    xo2= cs->GetY();
# Line 255  void CaloEnergy::DefineGeometry(){ Line 190  void CaloEnergy::DefineGeometry(){
190    xo5= cs->GetY();    xo5= cs->GetY();
191    // view y plane 1 strip 95    // view y plane 1 strip 95
192    cs->Set(1,1,95);    cs->Set(1,1,95);
193    xo6= cs->GetY();    xo6= cs->GetY();  
     
194    // view x plane 1 strip 0    // view x plane 1 strip 0
195    cs->Set(0,1,0);    cs->Set(0,1,0);
196    ye1= cs->GetX();    ye1= cs->GetX();
   //  z4= cs->GetZ();  
197    // view x plane 1 strip 31    // view x plane 1 strip 31
198    cs->Set(0,1,31);    cs->Set(0,1,31);
199    ye2= cs->GetX();    ye2= cs->GetX();
# Line 276  void CaloEnergy::DefineGeometry(){ Line 209  void CaloEnergy::DefineGeometry(){
209    // view x plane 1 strip 95    // view x plane 1 strip 95
210    cs->Set(0,1,95);    cs->Set(0,1,95);
211    ye6= cs->GetX();    ye6= cs->GetX();
   
212    //    //
213    for (Int_t p = 0; p<22; p ++){    for (Int_t p = 0; p<22; p ++){
214      for (Int_t v = 0; v<2; v++ ){        for (Int_t v = 0; v<2; v++ ){  
# Line 315  void CaloEnergy::Clear(){ Line 247  void CaloEnergy::Clear(){
247    yomax_en= 0.;    yomax_en= 0.;
248    yemax_en= 0.;    yemax_en= 0.;
249    //    //
250   memset(enstrip,0,2*22*96*(sizeof(Float_t)));      memset(enstrip,0,2*22*96*(sizeof(Float_t)));  
251    en = 0.;    en = 0.;
252    view = 0;    view = 0;
253    plane = 0;    plane = 0;
# Line 337  void CaloEnergy::Clear(){ Line 269  void CaloEnergy::Clear(){
269    memset(en_xop,0,11*sizeof(Float_t));    memset(en_xop,0,11*sizeof(Float_t));
270    memset(en_yop,0,11*sizeof(Float_t));    memset(en_yop,0,11*sizeof(Float_t));
271    //    //
272      fColumn = -1;
273      fColXE = -1;
274      fColXO = -1;
275      fColYE = -1;
276      fColYO = -1;
277      memset(encol,0,2*3*sizeof(Float_t));
278      entot[0] = 0.;
279      entot[1] = 0.;
280      //
281      X0pl = 0.76;
282      //
283  }  }
284    
285  void CaloEnergy::Print(){  void CaloEnergy::Print(){
# Line 369  void CaloEnergy::Print(){ Line 312  void CaloEnergy::Print(){
312    printf(" fM1     :.............. %f \n",fM1);    printf(" fM1     :.............. %f \n",fM1);
313    printf(" fRad    :.............. %i \n",fRad);    printf(" fRad    :.............. %i \n",fRad);
314    printf(" fPl     :.............. %i \n",fPl);    printf(" fPl     :.............. %i \n",fPl);
315      printf(" fColumn :.............. %i \n",fColumn);
316      printf(" multicol:.............. %i \n",multicol);
317      printf(" encol x :.............. %f \n",this->GetEncol(0));
318      printf(" encol y :.............. %f \n",this->GetEncol(1));
319      printf(" entot x :.............. %f \n",this->GetEntot(0));
320      printf(" entot y :.............. %f \n",this->GetEntot(1));
321      printf(" fColXE  :.............. %i \n",fColXE);
322      printf(" fColXO  :.............. %i \n",fColXO);
323      printf(" fColYE  :.............. %i \n",fColYE);
324      printf(" fColYO  :.............. %i \n",fColYO);
325    printf(" fConv_rxe ............. %f \n",fConv_rxe);    printf(" fConv_rxe ............. %f \n",fConv_rxe);
326    printf(" fConv_rxo ............. %f \n",fConv_rxo);    printf(" fConv_rxo ............. %f \n",fConv_rxo);
327    printf(" fConv_ryo ............. %f \n",fConv_ryo);    printf(" fConv_ryo ............. %f \n",fConv_ryo);
# Line 383  void CaloEnergy::Print(){ Line 336  void CaloEnergy::Print(){
336    printf(" fYEen_maxplane:........ %f \n",fYEen_maxplane);    printf(" fYEen_maxplane:........ %f \n",fYEen_maxplane);
337    printf(" fYOen_maxplane:........ %f \n",fYOen_maxplane);    printf(" fYOen_maxplane:........ %f \n",fYOen_maxplane);
338    printf(" x0max   :.............. %f \n",x0max);    printf(" x0max   :.............. %f \n",x0max);
339      printf(" X0pl    :.............. %f \n",X0pl);
340    printf(" debug   :.............. %i \n",debug);    printf(" debug   :.............. %i \n",debug);
341    
342    printf("========================================================================\n");    printf("========================================================================\n");
343    //    //
344  }  }
345    
346  void CaloEnergy::Delete(){  void CaloEnergy::SetMinimumContainment(TString section, Int_t plane){
347    Clear();    section.ToUpper();
348    delete this;    if ( section.Contains("XO") ) fXomin = plane;
349      if ( section.Contains("XE") ) fXemin = plane;
350      if ( section.Contains("YO") ) fYomin = plane;
351      if ( section.Contains("YE") ) fYemin = plane;
352    }
353    
354    void CaloEnergy::SetMinimumContainment(Int_t plane){
355      this->SetMinimumContainment("XEXOYEYO",plane);
356    }
357    
358    void CaloEnergy::SetConversionFactor(TString section, Float_t conv){
359      section.ToUpper();
360      if ( section.Contains("XO") ) fConv_rxo = conv;
361      if ( section.Contains("XE") ) fConv_rxe = conv;
362      if ( section.Contains("YO") ) fConv_ryo = conv;
363      if ( section.Contains("YE") ) fConv_rye = conv;
364    }
365    
366    void CaloEnergy::SetConversionFactor(Float_t conv){
367      this->SetConversionFactor("XEXOYEYO",conv);
368    }
369    
370    Int_t CaloEnergy::GetMinimumContainment(TString section){
371      section.ToUpper();
372      if ( section.Contains("XO") ) return(fXomin);
373      if ( section.Contains("XE") ) return(fXemin);
374      if ( section.Contains("YE") ) return(fYemin);
375      if ( section.Contains("YO") ) return(fYomin);
376      printf(" ERROR: section not recognized \n");
377      return(-1000);
378    }
379    
380    Float_t CaloEnergy::GetConversionFactor(TString section){
381      section.ToUpper();
382      if ( section.Contains("XO") ) return(fConv_rxo);
383      if ( section.Contains("XE") ) return(fConv_rxe);
384      if ( section.Contains("YO") ) return(fConv_ryo);
385      if ( section.Contains("YE") ) return(fConv_rye);
386      printf(" ERROR: section not recognized \n");
387      return(-1000.);
388    }
389    
390    Int_t CaloEnergy::GetMaxplane(TString section){
391      section.ToUpper();
392      if ( section.Contains("XO") ) return fMax_planexo;
393      if ( section.Contains("XE") ) return fMax_planexe;
394      if ( section.Contains("YO") ) return fMax_planeyo;
395      if ( section.Contains("YE") ) return fMax_planeye;
396      return(-1);
397    }
398    
399    Int_t CaloEnergy::GetColumn(TString section){
400      section.ToUpper();
401      if ( section.Contains("XO") ) return fColXO;
402      if ( section.Contains("XE") ) return fColXE;
403      if ( section.Contains("YO") ) return fColYO;
404      if ( section.Contains("YE") ) return fColYE;
405      return(-1);
406    }
407    
408    Float_t CaloEnergy::GetMipEnergyAtMaxplane(TString section){
409      printf(" WARNING: OBSOLETE METHOD, use GetMipEnergyAtMaxplane(TString) instead! \n");
410      return (this->GetEnergyAtMaxplane(section));
411    }
412    
413    Float_t CaloEnergy::GetEnergyAtMaxplane(TString section){
414      section.ToUpper();
415      if ( section.Contains("XO") ) return xomax_en;
416      if ( section.Contains("XE") ) return xemax_en;
417      if ( section.Contains("YO") ) return yomax_en;
418      if ( section.Contains("YE") ) return yemax_en;
419      return(-1);
420    }
421    
422    Float_t CaloEnergy::GetMaxEnergy(TString section){
423      printf(" WARNING: OBSOLETE METHOD, use GetMipEnergy(TString) instead! \n");
424      return (this->GetMipEnergy(section));
425    }
426    
427    Float_t CaloEnergy::GetMipEnergy(TString section){
428      section.ToUpper();
429      if ( fLong ){
430        this->Process(section);    
431        return fXOen_maxplane;
432      } else {
433        if ( section.Contains("XO") ) return fXOen_maxplane;
434        if ( section.Contains("XE") ) return fXEen_maxplane;
435        if ( section.Contains("YO") ) return fYOen_maxplane;
436        if ( section.Contains("YE") ) return fYEen_maxplane;
437      };
438      return(-1);
439    }
440    
441    Float_t CaloEnergy::GetEncol(Int_t i){  
442      if ( fColumn > -1 && (((fXesel || fXosel)&&i==1) || ((fYesel || fYosel)&&i==0)) ){
443        Int_t t = -1;
444        if ( i == 0 ){
445          if ( fColumn == 0 || fColumn == 3 || fColumn == 6 ) t = 0;
446          if ( fColumn == 1 || fColumn == 4 || fColumn == 7 ) t = 1;
447          if ( fColumn == 2 || fColumn == 5 || fColumn == 8 ) t = 2;
448        } else {
449          if ( fColumn == 0 || fColumn == 1 || fColumn == 2 ) t = 0;
450          if ( fColumn == 3 || fColumn == 4 || fColumn == 5 ) t = 1;
451          if ( fColumn == 6 || fColumn == 7 || fColumn == 8 ) t = 2;
452        };
453        if ( debug ) printf(" encol: i %i t %i encol %f \n",i,t,encol[i][t]);
454        return encol[i][t];
455      };
456      return(-1.);
457    }
458    
459    Float_t CaloEnergy::GetMaxEnergy(){  
460      printf(" WARNING: OBSOLETE METHOD, use GetMipEnergy() instead! \n");
461      return (this->GetMipEnergy());
462    }
463    
464    Float_t CaloEnergy::GetMipEnergy(){  
465      if ( fLong ){
466        if ( debug ) printf(" oh! call process! with asntr %s and sntr %s \n",asntr.Data(),sntr.Data());
467        this->Process(asntr);
468      };
469      return((fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane));
470  }  }
471    
472    
# Line 458  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 533  Bool_t CaloEnergy::IsInsideAcceptance(TS
533      enstrip[view][plane][strip]=en;      enstrip[view][plane][strip]=en;
534    };    };
535    //    //
536      if ( debug && ((fM1+0.122-0.244*(Float_t)fRad) < 0.) ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad));
537      //
538      // inclination factor (stolen from Daniele's code)
539      //
540      Float_t ytgx = 0;
541      Float_t ytgy = 0;
542      ytgx = 0.76 * cl2->tanx[0];
543      ytgy = 0.76 * cl2->tany[0];  
544      X0pl = sqrt( pow(0.76,2.) + pow(ytgx,2.) + pow(ytgy,2.) );
545      //
546    // sum energy plane by plane for each sections    // sum energy plane by plane for each sections
547    //    //
548      Float_t fen_xep[11];
549      Float_t fen_xop[11];
550      Float_t fen_yep[11];
551      Float_t fen_yop[11];
552      memset(fen_xep,0,11*sizeof(Float_t));
553      memset(fen_xop,0,11*sizeof(Float_t));
554      memset(fen_yep,0,11*sizeof(Float_t));
555      memset(fen_yop,0,11*sizeof(Float_t));
556      //
557    for (Int_t i=0;i<11;i++){    for (Int_t i=0;i<11;i++){
558      for(strip=0; strip<96; strip++) {      for(strip=0; strip<96; strip++) {
559          fen_xep[i] += enstrip[1][2*i][strip];
560          fen_yop[i] += enstrip[0][2*i][strip];
561          fen_xop[i] += enstrip[1][(2*i)+1][strip];
562          fen_yep[i] += enstrip[0][(2*i)+1][strip];
563        if ( fRad < 0 ){        if ( fRad < 0 ){
564          //          //
565          // run over all the strips of the plane          // run over all the strips of the plane
566          //          //
567          en_xep[i] += enstrip[1][2*i][strip];          en_xep[i] += enstrip[1][2*i][strip];
568          en_yop[i] += enstrip[0][2*i][strip];          en_yop[i] += enstrip[0][2*i][strip];
# Line 474  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 572  Bool_t CaloEnergy::IsInsideAcceptance(TS
572          //          //
573          // use only the strips inside a cylinder of given radius fRad          // use only the strips inside a cylinder of given radius fRad
574          //          //
575          if ( strip >= cl2->cibar[2*i][1]-1-fRad &&  strip <= cl2->cibar[2*i][1]-1+fRad ) en_xep[i] += enstrip[1][2*i][strip];          if ( cl2->cibar[2*i][1] >= 1 && cl2->cibar[2*i][1] <= 96 &&
576          if ( strip >= cl2->cibar[2*i][0]-1-fRad &&  strip <= cl2->cibar[2*i][0]-1+fRad ) en_yop[i] += enstrip[0][2*i][strip];               (strip >= (cl2->cibar[2*i][1]-1-fRad)) &&  (strip <= (cl2->cibar[2*i][1]-1+fRad)) ) en_xep[i] += enstrip[1][2*i][strip];
577          if ( strip >= cl2->cibar[(2*i)+1][1]-1-fRad &&  strip <= cl2->cibar[(2*i)+1][1]-1+fRad ) en_xop[i] += enstrip[1][(2*i)+1][strip];  
578          if ( strip >= cl2->cibar[(2*i)+1][0]-1-fRad &&  strip <= cl2->cibar[(2*i)+1][0]-1+fRad ) en_yep[i] += enstrip[0][(2*i)+1][strip];          if ( cl2->cibar[2*i][0] >= 1 && cl2->cibar[2*i][0] <= 96 &&
579                 (strip >= (cl2->cibar[2*i][0]-1-fRad)) &&  (strip <= (cl2->cibar[2*i][0]-1+fRad)) ) en_yop[i] += enstrip[0][2*i][strip];
580    
581            if ( cl2->cibar[(2*i)+1][1] >= 1 && cl2->cibar[(2*i)+1][1] <= 96 &&
582                 (strip >= (cl2->cibar[(2*i)+1][1]-1-fRad)) &&  (strip <= (cl2->cibar[(2*i)+1][1]-1+fRad)) ) en_xop[i] += enstrip[1][(2*i)+1][strip];
583    
584            if ( cl2->cibar[(2*i)+1][0] >= 1 && cl2->cibar[(2*i)+1][0] <= 96 &&
585                 (strip >= (cl2->cibar[(2*i)+1][0]-1-fRad)) &&  (strip <= (cl2->cibar[(2*i)+1][0]-1+fRad)) ) en_yep[i] += enstrip[0][(2*i)+1][strip];
586        };        };
587      };      };
588        if ( debug ) printf(" ex_xep[%i] %f cibar %i \n",i,en_xep[i],cl2->cibar[2*i][1]);
589        if ( debug ) printf(" ex_xop[%i] %f cibar %i \n",i,en_xop[i],cl2->cibar[(2*i)+1][1]);
590        if ( debug ) printf(" ex_yep[%i] %f cibar %i \n",i,en_yep[i],cl2->cibar[(2*i)+1][0]);
591        if ( debug ) printf(" ex_yop[%i] %f cibar %i \n",i,en_yop[i],cl2->cibar[2*i][0]);
592      energyxe += en_xep[i];      energyxe += en_xep[i];
593      energyyo += en_yop[i];      energyyo += en_yop[i];
594      energyxo += en_xop[i];      energyxo += en_xop[i];
# Line 493  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 602  Bool_t CaloEnergy::IsInsideAcceptance(TS
602    Int_t yon = 0;    Int_t yon = 0;
603    Int_t xon = 0;    Int_t xon = 0;
604    Int_t yen = 0;    Int_t yen = 0;
605      Float_t en = 0.;
606    //    //
607    if ( section.Contains("XE") ){    if ( section.Contains("XE") ){
608      yon++;      yon++;
609      xon++;      xon++;
610      yen++;      yen++;
611      for (Int_t ipl =0; ipl < 11; ipl ++) {      for (Int_t ipl =0; ipl < 11; ipl ++) {
612        if(en_xep[ipl] > xemax_en) {        en = fen_xep[ipl];
613          xemax_en = en_xep[ipl];        if ( !fAllpl ) en = en_xep[ipl];
614          if(en > xemax_en) {
615            xemax_en = en;
616          fMax_planexe = ipl;          fMax_planexe = ipl;
617        };        };
618      };      };
# Line 510  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 622  Bool_t CaloEnergy::IsInsideAcceptance(TS
622      xon++;      xon++;
623      yen++;      yen++;
624      for (Int_t ipl =0; ipl < 11; ipl ++) {      for (Int_t ipl =0; ipl < 11; ipl ++) {
625        if(en_yop[ipl] > yomax_en) {        en = fen_yop[ipl];
626          yomax_en = en_yop[ipl];        if ( !fAllpl ) en = en_yop[ipl];
627          if(en > yomax_en) {
628            yomax_en = en;
629          fMax_planeyo = ipl;          fMax_planeyo = ipl;
630        };        };
631      };      };
# Line 520  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 634  Bool_t CaloEnergy::IsInsideAcceptance(TS
634    if ( section.Contains("XO") ){    if ( section.Contains("XO") ){
635      yen++;      yen++;
636      for (Int_t ipl =0; ipl < 11; ipl ++) {      for (Int_t ipl =0; ipl < 11; ipl ++) {
637        if(en_xop[ipl] > xomax_en) {        en = fen_xop[ipl];
638          xomax_en = en_xop[ipl];        if ( !fAllpl ) en = en_xop[ipl];
639          if(en > xomax_en) {
640            xomax_en = en;
641          fMax_planexo = ipl;          fMax_planexo = ipl;
642        };        };
643      };      };
# Line 529  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 645  Bool_t CaloEnergy::IsInsideAcceptance(TS
645    //    //
646    if ( section.Contains("YE") ){    if ( section.Contains("YE") ){
647      for (Int_t ipl =0; ipl < 11; ipl ++) {      for (Int_t ipl =0; ipl < 11; ipl ++) {
648        if(en_yep[ipl] > yemax_en) {        en = fen_yep[ipl];
649          yemax_en = en_yep[ipl];        if ( !fAllpl ) en = en_yep[ipl];
650          if(en > yemax_en) {
651            yemax_en = en;
652          fMax_planeye = ipl;          fMax_planeye = ipl;
653        };        };
654      };      };
# Line 540  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 658  Bool_t CaloEnergy::IsInsideAcceptance(TS
658    //    //
659    if ( x0max > 0. ){    if ( x0max > 0. ){
660      if ( debug ) printf(" CALCULATE MAX PLANE GIVEN X0 ASSUMING PERPENDICULAR TRACK \n");      if ( debug ) printf(" CALCULATE MAX PLANE GIVEN X0 ASSUMING PERPENDICULAR TRACK \n");
661      Int_t wpl = (Int_t)roundf(x0max/0.76);      //    Int_t wpl = (Int_t)roundf(x0max/0.76);
662        Int_t wpl = (Int_t)roundf(x0max/X0pl);
663      Bool_t isY = false;      Bool_t isY = false;
664      if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true;      //    if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true;
665        if ( ((x0max/X0pl)-(Float_t)wpl) > 0. ) isY = true;
666      xomax_en = 0.;      xomax_en = 0.;
667      yemax_en = 0.;      yemax_en = 0.;
668      xemax_en = 0.;      xemax_en = 0.;
# Line 583  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 703  Bool_t CaloEnergy::IsInsideAcceptance(TS
703          if ( section.Contains("YO") ) yomax_en = 10.;          if ( section.Contains("YO") ) yomax_en = 10.;
704        };        };
705      };      };
706      if ( debug ) printf(" x0max %f x0max/0.76 %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/0.76),wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo);      //    if ( debug ) printf(" x0max %f x0max/0.76 %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/0.76),wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo);
707        if ( debug ) printf(" x0max %f x0max/X0pl %f X0pl %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/X0pl),X0pl,wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo);
708    };    };
709    //    //
710    Int_t nPl = fPl;    Int_t nPl = fPl;
# Line 655  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 776  Bool_t CaloEnergy::IsInsideAcceptance(TS
776    //    //
777    // for each plane of the calorimeter find the position of the track in the direction along the strip (where we do not have a measurement from the selected plane) by looking at the plane above/below of the other view and extrapolating the trajectory to the given plane    // for each plane of the calorimeter find the position of the track in the direction along the strip (where we do not have a measurement from the selected plane) by looking at the plane above/below of the other view and extrapolating the trajectory to the given plane
778    //    //
   Float_t  track_coordx[22][2];  
   Float_t  track_coordy[22][2];  
779    //    //
780    Float_t tgx_cl2;    Float_t tgx_cl2;
781    Float_t tgy_cl2;      Float_t tgy_cl2;  
# Line 666  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 785  Bool_t CaloEnergy::IsInsideAcceptance(TS
785    for (Int_t p=0; p<22; p++){        for (Int_t p=0; p<22; p++){    
786      track_coordy[p][1] = cl2->cbar[p][1];      track_coordy[p][1] = cl2->cbar[p][1];
787      track_coordx[p][1] = cl2->cbar[p][0] - fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2;      track_coordx[p][1] = cl2->cbar[p][0] - fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2;
788        //    track_coordx[p][1] = cl2->cbar[p][0] + fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2;
789      track_coordx[p][0] = cl2->cbar[p][0];      track_coordx[p][0] = cl2->cbar[p][0];
790      track_coordy[p][0] = cl2->cbar[p][1] - fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2;      track_coordy[p][0] = cl2->cbar[p][1] - fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2;
791        //    track_coordy[p][0] = cl2->cbar[p][1] + fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2;
792      if ( debug ) printf(" p %i track_coordy[p][1] %f track_coordx[p][1] %f track_coordx[p][0] %f  track_coordy[p][0] %f \n",p,track_coordy[p][1],track_coordx[p][1],track_coordx[p][0],track_coordy[p][0]);      if ( debug ) printf(" p %i track_coordy[p][1] %f track_coordx[p][1] %f track_coordx[p][0] %f  track_coordy[p][0] %f \n",p,track_coordy[p][1],track_coordx[p][1],track_coordx[p][0],track_coordy[p][0]);
793    };    };
794    //    //
# Line 675  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 796  Bool_t CaloEnergy::IsInsideAcceptance(TS
796    //    //
797    if ( section.Contains("XO") ){    if ( section.Contains("XO") ){
798      //      //
799      // check event is inside XO acceptance      // find the column hit in the first plane
800        //
801        Int_t ix = -1;
802        Int_t iy = -1;
803        if ( track_coordx[(2*0)+1][1] >= (-12.054+fM) && track_coordx[(2*0)+1][1] <= (-4.246-fM) ) ix = 0;
804        if ( track_coordx[(2*0)+1][1] >= ( -4.004+fM) && track_coordx[(2*0)+1][1] <= ( 3.804-fM) ) ix = 1;
805        if ( track_coordx[(2*0)+1][1] >= (  4.046+fM) && track_coordx[(2*0)+1][1] <= (11.854-fM) ) ix = 2;
806        if ( cl2->cbar[(2*0)+1][1] >= (xo1 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo2 - fM1) ) iy = 0;
807        if ( cl2->cbar[(2*0)+1][1] >= (xo3 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo4 - fM1) ) iy = 1;
808        if ( cl2->cbar[(2*0)+1][1] >= (xo5 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo6 - fM1) ) iy = 2;
809        if ( ix > -1 && iy > -1 ) fColXO = ix + iy*3;
810        //
811        // check event is inside XO acceptance, if multicol is false (SingleColumn mode) then the track must be contained in a column.
812      //      //
813      for (Int_t i=0; i<11; i++) {      for (Int_t i=0; i<11; i++) {
814        if        if ((
815          (             ( track_coordx[(2*i)+1][1] >= (-12.054+fM) && track_coordx[(2*i)+1][1] <= (-4.246-fM) && (ix == 0 || multicol) ) ||
816           (((track_coordx[(2*i)+1][1]>=(-12.054+fM))&&             ( track_coordx[(2*i)+1][1] >= ( -4.004+fM) && track_coordx[(2*i)+1][1] <= ( 3.804-fM) && (ix == 1 || multicol) ) ||
817             (track_coordx[(2*i)+1][1]<=(-4.246-fM)))&&             ( track_coordx[(2*i)+1][1] >= (  4.046+fM) && track_coordx[(2*i)+1][1] <= (11.854-fM) && (ix == 2 || multicol) )
818            (((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&&             ) && (
819              (cl2->cbar[(2*i)+1][1]<=xo2 - fM1 ))||                   ( cl2->cbar[(2*i)+1][1] >= (xo1 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo2 - fM1) && (iy == 0 || multicol) ) ||
820             ((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&&                   ( cl2->cbar[(2*i)+1][1] >= (xo3 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo4 - fM1) && (iy == 1 || multicol) ) ||      
821              (cl2->cbar[(2*i)+1][1]<=xo4 - fM1 ))||                   ( cl2->cbar[(2*i)+1][1] >= (xo5 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo6 - fM1) && (iy == 2 || multicol) )
822             ((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&&                   )){
             (cl2->cbar[(2*i)+1][1]<=xo6 - fM1 ))))||  
           
          (((track_coordx[(2*i)+1][1]>=(-4.004+fM))&&  
            (track_coordx[(2*i)+1][1]<=(3.804-fM)))&&  
           (((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo2 - fM1 ))||  
            ((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo4 - fM1))||  
            ((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo6 - fM1 ))))||  
           
          (((track_coordx[(2*i)+1][1]>=(4.046+fM))&&  
            (track_coordx[(2*i)+1][1]<=(11.854-fM)))&&  
           (((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo2 - fM1))||  
            ((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo4 - fM1 ))||  
            ((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&&  
             (cl2->cbar[(2*i)+1][1]<=xo6 - fM1 ))))  
          ){  
823          fXosel = true;          fXosel = true;
824          fXoout = i;          fXoout = i;
825        } else {        } else {
# Line 726  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 839  Bool_t CaloEnergy::IsInsideAcceptance(TS
839      // event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section      // event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section
840      //      //
841      if ( fXosel ){      if ( fXosel ){
842        for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)) ;iplm++){                for (Int_t iplm=0; iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)); iplm++){      
843          fXOen_maxplane += en_xop[iplm];          fXOen_maxplane += en_xop[iplm];
844          if ( debug ) printf(" XO iplm %i fXOen_maxplane %f  en_xop[iplm] %f\n",iplm,fXOen_maxplane,en_xop[iplm]);          if ( debug ) printf(" XO iplm %i fXOen_maxplane %f  en_xop[iplm] %f\n",iplm,fXOen_maxplane,en_xop[iplm]);
845        };        };
846        fEnergyxo = fXOen_maxplane/fConv_rxo;        fEnergyxo = fXOen_maxplane/fConv_rxo;
847          //
848          for (Int_t i=0;i<11;i++){
849            for(strip=0; strip<96; strip++) {
850              //
851              // run over all the strips of the plane
852              //
853              if ( strip >=  0 && strip < 32 ) encol[1][0] += enstrip[1][(2*i)+1][strip];
854              if ( strip >= 32 && strip < 64 ) encol[1][1] += enstrip[1][(2*i)+1][strip];
855              if ( strip >= 64 && strip < 96 ) encol[1][2] += enstrip[1][(2*i)+1][strip];
856              entot[1] += enstrip[1][(2*i)+1][strip];
857              //
858            };
859          };
860      };      };
861    };    };
862    //    //
863    if ( section.Contains("XE") ){    if ( section.Contains("XE") ){
864        //
865        // find the column hit in the first plane
866        //
867        Int_t ix = -1;
868        Int_t iy = -1;
869        if ( track_coordx[(2*0)][1] >= (-11.854+fM) && track_coordx[(2*0)][1] <= (-4.046-fM) ) ix = 0;
870        if ( track_coordx[(2*0)][1] >= ( -3.804+fM) && track_coordx[(2*0)][1] <= ( 4.004-fM) ) ix = 1;
871        if ( track_coordx[(2*0)][1] >= (  4.246+fM) && track_coordx[(2*0)][1] <= (12.054-fM) ) ix = 2;
872        if ( cl2->cbar[(2*0)][1] >= (xe1 + fM1) && cl2->cbar[(2*0)][1] <= (xe2 - fM1) ) iy = 0;
873        if ( cl2->cbar[(2*0)][1] >= (xe3 + fM1) && cl2->cbar[(2*0)][1] <= (xe4 - fM1) ) iy = 1;
874        if ( cl2->cbar[(2*0)][1] >= (xe5 + fM1) && cl2->cbar[(2*0)][1] <= (xe6 - fM1) ) iy = 2;
875        if ( ix > -1 && iy > -1 ) fColXE = ix + iy*3;
876        //
877        // check event is inside XO acceptance
878        //
879      for (Int_t i=0; i<11; i++) {      for (Int_t i=0; i<11; i++) {
880        if        if ((
881          (             ( track_coordx[(2*i)][1] >= (-11.854+fM) && track_coordx[(2*i)][1] <= (-4.046-fM) && (ix == 0 || multicol) ) ||
882           (((track_coordx[2*i][1]>=(-11.854+fM))&&             ( track_coordx[(2*i)][1] >= ( -3.804+fM) && track_coordx[(2*i)][1] <= ( 4.004-fM) && (ix == 1 || multicol) ) ||
883             (track_coordx[2*i][1]<=(-4.046-fM)))&&             ( track_coordx[(2*i)][1] >= (  4.246+fM) && track_coordx[(2*i)][1] <= (12.054-fM) && (ix == 2 || multicol) )
884            (((cl2->cbar[2*i][1]>=xe1 + fM1)&&             ) && (
885              (cl2->cbar[2*i][1]<=xe2 - fM1 ))||                   ( cl2->cbar[(2*i)][1] >= (xe1 + fM1) && cl2->cbar[(2*i)][1] <= (xe2 - fM1) && (iy == 0 || multicol) ) ||
886             ((cl2->cbar[2*i][1]>=xe3 + fM1)&&                   ( cl2->cbar[(2*i)][1] >= (xe3 + fM1) && cl2->cbar[(2*i)][1] <= (xe4 - fM1) && (iy == 1 || multicol) ) ||          
887              (cl2->cbar[2*i][1]<=xe4 - fM1 ))||                   ( cl2->cbar[(2*i)][1] >= (xe5 + fM1) && cl2->cbar[(2*i)][1] <= (xe6 - fM1) && (iy == 2 || multicol) )
888             ((cl2->cbar[2*i][1]>=xe5 + fM1)&&                   )){
             (cl2->cbar[2*i][1]<=xe6 - fM1 ))))||  
           
          (((track_coordx[2*i][1]>=(-3.804+fM))&&  
            (track_coordx[2*i][1]<=(4.004-fM)))&&  
           (((cl2->cbar[2*i][1]>=xe1 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe2 - fM1 ))||  
            ((cl2->cbar[2*i][1]>=xe3 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe4 - fM1))||  
            ((cl2->cbar[2*i][1]>=xe5 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe6 - fM1 ))))||  
           
          (((track_coordx[2*i][1]>=(4.246+fM))&&  
            (track_coordx[2*i][1]<=(12.054-fM)))&&  
           (((cl2->cbar[2*i][1]>=xe1 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe2 - fM1))||  
            ((cl2->cbar[2*i][1]>=xe3 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe4 - fM1 ))||  
            ((cl2->cbar[2*i][1]>=xe5 + fM1)&&  
             (cl2->cbar[2*i][1]<=xe6 - fM1 ))))  
          ){  
889          fXesel = true;          fXesel = true;
890          fXeout = i;          fXeout = i;
891        } else {        } else {
# Line 772  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 893  Bool_t CaloEnergy::IsInsideAcceptance(TS
893          break;          break;
894        };        };
895      };      };
896        //
897      if ( !fXesel && fXeout >= fXemin && fXeout >= (Int_t)(fMax_planexe+nPl) ){      if ( !fXesel && fXeout >= fXemin && fXeout >= (Int_t)(fMax_planexe+nPl) ){
898        if ( debug ) printf(" XE: this event is only partially contained: fXeout %i fXemin %i fMax_planexe + nPl %i \n",fXeout,fXemin,(Int_t)(fMax_planexe+nPl));        if ( debug ) printf(" XE: this event is only partially contained: fXeout %i fXemin %i fMax_planexe + nPl %i \n",fXeout,fXemin,(Int_t)(fMax_planexe+nPl));
899        fPartsel = true;        fPartsel = true;
# Line 783  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 905  Bool_t CaloEnergy::IsInsideAcceptance(TS
905          if ( debug ) printf(" XE iplm %i fXOen_maxplane %f  en_xop[iplm] %f\n",iplm,fXEen_maxplane,en_xep[iplm]);          if ( debug ) printf(" XE iplm %i fXOen_maxplane %f  en_xop[iplm] %f\n",iplm,fXEen_maxplane,en_xep[iplm]);
906        };        };
907        fEnergyxe = fXEen_maxplane/fConv_rxe;        fEnergyxe = fXEen_maxplane/fConv_rxe;
908          //
909          for (Int_t i=0;i<11;i++){
910            for(strip=0; strip<96; strip++) {
911              //
912              // run over all the strips of the plane
913              //
914              if ( strip >=  0 && strip < 32 ) encol[1][0] += enstrip[1][(2*i)][strip];
915              if ( strip >= 32 && strip < 64 ) encol[1][1] += enstrip[1][(2*i)][strip];
916              if ( strip >= 64 && strip < 96 ) encol[1][2] += enstrip[1][(2*i)][strip];
917              entot[1] += enstrip[1][(2*i)][strip];
918              //
919            };
920          };
921      };      };
922    };      };  
923    //    //
924    if ( section.Contains("YE") ){    if ( section.Contains("YE") ){
925        //
926        // find the column hit in the first plane
927        //
928        Int_t ix = -1;
929        Int_t iy = -1;
930        if ( track_coordy[(2*0)+1][0] >= (-12.154+fM) && track_coordy[(2*0)+1][0] <= (-4.346-fM) ) iy = 0;
931        if ( track_coordy[(2*0)+1][0] >= ( -4.104+fM) && track_coordy[(2*0)+1][0] <= ( 3.704-fM) ) iy = 1;
932        if ( track_coordy[(2*0)+1][0] >= (  3.946+fM) && track_coordy[(2*0)+1][0] <= (11.754-fM) ) iy = 2;
933        if ( cl2->cbar[(2*0)+1][0] >= (ye1 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye2 - fM1) ) ix = 0;
934        if ( cl2->cbar[(2*0)+1][0] >= (ye3 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye4 - fM1) ) ix = 1;
935        if ( cl2->cbar[(2*0)+1][0] >= (ye5 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye6 - fM1) ) ix = 2;
936        if ( ix > -1 && iy > -1 ) fColYE = ix + iy*3;
937        //
938        // check event is inside XO acceptance
939        //
940      for (Int_t i=0; i<11; i++) {      for (Int_t i=0; i<11; i++) {
941        if        if ((
942          (             ( track_coordy[(2*i)+1][0] >= (-12.154+fM) && track_coordy[(2*i)+1][0] <= (-4.346-fM) && (iy == 0 || multicol) ) ||
943           (((track_coordy[(2*i)+1][0]>=(-12.154+fM))&&             ( track_coordy[(2*i)+1][0] >= ( -4.104+fM) && track_coordy[(2*i)+1][0] <= ( 3.704-fM) && (iy == 1 || multicol) ) ||
944             (track_coordy[(2*i)+1][0]<=(-4.346-fM)))&&             ( track_coordy[(2*i)+1][0] >= (  3.946+fM) && track_coordy[(2*i)+1][0] <= (11.754-fM) && (iy == 2 || multicol) )
945            (((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&&             ) && (
946              (cl2->cbar[(2*i)+1][0]<=ye2 - fM1 ))||                   ( cl2->cbar[(2*i)+1][0] >= (ye1 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye2 - fM1) && (ix == 0 || multicol) ) ||
947             ((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&&                   ( cl2->cbar[(2*i)+1][0] >= (ye3 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye4 - fM1) && (ix == 1 || multicol) ) ||      
948              (cl2->cbar[(2*i)+1][0]<=ye4 - fM1 ))||                   ( cl2->cbar[(2*i)+1][0] >= (ye5 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye6 - fM1) && (ix == 2 || multicol) )
949             ((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&&                   )){
             (cl2->cbar[(2*i)+1][0]<=ye6 - fM1 ))))||  
           
          (((track_coordy[(2*i)+1][0]>=(-4.104+fM))&&  
            (track_coordy[(2*i)+1][0]<=(3.704-fM)))&&  
           (((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye2 - fM1 ))||  
            ((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye4 - fM1))||  
            ((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye6 - fM1 ))))||  
           
          (((track_coordy[(2*i)+1][0]>=(3.946+fM))&&  
            (track_coordy[(2*i)+1][0]<=(11.754-fM)))&&  
           (((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye2 - fM1))||  
            ((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye4 - fM1 ))||  
            ((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&&  
             (cl2->cbar[(2*i)+1][0]<=ye6 - fM1 ))))  
          ){  
950          fYesel = true;          fYesel = true;
951          fYeout = i;          fYeout = i;
952        } else {        } else {
# Line 824  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 954  Bool_t CaloEnergy::IsInsideAcceptance(TS
954          break;          break;
955        };        };
956      };      };
957        //
958      if ( !fYesel && fYeout >= fYemin && fYeout >= (Int_t)(fMax_planeye+nPl) ){      if ( !fYesel && fYeout >= fYemin && fYeout >= (Int_t)(fMax_planeye+nPl) ){
959        if ( debug ) printf(" YE: this event is only partially contained: fYeout %i fYemin %i fMax_planeye + nPl %i \n",fYeout,fYemin,(Int_t)(fMax_planeye+nPl));        if ( debug ) printf(" YE: this event is only partially contained: fYeout %i fYemin %i fMax_planeye + nPl %i \n",fYeout,fYemin,(Int_t)(fMax_planeye+nPl));
960        fPartsel = true;        fPartsel = true;
# Line 832  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 963  Bool_t CaloEnergy::IsInsideAcceptance(TS
963      if ( fYesel ){      if ( fYesel ){
964        for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeye+nPl)) ;iplm++) fYEen_maxplane += en_yep[iplm];        for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeye+nPl)) ;iplm++) fYEen_maxplane += en_yep[iplm];
965        fEnergyye = fYEen_maxplane/fConv_rye;        fEnergyye = fYEen_maxplane/fConv_rye;
966          //
967          for (Int_t i=0;i<11;i++){
968            for(strip=0; strip<96; strip++) {
969              //
970              // run over all the strips of the plane
971              //
972              if ( strip >=  0 && strip < 32 ) encol[0][0] += enstrip[0][(2*i)+1][strip];
973              if ( strip >= 32 && strip < 64 ) encol[0][1] += enstrip[0][(2*i)+1][strip];
974              if ( strip >= 64 && strip < 96 ) encol[0][2] += enstrip[0][(2*i)+1][strip];
975              entot[0] += enstrip[0][(2*i)+1][strip];
976              //
977            };
978          };
979          //
980      };      };
981    };      };  
982    //    //
983    if ( section.Contains("YO") ){    if ( section.Contains("YO") ){
984        //
985        // find the column hit in the first plane
986        //
987        Int_t ix = -1;
988        Int_t iy = -1;
989        if ( track_coordy[(2*0)][0] >= (-11.954+fM) && track_coordy[(2*0)][0] <= (-4.146-fM) ) iy = 0;
990        if ( track_coordy[(2*0)][0] >= ( -3.904+fM) && track_coordy[(2*0)][0] <= ( 3.904-fM) ) iy = 1;
991        if ( track_coordy[(2*0)][0] >= (  4.146+fM) && track_coordy[(2*0)][0] <= (11.954-fM) ) iy = 2;
992        if ( cl2->cbar[(2*0)][0] >= (yo1 + fM1) && cl2->cbar[(2*0)][0] <= (yo2 - fM1) ) ix = 0;
993        if ( cl2->cbar[(2*0)][0] >= (yo3 + fM1) && cl2->cbar[(2*0)][0] <= (yo4 - fM1) ) ix = 1;
994        if ( cl2->cbar[(2*0)][0] >= (yo5 + fM1) && cl2->cbar[(2*0)][0] <= (yo6 - fM1) ) ix = 2;
995        if ( ix > -1 && iy > -1 ) fColYO = ix + iy*3;
996        //
997        // check event is inside XO acceptance
998        //
999      for (Int_t i=0; i<11; i++) {      for (Int_t i=0; i<11; i++) {
1000        if ( debug ) printf(" i %i track_coordy[2*i][0] %f cl2->cbar[2*i][0] %f \n",i,track_coordy[2*i][0],cl2->cbar[2*i][0]);        if ((
1001        if ( debug ) printf(" i %i fm %f fm1 %f yo1 %g yo2 %f yo3 %f yo4 %f yo5 %f yo6 %f \n",i,fM,fM1,yo1,yo2,yo3,yo4,yo5,yo6);             ( track_coordy[(2*i)][0] >= (-11.954+fM) && track_coordy[(2*i)][0] <= (-4.146-fM) && (iy == 0 || multicol) ) ||
1002        if             ( track_coordy[(2*i)][0] >= ( -3.904+fM) && track_coordy[(2*i)][0] <= ( 3.904-fM) && (iy == 1 || multicol) ) ||
1003          (             ( track_coordy[(2*i)][0] >= (  4.146+fM) && track_coordy[(2*i)][0] <= (11.954-fM) && (iy == 2 || multicol) )
1004           (((track_coordy[2*i][0]>=(-11.954+fM))&&             ) && (
1005             (track_coordy[2*i][0]<=(-4.146-fM)))&&                   ( cl2->cbar[(2*i)][0] >= (yo1 + fM1) && cl2->cbar[(2*i)][0] <= (yo2 - fM1) && (ix == 0 || multicol) ) ||
1006            (((cl2->cbar[2*i][0]>=yo1 + fM1)&&                   ( cl2->cbar[(2*i)][0] >= (yo3 + fM1) && cl2->cbar[(2*i)][0] <= (yo4 - fM1) && (ix == 1 || multicol) ) ||          
1007              (cl2->cbar[2*i][0]<=yo2 - fM1 ))||                   ( cl2->cbar[(2*i)][0] >= (yo5 + fM1) && cl2->cbar[(2*i)][0] <= (yo6 - fM1) && (ix == 2 || multicol) )
1008             ((cl2->cbar[2*i][0]>=yo3 + fM1)&&                   )){
             (cl2->cbar[2*i][0]<=yo4 - fM1 ))||  
            ((cl2->cbar[2*i][0]>=yo5 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo6 - fM1 ))))||  
           
          (((track_coordy[2*i][0]>=(-3.904+fM))&&  
            (track_coordy[2*i][0]<=(+3.904-fM)))&&  
           (((cl2->cbar[2*i][0]>=yo1 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo2 - fM1 ))||  
            ((cl2->cbar[2*i][0]>=yo3 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo4 - fM1))||  
            ((cl2->cbar[2*i][0]>=yo5 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo6 - fM1 ))))||  
           
          (((track_coordy[2*i][0]>=(4.146+fM))&&  
            (track_coordy[2*i][0]<=(11.954-fM)))&&  
           (((cl2->cbar[2*i][0]>=yo1 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo2 - fM1))||  
            ((cl2->cbar[2*i][0]>=yo3 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo4 - fM1 ))||  
            ((cl2->cbar[2*i][0]>=yo5 + fM1)&&  
             (cl2->cbar[2*i][0]<=yo6 - fM1 ))))  
          ){  
1009          fYosel = true;          fYosel = true;
1010          fYoout = i;          fYoout = i;
1011        } else {        } else {
# Line 875  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 1013  Bool_t CaloEnergy::IsInsideAcceptance(TS
1013          break;          break;
1014        };        };
1015      };      };
1016        //
1017      if ( !fYosel && fYoout >= fYomin && fYoout >= (Int_t)(fMax_planeyo+nPl) ){      if ( !fYosel && fYoout >= fYomin && fYoout >= (Int_t)(fMax_planeyo+nPl) ){
1018        if ( debug ) printf(" YO: this event is only partially contained: fYoout %i fYomin %i fMax_planeyo + nPl %i \n",fYoout,fYomin,(Int_t)(fMax_planeyo+nPl));        if ( debug ) printf(" YO: this event is only partially contained: fYoout %i fYomin %i fMax_planeyo + nPl %i \n",fYoout,fYomin,(Int_t)(fMax_planeyo+nPl));
1019        fPartsel = true;        fPartsel = true;
# Line 883  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 1022  Bool_t CaloEnergy::IsInsideAcceptance(TS
1022      if ( fYosel ){      if ( fYosel ){
1023        for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeyo+nPl)) ;iplm++) fYOen_maxplane += en_yop[iplm];        for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeyo+nPl)) ;iplm++) fYOen_maxplane += en_yop[iplm];
1024        fEnergyyo = fYOen_maxplane/fConv_ryo;        fEnergyyo = fYOen_maxplane/fConv_ryo;
1025          //
1026          for (Int_t i=0;i<11;i++){
1027            for(strip=0; strip<96; strip++) {
1028              //
1029              // run over all the strips of the plane
1030              //
1031              if ( strip >=  0 && strip < 32 ) encol[0][0] += enstrip[0][(2*i)][strip];
1032              if ( strip >= 32 && strip < 64 ) encol[0][1] += enstrip[0][(2*i)][strip];
1033              if ( strip >= 64 && strip < 96 ) encol[0][2] += enstrip[0][(2*i)][strip];
1034              entot[0] += enstrip[0][(2*i)][strip];
1035              //
1036            };
1037          };
1038      };      };
1039    };      };  
1040    //    //
# Line 908  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 1060  Bool_t CaloEnergy::IsInsideAcceptance(TS
1060        fSel = true;        fSel = true;
1061      };      };
1062      fEnergy = (fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)/fConv_rxo;      fEnergy = (fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)/fConv_rxo;
1063        if ( fSel ){
1064          if ( fXesel ) fColumn = fColXE;
1065          if ( fXosel ){
1066            if ( fColXO != fColumn && fColumn > -1 ){
1067              printf(" ERROR! mismatch in column number between different sections! fColumn %i fColXO %i \n",fColumn,fColXO);
1068            } else {
1069              fColumn = fColXO;
1070            };
1071          };
1072          if ( fYesel ){
1073            if ( fColYE != fColumn && fColumn > -1 ){
1074              printf(" ERROR! mismatch in column number between different sections! fColumn %i fColYE %i \n",fColumn,fColYE);
1075            } else {
1076              fColumn = fColYE;
1077            };
1078          };
1079          if ( fYosel ){
1080            if ( fColYO != fColumn && fColumn > -1 ){
1081              printf(" ERROR! mismatch in column number between different sections! fColumn %i fColYO %i \n",fColumn,fColYO);
1082            } else {
1083              fColumn = fColYO;
1084            };
1085          };
1086        };
1087    };    };
1088    //    //
1089    if ( debug ) printf("sel %i indep %i fMax_plane %f conv_r %f en_maxplane %f encalo %f  \n",fSel,indep,fMax_plane,fConv_rxo,fXOen_maxplane,fEnergy);    if ( debug ) printf("sel %i indep %i fMax_plane %f conv_r %f en_maxplane %f encalo %f  \n",fSel,indep,fMax_plane,fConv_rxo,fXOen_maxplane,fEnergy);
# Line 960  void CaloEnergy::Process(TString section Line 1136  void CaloEnergy::Process(TString section
1136    //    //
1137    // check if the cylinder of integration can go out of the sensor given the frame which has been set (if we use all the calorimeter fRad is < 0 and the printout is suppressed)    // check if the cylinder of integration can go out of the sensor given the frame which has been set (if we use all the calorimeter fRad is < 0 and the printout is suppressed)
1138    //    //
1139    if ( (fM1+0.122-0.244*(Float_t)fRad) < 0. ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad));    if ( (fM1+0.122-0.244*(Float_t)fRad) < 0.001 ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad));
1140    //    //
1141    if ( fLong ){    if ( fLong ){
1142      if ( debug ) printf(" ==================================================================> LONGITUDINAL FIT! \n");      if ( debug ) printf(" ==================================================================> LONGITUDINAL FIT! \n");
# Line 996  void CaloEnergy::Process(TString section Line 1172  void CaloEnergy::Process(TString section
1172        } else {        } else {
1173          fXOen_maxplane = clong->Get_E0();          fXOen_maxplane = clong->Get_E0();
1174        };        };
1175          fMax_plane = clong->Get_tmax();
1176        fYOen_maxplane  = 0.;        fYOen_maxplane  = 0.;
1177        fYEen_maxplane = 0.;        fYEen_maxplane = 0.;
1178        fXEen_maxplane = 0.;        fXEen_maxplane = 0.;
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.23