/[PAMELA software]/DarthVader/CalorimeterLevel2/src/CaloLevel2.cpp
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/CaloLevel2.cpp

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

revision 1.7 by mocchiut, Thu Nov 9 17:05:44 2006 UTC revision 1.14 by mocchiut, Mon Sep 3 08:42:13 2007 UTC
# Line 5  Line 5 
5  **/  **/
6  #include <TObject.h>  #include <TObject.h>
7  #include <CaloLevel2.h>  #include <CaloLevel2.h>
8    
9    //
10  ClassImp(CaloTrkVar);  ClassImp(CaloTrkVar);
11  ClassImp(CaloLevel2);  ClassImp(CaloLevel2);
12    
# Line 84  CaloTrkVar::CaloTrkVar(const CaloTrkVar Line 86  CaloTrkVar::CaloTrkVar(const CaloTrkVar
86  **/  **/
87  CaloLevel2::CaloLevel2() {      CaloLevel2::CaloLevel2() {    
88    //    //
89    CaloTrk = new TClonesArray("CaloTrkVar",1);  //  CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
90    estrip = TArrayF(0,NULL);    CaloTrk = 0; //ELENA
91    //    //
92    this->Clear();    this->Clear();
93    //    //
94  };  };
95    /**
96     * Create the TClonesArray
97    **/
98    void CaloLevel2::Set(){//ELENA
99        if(!CaloTrk)CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
100    }//ELENA
101    
102    /**
103     * Clear the CaloLevel2 object
104     **/
105  void CaloLevel2::Clear() {      void CaloLevel2::Clear() {    
106    //    //
107    CaloTrk->Clear();  //  CaloTrk->Clear(); //ELENA
108      if(CaloTrk)CaloTrk->Delete(); //ELENA
109    //    //
110    nstrip = 0;    nstrip = 0;
111    qtot = 0.;    qtot = 0.;
112    impx = 0.;    //  impx = 0.;
113    impy = 0.;    //  impy = 0.;
   tanx = 0.;  
   tany = 0.;  
114    qmax = 0.;    qmax = 0.;
115    nx22 = 0;    nx22 = 0;
116    qx22 = 0.;    qx22 = 0.;
# Line 110  void CaloLevel2::Clear() {     Line 120  void CaloLevel2::Clear() {    
120    memset(swerr, 0, 4*sizeof(Int_t));    memset(swerr, 0, 4*sizeof(Int_t));
121    memset(crc, 0, 4*sizeof(Int_t));    memset(crc, 0, 4*sizeof(Int_t));
122    memset(qq, 0, 4*sizeof(Int_t));    memset(qq, 0, 4*sizeof(Int_t));
123    memset(varcfit, 0, 2*sizeof(Float_t));    memset(varcfit, 0, 4*sizeof(Float_t));
124    memset(npcfit, 0, 2*sizeof(Int_t));    memset(npcfit, 0, 4*sizeof(Int_t));
125      memset(tanx, 0, 2*sizeof(Int_t));
126      memset(tany, 0, 2*sizeof(Int_t));
127      memset(fitmode, 0, 2*sizeof(Int_t));
128    memset(planemax, 0, 2*sizeof(Int_t));    memset(planemax, 0, 2*sizeof(Int_t));
129      memset(selfdelay, 0, 4*7*sizeof(Int_t));
130    memset(cibar, 0, 2*22*sizeof(Int_t));    memset(cibar, 0, 2*22*sizeof(Int_t));
131    memset(cbar, 0, 2*22*sizeof(Float_t));    memset(cbar, 0, 2*22*sizeof(Float_t));
132    good = 0;    good = 0;
133    selftrigger = 0;    selftrigger = 0;
134    estrip.Reset();    //
135  };  };
136    
137    /**
138     * Delete the CaloLevel2 object
139     **/
140    void CaloLevel2::Delete() {     //ELENA
141      if(CaloTrk){ //ELENA
142          CaloTrk->Delete(); //ELENA
143          delete CaloTrk; //ELENA
144      }     //ELENA
145    } //ELENA
146    
147  /**  /**
148   * Fills a struct cCaloLevel2 with values from a CaloLevel2 object (to put data into a F77 common).   * Fills a struct cCaloLevel2 with values from a CaloLevel2 object (to put data into a F77 common).
# Line 133  void CaloLevel2::GetLevel2Struct(cCaloLe Line 156  void CaloLevel2::GetLevel2Struct(cCaloLe
156    l2->qtot = qtot;        l2->qtot = qtot;    
157    l2->qx22 = qx22;        l2->qx22 = qx22;    
158    l2->qmax = qmax;        l2->qmax = qmax;    
159    l2->impx = impx;        //  l2->impx = impx;    
160    l2->impy = impy;        //  l2->impy = impy;    
161    l2->tanx = tanx;        //  l2->tanx = tanx;    
162    l2->tany = tany;      //  l2->tany = tany;  
163    l2->elen = elen;          l2->elen = elen;      
164    l2->selen = selen;    l2->selen = selen;
165                    
# Line 152  void CaloLevel2::GetLevel2Struct(cCaloLe Line 175  void CaloLevel2::GetLevel2Struct(cCaloLe
175      l2->qq[i] = qq[i];        l2->qq[i] = qq[i];  
176    }    }
177    
178    l2->calntrk = CaloTrk->GetEntries();    if(CaloTrk){ //ELENA
179            l2->calntrk = CaloTrk->GetEntries();  
180    for(Int_t i=0;i<l2->calntrk;i++){        for(Int_t i=0;i<l2->calntrk;i++){
181      l2->caltrkseqno[i] = ((CaloTrkVar *)CaloTrk->At(i))->trkseqno;              l2->caltrkseqno[i] = ((CaloTrkVar *)CaloTrk->At(i))->trkseqno;  
182      l2->ncore[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncore;                l2->ncore[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncore;    
183      l2->noint[i] = ((CaloTrkVar *)CaloTrk->At(i))->noint;                l2->noint[i] = ((CaloTrkVar *)CaloTrk->At(i))->noint;    
184      l2->ncyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncyl;                  l2->ncyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncyl;      
185      l2->nlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlast;                l2->nlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlast;    
186      l2->npre[i] = ((CaloTrkVar *)CaloTrk->At(i))->npre;                  l2->npre[i] = ((CaloTrkVar *)CaloTrk->At(i))->npre;      
187      l2->npresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->npresh;                l2->npresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->npresh;    
188      l2->ntr[i] = ((CaloTrkVar *)CaloTrk->At(i))->ntr;                  l2->ntr[i] = ((CaloTrkVar *)CaloTrk->At(i))->ntr;      
189      l2->planetot[i] = ((CaloTrkVar *)CaloTrk->At(i))->planetot;              l2->planetot[i] = ((CaloTrkVar *)CaloTrk->At(i))->planetot;  
190      l2->nlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlow;                  l2->nlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlow;      
191      l2->qcore[i] =((CaloTrkVar *)CaloTrk->At(i))->qcore ;              l2->qcore[i] =((CaloTrkVar *)CaloTrk->At(i))->qcore ;  
192      l2->qcyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->qcyl;                l2->qcyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->qcyl;    
193      l2->qlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlast;              l2->qlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlast;  
194      l2->qpre[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpre;                l2->qpre[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpre;    
195      l2->qpresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpresh;              l2->qpresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpresh;  
196      l2->qtr[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtr;                l2->qtr[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtr;    
197      l2->qtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrack;              l2->qtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrack;  
198      l2->qtrackx[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrackx;            l2->qtrackx[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrackx;
199      l2->qtracky[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtracky;            l2->qtracky[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtracky;
200      l2->dxtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dxtrack;            l2->dxtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dxtrack;
201      l2->dytrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dytrack;            l2->dytrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dytrack;
202      l2->qmean[i] = ((CaloTrkVar *)CaloTrk->At(i))->qmean;              l2->qmean[i] = ((CaloTrkVar *)CaloTrk->At(i))->qmean;  
203      l2->qlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlow;                l2->qlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlow;    
204      l2->dX0l[i] = ((CaloTrkVar *)CaloTrk->At(i))->dX0l;              l2->dX0l[i] = ((CaloTrkVar *)CaloTrk->At(i))->dX0l;  
205      for (Int_t j=0; j<2; j++){            for (Int_t j=0; j<2; j++){
206        for (Int_t k=0; k<22; k++){                for (Int_t k=0; k<22; k++){
207          l2->tbar[i][k][j] = ((CaloTrkVar *)CaloTrk->At(i))->tbar[k][j];                    l2->tbar[i][k][j] = ((CaloTrkVar *)CaloTrk->At(i))->tbar[k][j];
208        };                };
209      };            };
210    }          }  
211      } //ELENA
212  }  }
213    
214  /**  /**
215   * Gives the detected energy for the given strip once loaded the event   * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
216     * @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
217  **/  **/
218  Float_t CaloLevel2::GetEstrip(Int_t view, Int_t plane, Int_t strip){  Float_t CaloLevel2::impx(Int_t tr){
219    Int_t splane = 0;    if ( tr == 0 ) return(cbar[0][0]);
220    Int_t sstrip = 0;    if ( tr == 1 ) {
221    //      if ( !CaloTrk ) return(-110.);
222    if ( nstrip == 0 ) return(0.);      TClonesArray &t = *(CaloTrk);
223    //      for (Int_t itrk=0; itrk<ntrk(); itrk++){
224    for (Int_t i = 0; i<nstrip; i++ ){        CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
225      if ( view == 0 ){        if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][0]);
       if ( estrip.At(i) > 0. ){  
         splane = (Int_t)trunc(estrip.At(i)/1000000.);  
         sstrip = (Int_t)trunc((estrip.At(i)-((Float_t)splane*1000000.))/10000.);  
         if ( splane == plane && sstrip == strip ) return(estrip.At(i)-(Float_t)splane*1000000.-(Float_t)sstrip*10000.);  
       };          
     } else {  
       if ( estrip.At(i) < 0. ){  
         splane = (Int_t)trunc(-estrip.At(i)/1000000.);  
         sstrip = (Int_t)trunc((-estrip.At(i)-((Float_t)splane*1000000.))/10000.);  
         if ( splane == plane && sstrip == strip ) return(-estrip.At(i)-(Float_t)splane*1000000.-(Float_t)sstrip*10000.);          
       };          
226      };      };
227    };    };
228    return(0.);    if ( tr !=0 && tr !=1 ){
229        printf(" Cannot get impx for other than calo or selftrigger tracks!\n");
230      } else {
231        printf(" Cannot find selftrigger block\n");
232      };
233      return(-100.);
234  };  };
235    
236  /**  /**
237   * Given estrip entry returns energy and strip   * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
238     * @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
239  **/  **/
240  Float_t CaloLevel2::DecodeEstrip(Int_t entry, Int_t &view, Int_t &plane, Int_t &strip){  Float_t CaloLevel2::impy(Int_t tr){
241    if ( entry>nstrip ) return(0.);    if ( tr == 0 ) return(cbar[0][1]);
242    //    if ( tr == 1 ) {
243    if ( estrip.At(entry) > 0. ){      if ( !CaloTrk ) return(-110.);
244      view = 0;      TClonesArray &t = *(CaloTrk);
245      plane = (Int_t)trunc(estrip.At(entry)/1000000.);      for (Int_t itrk=0; itrk<ntrk(); itrk++){
246      strip = (Int_t)trunc((estrip.At(entry)-((Float_t)plane*1000000.))/10000.);        CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
247      return(estrip.At(entry)-(Float_t)plane*1000000.-(Float_t)strip*10000.);        if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][1]);
248        };
249    };    };
250    if ( estrip.At(entry) < 0. ){    if ( tr !=0 && tr !=1 ){
251      view = 1;      printf(" Cannot get impy for other than calo or selftrigger tracks!\n");
252      plane = (Int_t)trunc(-estrip.At(entry)/1000000.);    } else {
253      strip = (Int_t)trunc((-estrip.At(entry)-((Float_t)plane*1000000.))/10000.);      printf(" Cannot find selftrigger block\n");
254      return(-estrip.At(entry)-(Float_t)plane*1000000.-(Float_t)strip*10000.);        };
255    };        return(-100.);
256    //  };
257    printf(" WARNING: problems decoding value %f at entry %i \n",estrip.At(entry),entry);  /**
258    //   * Should return the energy in GeV if the particle would be an electron
259    view = -1;   * using a parametrization taken from Monte Carlo simulation
260    plane = -1;  **/
   strip = -1;  
   return(0.);    
 }  
   
