/[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.14 by mocchiut, Mon Sep 3 08:42:13 2007 UTC revision 1.18 by mocchiut, Thu Dec 4 13:50:43 2008 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  **/  **/
# Line 102  void CaloLevel2::Set(){//ELENA Line 103  void CaloLevel2::Set(){//ELENA
103  /**  /**
104   * Clear the CaloLevel2 object   * Clear the CaloLevel2 object
105   **/   **/
106  void CaloLevel2::Clear() {      void CaloLevel2::Clear(Option_t *t ) {    
107    //    //
108  //  CaloTrk->Clear(); //ELENA    //  CaloTrk->Clear(); //ELENA
109    if(CaloTrk)CaloTrk->Delete(); //ELENA    if(CaloTrk)CaloTrk->Delete(); //ELENA
110    //    //
111    nstrip = 0;    nstrip = 0;
112      nsatstrip = 0;
113    qtot = 0.;    qtot = 0.;
114    //  impx = 0.;    //  impx = 0.;
115    //  impy = 0.;    //  impy = 0.;
# Line 132  void CaloLevel2::Clear() {     Line 134  void CaloLevel2::Clear() {    
134    good = 0;    good = 0;
135    selftrigger = 0;    selftrigger = 0;
136    //    //
137  };  }
138    
139  /**  /**
140   * Delete the CaloLevel2 object   * Delete the CaloLevel2 object
141   **/   **/
142  void CaloLevel2::Delete() {     //ELENA  void CaloLevel2::Delete(Option_t *t) {     //ELENA
143    if(CaloTrk){ //ELENA    if(CaloTrk){ //ELENA
144        CaloTrk->Delete(); //ELENA        CaloTrk->Delete(); //ELENA
145        delete CaloTrk; //ELENA        delete CaloTrk; //ELENA
# Line 145  void CaloLevel2::Delete() {     //ELENA Line 147  void CaloLevel2::Delete() {     //ELENA
147  } //ELENA  } //ELENA
148    
149  /**  /**
150     * CaloLevel2
151    **/
152    Bool_t CaloLevel2::IsGood(Bool_t strict) {  
153      //
154      if ( strict ){
155        if ( !good ) return(false);
156        if ( !perr[0] ) return(false);
157        if ( !perr[1] ) return(false);
158        if ( !perr[2] ) return(false);
159        if ( !perr[3] ) return(false);
160        if ( !swerr[0] ) return(false);
161        if ( !swerr[1] ) return(false);
162        if ( !swerr[2] ) return(false);
163        if ( !swerr[3] ) return(false);
164        if ( !crc[0] ) return(false);
165        if ( !crc[1] ) return(false);
166        if ( !crc[2] ) return(false);
167        if ( !crc[3] ) return(false);
168      } else {
169        if ( perr[0] == 129 ||  perr[0] == 136 || perr[0] == 142 || perr[0] == 143 ) return(false);
170        if ( perr[1] == 129 ||  perr[1] == 136 || perr[1] == 142 || perr[1] == 143 ) return(false);
171        if ( perr[2] == 129 ||  perr[2] == 136 || perr[2] == 142 || perr[2] == 143 ) return(false);
172        if ( perr[3] == 129 ||  perr[3] == 136 || perr[3] == 142 || perr[3] == 143 ) return(false);
173      };
174      //
175      return(true);
176    }
177    
178    /**
179   * 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).
180   */   */
181  void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {  void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {
# Line 231  Float_t CaloLevel2::impx(Int_t tr){ Line 262  Float_t CaloLevel2::impx(Int_t tr){
262      printf(" Cannot find selftrigger block\n");      printf(" Cannot find selftrigger block\n");
263    };    };
264    return(-100.);    return(-100.);
265  };  }
266    
267  /**  /**
268   * 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 253  Float_t CaloLevel2::impy(Int_t tr){ Line 284  Float_t CaloLevel2::impy(Int_t tr){
284      printf(" Cannot find selftrigger block\n");      printf(" Cannot find selftrigger block\n");
285    };    };
286    return(-100.);    return(-100.);
287  };  }
288  /**  /**
289   * 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
290   * using a parametrization taken from Monte Carlo simulation   * using a parametrization taken from Monte Carlo simulation
291  **/  **/
292  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){  void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){
293    if ( nstrip == 0 ) return;    if ( nstrip == 0 ) return;
294    energy = qtot * 40.82 * 0.000106;    energy = qtot / 260.;
295      //  energy = qtot * 40.82 * 0.000106;
296    sigma = 0.;    sigma = 0.;
297    if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));    if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));
298    return;    return;
299  };  }
300    
301  /**  /**
302   * Returns pointer to the set of track-related variables "itrk"   * Returns pointer to the set of track-related variables "itrk"
# Line 307  CaloTrkVar *CaloLevel2::GetCaloStoredTra Line 339  CaloTrkVar *CaloLevel2::GetCaloStoredTra
339    };    };
340    return c;    return c;
341            
342  };  }
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23