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

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

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

revision 1.19 by mocchiut, Mon Dec 14 14:42:53 2009 UTC revision 1.20 by pam-fi, Mon Apr 20 09:12:05 2015 UTC
# Line 7  Line 7 
7  /**  /**
8   * Default constructor   * Default constructor
9   */   */
10  CaloNuclei::CaloNuclei(){  // CaloNuclei::CaloNuclei(){
11    Clear();  //   Clear();
12  };  // };
13    
14  CaloNuclei::CaloNuclei(PamLevel2 *l2p){    CaloNuclei::CaloNuclei(PamLevel2 *l2p,const char* alg){  
15    //    //
16    Clear();    Clear();
17    //    //
# Line 29  CaloNuclei::CaloNuclei(PamLevel2 *l2p){ Line 29  CaloNuclei::CaloNuclei(PamLevel2 *l2p){
29    // debug = true;    // debug = true;
30    usetrack = true;    usetrack = true;
31    usepl18x = false;    usepl18x = false;
32      trkAlg = alg;
33    //    //
34  };  };
35    
# Line 190  void CaloNuclei::Process(Int_t ntr){ Line 191  void CaloNuclei::Process(Int_t ntr){
191    //    //
192    if ( usetrack ){    if ( usetrack ){
193      if ( ntr >= 0 ){      if ( ntr >= 0 ){
194        ptrack = L2->GetTrack(ntr);          ptrack = L2->GetTrack(ntr,trkAlg);
195        if ( ptrack ) track = ptrack->GetCaloTrack();        if ( ptrack ) track = ptrack->GetCaloTrack();
196      } else {      } else {
197        track = L2->GetCaloStoredTrack(ntr);        track = L2->GetCaloStoredTrack(ntr);
198      };      };
199      //      //
200      if ( !track && ntr >= 0 ){      if ( !track && ntr >= 0 ){
201        printf(" ERROR: cannot find any track!\n");        printf(" ERROR: cannot find any track! \n");      
202          cout << "ERROR: trk.algorythm --> "<<trkAlg<<endl;
203        printf(" ERROR: CaloNuclei variables not completely filled \n");        printf(" ERROR: CaloNuclei variables not completely filled \n");
204        return;          return;  
205      };      };
# Line 632  Float_t charge = 1000.; Line 634  Float_t charge = 1000.;
634  Float_t beta = 100.;  Float_t beta = 100.;
635    
636  //------- First try track dependent beta  //------- First try track dependent beta
637      if( L2->GetTrkLevel2()->GetNTracks()>=1 ){      if( L2->GetNTracks(trkAlg)>=1 ){
638      PamTrack *TRKtrack = L2->GetTrack(0);          PamTrack *TRKtrack = L2->GetTrack(0,trkAlg);
639      if (fabs(TRKtrack->GetToFTrack()->beta[12]) < 100.) beta = fabs(TRKtrack->GetToFTrack()->beta[12]);      if (fabs(TRKtrack->GetToFTrack()->beta[12]) < 100.) beta = fabs(TRKtrack->GetToFTrack()->beta[12]);
640                                                  }                                                  }
641  //------- If no beta found, try standalone beta  //------- If no beta found, try standalone beta
# Line 741  Float_t beta = 100.; Line 743  Float_t beta = 100.;
743      if (beta<2.) { // it makes no sense to allow beta=5 or so...      if (beta<2.) { // it makes no sense to allow beta=5 or so...
744                
745            
746        if( L2->GetTrkLevel2()->GetNTracks()>=1 ){        if( L2->GetNTracks(trkAlg)>=1 ){
747          mip=dedx1;          mip=dedx1;
748        }        }
749        if (mip==0) mip=stdedx1;        if (mip==0) mip=stdedx1;
# Line 790  Float_t beta = 100.; Line 792  Float_t beta = 100.;
792    
793      if (beta<2.) { // it makes no sense to allow beta=5 or so...      if (beta<2.) { // it makes no sense to allow beta=5 or so...
794    
795        if( L2->GetTrkLevel2()->GetNTracks()>=1 ){        if( L2->GetNTracks(trkAlg)>=1 ){
796          PamTrack *TRKtrack = L2->GetTrack(0);            PamTrack *TRKtrack = L2->GetTrack(0,trkAlg);
797          mip=dedx1;          mip=dedx1;
798          if (mip==0) mip=stdedx1;          if (mip==0) mip=stdedx1;
799          defl=TRKtrack->GetTrkTrack()->al[4];          defl=TRKtrack->GetTrkTrack()->al[4];
# Line 839  Float_t beta = 100.; Line 841  Float_t beta = 100.;
841    
842      if (beta<2.) { // it makes no sense to allow beta=5 or so...      if (beta<2.) { // it makes no sense to allow beta=5 or so...
843    
844        if( L2->GetTrkLevel2()->GetNTracks()>=1 ){        if( L2->GetNTracks(trkAlg)>=1 ){
845          mip=qNmin1;          mip=qNmin1;
846                    
847          if (mip>0 && defl<0.7 && defl>0)   {          if (mip>0 && defl<0.7 && defl>0)   {
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.23