CalorimeterLevel2/src/CaloLevel2.cpp

/**
 * \file src/CaloLevel2.cpp
 * \author Emiliano Mocchiutti
 *
**/
#include <TObject.h>
#include <CaloLevel2.h>

//
ClassImp(CaloTrkVar);
ClassImp(CaloLevel2);

/**
 * CaloTrkVar constructor
**/
CaloTrkVar::CaloTrkVar() {
  this->Clear();
}

/**
 * Clear variables
**/
void CaloTrkVar::Clear(Option_t *t) {
  //void CaloTrkVar::Clear() {
  trkseqno = 0;
  noint = 0;
  ncore = 0;
  qcore = 0.;
  ncyl = 0;
  qcyl = 0.;
  qtrack = 0.;
  qtrackx = 0.;
  qtracky = 0.;
  dxtrack = 0.;
  dytrack = 0.;
  qlast = 0.;
  nlast = 0;
  qpre = 0.;
  npre = 0;
  qpresh = 0.;
  npresh = 0;
  qtr = 0.;
  ntr = 0;
  planetot = 0;
  qmean = 0.;
  qlow = 0.;
  nlow = 0;
  dX0l = 0.;
  memset(tbar, 0, 2*22*sizeof(Float_t));
  memset(tibar, 0, 2*22*sizeof(Int_t));
}

/**
 * Copies from t to this
**/
CaloTrkVar::CaloTrkVar(const CaloTrkVar &t){
  trkseqno = t.trkseqno;
  noint = t.noint;
  ncore = t.ncore;
  qcore = t.qcore;
  ncyl = t.ncyl;
  qcyl = t.qcyl;
  qtrack = t.qtrack;
  qtrackx = t.qtrackx;
  qtracky = t.qtracky;
  dxtrack = t.dxtrack;
  dytrack = t.dytrack;
  qlast = t.qlast;
  nlast = t.nlast;
  qpre = t.qpre;
  npre = t.npre;
  qpresh = t.qpresh;
  npresh = t.npresh;
  qtr = t.qtr;
  ntr = t.ntr;  
  planetot = t.planetot;
  qmean = t.qmean;
  dX0l = t.dX0l;
  qlow = t.qlow;
  nlow = t.nlow;
  memcpy(tibar,t.tibar,sizeof(tibar));
  memcpy(tbar,t.tbar,sizeof(tbar));
}

/**
 * CaloLevel2 constructor
**/
CaloLevel2::CaloLevel2() {    
  //
//  CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
  CaloTrk = 0; //ELENA
  //
  this->Clear();
  //
}
/**
 * Create the TClonesArray
**/
void CaloLevel2::Set(){//ELENA
    if(!CaloTrk)CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
}//ELENA

/**
 * Clear the CaloLevel2 object
 **/
void CaloLevel2::Clear(Option_t *t ) {    
  //
  //  CaloTrk->Clear(); //ELENA
  if(CaloTrk)CaloTrk->Delete(); //ELENA
  //
  nstrip = 0;
  nsatstrip = 0;
  qtot = 0.;
  //  impx = 0.;
  //  impy = 0.;
  qmax = 0.;
  nx22 = 0;
  qx22 = 0.;
  elen = 0.;
  selen = 0.;
  memset(perr, 0, 4*sizeof(Int_t));
  memset(swerr, 0, 4*sizeof(Int_t));
  memset(crc, 0, 4*sizeof(Int_t));
  memset(qq, 0, 4*sizeof(Int_t));
  memset(varcfit, 0, 4*sizeof(Float_t));
  memset(npcfit, 0, 4*sizeof(Int_t));
  memset(tanx, 0, 2*sizeof(Int_t));
  memset(tany, 0, 2*sizeof(Int_t));
  memset(fitmode, 0, 2*sizeof(Int_t));
  memset(planemax, 0, 2*sizeof(Int_t));
  memset(selfdelay, 0, 4*7*sizeof(Int_t));
  memset(cibar, 0, 2*22*sizeof(Int_t));
  memset(cbar, 0, 2*22*sizeof(Float_t));
  good = 0;
  selftrigger = 0;
  //
}

/**
 * Delete the CaloLevel2 object
 **/
void CaloLevel2::Delete(Option_t *t) {     //ELENA
  if(CaloTrk){ //ELENA
      CaloTrk->Delete(); //ELENA
      delete CaloTrk; //ELENA
  }     //ELENA
} //ELENA

/**
 * Fills a struct cCaloLevel2 with values from a CaloLevel2 object (to put data into a F77 common).
 */