261  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){
262    if ( nstrip == 0 ) return;    if ( nstrip == 0 ) return;
263    energy = qtot * 40.82 * 0.000106;    energy = qtot * 40.82 * 0.000106;
# Line 260  CaloTrkVar *CaloLevel2::GetCaloTrkVar(In Line 276  CaloTrkVar *CaloLevel2::GetCaloTrkVar(In
276      printf("                   stored track related variables = %i \n",ntrk());      printf("                   stored track related variables = %i \n",ntrk());
277      return(NULL);      return(NULL);
278    }    }
279      if(!CaloTrk)return 0; //ELENA
280    TClonesArray &t = *(CaloTrk);    TClonesArray &t = *(CaloTrk);
281    CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];    CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
282    return calotrack;    return calotrack;
283  }  }
284    
285    /**
286     * Retrieves the calorimeter track matching the seqno-th tracker stored track.
287     * (If seqno = -1 retrieves the self-trigger calorimeter track)
288     */
289    CaloTrkVar *CaloLevel2::GetCaloStoredTrack(int seqno){
290    
291      if( ntrk()==0 ){
292        printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i but no Calorimeter tracks are stored\n",seqno);
293        return NULL;
294      };
295      
296      CaloTrkVar *c = 0;
297      Int_t it_calo=0;
298        
299      do {
300        c = GetCaloTrkVar(it_calo);
301        it_calo++;
302      } while( c && seqno != c->trkseqno && it_calo < ntrk());      
303      
304      if(!c || seqno != c->trkseqno){
305        c = 0;
306        if(seqno!=-1 && seqno !=-2 && seqno!=-3 ) printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i does not match Calorimeter stored tracks\n",seqno);
307      };
308      return c;
309        
310    };
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23