/[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.12 by mocchiut, Mon Jan 22 09:16:59 2007 UTC revision 1.20 by pam-fi, Wed Oct 15 12:06:21 2014 UTC
# Line 15  ClassImp(CaloLevel2); Line 15  ClassImp(CaloLevel2);
15  **/  **/
16  CaloTrkVar::CaloTrkVar() {  CaloTrkVar::CaloTrkVar() {
17    this->Clear();    this->Clear();
18  };  }
19    
20  /**  /**
21   * Clear variables   * Clear variables
22  **/  **/
23  void CaloTrkVar::Clear() {  void CaloTrkVar::Clear(Option_t *t) {
24      //void CaloTrkVar::Clear() {
25    trkseqno = 0;    trkseqno = 0;
26    noint = 0;    noint = 0;
27    ncore = 0;    ncore = 0;
# Line 91  CaloLevel2::CaloLevel2() {     Line 92  CaloLevel2::CaloLevel2() {    
92    //    //
93    this->Clear();    this->Clear();
94    //    //
95  };  }
96  /**  /**
97   * Create the TClonesArray   * Create the TClonesArray
98  **/  **/
99  void CaloLevel2::Set(){//ELENA  void CaloLevel2::Set(){//ELENA
100      if(!CaloTrk)CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA      if(!CaloTrk)CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
101  }//ELENA  }//ELENA
102    //--------------------------------------
103    //
104    //
105    //--------------------------------------
106    void CaloLevel2::SetTrackArray(TClonesArray *track){//ELENA
107    //    if(track && strcmp(track->GetClass()->GetName(),"CaloTrkVar")==0){
108        if(track ){
109            if(CaloTrk)CaloTrk->Clear("C");
110            CaloTrk = track;
111        }
112    }
113    
114  /**  /**
115   * Clear the CaloLevel2 object   * Clear the CaloLevel2 object
116   **/   **/
117  void CaloLevel2::Clear() {      void CaloLevel2::Clear(Option_t *t ) {    
118    //    //
119  //  CaloTrk->Clear(); //ELENA    //  CaloTrk->Clear(); //ELENA
120    if(CaloTrk)CaloTrk->Delete(); //ELENA    if(CaloTrk)CaloTrk->Delete(); //ELENA
121    //    //
122    nstrip = 0;    nstrip = 0;
123      nsatstrip = 0;
124    qtot = 0.;    qtot = 0.;
125    //  impx = 0.;    //  impx = 0.;
126    //  impy = 0.;    //  impy = 0.;
# Line 126  void CaloLevel2::Clear() {     Line 139  void CaloLevel2::Clear() {    
139    memset(tany, 0, 2*sizeof(Int_t));    memset(tany, 0, 2*sizeof(Int_t));
140    memset(fitmode, 0, 2*sizeof(Int_t));    memset(fitmode, 0, 2*sizeof(Int_t));
141    memset(planemax, 0, 2*sizeof(Int_t));    memset(planemax, 0, 2*sizeof(Int_t));
142      memset(selfdelay, 0, 4*7*sizeof(Int_t));
143    memset(cibar, 0, 2*22*sizeof(Int_t));    memset(cibar, 0, 2*22*sizeof(Int_t));
144    memset(cbar, 0, 2*22*sizeof(Float_t));    memset(cbar, 0, 2*22*sizeof(Float_t));
145    good = 0;    good = 0;
146    selftrigger = 0;    selftrigger = 0;
147    //    //
148  };  }
149    
150  /**  /**
151   * Delete the CaloLevel2 object   * Delete the CaloLevel2 object
152   **/   **/
153  void CaloLevel2::Delete() {     //ELENA  void CaloLevel2::Delete(Option_t *t) {     //ELENA
154    if(CaloTrk){ //ELENA    if(CaloTrk){ //ELENA
155        CaloTrk->Delete(); //ELENA        CaloTrk->Delete(); //ELENA
156        delete CaloTrk; //ELENA        delete CaloTrk; //ELENA
# Line 144  void CaloLevel2::Delete() {     //ELENA Line 158  void CaloLevel2::Delete() {     //ELENA
158  } //ELENA  } //ELENA
159    
160  /**  /**
161     * CaloLevel2
162    **/
163    Bool_t CaloLevel2::IsGood(Bool_t strict) {  
164      //
165      if ( strict ){
166        if ( !good ) return(false);
167        if ( perr[0] ) return(false);
168        if ( perr[1] ) return(false);
169        if ( perr[2] ) return(false);
170        if ( perr[3] ) return(false);
171        if ( swerr[0] ) return(false);
172        if ( swerr[1] ) return(false);
173        if ( swerr[2] ) return(false);
174        if ( swerr[3] ) return(false);
175        if ( crc[0] ) return(false);
176        if ( crc[1] ) return(false);
177        if ( crc[2] ) return(false);
178        if ( crc[3] ) return(false);
179      } else {
180        if ( perr[0] == 129 ||  perr[0] == 136 || perr[0] == 142 || perr[0] == 143 ) return(false);
181        if ( perr[1] == 129 ||  perr[1] == 136 || perr[1] == 142 || perr[1] == 143 ) return(false);
182        if ( perr[2] == 129 ||  perr[2] == 136 || perr[2] == 142 || perr[2] == 143 ) return(false);
183        if ( perr[3] == 129 ||  perr[3] == 136 || perr[3] == 142 || perr[3] == 143 ) return(false);
184      };
185      //
186      return(true);
187    }
188    
189    /**
190   * 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).
191   */   */
192  void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {  void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {
# Line 230  Float_t CaloLevel2::impx(Int_t tr){ Line 273  Float_t CaloLevel2::impx(Int_t tr){
273      printf(" Cannot find selftrigger block\n");      printf(" Cannot find selftrigger block\n");
274    };    };
275    return(-100.);    return(-100.);
276  };  }
277    
278  /**  /**
279   * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.   * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
# Line 252  Float_t CaloLevel2::impy(Int_t tr){ Line 295  Float_t CaloLevel2::impy(Int_t tr){
295      printf(" Cannot find selftrigger block\n");      printf(" Cannot find selftrigger block\n");
296    };    };
297    return(-100.);    return(-100.);
298  };  }
299  /**  /**
300   * Should return the energy in GeV if the particle would be an electron   * Should return the energy in GeV if the particle would be an electron
301   * using a parametrization taken from Monte Carlo simulation   * using a parametrization taken from Monte Carlo simulation
302  **/  **/
303  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){
304    if ( nstrip == 0 ) return;    if ( nstrip == 0 ) return;
305    energy = qtot * 40.82 * 0.000106;    energy = qtot / 260.;
306      //  energy = qtot * 40.82 * 0.000106;
307    sigma = 0.;    sigma = 0.;
308    if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));    if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));
309    return;    return;
310  };  }
311    
312  /**  /**
313   * Returns pointer to the set of track-related variables "itrk"   * Returns pointer to the set of track-related variables "itrk"
# Line 280  CaloTrkVar *CaloLevel2::GetCaloTrkVar(In Line 324  CaloTrkVar *CaloLevel2::GetCaloTrkVar(In
324    CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];    CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
325    return calotrack;    return calotrack;
326  }  }
327    
328    /**
329     * Retrieves the calorimeter track matching the seqno-th tracker stored track.
330     * (If seqno = -1 retrieves the self-trigger calorimeter track)
331     */
332    CaloTrkVar *CaloLevel2::GetCaloStoredTrack(int seqno){
333    
334      if( ntrk()==0 ){
335        printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i but no Calorimeter tracks are stored\n",seqno);
336        return NULL;
337      };
338      
339      CaloTrkVar *c = 0;
340      Int_t it_calo=0;
341        
342      do {
343        c = GetCaloTrkVar(it_calo);
344        it_calo++;
345      } while( c && seqno != c->trkseqno && it_calo < ntrk());      
346      
347      if(!c || seqno != c->trkseqno){
348        c = 0;
349        if(seqno!=-1 && seqno !=-2 && seqno!=-3 ) printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i does not match Calorimeter stored tracks\n",seqno);
350      };
351      return c;
352        
353    }
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.23