void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {
 
  l2->good = good;      
  l2->selftrigger = selftrigger;
  l2->nstrip = nstrip;     
  l2->nx22 = nx22;  
  l2->qtot = qtot;     
  l2->qx22 = qx22;     
  l2->qmax = qmax;     
  //  l2->impx = impx;     
  //  l2->impy = impy;     
  //  l2->tanx = tanx;     
  //  l2->tany = tany;  
  l2->elen = elen;      
  l2->selen = selen; 
        
  for(Int_t i=0;i<2;i++){
    l2->planemax[i] = planemax[i];       
    l2->varcfit[i] = varcfit[i];
    l2->npcfit[i] = npcfit[i];
  }
  for(Int_t i=0;i<4;i++){
    l2->perr[i] = perr[i];   
    l2->swerr[i] = swerr[i];  
    l2->calcrc[i] = crc[i];     
    l2->qq[i] = qq[i];  
  }

  if(CaloTrk){ //ELENA
      l2->calntrk = CaloTrk->GetEntries();  
      for(Int_t i=0;i<l2->calntrk;i++){
          l2->caltrkseqno[i] = ((CaloTrkVar *)CaloTrk->At(i))->trkseqno;  
          l2->ncore[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncore;     
          l2->noint[i] = ((CaloTrkVar *)CaloTrk->At(i))->noint;     
          l2->ncyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncyl;      
          l2->nlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlast;     
          l2->npre[i] = ((CaloTrkVar *)CaloTrk->At(i))->npre;      
          l2->npresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->npresh;    
          l2->ntr[i] = ((CaloTrkVar *)CaloTrk->At(i))->ntr;       
          l2->planetot[i] = ((CaloTrkVar *)CaloTrk->At(i))->planetot;  
          l2->nlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlow;      
          l2->qcore[i] =((CaloTrkVar *)CaloTrk->At(i))->qcore ;   
          l2->qcyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->qcyl;    
          l2->qlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlast;   
          l2->qpre[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpre;    
          l2->qpresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpresh;  
          l2->qtr[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtr;     
          l2->qtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrack;  
          l2->qtrackx[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrackx; 
          l2->qtracky[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtracky; 
          l2->dxtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dxtrack; 
          l2->dytrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dytrack; 
          l2->qmean[i] = ((CaloTrkVar *)CaloTrk->At(i))->qmean;   
          l2->qlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlow;    
          l2->dX0l[i] = ((CaloTrkVar *)CaloTrk->At(i))->dX0l;  
          for (Int_t j=0; j<2; j++){
              for (Int_t k=0; k<22; k++){
                  l2->tbar[i][k][j] = ((CaloTrkVar *)CaloTrk->At(i))->tbar[k][j];
              };
          };
      }  
  } //ELENA
}

/**
 * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
 * @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
**/
Float_t CaloLevel2::impx(Int_t tr){
  if ( tr == 0 ) return(cbar[0][0]);
  if ( tr == 1 ) {
    if ( !CaloTrk ) return(-110.);
    TClonesArray &t = *(CaloTrk);
    for (Int_t itrk=0; itrk<ntrk(); itrk++){
      CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
      if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][0]);
    };
  };
  if ( tr !=0 && tr !=1 ){
    printf(" Cannot get impx for other than calo or selftrigger tracks!\n");
  } else {
    printf(" Cannot find selftrigger block\n");
  };
  return(-100.);
}

/**
 * Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
 * @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
**/
Float_t CaloLevel2::impy(Int_t tr){
  if ( tr == 0 ) return(cbar[0][1]);
  if ( tr == 1 ) {
    if ( !CaloTrk ) return(-110.);
    TClonesArray &t = *(CaloTrk);
    for (Int_t itrk=0; itrk<ntrk(); itrk++){
      CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
      if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][1]);
    };
  };
  if ( tr !=0 && tr !=1 ){
    printf(" Cannot get impy for other than calo or selftrigger tracks!\n");
  } else {
    printf(" Cannot find selftrigger block\n");
  };
  return(-100.);
}
/**
 * Should return the energy in GeV if the particle would be an electron
 * using a parametrization taken from Monte Carlo simulation
**/
void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){
  if ( nstrip == 0 ) return;
  energy = qtot / 260.;
  //  energy = qtot * 40.82 * 0.000106;
  sigma = 0.;
  if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));
  return;
}

