--- calo/flight/CaloPreSampler/src/CaloPreSampler.cpp	2007/12/03 13:57:11	1.7
+++ calo/flight/CaloPreSampler/src/CaloPreSampler.cpp	2009/09/10 12:52:51	1.14
@@ -41,6 +41,9 @@
   rigdefault = 50.;
   nox = false;
   noy = false;
+  forcecalo = false;
+  forcefitmode = -1;
+  memset(mask,0,2*22*sizeof(Int_t));
   //
   Clear();
   //
@@ -93,6 +96,18 @@
   };
 }
 
+void CaloPreSampler::SplitInto(Int_t NoWpreSampler, Int_t NoWcalo){
+  this->SetNoWcalo(0);
+  this->SetNoWpreSampler(0);
+  if ( NoWpreSampler < NoWcalo ){
+	this->SetNoWpreSampler(NoWpreSampler);
+	this->SetNoWcalo(NoWcalo);
+  } else {
+	this->SetNoWcalo(NoWcalo);
+	this->SetNoWpreSampler(NoWpreSampler);
+  };
+}
+
 void CaloPreSampler::Clear(){
   //
   pcalo->Clear();
@@ -136,9 +151,23 @@
     return;
   };
   //
+  // Clear structures used to communicate with fortran
+  //
+  event->ClearStructs();//ELENA
+  if ( forcefitmode > 0 ){
+    if ( forcefitmode != 1000 && forcefitmode != 1001 && forcefitmode != 1002 ){
+      printf(" ERROR! forcefitmode=%i \n Use forcefitmode = 1000 for fit mode 0, 1001 fit mode 1, 1002 fit mode 3\n",forcefitmode);
+      forcefitmode = -1;
+    } else {
+      event->clevel2->fmode[0] = forcefitmode;
+      event->clevel2->fmode[1] = forcefitmode;
+    };
+  };
+  //
   Bool_t newentry = false;
   //
   if ( L2->IsORB() ){
+    if ( debug ) printf(" I am here, we have orbital infos \n");
     if ( L2->GetOrbitalInfo()->pkt_num != PKT || L2->GetOrbitalInfo()->OBT != OBT || L2->GetOrbitalInfo()->absTime != atime || sel != ssel ){
       newentry = true;
       OBT = L2->GetOrbitalInfo()->OBT;
@@ -199,16 +228,19 @@
     //
     // Mask x or y view if nox and/or noy are true (default false)
     //
+    if ( mask[view][plane] ) mip = 0.;
     if ( nox && view == 0 ) mip = 0.; 
     if ( noy && view == 1 ) mip = 0.; 
     //
+    if ( emulate18 && view == 0 && plane == 18 ) mip = 0.;
+    //
     // Selection mode: fill the matrix only for plane < (22 - N) REMEMBER N = number of W planes to be used as presampler, ie if N = 2 then we want to use planes from 0 to 19
     //                 included so plane < (22 - 2)
     //
     if ( sel ){
       if ( plane < (22 - N) ){
 	//
-	if ( emulate18 && plane == (18 - N) ) mip = 0.;
+	if ( emulate18 && plane == (18 - N) && view == 0 ) mip = 0.;
 	if ( plane >= NC ) mip = 0.;
 	event->clevel1->estrip[strip][plane][view] = mip;
 	//
@@ -220,7 +252,7 @@
     if ( cont ){
       if ( plane >= N ){
 	//
-	if ( emulate18 && plane == (18 + N) ) mip = 0.;
+	if ( emulate18 && plane == (18 + N) && view == 0 ) mip = 0.;
 	if ( (plane-N) >= NC ) mip = 0.;
 	event->clevel1->estrip[strip][(plane-N)][view] = mip;
 	//
@@ -357,7 +389,7 @@
   // conditions are: 0) no track from the tracker 1) we have a track fit both in x and y 2) no problems with calo for this event 3) no selftrigger event
   //
   //  if ( trackanyway && !filled && event->clevel2->npcfit[0] >= 2 && event->clevel2->npcfit[1] >= 2 && event->clevel2->good != 0 && event->clevel2->trigty < 2. ){
-  if ( trackanyway && m3 ){
+  if ( (trackanyway && m3) || forcecalo ){
     if ( debug ) printf(" Event with a track not fitted by the tracker \n");
     //
     // Disable "track mode" in the fortran routine
@@ -489,10 +521,11 @@
       if ( debug ) printf(" Selftrigger: problems with event \n");
     };
   };
-  //
-  // Clear structures used to communicate with fortran
-  //
-  event->ClearStructs();
+//   //
+//   // Clear structures used to communicate with fortran
+//   //
+//   event->ClearStructs();
+// ELENA: moved @ beginning
   //
   //
   //
@@ -500,3 +533,50 @@
   if ( debug ) printf(" exit \n");
   //
 }
+
+//
+// Method to add a calorimeter track, evaluated around a tracker track defined by a status vector.
+// (can be used to evaluate the calorimeter track around an arbitrary axis, by setting the status vector with zero deflection )
+// 
+//
+CaloTrkVar* CaloPreSampler::AddCaloTrkVar(float *al,int trktag){
+
+    int ntrkentry = pcalo->ntrk();
+    //
+    for (Int_t nt=0; nt < ntrkentry; nt++){  
+	if( pcalo->GetCaloTrkVar(nt)->trkseqno == trktag){
+	    cout << " CaloTrkVar* CaloPreSampler::AddCaloTrkVar(float *al,int trktag)"<<endl;
+	    cout << " --> trktag = "<<trktag<<" already defined "<<endl;
+	    return NULL;
+	}
+    }
+    //
+    event->clevel1->good2 = 1; //is a trk track
+    event->clevel1->trkchi2 = 0;
+    event->clevel1->hzn = 0; 
+    //
+    // Copy the alpha vector in the input structure
+    //
+    for (Int_t e = 0; e < 5 ; e++){
+	event->clevel1->al_p[e][0] = al[e];
+    };	  
+    //
+    // Get tracker related variables for this track
+    //
+    if ( debug ) printf("track %i Call GetTrkVar() \n",trktag);
+    event->GetTrkVar();
+    if ( debug ) printf(" event->clevel2->dX0l %f \n",event->clevel2->dX0l);
+    //
+    // Save tracker track sequence number
+    //	
+    event->trkseqno = trktag;
+    //
+    // Copy values in the class ca from the structure clevel2
+    //
+    if ( debug ) printf("track %i Call FillTrkVar() \n",trktag);
+    event->FillTrkVar(pcalo,ntrkentry);
+
+    return pcalo->GetCaloTrkVar(ntrkentry);
+    
+
+};//ELENA