/**
 * Returns pointer to the set of track-related variables "itrk"
**/
CaloTrkVar *CaloLevel2::GetCaloTrkVar(Int_t itrk){
  //    
  if(itrk >= ntrk()){
    printf(" CaloLevel2 ERROR: track related variables set %i does not exists! \n",itrk);
    printf("                   stored track related variables = %i \n",ntrk());
    return(NULL);
  }
  if(!CaloTrk)return 0; //ELENA
  TClonesArray &t = *(CaloTrk);
  CaloTrkVar *calotrack = (CaloTrkVar*)t[itrk];
  return calotrack;
}

/**
 * Retrieves the calorimeter track matching the seqno-th tracker stored track. 
 * (If seqno = -1 retrieves the self-trigger calorimeter track)
 */
CaloTrkVar *CaloLevel2::GetCaloStoredTrack(int seqno){

  if( ntrk()==0 ){
    printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i but no Calorimeter tracks are stored\n",seqno);
    return NULL;
  };
  
  CaloTrkVar *c = 0;
  Int_t it_calo=0;
     
  do { 
    c = GetCaloTrkVar(it_calo);
    it_calo++; 
  } while( c && seqno != c->trkseqno && it_calo < ntrk());      
  
  if(!c || seqno != c->trkseqno){
    c = 0;
    if(seqno!=-1 && seqno !=-2 && seqno!=-3 ) printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i does not match Calorimeter stored tracks\n",seqno);
  };
  return c;
     
}
1	/**
2	* \file src/CaloLevel2.cpp
3	* \author Emiliano Mocchiutti
4	*
5	**/
6	#include <TObject.h>
7	#include <CaloLevel2.h>
8
9	//
10	ClassImp(CaloTrkVar);
11	ClassImp(CaloLevel2);
12
13	/**
14	* CaloTrkVar constructor
15	**/
16	CaloTrkVar::CaloTrkVar() {
17	this->Clear();
18	}
19
20	/**
21	* Clear variables
22	**/
23	void CaloTrkVar::Clear(Option_t *t) {
24	//void CaloTrkVar::Clear() {
25	trkseqno = 0;
26	noint = 0;
27	ncore = 0;
28	qcore = 0.;
29	ncyl = 0;
30	qcyl = 0.;
31	qtrack = 0.;
32	qtrackx = 0.;
33	qtracky = 0.;
34	dxtrack = 0.;
35	dytrack = 0.;
36	qlast = 0.;
37	nlast = 0;
38	qpre = 0.;
39	npre = 0;
40	qpresh = 0.;
41	npresh = 0;
42	qtr = 0.;
43	ntr = 0;
44	planetot = 0;
45	qmean = 0.;
46	qlow = 0.;
47	nlow = 0;
48	dX0l = 0.;
49	memset(tbar, 0, 222sizeof(Float_t));
50	memset(tibar, 0, 222sizeof(Int_t));
51	}
52
53	/**
54	* Copies from t to this
55	**/
56	CaloTrkVar::CaloTrkVar(const CaloTrkVar &t){
57	trkseqno = t.trkseqno;
58	noint = t.noint;
59	ncore = t.ncore;
60	qcore = t.qcore;
61	ncyl = t.ncyl;
62	qcyl = t.qcyl;
63	qtrack = t.qtrack;
64	qtrackx = t.qtrackx;
65	qtracky = t.qtracky;
66	dxtrack = t.dxtrack;
67	dytrack = t.dytrack;
68	qlast = t.qlast;
69	nlast = t.nlast;
70	qpre = t.qpre;
71	npre = t.npre;
72	qpresh = t.qpresh;
73	npresh = t.npresh;
74	qtr = t.qtr;
75	ntr = t.ntr;
76	planetot = t.planetot;
77	qmean = t.qmean;
78	dX0l = t.dX0l;
79	qlow = t.qlow;
80	nlow = t.nlow;
81	memcpy(tibar,t.tibar,sizeof(tibar));
82	memcpy(tbar,t.tbar,sizeof(tbar));
83	}
84
85	/**
86	* CaloLevel2 constructor
87	**/
88	CaloLevel2::CaloLevel2() {
89	//
90	// CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
91	CaloTrk = 0; //ELENA
92	//
93	this->Clear();
94	//
95	}
96	/**
97	* Create the TClonesArray
98	**/
99	void CaloLevel2::Set(){//ELENA
100	if(!CaloTrk)CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
101	}//ELENA
102
103	/**
104	* Clear the CaloLevel2 object
105	**/
106	void CaloLevel2::Clear(Option_t *t ) {
107	//
108	// CaloTrk->Clear(); //ELENA
109	if(CaloTrk)CaloTrk->Delete(); //ELENA
110	//
111	nstrip = 0;
112	nsatstrip = 0;
113	qtot = 0.;
114	// impx = 0.;
115	// impy = 0.;
116	qmax = 0.;
117	nx22 = 0;
118	qx22 = 0.;
119	elen = 0.;
120	selen = 0.;
121	memset(perr, 0, 4*sizeof(Int_t));
122	memset(swerr, 0, 4*sizeof(Int_t));
123	memset(crc, 0, 4*sizeof(Int_t));
124	memset(qq, 0, 4*sizeof(Int_t));
125	memset(varcfit, 0, 4*sizeof(Float_t));
126	memset(npcfit, 0, 4*sizeof(Int_t));
127	memset(tanx, 0, 2*sizeof(Int_t));
128	memset(tany, 0, 2*sizeof(Int_t));
129	memset(fitmode, 0, 2*sizeof(Int_t));
130	memset(planemax, 0, 2*sizeof(Int_t));
131	memset(selfdelay, 0, 47sizeof(Int_t));
132	memset(cibar, 0, 222sizeof(Int_t));
133	memset(cbar, 0, 222sizeof(Float_t));
134	good = 0;
135	selftrigger = 0;
136	//
137	}
138
139	/**
140	* Delete the CaloLevel2 object
141	**/
142	void CaloLevel2::Delete(Option_t *t) { //ELENA
143	if(CaloTrk){ //ELENA
144	CaloTrk->Delete(); //ELENA
145	delete CaloTrk; //ELENA
146	} //ELENA
147	} //ELENA
148
149	/**
150	* Fills a struct cCaloLevel2 with values from a CaloLevel2 object (to put data into a F77 common).
151	*/
152	void CaloLevel2::GetLevel2Struct(cCaloLevel2 *l2) const {
153
154	l2->good = good;
155	l2->selftrigger = selftrigger;
156	l2->nstrip = nstrip;
157	l2->nx22 = nx22;
158	l2->qtot = qtot;
159	l2->qx22 = qx22;
160	l2->qmax = qmax;
161	// l2->impx = impx;
162	// l2->impy = impy;
163	// l2->tanx = tanx;
164	// l2->tany = tany;
165	l2->elen = elen;
166	l2->selen = selen;
167
168	for(Int_t i=0;i<2;i++){
169	l2->planemax[i] = planemax[i];
170	l2->varcfit[i] = varcfit[i];
171	l2->npcfit[i] = npcfit[i];
172	}
173	for(Int_t i=0;i<4;i++){
174	l2->perr[i] = perr[i];
175	l2->swerr[i] = swerr[i];
176	l2->calcrc[i] = crc[i];
177	l2->qq[i] = qq[i];
178	}
179
180	if(CaloTrk){ //ELENA
181	l2->calntrk = CaloTrk->GetEntries();
182	for(Int_t i=0;i<l2->calntrk;i++){
183	l2->caltrkseqno[i] = ((CaloTrkVar *)CaloTrk->At(i))->trkseqno;
184	l2->ncore[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncore;
185	l2->noint[i] = ((CaloTrkVar *)CaloTrk->At(i))->noint;
186	l2->ncyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->ncyl;
187	l2->nlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlast;
188	l2->npre[i] = ((CaloTrkVar *)CaloTrk->At(i))->npre;
189	l2->npresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->npresh;
190	l2->ntr[i] = ((CaloTrkVar *)CaloTrk->At(i))->ntr;
191	l2->planetot[i] = ((CaloTrkVar *)CaloTrk->At(i))->planetot;
192	l2->nlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->nlow;
193	l2->qcore[i] =((CaloTrkVar *)CaloTrk->At(i))->qcore ;
194	l2->qcyl[i] = ((CaloTrkVar *)CaloTrk->At(i))->qcyl;
195	l2->qlast[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlast;
196	l2->qpre[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpre;
197	l2->qpresh[i] = ((CaloTrkVar *)CaloTrk->At(i))->qpresh;
198	l2->qtr[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtr;
199	l2->qtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrack;
200	l2->qtrackx[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtrackx;
201	l2->qtracky[i] = ((CaloTrkVar *)CaloTrk->At(i))->qtracky;
202	l2->dxtrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dxtrack;
203	l2->dytrack[i] = ((CaloTrkVar *)CaloTrk->At(i))->dytrack;
204	l2->qmean[i] = ((CaloTrkVar *)CaloTrk->At(i))->qmean;
205	l2->qlow[i] = ((CaloTrkVar *)CaloTrk->At(i))->qlow;
206	l2->dX0l[i] = ((CaloTrkVar *)CaloTrk->At(i))->dX0l;
207	for (Int_t j=0; j<2; j++){
208	for (Int_t k=0; k<22; k++){
209	l2->tbar[i][k][j] = ((CaloTrkVar *)CaloTrk->At(i))->tbar[k][j];
210	};
211	};
212	}
213	} //ELENA
214	}
215
216	/**
217	* Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
218	* @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
219	**/
220	Float_t CaloLevel2::impx(Int_t tr){
221	if ( tr == 0 ) return(cbar[0][0]);
222	if ( tr == 1 ) {
223	if ( !CaloTrk ) return(-110.);
224	TClonesArray &t = *(CaloTrk);
225	for (Int_t itrk=0; itrk<ntrk(); itrk++){
226	CaloTrkVar calotrack = (CaloTrkVar)t[itrk];
227	if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][0]);
228	};
229	};
230	if ( tr !=0 && tr !=1 ){
231	printf(" Cannot get impx for other than calo or selftrigger tracks!\n");
232	} else {
233	printf(" Cannot find selftrigger block\n");
234	};
235	return(-100.);
236	}
237
238	/**
239	* Returns the impact position on the top of the calorimeter as determined by the calorimeter itself.
240	* @param tr : if tr = 0 use the calorimeter "normal" fit, if 1 use the calorimeter "selftrigger" fit (if any!)
241	**/
242	Float_t CaloLevel2::impy(Int_t tr){
243	if ( tr == 0 ) return(cbar[0][1]);
244	if ( tr == 1 ) {
245	if ( !CaloTrk ) return(-110.);
246	TClonesArray &t = *(CaloTrk);
247	for (Int_t itrk=0; itrk<ntrk(); itrk++){
248	CaloTrkVar calotrack = (CaloTrkVar)t[itrk];
249	if ( calotrack->trkseqno == -1 ) return(calotrack->tbar[0][1]);
250	};
251	};
252	if ( tr !=0 && tr !=1 ){
253	printf(" Cannot get impy for other than calo or selftrigger tracks!\n");
254	} else {
255	printf(" Cannot find selftrigger block\n");
256	};
257	return(-100.);
258	}
259	/**
260	* Should return the energy in GeV if the particle would be an electron
261	* using a parametrization taken from Monte Carlo simulation
262	**/
263	void CaloLevel2::GetElectronEnergy(Float_t &energy, Float_t &sigma){
264	if ( nstrip == 0 ) return;
265	energy = qtot / 260.;
266	// energy = qtot * 40.82 * 0.000106;
267	sigma = 0.;
268	if ( energy > 0. ) sigma = energy * (0.01183 + 0.121/sqrt(energy));
269	return;
270	}
271
272	/**
273	* Returns pointer to the set of track-related variables "itrk"
274	**/
275	CaloTrkVar *CaloLevel2::GetCaloTrkVar(Int_t itrk){
276	//
277	if(itrk >= ntrk()){
278	printf(" CaloLevel2 ERROR: track related variables set %i does not exists! \n",itrk);
279	printf(" stored track related variables = %i \n",ntrk());
280	return(NULL);
281	}
282	if(!CaloTrk)return 0; //ELENA
283	TClonesArray &t = *(CaloTrk);
284	CaloTrkVar calotrack = (CaloTrkVar)t[itrk];
285	return calotrack;
286	}
287
288	/**
289	* Retrieves the calorimeter track matching the seqno-th tracker stored track.
290	* (If seqno = -1 retrieves the self-trigger calorimeter track)
291	*/
292	CaloTrkVar *CaloLevel2::GetCaloStoredTrack(int seqno){
293
294	if( ntrk()==0 ){
295	printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i but no Calorimeter tracks are stored\n",seqno);
296	return NULL;
297	};
298
299	CaloTrkVar *c = 0;
300	Int_t it_calo=0;
301
302	do {
303	c = GetCaloTrkVar(it_calo);
304	it_calo++;
305	} while( c && seqno != c->trkseqno && it_calo < ntrk());
306
307	if(!c \|\| seqno != c->trkseqno){
308	c = 0;
309	if(seqno!=-1 && seqno !=-2 && seqno!=-3 ) printf("CaloLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo %i does not match Calorimeter stored tracks\n",seqno);
310	};
311	return c;
312
313	}